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}
(59) 
_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 




1
(1) 
2

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

12
(3) 
13

14
(3) 
15
(1) 
16

17
(1) 
18
(4) 
19
(1) 
20
(1) 
21

22

23

24
(1) 
25
(1) 
26
(3) 
27
(3) 
28
(1) 
29
(8) 
30
(4) 


From: SourceForge.net <noreply@so...>  20061109 23:45:35

Bugs item #635606, was opened at 20021108 12:47 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=635606&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(abs(log(x))) internal error, UND Initial Comment: Maxima 5.5/Windows/gcl limit(abs(log(x)),x,0) Error: The tag LIMIT is undefined. Should of course be INF. More controversially, perhaps, limit(log(x),x,0) gives UND  I believe it should give INFINITY; after all, limit(log (x),x,0,minus) gives INFINITY.  >Comment By: Stavros Macrakis (macrakis) Date: 20061109 18:45 Message: Logged In: YES user_id=588346 re limit(abs(log(x)),x,0) If we're working over real x and the complex log function, then if x>0, it should clearly be inf. If x<0, log(x) is not real, but abs(log(x)) is, and sure enough limit(abs(log(x)),x,0,minus) gives inf. This is I believe correct. And it's all consistent with limit(abs(log(x)),x,0) => inf. If on the other hand you take the position that limit operates on real functions, then negative x are not part of the domain of log(x), so we should look only at the limit for x>0, which also gives us the result inf. But limit is actually happy to return imaginary results, e.g. limit(sqrt(x),x,1) => %i. By the way, currently, limit(log(x),x,0,minus) gives minf+%i*%pi, which makes some sort of intuitive sense, but isn't really a valid expression  it should be infinity.  Comment By: Raymond Toy (rtoy) Date: 20061108 21:49 Message: Logged In: YES user_id=28849 The error no longer occurs in current CVS. It returns UND, after asking if x is positive or negative. Why should the answer be INF? log(x) is undefined for negative real x.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=635606&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 17:50:15

Bugs item #965926, was opened at 20040603 13:21 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=965926&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  Trigonometry Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: trigsimp exponentially slow on lists Initial Comment: trigsimp of a [] list can sometimes take time exponential in the length of the list. For example: trigsimp(makelist(sin(i)^2+cos(i)^2,i,1,N)) for N=4,5,... takes 0.12, 0.52, 1.50, 4.45, 13.97 secs. Also for sin(i)^2. Of course, it should take linear time. This doesn't happen  >Comment By: Raymond Toy (rtoy) Date: 20061109 12:50 Message: Logged In: YES user_id=28849 This is probably due to the way trigsimp1 and improve work (share/trigonometry/trgsmp.mac). It looks like caused by trigsimp1 and improve, which cause quadratic behavior, I think. I think if trigsimp3 is modified to map(trigsimp1, num(expn))/map(trigsimp1,denom(expn), things will work much faster. Some care must be taken in case expn is not a list, but that's not too difficult.  Comment By: Stavros Macrakis (macrakis) Date: 20040603 13:22 Message: Logged In: YES user_id=588346 Oops, "this doesn't happen" should continue... in other cases, like sin(1)^2, sin(x)^2+cos(x)^2, etc.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=965926&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 17:32:14

