You can subscribe to this list here.
2002 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}

_{Jun}
(67) 
_{Jul}
(61) 
_{Aug}
(49) 
_{Sep}
(43) 
_{Oct}
(59) 
_{Nov}
(24) 
_{Dec}
(18) 

2003 
_{Jan}
(34) 
_{Feb}
(35) 
_{Mar}
(72) 
_{Apr}
(42) 
_{May}
(46) 
_{Jun}
(15) 
_{Jul}
(64) 
_{Aug}
(62) 
_{Sep}
(22) 
_{Oct}
(41) 
_{Nov}
(57) 
_{Dec}
(56) 
2004 
_{Jan}
(48) 
_{Feb}
(47) 
_{Mar}
(33) 
_{Apr}
(39) 
_{May}
(6) 
_{Jun}
(17) 
_{Jul}
(19) 
_{Aug}
(10) 
_{Sep}
(14) 
_{Oct}
(74) 
_{Nov}
(80) 
_{Dec}
(22) 
2005 
_{Jan}
(43) 
_{Feb}
(33) 
_{Mar}
(52) 
_{Apr}
(74) 
_{May}
(32) 
_{Jun}
(58) 
_{Jul}
(18) 
_{Aug}
(41) 
_{Sep}
(71) 
_{Oct}
(28) 
_{Nov}
(65) 
_{Dec}
(68) 
2006 
_{Jan}
(54) 
_{Feb}
(37) 
_{Mar}
(82) 
_{Apr}
(211) 
_{May}
(69) 
_{Jun}
(75) 
_{Jul}
(279) 
_{Aug}
(139) 
_{Sep}
(135) 
_{Oct}
(58) 
_{Nov}
(81) 
_{Dec}
(78) 
2007 
_{Jan}
(141) 
_{Feb}
(134) 
_{Mar}
(65) 
_{Apr}
(49) 
_{May}
(61) 
_{Jun}
(90) 
_{Jul}
(72) 
_{Aug}
(53) 
_{Sep}
(86) 
_{Oct}
(61) 
_{Nov}
(62) 
_{Dec}
(101) 
2008 
_{Jan}
(100) 
_{Feb}
(66) 
_{Mar}
(76) 
_{Apr}
(95) 
_{May}
(77) 
_{Jun}
(93) 
_{Jul}
(103) 
_{Aug}
(76) 
_{Sep}
(42) 
_{Oct}
(55) 
_{Nov}
(44) 
_{Dec}
(75) 
2009 
_{Jan}
(103) 
_{Feb}
(105) 
_{Mar}
(121) 
_{Apr}
(59) 
_{May}
(103) 
_{Jun}
(82) 
_{Jul}
(67) 
_{Aug}
(76) 
_{Sep}
(85) 
_{Oct}
(75) 
_{Nov}
(181) 
_{Dec}
(133) 
2010 
_{Jan}
(107) 
_{Feb}
(116) 
_{Mar}
(145) 
_{Apr}
(89) 
_{May}
(138) 
_{Jun}
(85) 
_{Jul}
(82) 
_{Aug}
(111) 
_{Sep}
(70) 
_{Oct}
(83) 
_{Nov}
(60) 
_{Dec}
(16) 
2011 
_{Jan}
(61) 
_{Feb}
(16) 
_{Mar}
(52) 
_{Apr}
(41) 
_{May}
(34) 
_{Jun}
(41) 
_{Jul}
(57) 
_{Aug}
(73) 
_{Sep}
(21) 
_{Oct}
(45) 
_{Nov}
(50) 
_{Dec}
(28) 
2012 
_{Jan}
(70) 
_{Feb}
(36) 
_{Mar}
(71) 
_{Apr}
(29) 
_{May}
(48) 
_{Jun}
(61) 
_{Jul}
(44) 
_{Aug}
(54) 
_{Sep}
(20) 
_{Oct}
(28) 
_{Nov}
(41) 
_{Dec}
(137) 
2013 
_{Jan}
(62) 
_{Feb}
(55) 
_{Mar}
(31) 
_{Apr}
(23) 
_{May}
(54) 
_{Jun}
(54) 
_{Jul}
(90) 
_{Aug}
(46) 
_{Sep}
(38) 
_{Oct}
(60) 
_{Nov}
(92) 
_{Dec}
(17) 
2014 
_{Jan}
(62) 
_{Feb}
(35) 
_{Mar}
(72) 
_{Apr}
(30) 
_{May}
(97) 
_{Jun}
(81) 
_{Jul}
(63) 
_{Aug}
(64) 
_{Sep}
(28) 
_{Oct}
(45) 
_{Nov}
(48) 
_{Dec}
(109) 
2015 
_{Jan}
(106) 
_{Feb}
(36) 
_{Mar}
(65) 
_{Apr}
(63) 
_{May}
(95) 
_{Jun}
(56) 
_{Jul}
(48) 
_{Aug}
(55) 
_{Sep}
(100) 
_{Oct}
(57) 
_{Nov}
(33) 
_{Dec}
(46) 
2016 
_{Jan}
(76) 
_{Feb}
(53) 
_{Mar}
(88) 
_{Apr}
(79) 
_{May}
(62) 
_{Jun}
(65) 
_{Jul}
(37) 
_{Aug}
(23) 
_{Sep}
(108) 
_{Oct}
(68) 
_{Nov}
(66) 
_{Dec}
(47) 
2017 
_{Jan}
(55) 
_{Feb}
(11) 
_{Mar}
(30) 
_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 





