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}
(28) 
_{Mar}

_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 






1
(3) 
2

3

4
(7) 
5
(1) 
6
(3) 
7
(1) 
8
(3) 
9
(2) 
10
(1) 
11
(3) 
12
(4) 
13
(8) 
14
(7) 
15
(1) 
16

17

18
(3) 
19
(5) 
20

21
(5) 
22
(3) 
23

24
(2) 
25
(1) 
26
(2) 
27
(3) 
28
(2) 
29
(2) 
30

31
(2) 






From: SourceForge.net <noreply@so...>  20041011 13:56:59

Bugs item #1044339, was opened at 20041011 00:19 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1044339&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: diff w.r.t. subscript wrong Initial Comment: diff(bessel_j[i+1](x),i) => 'diff(bessel_j[i+1](x), i+1) diff(li[i+1](x),i) => 'diff(..., i+1) But diff(..., i+1) is nonsense. Indeed, if you ev(..., diff) the above results, you get the error "nonvariable 2nd argument to diff".  >Comment By: Barton Willis (willisbl) Date: 20041011 08:56 Message: Logged In: YES user_id=895922 Maybe this has been fixed? Seems to be okay with 5.9.1.1cvs. (%i1) display2d : false$ (%i2) diff(bessel_j[i+1](x),i); (%o2) 'DIFF(BESSEL_J[i+1](x),i,1) (%i3) build_info(); Maxima version: 5.9.1.1cvs Maxima build date: 9:34 10/4/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.5 Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1044339&group_id=4933 
From: SourceForge.net <noreply@so...>  20041011 05:19:11

Bugs item #1044339, was opened at 20041011 01:19 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=1044339&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: diff w.r.t. subscript wrong Initial Comment: diff(bessel_j[i+1](x),i) => 'diff(bessel_j[i+1](x), i+1) diff(li[i+1](x),i) => 'diff(..., i+1) But diff(..., i+1) is nonsense. Indeed, if you ev(..., diff) the above results, you get the error "nonvariable 2nd argument to diff".  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1044339&group_id=4933 
From: SourceForge.net <noreply@so...>  20041011 03:40:56

Bugs item #1044318, was opened at 20041010 23:40 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=1044318&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: defint(1/(sin(x)^2+1),x,0,3*%pi) wrong Initial Comment: INTEGRATE(1/(SIN(x)^2+1),x,0,3*%PI) => 0 Since the integrand is everywhere >= 1/2, the integral cannot be zero  in fact integrate(1/(sin(x)^2+1),x,q,q+k*%pi) = k*pi/sqrt(2) (real q, integral k) Presumably defint is using the indefinite integral atan(2*tan(x)/sqrt(2))/sqrt(2) inappropriately.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1044318&group_id=4933 
From: SourceForge.net <noreply@so...>  20041010 06:30:20

Bugs item #808280, was opened at 20030917 19:16 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=808280&group_id=4933 >Category: Documentation Group: None Status: Open Resolution: None Priority: 5 Submitted By: Robert Dodier (robert_dodier) Assigned to: Nobody/Anonymous (nobody) Summary: describe(psi); claims Maxima can't evaluate psi Initial Comment: describe(psi); claims that Maxima can't evaluate the psi (the logarithmic derivative of gamma). Indeed psi[n](x); returns a simple expression involving zeta(n+1) for integer arguments x, and just psi[n](x) for other arguments. But there is another implementation of psi in the bffac module, so load(bffac); bpsi(n, x, precision); does yield numerical results for general arguments. It seems like the right thing to do is to have psi call bfpsi for arguments it can't evaluate. Another slightly less right thing is to change the description of psi to say that bfpsi in bffac can handle general arguments. For what it's worth, Robert Dodier  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=808280&group_id=4933 
From: SourceForge.net <noreply@so...>  20041009 04:17:25

Bugs item #1042244, was opened at 20041007 09:47 Message generated for change (Settings changed) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042244&group_id=4933 Category: Tests Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) >Summary: gfactor(a^4+b^4) fatal error Initial Comment: The command : factorsum(a**4+b**4) works, while the command gfactorsum(a**4+b**4) does not work. I use the version 5.5 of maxima  Comment By: Stavros Macrakis (macrakis) Date: 20041007 21:16 Message: Logged In: YES user_id=588346 Error is in gfactor(a^4+b^4)  result should be (b^2%I*a^2)*(b^2+%I*a^2)  Comment By: Stavros Macrakis (macrakis) Date: 20041007 21:11 Message: Logged In: YES user_id=588346 Specifically: gfactorsum(a^4+b^4) gives a fatal error in PSHIFT  tested in 5.9.0/GCL2.5.0/Mingw32/W2k  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042244&group_id=4933 
From: SourceForge.net <noreply@so...>  20041009 00:17:40