Bugs item #1370433, was opened at 20051130 17:38 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1370433&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  Trigonometry Group: Fix for 5.9.2 Status: Open Resolution: None Priority: 5 Private: No Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: trigsimp(sqrt(%i2)) != sqrt(trigsimp(%i2)) Initial Comment:  Maxima version: 5.9.2 Maxima build date: 9:5 10/12/2005 host type: i686pcmingw32 lispimplementationtype: GNU Common Lisp (GCL) lispimplementationversion: GCL 2.6.7  ################################################## # Start problem with sqrt and trigsimp: # # # # (%i2) 2*(cos(x)^2sin(x)^2)+2; # # 2 2 # # (%o2) 2 (cos (x)  sin (x)) + 2 # # (%i3) trigsimp(sqrt(%i2)); # # (%o3)  2 abs(cos(x)) # # (%i4) sqrt(trigsimp(%i2)); # # (%o4) 2 abs(cos(x)) # # # # End problem with sqrt and trigsimp. # ##################################################  >Comment By: Raymond Toy (rtoy) Date: 20061109 12:32 Message: Logged In: YES user_id=28849 This bug is caused by the call to radcan in trigsimp (share/trigonometry/trgsmp.mac). For the case trigsimp(sqrt(%i2)), radcan converts sqrt(2*(cos(x)^2sin(x)^2)+2) int sqrt(2)*%i*sqrt(sin(x)^2cos(x)^21). I don't think this is what we want. Perhaps replacing radcan with ratsimp would be better.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1370433&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 15:34:46

Bugs item #1448605, was opened at 20060312 21:26 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1448605&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  Trigonometry Group: Fix for 5.9.2 Status: Open Resolution: None Priority: 5 Private: No Submitted By: Jeffrey Pikul (jpikul) Assigned to: Nobody/Anonymous (nobody) Summary: atan returns illegal value Initial Comment: (%i1) atan(tan(4)); (%o1) 4 (should be 4  %pi, or 0.858407346 as a float) (%i2) declare(z,complex); (%o2) done (%i3) atan(tan(z)); (%o3) z (see below) atan(tan(z)) ==> z is only true if %pi/2<z<%pi/2, so this should return atan(tan(z)) unless qualified with: assume(z<%pi/2,z>%pi/2); Maxima version: 5.9.2 Maxima build date: 9:5 10/12/2005 host type: i686pcmingw32 lispimplementationtype: GNU Common Lisp (GCL) lispimplementationversion: GCL 2.6.7  >Comment By: Raymond Toy (rtoy) Date: 20061109 10:34 Message: Logged In: YES user_id=28849 This simplification is controlled by the variable triginverses. It defaults to all, which is documented to convert atan(tan(x)) to x. I think this is not a bug.  Comment By: Nobody/Anonymous (nobody) Date: 20060313 10:38 Message: Logged In: NO Looks like SIMP%ASIN etc in src/trigo.lisp make the same type of simplification  afoo(foo(x)) => x for foo in {sin, cos, ...}.  Comment By: Nobody/Anonymous (nobody) Date: 20060313 00:46 Message: Logged In: NO Source of this bug seems to be SIMP%ATAN in src/trigi.lisp, in particular this line: (if (eq (caar y) '%tan) (cadr y)) where y is the argument of atan. It seems likely that the other simplification functions in the same file might suffer from similar naive assertions about function inverses. Robert Dodier  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1448605&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 14:00:30

Bugs item #626721, was opened at 20021022 03:15 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626721&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  Trigonometry Group: None >Status: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: logarc of atan2 wrong Initial Comment: res : logarc(atan2(y,x))$ rectform(res),y=1,x=1; => %pi/4 BUT atan2(1,1) => 3*%pi/4 The fix is to change the formula in $logarc and in simpatan2. Currently, logarc(atan2(y,x)) => logarc(atan (y/x)), which gives incorrect results as above. This formula should be replaced by %i*log((y+%i*x)/sqrt(x^2+y^2))  >Comment By: Raymond Toy (rtoy) Date: 20061109 09:00 Message: Logged In: YES user_id=28849 Fixed as suggested, with slight correction.  Comment By: Raymond Toy (rtoy) Date: 20061108 23:26 Message: Logged In: YES user_id=28849 Oops. I think the formula was intended to be %i*log((x+%i*y)/sqrt(x^2+y^2)). That produces the correct values. At least factor(rectform(<formula>)) produces atan2(y/r,x/r) with r = sqrt(x^2+y^2). It would be nice if atan2 simplified out the denominator to leave only atan2(y,x).  Comment By: Raymond Toy (rtoy) Date: 20061108 23:13 Message: Logged In: YES user_id=28849 Isn't this formula also wrong? For y=1, x=1, the formula becomes %i*log((1%i)/sqrt(2)). But log is the principal log so we still get %pi/4. I think it would be better to make atan2 be in the range (%pi,%pi].  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626721&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 04:26:30