1
(4) 
2
(2) 
3
(3) 
4
(10) 
5
(2) 
6
(4) 
7
(4) 
8
(5) 
9
(1) 
10

11
(11) 
12
(2) 
13
(1) 
14
(8) 
15
(1) 
16
(3) 
17
(4) 
18
(5) 
19
(14) 
20
(5) 
21
(5) 
22
(15) 
23
(16) 
24
(1) 
25
(3) 
26
(3) 
27

28
(2) 



From: SourceForge.net <noreply@so...>  20070222 19:43:04

Bugs item #1654602, was opened at 20070207 16:43 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1654602&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Limit Group: None Status: Open Resolution: None Priority: 3 Private: No Submitted By: Kirill Smelkov (kirr) Assigned to: Nobody/Anonymous (nobody) Summary: limit(%e^b*erf(x), x, inf) asks sign(b) Initial Comment: (%i10) limit(%e^b*erf(x), x, inf); Is b positive, negative, or zero? pos; b (%o10) %e Maxima version: 5.11.0cvs Maxima build date: 21:36 2/7/2007 host type: i686pclinuxgnu lispimplementationtype: GNU Common Lisp (GCL) lispimplementationversion: GCL 2.6.6  >Comment By: Raymond Toy (rtoy) Date: 20070222 14:40 Message: Logged In: YES user_id=28849 Originator: NO For the record, the question comes from posnegp in forq. Perhaps limit1 should factor out terms not involving the limit variable?  Comment By: Stavros Macrakis (macrakis) Date: 20070220 17:49 Message: Logged In: YES user_id=588346 Originator: NO Reduced priority to 3. Result is correct, even if question is unnecessary. Maxima alas asks lots of unnecessary questions.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1654602&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 19:26:39