Bugs item #857171, was opened at 20031209 14:38 Message generated for change (Comment added) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=857171&group_id=4933 Category: Documentation Group: None >Status: Closed Resolution: None Priority: 5 Submitted By: Robert Dodier (robert_dodier) Assigned to: Nobody/Anonymous (nobody) Summary: variable needs to be quoted for DESCRIBE to work Initial Comment: I'm running Maxima 5.9.0 cvs version of 20031128, with Clisp 2.31, on linux kernel 2.4.2. I haven't tried this on any other system. I find that DESCRIBE(FOO) doesn't work as intended if FOO is a variable  it has to be DESCRIBE('FOO). (If FOO is a function then DESCRIBE works as intended.) I don't know what's the best solution  (1) change DESCRIBE so that it doesn't evaluate its argument? (Is that possible? Desirable?) (2) Fix up the documentation so that references to FOO say "Do DESCRIBE('FOO)" ? There are many places in the texi documents which say "Do DESCRIBE(FOO)" or "See FOO". Here is a partial list of those: Contexts.texi: CONTEXTS Expressions.texi, Runtime.texi: FEATURES, features Expressions.texi: PRODHACK Expressions.texi: SUMHACK Function.texi, xrefs.texi: MACROS Input.texi: INFOLISTS, GRIND, %% Integration.texi: ROMBERGABS, ROMBERGMIN Matrices.texi: LMXCHAR, RMXCHAR, RATMX, LISTARITH, DETOUT, DOALLMXOPS DOMXEXPT, DOMXMXOPS, DOSCMXOPS, DOSCMXPLUS, SCALARMATRIX, and SPARSE near "Do DESCRIBE(option)" Operators.texi: MODULUS Runtime.texi: RATEXPAND Trigonometric.texi: TRIGEXPAND, TRIGSIGN Debugging.texi: SETCHECKBREAK, DEBUGMODE Equations.texi: IEQNPRINT, SOLVE_INCONSISTENT_ERROR Miscellaneous.texi: INF, MINF Series.texi: TAYLOR_LOGEXPAND Simplification.texi: MAXPOSEX, MAXNEGEX, CAUCHYSUM Special.texi: GAMMALIM  >Comment By: Robert Dodier (robert_dodier) Date: 20041008 18:17 Message: Logged In: YES user_id=501686 I'm closing this bug report since the change suggested (to make DESCRIBE an MSPEC) was implemented in rev 1.13 of src/macdes.lisp.  Comment By: Stavros Macrakis (macrakis) Date: 20031209 23:44 Message: Logged In: YES user_id=588346 Describe could certainly be made into an MSPEC (like an fexpr in Lisp). In general, I am against making things quote their arguments for convenience (because it mucks up programmatic use), but in this case I think it makes sense, since this really is a commandline, useroriented command. Even after N years of using Maxima, I still commit the error of writing describe(gcd), for example. It seems very unlikely that anyone would want to have a variable argument to describe. Well, I suppose someone imaginative or misguided might try for i in [first,last] do describe(i), but given that describe is usually interactive (asking for the particular case to describe), this is not a great idea. So let's make it into an mspec.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=857171&group_id=4933 
From: SourceForge.net <noreply@so...>  20041008 01:29:22

Bugs item #1042702, was opened at 20041007 21:29 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=1042702&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor TAN(a*ATAN(q/a)) stack overflow Initial Comment: taylor( TAN(a*ATAN(q/a)), a, 1, 1) > stack overflow No stack overflow with q=1  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042702&group_id=4933 
From: SourceForge.net <noreply@so...>  20041008 01:16:39

Bugs item #1042244, was opened at 20041007 09:47 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042244&group_id=4933 Category: Tests Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: factorsum and gfactorsum functions Initial Comment: The command : factorsum(a**4+b**4) works, while the command gfactorsum(a**4+b**4) does not work. I use the version 5.5 of maxima  >Comment By: Stavros Macrakis (macrakis) Date: 20041007 21:16 Message: Logged In: YES user_id=588346 Error is in gfactor(a^4+b^4)  result should be (b^2%I*a^2)*(b^2+%I*a^2)  Comment By: Stavros Macrakis (macrakis) Date: 20041007 21:11 Message: Logged In: YES user_id=588346 Specifically: gfactorsum(a^4+b^4) gives a fatal error in PSHIFT  tested in 5.9.0/GCL2.5.0/Mingw32/W2k  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042244&group_id=4933 
From: SourceForge.net <noreply@so...>  20041008 01:11:42