Bugs item #626721, was opened at 20021022 03:15 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626721&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  Trigonometry Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: logarc of atan2 wrong Initial Comment: res : logarc(atan2(y,x))$ rectform(res),y=1,x=1; => %pi/4 BUT atan2(1,1) => 3*%pi/4 The fix is to change the formula in $logarc and in simpatan2. Currently, logarc(atan2(y,x)) => logarc(atan (y/x)), which gives incorrect results as above. This formula should be replaced by %i*log((y+%i*x)/sqrt(x^2+y^2))  >Comment By: Raymond Toy (rtoy) Date: 20061108 23:26 Message: Logged In: YES user_id=28849 Oops. I think the formula was intended to be %i*log((x+%i*y)/sqrt(x^2+y^2)). That produces the correct values. At least factor(rectform(<formula>)) produces atan2(y/r,x/r) with r = sqrt(x^2+y^2). It would be nice if atan2 simplified out the denominator to leave only atan2(y,x).  Comment By: Raymond Toy (rtoy) Date: 20061108 23:13 Message: Logged In: YES user_id=28849 Isn't this formula also wrong? For y=1, x=1, the formula becomes %i*log((1%i)/sqrt(2)). But log is the principal log so we still get %pi/4. I think it would be better to make atan2 be in the range (%pi,%pi].  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626721&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 04:13:04

Bugs item #626721, was opened at 20021022 03:15 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626721&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  Trigonometry Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: logarc of atan2 wrong Initial Comment: res : logarc(atan2(y,x))$ rectform(res),y=1,x=1; => %pi/4 BUT atan2(1,1) => 3*%pi/4 The fix is to change the formula in $logarc and in simpatan2. Currently, logarc(atan2(y,x)) => logarc(atan (y/x)), which gives incorrect results as above. This formula should be replaced by %i*log((y+%i*x)/sqrt(x^2+y^2))  >Comment By: Raymond Toy (rtoy) Date: 20061108 23:13 Message: Logged In: YES user_id=28849 Isn't this formula also wrong? For y=1, x=1, the formula becomes %i*log((1%i)/sqrt(2)). But log is the principal log so we still get %pi/4. I think it would be better to make atan2 be in the range (%pi,%pi].  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626721&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 03:14:26

Bugs item #593344, was opened at 20020809 22:14 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593344&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: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(abs(infinity)) strange result (was: infinite loop) Initial Comment: limit(abs(infinity)) appears to get into an infinite loop.  >Comment By: Raymond Toy (rtoy) Date: 20061108 22:14 Message: Logged In: YES user_id=28849 Current CVS returns infinity. That seems right. Closing this bug.  Comment By: Robert Dodier (robert_dodier) Date: 20050802 01:48 Message: Logged In: YES user_id=501686 Well, I don't see an infinite loop, but I do see this: (%i3) limit(abs(infinity)); (%o3) 'limit(abs(?foo),?foo,infinity) Maybe this is actually OK ?? Not sure. Maxima version: 5.9.1.1cvs Maxima build date: 10:5 7/28/2005 host type: i686redhatlinuxgnu lispimplementationtype: CLISP lispimplementationversion: 2.33.2 (20040602)  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593344&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 03:13:10