Bugs item #1664729, was opened at 20070220 18:00 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664729&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sin(x)cos(x),x,0,inf) wrong Initial Comment: integrate(sin(x)cos(x),x,0,inf) => 0 (!!!) This seems to be caused by: integrate(sin(x),x,0,inf) => IND I suppose this could be interpreted as correct, but integrate is supposed to give an "Integral is divergent" error, as it still does in other cases, e.g. integrate(1/x,x,0,inf). And then, IND doesn't simplify properly, so you get nonsense like integrate(sin(x)cos(x),x,0,inf) => 0 (from INDIND) Along the same lines, integrate(sin(x^2)*x,x,0,inf) charmingly gives ind+1/2, which is sort of right, but didn't simplify to ind. 5.11.0 GCL 2.6.8 Windows 2k  >Comment By: Raymond Toy (rtoy) Date: 20070222 14:26 Message: Logged In: YES user_id=28849 Originator: NO Don't know about the conventions, but rev 1.1 of defint.lisp shows that scaxn returns ind. However, I don't know if the original version computes this integral via scaxn. I assume it did. (We really need to make ind, inf, minf, und, and friends contagious.)  Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:45 Message: Logged In: YES user_id=588346 Originator: YES Re > integrate(sin(x),x,0,inf) returns ind because scaxn returns ind. A quick look at the code seems to show that it's prepared to handle IND when returned by LIMIT (takeprincipal etc.), but it doesn't generate it internally. Perhaps this is new code that doesn't follow defint's conventions?  Comment By: Raymond Toy (rtoy) Date: 20070222 12:26 Message: Logged In: YES user_id=28849 Originator: NO integrate(sin(x),x,0,inf) returns ind because scaxn returns ind.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664729&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 19:20:08

Bugs item #1664726, was opened at 20070220 17:56 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Closed Resolution: Invalid Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sin(x^2),x,0,inf) wrong Initial Comment: integrate(sin(x^2),x,0,inf) => sqrt(%pi)/(2*sqrt(2)) but the integral is not convergent  after all, integrate(sin(x^2),x,0,a) oscillates by +/1 all the way to infinity (as it were). 5.11.0 GCL 2.6.8 Windows 2k  >Comment By: Raymond Toy (rtoy) Date: 20070222 14:20 Message: Logged In: YES user_id=28849 Originator: NO My typo. The thesis says n > 1, as does the code.  Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:41 Message: Logged In: YES user_id=588346 Originator: YES I was wrong about the original case, sin(x^2). That *DOES* converge. My mistake. I will close the bug. But it doesn't converge for n=1. Maybe n>1?  Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:38 Message: Logged In: YES user_id=588346 Originator: YES > Wang's thesis on page 81 considers integrals of the form sin(k*x^n) > for n >= 1. He proves that these integrals do exist via contour > integration over a sector of a circle. integrate(Q,x,0,inf) = limit( integrate(Q,x,0,a), a, inf ). If Q=sin(x) (sin(k*x^n) with k=1,n=1), then integrate(Q,x,0,a) = 1cos(a) limit(1cos(a),a,inf) is undefined (ind). Therefore the integral is divergent. A quick look at the code shows that it is set up to handle IND when it comes back from Limit, but not to use it internally, as in  Comment By: Raymond Toy (rtoy) Date: 20070222 11:57 Message: Logged In: YES user_id=28849 Originator: NO Hmm. Wang's thesis on page 81 considers integrals of the form sin(k*x^n) for n >= 1. He proves that these integrals do exist via contour integration over a sector of a circle.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 18:45:33

Bugs item #1664729, was opened at 20070220 18:00 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664729&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sin(x)cos(x),x,0,inf) wrong Initial Comment: integrate(sin(x)cos(x),x,0,inf) => 0 (!!!) This seems to be caused by: integrate(sin(x),x,0,inf) => IND I suppose this could be interpreted as correct, but integrate is supposed to give an "Integral is divergent" error, as it still does in other cases, e.g. integrate(1/x,x,0,inf). And then, IND doesn't simplify properly, so you get nonsense like integrate(sin(x)cos(x),x,0,inf) => 0 (from INDIND) Along the same lines, integrate(sin(x^2)*x,x,0,inf) charmingly gives ind+1/2, which is sort of right, but didn't simplify to ind. 5.11.0 GCL 2.6.8 Windows 2k  >Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:45 Message: Logged In: YES user_id=588346 Originator: YES Re > integrate(sin(x),x,0,inf) returns ind because scaxn returns ind. A quick look at the code seems to show that it's prepared to handle IND when returned by LIMIT (takeprincipal etc.), but it doesn't generate it internally. Perhaps this is new code that doesn't follow defint's conventions?  Comment By: Raymond Toy (rtoy) Date: 20070222 12:26 Message: Logged In: YES user_id=28849 Originator: NO integrate(sin(x),x,0,inf) returns ind because scaxn returns ind.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664729&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 18:41:47