Bugs item #1042244, was opened at 20041007 09:47 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042244&group_id=4933 Category: Tests Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: factorsum and gfactorsum functions Initial Comment: The command : factorsum(a**4+b**4) works, while the command gfactorsum(a**4+b**4) does not work. I use the version 5.5 of maxima  >Comment By: Stavros Macrakis (macrakis) Date: 20041007 21:11 Message: Logged In: YES user_id=588346 Specifically: gfactorsum(a^4+b^4) gives a fatal error in PSHIFT  tested in 5.9.0/GCL2.5.0/Mingw32/W2k  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042244&group_id=4933 
From: SourceForge.net <noreply@so...>  20041007 13:47:26

Bugs item #1042244, was opened at 20041007 06:47 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=1042244&group_id=4933 Category: Tests Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: factorsum and gfactorsum functions Initial Comment: The command : factorsum(a**4+b**4) works, while the command gfactorsum(a**4+b**4) does not work. I use the version 5.5 of maxima  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042244&group_id=4933 
From: SourceForge.net <noreply@so...>  20041006 16:59:38

Bugs item #1041148, was opened at 20041005 21:46 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1041148&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 3 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(x^a,[x],0,1) unsimplified Initial Comment: taylor(x^a,[x],0,1) => 1^2*x^a+... Note unsimplified 1^2.  >Comment By: Barton Willis (willisbl) Date: 20041006 11:59 Message: Logged In: YES user_id=895922 Maybe there is a connection between this bug and bug 974734. Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1041148&group_id=4933 
From: SourceForge.net <noreply@so...>  20041006 16:05:40

Bugs item #1041570, was opened at 20041006 12:05 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=1041570&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: assume(abs(x)<1) should imply x<1 and x>1 Initial Comment: assume(abs(x)<1)$ sign(x1) => doesn't know (should be NEG) sign(1x) => doesn't know (should be POS) assume(x>1) => doesn't notice contradiction The easy fix is of course to make the original assumption add x<1 and x>1 to the database, but how do you ensure that they get *removed* if you remove (forget) the original assertion. And what if x<1 has been explicitly added in the meantime? Same problem in the reverse direction: assume(y<1,y>1)$ sign(abs(y)1) => PNZ (should be NEG)  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1041570&group_id=4933 
From: SourceForge.net <noreply@so...>  20041006 02:46:55

Bugs item #1041148, was opened at 20041005 22:46 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=1041148&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 3 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(x^a,[x],0,1) unsimplified Initial Comment: taylor(x^a,[x],0,1) => 1^2*x^a+... Note unsimplified 1^2.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1041148&group_id=4933 
From: SourceForge.net <noreply@so...>  20041005 18:58:48

Bugs item #956730, was opened at 20040519 09:04 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=956730&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^n,x,0,inf) => sinrx error Initial Comment: integrate(x^n,x,0,inf); Is n positive, negative, or zero? neg; Is n + 1 positive, negative, or zero? neg; Is n an integer? y; Error: $x is not of type LIST. Error signalled by SINRX. The integral is of course divergent, but it should not be getting an internal error. Why does it even ask whether n is integral? Strangely, with n noninteger, it correctly diverges.  >Comment By: Barton Willis (willisbl) Date: 20041005 13:58 Message: Logged In: YES user_id=895922 It seems that this bug has been fixed. (%i2) integrate(x^n,x,0,inf); Is n positive, negative, or zero? neg;Is n + 1 positive, negative, or zero? neg;Is n an integer? yes;Integral is divergent  an error. Quitting. To debug this try DEBUGMODE(TRUE); (%i3) build_info(); Maxima version: 5.9.1.1cvs Maxima build date: 9:34 10/4/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.5  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=956730&group_id=4933 
From: SourceForge.net <noreply@so...>  20041004 21:16:03

Bugs item #924868, was opened at 20040328 07:53 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=924868&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: defint log(sqrt(q^21)+1) asks about YX Initial Comment: T: log(sqrt(q^21)+1); integrate(T,q,0,1) asks the question: Is YX  1 positive, negative, or zero? Note that the variable YX does not appear in the integrand! By the way, if you answer Pos, you get the error SIGN called on an imaginary argument: %I and if you answer Neg, you get the error Quotient by a polynomial of higher degree (with default GCD) Maxima 5.9.0 gcl 2.5.0 W2k  >Comment By: Raymond Toy (rtoy) Date: 20041004 17:15 Message: Logged In: YES user_id=28849 This works much better in maxima 5.9.1, but it still asks about YX. It's trying to solve yx  sqrt(q^1+1)  1 for q. Perhaps for some kind of substitution.  Comment By: Stavros Macrakis (macrakis) Date: 20041001 12:28 Message: Logged In: YES user_id=588346 integrate(exp(x^(1/2)),x,0,1) also asks about YX, but gets the correct answer (2) regardless of what you answer. Maxima version: 5.9.0.9beta2 Maxima build date: 10:50 7/27/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.3  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=924868&group_id=4933 
From: SourceForge.net <noreply@so...>  20041004 21:03:23