Bugs item #535363, was opened at 20020326 14:21 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=535363&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: Daniel Lemire (lemire) Assigned to: Nobody/Anonymous (nobody) Summary: Hospital gives up without warning Initial Comment: Behavior: (C1) limit(exp(1/x)/x^5,x,0,PLUS); (D1) 0 (C2) limit(exp(1/x)/x^6,x,0,PLUS);  1/x %E (D2) limit  x > 0+ 6 x Explanation (Richard Fateman): lhospitallim is set to 5. Yes it is arbitrary, but you can change it. Expected behavior: Wouldn't it be nicer if there was some warning message telling us "I stopped, but by changing this variable you could maybe get me to evaluate the limit in full"...? This gives the wrong impression... like "oh! Maxima can't solve that!". Further comments (Richard Fateman): I suppose this could be done, but it might also be the case that after giving up on L'Hopital's rule it tries something else that might succeed. The source code for limit is available. Further comments (James Amundson): I agree. It also might be a good idea to increase the default value. Computers are bigger and stronger now.  >Comment By: Raymond Toy (rtoy) Date: 20061108 22:13 Message: Logged In: YES user_id=28849 I think we should just increase the default to 8 (just as arbitrary as 5) and close this bug. Alternatively, we could make L'Hopital check to see if the numerator or denominator is a polynomial and change the limit appropriately.  Comment By: Robert Dodier (robert_dodier) Date: 20060326 18:21 Message: Logged In: YES user_id=501686 For the record, same behavior in Maxima 5.9.3.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=535363&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 03:05:21

Bugs item #593357, was opened at 20020809 23:36 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593357&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: Various problems with Limit and Atan2 Initial Comment: limit(atan2(sin(x),1/x),x,0) returns  FALSE Note that this only happens when the argument expression is negated (!).  limit(atan2(x,2*x),x,0) and limit(atan2(2*x^2,x^2),x,0) give the error Atan2(0,0) has been generated. But the first should give IND, and the second should give atan(2).  limit(atan2(x*abs(a),x),inf) correctly gives atan(abs(a)), but limit(atan2(x,abs(a)*x),x,inf) gives the noun form   >Comment By: Raymond Toy (rtoy) Date: 20061108 22:05 Message: Logged In: YES user_id=28849 Current CVS (including the fix for bug 626697 does this: limit(atan2(sin(x),1/x),x,0) > nounform limit(atan2(x*abs(a),x),inf) > atan(a) limit(atan2(x,abs(a)*x),x,inf) > atan(1/a) limit(atan2(x,2*x),x,0) and limit(atan2(2*x^2,x^2),x,0) both return atan2(0,0) generated.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593357&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 02:59:59

Bugs item #626697, was opened at 20021022 01:27 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626697&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: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(atan2(y,x),y,minf) => FALSE Initial Comment: limit(atan2(y,x),y,minf) => FALSE The fix is in the very last clause of SIMPLIMIT. Currently, it is (if $limsubst <stuff>) It should be (if $limsubst <stuff> (nounlimit exp var val))  >Comment By: Raymond Toy (rtoy) Date: 20061108 21:59 Message: Logged In: YES user_id=28849 Fixed as suggested. Closing.  Comment By: Stavros Macrakis (macrakis) Date: 20031009 15:38 Message: Logged In: YES user_id=588346 Same problem, same solution for limit(BETA((a+1)/b,(ba 1)/b)/b,a,b1);  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=626697&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 02:49:50

Bugs item #635606, was opened at 20021108 12:47 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=635606&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(abs(log(x))) internal error, UND Initial Comment: Maxima 5.5/Windows/gcl limit(abs(log(x)),x,0) Error: The tag LIMIT is undefined. Should of course be INF. More controversially, perhaps, limit(log(x),x,0) gives UND  I believe it should give INFINITY; after all, limit(log (x),x,0,minus) gives INFINITY.  >Comment By: Raymond Toy (rtoy) Date: 20061108 21:49 Message: Logged In: YES user_id=28849 The error no longer occurs in current CVS. It returns UND, after asking if x is positive or negative. Why should the answer be INF? log(x) is undefined for negative real x.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=635606&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 02:19:05