Bugs item #1664726, was opened at 20070220 17:56 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None >Status: Closed >Resolution: Invalid Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sin(x^2),x,0,inf) wrong Initial Comment: integrate(sin(x^2),x,0,inf) => sqrt(%pi)/(2*sqrt(2)) but the integral is not convergent  after all, integrate(sin(x^2),x,0,a) oscillates by +/1 all the way to infinity (as it were). 5.11.0 GCL 2.6.8 Windows 2k  >Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:41 Message: Logged In: YES user_id=588346 Originator: YES I was wrong about the original case, sin(x^2). That *DOES* converge. My mistake. I will close the bug. But it doesn't converge for n=1. Maybe n>1?  Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:38 Message: Logged In: YES user_id=588346 Originator: YES > Wang's thesis on page 81 considers integrals of the form sin(k*x^n) > for n >= 1. He proves that these integrals do exist via contour > integration over a sector of a circle. integrate(Q,x,0,inf) = limit( integrate(Q,x,0,a), a, inf ). If Q=sin(x) (sin(k*x^n) with k=1,n=1), then integrate(Q,x,0,a) = 1cos(a) limit(1cos(a),a,inf) is undefined (ind). Therefore the integral is divergent. A quick look at the code shows that it is set up to handle IND when it comes back from Limit, but not to use it internally, as in  Comment By: Raymond Toy (rtoy) Date: 20070222 11:57 Message: Logged In: YES user_id=28849 Originator: NO Hmm. Wang's thesis on page 81 considers integrals of the form sin(k*x^n) for n >= 1. He proves that these integrals do exist via contour integration over a sector of a circle.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 18:40:04

Bugs item #1665657, was opened at 20070221 17:30 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1665657&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Limit Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Robert Dodier (robert_dodier) Assigned to: Nobody/Anonymous (nobody) Summary: limit fails to find easy limit Initial Comment: limit(x/(x1)1/log(x),x,1,plus); => 'limit(x1/log(1/x+1)+1,x,inf) But Maxima can handle this with a little encouragement: tlimit(x/(x1)1/log(x),x,1,plus) => 1/2 limit(ratsimp(x/(x1)1/log(x)),x,1,plus) => 1/2 Reported on mailing list 2007/02/21.  >Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:40 Message: Logged In: YES user_id=588346 Originator: NO My reporting error. Actually: limit(ratsimp(x/(x1)1/log(x)),x,1,plus) fails though limit(ratsimp(x/(x1)1/log(x)),x,1) => 1/2 Clearly if the twosided limit exists, the onesided limits must be equal to it....  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1665657&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 18:38:17

Bugs item #1664726, was opened at 20070220 17:56 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sin(x^2),x,0,inf) wrong Initial Comment: integrate(sin(x^2),x,0,inf) => sqrt(%pi)/(2*sqrt(2)) but the integral is not convergent  after all, integrate(sin(x^2),x,0,a) oscillates by +/1 all the way to infinity (as it were). 5.11.0 GCL 2.6.8 Windows 2k  >Comment By: Stavros Macrakis (macrakis) Date: 20070222 13:38 Message: Logged In: YES user_id=588346 Originator: YES > Wang's thesis on page 81 considers integrals of the form sin(k*x^n) > for n >= 1. He proves that these integrals do exist via contour > integration over a sector of a circle. integrate(Q,x,0,inf) = limit( integrate(Q,x,0,a), a, inf ). If Q=sin(x) (sin(k*x^n) with k=1,n=1), then integrate(Q,x,0,a) = 1cos(a) limit(1cos(a),a,inf) is undefined (ind). Therefore the integral is divergent. A quick look at the code shows that it is set up to handle IND when it comes back from Limit, but not to use it internally, as in  Comment By: Raymond Toy (rtoy) Date: 20070222 11:57 Message: Logged In: YES user_id=28849 Originator: NO Hmm. Wang's thesis on page 81 considers integrals of the form sin(k*x^n) for n >= 1. He proves that these integrals do exist via contour integration over a sector of a circle.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 17:31:20