Bugs item #620070, was opened at 20021008 00:30 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=620070&group_id=4933 Category: None Group: None >Status: Closed Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: rat>float loses 9 digits of precision! Initial Comment: Maxima 5.5 / gcl 23.6 / Windows 2000 / Athlon num: 1099511627775 $ den: 1099511627776 $ fnfd : float(num)/float(den) => 0.9999999999990905 (correct) but fnd: float(num/den) => 0.9999999403953552 That is, we're only getting 7 digits of precision when our floating point actually supports 16! The problem is in Fpcofrat1, with the value of machine mantissaprecision, which is substituted at readtime, and is defined as: (evalwhen (compile load) (defconstant machinemantissaprecision (floatdigits 1d0)) ) Apparently when Maxima was built, this evaluated to 24, not to 53 as it should have. Correcting machinemantissaprecision to 53 and redefining fpcofrat1 with the correct value gives the correct result. Why does machinemantissaprecision have the wrong value? Is Maxima crosscompiled for Windows? Anyway, while looking at the fpcofrat1 code, I realized that it could be simplified and optimized at the same time. Here is my code. Does someone want to look it over and check it in?  The following code takes advantage of modern floating point's extended range for the common case, and Maxima's bfloat routines to simplify the (very rare) case where a numerator or denominator has more than 300 digits. ;; Convert rational number NU/D to float (defun fpcofrat1 (nu d) (if (and (< (abs nu) mostpositivedoublefloat) (< (abs d) mostpositivedoublefloat)) ;; Machine arithmetic is excellent (/ (float nu) d) ;; ...except in the rare case where nu or d are huge (let ((fpprec #.(+ machinemantissaprecision 3))) (fp2flo ($bfloat (list '(rat simp) nu d))))))  >Comment By: Stavros Macrakis (macrakis) Date: 20041004 10:51 Message: Logged In: YES user_id=588346 Works correctly in current version: 5.9.0.9beta2 Maxima build date: 10:50 7/27/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.3  Comment By: Raymond Toy (rtoy) Date: 20021015 12:51 Message: Logged In: YES user_id=28849 This appears to be fixed in the CVS sources.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=620070&group_id=4933 
From: SourceForge.net <noreply@so...>  20041004 20:39:24

Bugs item #1037903, was opened at 20040930 13:17 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1037903&group_id=4933 Category: Xmaxima Group: Fix for 5.9.0 Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: Prolems with the output "stringout()" Initial Comment: When we use the output "stringout()" some problems with anteriority occurs. e.g: the answer: (x^2)/3 is showed x^2/3, and can mean x^(2/3) So, can we specify for "stringout()" to break the sentences of the equation? Thank you for your attention. We're waiting for an answer. EMat developers team, http://www.gritee.com info@...  >Comment By: Raymond Toy (rtoy) Date: 20041004 16:39 Message: Logged In: YES user_id=28849 Why is this a problem? x^2/3 is supposed to mean (x^2)/3. I think it is wrong to interpret x^2/3 as x^(2/3).  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1037903&group_id=4933 
From: SourceForge.net <noreply@so...>  20041004 20:35:41

Bugs item #1036900, was opened at 20040929 05:51 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1036900&group_id=4933 Category: None Group: None >Status: Closed >Resolution: Fixed Priority: 5 Submitted By: Jaroslaw Piskorski (jaropis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(7^n/8^n,n,inf); unevaluated Initial Comment: Maxima returns this limit unevaluated, while it should return 0. Another example is limit(4^n/2^(2*n),n,inf); Perhaps the same problem causes Maxima to return limit((10^n+9^n+8^n)^(1/n),n,inf); unevaluated (well, in a way), while it should return 10.  >Comment By: Raymond Toy (rtoy) Date: 20041004 16:35 Message: Logged In: YES user_id=28849 Closed. See Bug 1039965.  Comment By: Raymond Toy (rtoy) Date: 20041004 10:28 Message: Logged In: YES user_id=28849 Partially fixed in CVS. The 4^n/2^(2*n) case still returns infinity instead of 1. I'm going to close this bug and add a new one about this case.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1036900&group_id=4933 
From: SourceForge.net <noreply@so...>  20041004 15:02:34