Bugs item #1036901, was opened at 20040929 05:53 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1036901&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: Jaroslaw Piskorski (jaropis) Assigned to: Nobody/Anonymous (nobody) Summary: tlimit(7^(n^2)/8^n,n,inf); wrong result Initial Comment: Maxima returns minf, while it should return inf here.  >Comment By: Raymond Toy (rtoy) Date: 20061108 21:19 Message: Logged In: YES user_id=28849 I think the bug is in taylor, which is called by taylim, via limit. taylor(7^(n^2)/8^n,n,inf,1) > log(8)*n+1 I don't think that that is right.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1036901&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 01:59:10

Bugs item #1100129, was opened at 20050111 08:05 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1100129&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: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: limit finds not a limit from defint, but finds it from topl Initial Comment: I found a problem Maxima hangs at, could someone take a look? I tried to solve a definite integration problem, and then I tried to track down the problem in maxima internals. I found that $limit fails on an expression which is easily solved from toplevel. Any ideas on what's happening? One more question: unless I explicilty do (trace sratsimp behaviorbydiff), :bt reporst a stop in $limit. It take quite a time to type in most of limit.lisp function names to get to the actual timeeater. Please, show me a technique to dig an actual call sequence?  Andrei Zorine Maxima 5.9.1 http://maxima.sourceforge.net Using Lisp Kyoto Common Lisp GCL 2.6.5 (aka GCL) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. This is a development version of Maxima. The function bug_report() provides bug reporting information. (%i1) (sqrt(sqrt(x^2+1)+x)sqrt(sqrt(x^2+1)x))/x/(x^2+1); 2 2 SQRT(SQRT(x + 1) + x)  SQRT(SQRT(x + 1)  x) (%o1)  2 x (x + 1) (%i2) limit(x*%,x,0); (%o2) 0 (%i3) defint(%o1,x,0,inf); Maxima encountered a Lisp error: Console interrupt. Automatically continuing. To reenable the Lisp debugger set *debuggerhook* to nil. (%i4) :lisp(trace sratsimp) (SRATSIMP) (%i4) :lisp(trace behaviorbydiff) (BEHAVIORBYDIFF) (%i5) :lisp(setq *debuggerhook* nil) NIL (%i5) defint(%o1,x,0,inf); console interrupt. Fast links are on: do (usefastlinks nil) for debugging Broken at SYSTEM::CLCSTERMINALINTERRUPT. Type :H for Help. 1 (Continue) Continues execution. 2 (Abort) Return to top level. dbl:MAXIMA>>:bt #0 CLCSTERMINALINTERRUPT {loc0=t} [ihs=17] #1 TRACECALL {tempname=nil,args=nil,cond=nil,entrycond=((system::arglist (#))),entry=nil,exi...} [ihs=16] #2 SRATSIMP {(#0=(mplus . #1=(simp)) (#2=(mtimes . #1#) (#3=# 1 2) (#4=# #5=# 1) ...) (#2#...} [ihs=15] #3 TRACECALL {tempname=nil,args=(#0=(mplus . #1=(simp)) (#2=(mtimes . #1#) (#3=# 1 2) (#4=#...} [ihs=14] #4 BEHAVIORBYDIFF {(#0=(mplus . #1=(simp)) (#2=(mtimes . #1#) 1 (#3=# # 1) ...) (#2# (#3# # 1) ...} [ihs=13] #5 $LIMIT {loc0=((mtimes . #0=(simp)) (#1=(mexpt . #0#) (#2=# 1 #) 1) (#2# (# 1 #) (#1# ...} [ihs=12] #6 LIMITNOERR {loc0=((mtimes . #0=(simp)) (#1=(mexpt . #0#) (#2=# 1 #) 1) (#2# (# 1 #) (#1# ...} [ihs=11] #7 MACSYMATOPLEVEL {inputstream=:internal,batchflag=$inf,loc2=nil,loc3=stringchar,loc4="0",loc...} [ihs=10] #8 RUN {} [ihs=6] #9 TOPLEVEL {loc0=nil,loc1=nil,loc2=nil,loc3=(lambdablock run nil ...)} [ihs=5] #10 FUNCALL {loc0=#<compiledfunction system:toplevel>,loc1=nil,loc2=0,loc3=0,loc4=nil,loc5...} [ihs=4] NIL  >Comment By: Raymond Toy (rtoy) Date: 20061108 20:59 Message: Logged In: YES user_id=28849 I'm not familiar with GCL's backtrace mechanism. CMUCL and SBCL produce much easier to understand backtraces. In any case, you can see there's a call to $limit. trace(limit) shows that it's getting stuck finding the limit of (sqrt(sqrt(1+eps^2)+eps)sqrt(sqrt(1+eps^2)eps))/(1+eps^2) for eps > 0 from above. Don't know why.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1100129&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 01:35:22