Bugs item #1666354, was opened at 20070222 12:31 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1666354&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: None Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: limit(gamma(1/x)*sin(%pi/(2*x))/x,x,inf) wrong Initial Comment: limit(gamma(1/x)*sin(%pi/(2*x))/x,x,inf) returns a noun form. tlimit returns 0, which is right.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1666354&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 17:26:42

Bugs item #1664729, was opened at 20070220 18:00 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664729&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sin(x)cos(x),x,0,inf) wrong Initial Comment: integrate(sin(x)cos(x),x,0,inf) => 0 (!!!) This seems to be caused by: integrate(sin(x),x,0,inf) => IND I suppose this could be interpreted as correct, but integrate is supposed to give an "Integral is divergent" error, as it still does in other cases, e.g. integrate(1/x,x,0,inf). And then, IND doesn't simplify properly, so you get nonsense like integrate(sin(x)cos(x),x,0,inf) => 0 (from INDIND) Along the same lines, integrate(sin(x^2)*x,x,0,inf) charmingly gives ind+1/2, which is sort of right, but didn't simplify to ind. 5.11.0 GCL 2.6.8 Windows 2k  >Comment By: Raymond Toy (rtoy) Date: 20070222 12:26 Message: Logged In: YES user_id=28849 Originator: NO integrate(sin(x),x,0,inf) returns ind because scaxn returns ind.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664729&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 16:57:24

Bugs item #1664726, was opened at 20070220 17:56 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sin(x^2),x,0,inf) wrong Initial Comment: integrate(sin(x^2),x,0,inf) => sqrt(%pi)/(2*sqrt(2)) but the integral is not convergent  after all, integrate(sin(x^2),x,0,a) oscillates by +/1 all the way to infinity (as it were). 5.11.0 GCL 2.6.8 Windows 2k  >Comment By: Raymond Toy (rtoy) Date: 20070222 11:57 Message: Logged In: YES user_id=28849 Originator: NO Hmm. Wang's thesis on page 81 considers integrals of the form sin(k*x^n) for n >= 1. He proves that these integrals do exist via contour integration over a sector of a circle.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1664726&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 11:24:09

Bugs item #1663536, was opened at 20070219 08:50 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1663536&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: None Group: None >Status: Closed >Resolution: Invalid Priority: 5 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: newton sometines very slow Initial Comment: Fast: (%i19) newton(x^2  %i, 1  %i/5); (%o19) 7.071067853080116b1*%i+7.071067813155526b1 So slow I don't know that it ever finishes: (%i20) newton(x^2  %i, 1  %i);  >Comment By: Barton Willis (willisbl) Date: 20070222 05:24 Message: Logged In: YES user_id=895922 Originator: YES You are correct. I was confused.  Comment By: Raymond Toy (rtoy) Date: 20070220 10:19 Message: Logged In: YES user_id=28849 Originator: NO Is newton required to converge from this starting point? I modified newton to print out each guess. It looks like newton is stuck, and keeps moving back and forth on some line, constantly producing guesses of the form aa*%i. Plugging this into the newton iteration shows that the next guess would be (1%i)*(2*a^21)/(4*a), so the form of the next answer never changes.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1663536&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 11:12:10

Bugs item #1666011, was opened at 20070222 05:12 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1666011&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Share Libraries Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: mnewton & Polynomial quotient is not exact Initial Comment: (%i1) load(mnewton); (%i2) mnewton([x^y  %i,x+y2],[x,y],[1,1.5]); Polynomial quotient is not exact I can't eliminate the error by setting gcd : spmod.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1666011&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 03:20:25