Bugs item #1036900, was opened at 20040929 05:51 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1036900&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Jaroslaw Piskorski (jaropis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(7^n/8^n,n,inf); unevaluated Initial Comment: Maxima returns this limit unevaluated, while it should return 0. Another example is limit(4^n/2^(2*n),n,inf); Perhaps the same problem causes Maxima to return limit((10^n+9^n+8^n)^(1/n),n,inf); unevaluated (well, in a way), while it should return 10.  >Comment By: Raymond Toy (rtoy) Date: 20041004 10:28 Message: Logged In: YES user_id=28849 Partially fixed in CVS. The 4^n/2^(2*n) case still returns infinity instead of 1. I'm going to close this bug and add a new one about this case.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1036900&group_id=4933 
From: SourceForge.net <noreply@so...>  20041004 15:02:27

Bugs item #1038584, was opened at 20041001 12:39 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1038584&group_id=4933 Category: None Group: None >Status: Closed >Resolution: Fixed Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(exp(x^3),x,0,1) => internal error Initial Comment: integrate(exp(x^a),x,0,1),a=3; Maxima encountered a Lisp error: Error in MACSYMATOPLEVEL [or a callee]: ((RAT SIMP) 10 3) is not of type NUMBER. Same problem for a=4, a=3/2 (after questions about YX cf bug report # 924868) Maxima version: 5.9.0.9beta2 Maxima build date: 10:50 7/27/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.3  >Comment By: Raymond Toy (rtoy) Date: 20041004 10:30 Message: Logged In: YES user_id=28849 Fixed. Needed to use LESSTHAN instead of <.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1038584&group_id=4933 
From: SourceForge.net <noreply@so...>  20041004 15:02:24

Bugs item #1039965, was opened at 20041004 10:34 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=1039965&group_id=4933 Category: Lisp Core Group: None Status: Open Resolution: None Priority: 5 Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: limit(4^n/2^(2*n),n,inf) is wrong Initial Comment: The limit is returned as infinity instead of 1. The limit routines are confused by this and it the exponent for 4^n is larger than the 2^(2*n). More investigation needed. Compare bug #1036900.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1039965&group_id=4933 
From: SourceForge.net <noreply@so...>  20041001 17:58:39

Bugs item #1038624, was opened at 20041001 13: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=1038624&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: askinteger ignores asksign database Initial Comment: assume(equal(zzz,0))$ askinteger(zzz) asks the user Last time I looked, zero was an integer. You can elicit this also in cases like integrate(x^r),x,0,1), which ask Is r pnz? answer=>ZERO Is r an integer? It is also very confusing (though duly documented) that askinteger affects the *global* database, while asksign only affects the *current computation* database. Fortunately, this is not true for internal uses of askinteger like the integrate example above.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1038624&group_id=4933 
From: SourceForge.net <noreply@so...>  20041001 16:39:13

Bugs item #1038584, was opened at 20041001 12:39 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=1038584&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(exp(x^3),x,0,1) => internal error Initial Comment: integrate(exp(x^a),x,0,1),a=3; Maxima encountered a Lisp error: Error in MACSYMATOPLEVEL [or a callee]: ((RAT SIMP) 10 3) is not of type NUMBER. Same problem for a=4, a=3/2 (after questions about YX cf bug report # 924868) Maxima version: 5.9.0.9beta2 Maxima build date: 10:50 7/27/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.3  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1038584&group_id=4933 
From: SourceForge.net <noreply@so...>  20041001 16:28:20

Bugs item #924868, was opened at 20040328 07:53 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=924868&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: defint log(sqrt(q^21)+1) asks about YX Initial Comment: T: log(sqrt(q^21)+1); integrate(T,q,0,1) asks the question: Is YX  1 positive, negative, or zero? Note that the variable YX does not appear in the integrand! By the way, if you answer Pos, you get the error SIGN called on an imaginary argument: %I and if you answer Neg, you get the error Quotient by a polynomial of higher degree (with default GCD) Maxima 5.9.0 gcl 2.5.0 W2k  >Comment By: Stavros Macrakis (macrakis) Date: 20041001 12:28 Message: Logged In: YES user_id=588346 integrate(exp(x^(1/2)),x,0,1) also asks about YX, but gets the correct answer (2) regardless of what you answer. Maxima version: 5.9.0.9beta2 Maxima build date: 10:50 7/27/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.3  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=924868&group_id=4933 