Bugs item #1106912, was opened at 20050121 14:07 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1106912&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) >Summary: limit(x/sin(x)^2,x,inf) Initial Comment: limit(x/sin(x)^2,x,inf) => UND actually = inf  >Comment By: Raymond Toy (rtoy) Date: 20061108 20:35 Message: Logged In: YES user_id=28849 Fix summary. Problem still exists in current CVS.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1106912&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 01:33:53

Bugs item #1281736, was opened at 20050904 15:26 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1281736&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: Vadim V. Zhytnikov (vvzhy) Assigned to: Nobody/Anonymous (nobody) Summary: limit((x/log(x))*(x^(1/x)1),x,inf)  wrong result Initial Comment: (%i1) limit((x/log(x))*(x^(1/x)1),x,inf); (%o1) inf Correct result should be 1  >Comment By: Raymond Toy (rtoy) Date: 20061108 20:33 Message: Logged In: YES user_id=28849 Current CVS returns the nounform for the limit. Better than inf, but could be better. tlimit actually returns 1 in this case.  Comment By: Stavros Macrakis (macrakis) Date: 20050928 10:25 Message: Logged In: YES user_id=588346 This is a bug. Even worse, tlimit also doesn't work, and for no good reason: tlimit(...) => nounform but limit ( taylor ( (x/log(x))*(x^(1/x)1) , x, inf, 0 ) ) => 0  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1281736&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 01:27:37

Bugs item #1593083, was opened at 20061108 20:27 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=1593083&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: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: tlimit(t^2*exp(4*t/38*exp(t)),t,inf) gives error Initial Comment: tlimit(t^2*exp(4*t/38*exp(t)),t,inf) produces an error: Invalid call to varexpand  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1593083&group_id=4933 
From: SourceForge.net <noreply@so...>  20061109 01:17:24

Bugs item #1548643, was opened at 20060829 11:09 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1548643&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: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: limit with abs: tag LIMIT is undefined Initial Comment: ex:abs(sqrt(11/x)1) limit(ex,x,0); Maxima encountered a Lisp error: Error in PROGN [or a callee]: The tag LIMIT is undefined.  >Comment By: Raymond Toy (rtoy) Date: 20061108 20:17 Message: Logged In: YES user_id=28849 Suggested fix implemented some time ago, so we don't get an error anymore. Closing.  Comment By: Raymond Toy (rtoy) Date: 20060830 21:40 Message: Logged In: YES user_id=28849 In mabssubst, near the very beginning, there's the line (equal ($imagpart (limit d var val 'think)) 0) Replace the call to limit with (let ((v (limitcatch d var val))) (unless v (throw 'mabs '$und)) v) This causes the given test to return "und". Not sure if that's what we want to return, but certainly better than an error. The testsuite passes with this change, FWIW.  Comment By: Raymond Toy (rtoy) Date: 20060830 19:49 Message: Logged In: YES user_id=28849 Ok. This is caused by mabssubst calling limit, but forgetting that limit might throw to the tag LIMIT. Not exactly sure what the right solution would be. Perhaps in that case mabssubst should, itself, throw 'mabs 'retn?  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1548643&group_id=4933 