Bugs item #1644378, was opened at 20070125 04:48 Message generated for change (Comment added) made by sfrobot You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1644378&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None >Status: Closed Resolution: Fixed Priority: 5 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: floor and ceiling change bigfloatone Initial Comment: Look at bigfloatone: (%i1) ?print(?bigfloatone); ((BIGFLOAT SIMP 56) 36028797018963968 1) Compute a ceiling of a constant: (%i2) ceiling(log(1.3b0)/log(2)); (%o2) 1 Look at bigfloatone: Yikes! It changed (%i3) ?print(?bigfloatone); ((BIGFLOAT SIMP 218) <junk>) Now compute log(2.0b0) (%i4) log(2.0b0); (%o4) 6.9314718055[lots of digits] 0095b1 But fpprec is still the default. (%i5) fpprec; (%o5) 16  >Comment By: SourceForge Robot (sfrobot) Date: 20070221 19:20 Message: Logged In: YES user_id=1312539 Originator: NO This Tracker item was closed automatically by the system. It was previously set to a Pending status, and the original submitter did not respond within 14 days (the time period specified by the administrator of this Tracker).  Comment By: Barton Willis (willisbl) Date: 20070207 11:18 Message: Logged In: YES user_id=895922 Originator: YES Fixed by nummod.lisp CVS r 1.4; also added test in rtest_allnumod.mac.  Comment By: Stavros Macrakis (macrakis) Date: 20070206 15:03 Message: Logged In: YES user_id=588346 Originator: NO This will do the trick: (defmacro bindfpprec (val &rest exprs) `(let ($fpprec fpprec bigfloatzero bigfloatone bfhalf bfmhalf) (fpprec1 nil ,val) ,@exprs)) I was actually just looking at this code....  Comment By: Barton Willis (willisbl) Date: 20070125 05:33 Message: Logged In: YES user_id=895922 Originator: YES I think maxmin.lisp needs function something like (defun protectedbigfloat (e n) (unwindprotect (progn (fpprec1 '$fpprec n) (simplify ($bfloat e))) (fpprec1 '$fpprec $fpprec))) Comments?  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1644378&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 03:20:25

Bugs item #1653672, was opened at 20070206 14:56 Message generated for change (Comment added) made by sfrobot You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1653672&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Simplification Group: None >Status: Closed Resolution: Fixed Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Barton Willis (willisbl) Summary: floor(...) gives internal error Initial Comment: ceiling(asin(sin(sqrt(2))^2+cos(sqrt(2))^21/7)) => Error ((MPLUS SIMP) ... $%I)))) is not of type NUMBER. I sort of understand how this might happen if it were 2^100 and not 1/7, but even in that case it should at worst return a noun form, not an error. Compare floor(1/2sin(%pi/4+q)^2). For q=2150 (not q=1/7 !!), it correctly gives 1, for q=^2160 it sensibly gives a noun form. Maxima 5.11.0 GCL 2.6.8  >Comment By: SourceForge Robot (sfrobot) Date: 20070221 19:20 Message: Logged In: YES user_id=1312539 Originator: NO This Tracker item was closed automatically by the system. It was previously set to a Pending status, and the original submitter did not respond within 14 days (the time period specified by the administrator of this Tracker).  Comment By: Barton Willis (willisbl) Date: 20070207 11:17 Message: Logged In: YES user_id=895922 Originator: NO Fixed by nummod.lisp CVS r 1.4; also added test in rtest_allnumod.mac.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1653672&group_id=4933 
From: SourceForge.net <noreply@so...>  20070222 02:58:51

Bugs item #1665777, was opened at 20070221 20:58 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1665777&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: None Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: determ broken Initial Comment: /share/matrix/determ is broken: (%i1) load(determ)$ (%i2) m : genmatrix(lambda([i,j],1/(1+i+j)),4,4)$ (%i3) det(m); (%o3) 105 < wrong (%i4) det(m); (%o4) 0 < wrong again, but different  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1665777&group_id=4933 