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}
(19) 
_{May}
(14) 
_{Jun}
(10) 
_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 


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

9
(2) 
10
(2) 
11
(3) 
12
(2) 
13
(1) 
14
(2) 
15
(1) 
16

17
(1) 
18
(1) 
19
(2) 
20
(1) 
21
(2) 
22
(5) 
23
(2) 
24
(7) 
25
(6) 
26
(1) 
27
(3) 
28
(1) 
29

30
(1) 
31




From: SourceForge.net <noreply@so...>  20110804 21:12:12

Bugs item #3004558, was opened at 20100520 10:12 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3004558&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: Pending Resolution: None Priority: 5 Private: No Submitted By: l_butler () Assigned to: Nobody/Anonymous (nobody) Summary: ordfna error Initial Comment: I encounter this Lisp error in the course of computations. I've provided the backtrace. Typeerror in KERNEL::OBJECTNOTLISTERRORHANDLER: #:(4/a2^2+a0^2+1)^(1/2)28262 is not of type LIST [Condition of type TYPEERROR] Restarts: 0: [MACSYMAQUIT] Maxima toplevel 1: [ABORT ] Skip remaining initializations. Debug (type H for help) (ORDFNA (#:(4/a2^2+a0^2+1)^(1/2)28262 2 (#:min((a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^24*sqrt(4/a2^2+a0^2+1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2),(a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^2+4*sqrt(4/a2^2+a0^2\ +1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2))28304 2 1 1 2 ...) 1 (#:min((a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^24*sqrt(4/a2^2+a0^2+1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2),(a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^2+4*sqrt(4/a2^2+a0^2\ +1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2))28304 2 2 1 6 ...) ...) $A0) Source: (CAAR E) 0] backtrace 30 0: (ORDFNA (#:(4/a2^2+a0^2+1)^(1/2)28262 2 (#:min((a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^24*sqrt(4/a2^2+a0^2+1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2),(a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^2+4*sqrt(4/a2^2+a\ 0^2+1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2))28304 2 1 1 2 ...) 1 (#:min((a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^24*sqrt(4/a2^2+a0^2+1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2),(a2^2*a0^2+(sqrt(4/a2^2+a0^2+1)^2+1)*a2^2+4*sqrt(4/a2^2+a\ 0^2+1)*a2+4)/(sqrt(4/a2^2+a0^2+1)^2*a2^2))28304 2 2 1 6 ...) ...) $A0) 1: (GREAT $A0 ((MPLUS SIMP) 26 (# 27 #))) 2: (TIMESIN ($A0 2) (1 (# 26 #)) 2) 3: (TMS ((MEXPT SIMP RATSIMP) $A0 2) 1 ((MTIMES) 1 (# 26 #))) 4: (SIMPTIMES ((# $A0 2)) ((MEXPT SIMP RATSIMP) $A0 2) NIL) 5: (SIMPLUS ((# # #) (# # $A0) (# 17 #) 1) 1 NIL) 6: (SIMPTIMES ((# # # # 1) (# # 1)) 1 NIL) 7: ($TOTALDISREP ((MRAT SIMP # #) (#:(4/a2^2+a0^2+1)^(1/2)28262 2 # 1 # ...) #:(4/a2^2+a0^2+1)^(1/2)28262 2 27 ...)) 8: (FULLRATSIMP ((MRAT SIMP # #) (#:(4/a2^2+a0^2+1)^(1/2)28262 2 # 1 # ...) #:(4/a2^2+a0^2+1)^(1/2)28262 2 27 ...)) 9: (MEVAL1 (($RATSIMP) $S)) 10: (MEVAL (($RATSIMP) $S)) 11: ("DEFMSPEC MSETQ" ((MSETQ) $SRAT (# $S))) 12: (MEVAL1 ((MSETQ) $SRAT (# $S))) 13: (MEVAL ((MSETQ) $SRAT (# $S))) 14: ("DEFMSPEC MPROG" ((# $A0) (# #) (# $S #) (# $SRAT #) (# $SRAT) ...)) 15: (MEVAL1 ((MPROG) (# $L $S $SRAT $A0 ...) (# $A0) (# #) (# $S #) ...)) 16: (MEVAL ((MPROG) (# $L $S $SRAT $A0 ...) (# $A0) (# #) (# $S #) ...)) 17: (MEVAL* ((MPROG) (# $L $S $SRAT $A0 ...) (# $A0) (# #) (# $S #) ...)) 18: (CONTINUE #<TwoWay Stream, Input = #<Synonym Stream to SYSTEM:*STDIN*>, Output = #<Synonym Stream to SYSTEM:*STDOUT*>> NIL) 19: (MACSYMATOPLEVEL #<TwoWay Stream, Input = #<Synonym Stream to SYSTEM:*STDIN*>, Output = #<Synonym Stream to SYSTEM:*STDOUT*>> NIL)  >Comment By: Dieter Kaiser (crategus) Date: 20110804 23:12 Message: I think we do not have enough information to work on the problem. Setting the status to pending and the resolution to "none". Dieter Kaiser  Comment By: Dieter Kaiser (crategus) Date: 20100528 19:40 Message: Sorry, perhaps I oversee something. But I do not see the initial problem which triggers the error. Please can you post it. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3004558&group_id=4933 
From: SourceForge.net <noreply@so...>  20110804 21:05:41

Bugs item #2974948, was opened at 20100323 00:52 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2974948&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  Floating point Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: rectform loses some accuracy Initial Comment: Consider: rectform(1e300/(1e300+3/2*%i)) > 1.0 But float(rectform(10^300/(10^300+3/2*%i))) > 1.0  1.5e300 * %i The inaccuracy is caused by Maxima computing the rectform in two parts. Maxima computes 1/(1e300+3/2*%i) and then multiplies by 1e300. Maxima correctly computes 1/(1e300+3/2*%i) as 1e300. But doing it in two steps causes the very tiny imaginary part to be lost.  >Comment By: Dieter Kaiser (crategus) Date: 20110804 23:05 Message: For the record: With SBCL and Maxima 5.25post I get for the example an overflow error: (%i1) rectform(1e300/(1e300+3/2*%i)); Maxima encountered a Lisp error: arithmetic error FLOATINGPOINTOVERFLOW signalled Automatically continuing. To reenable the Lisp debugger set *debuggerhook* to nil. I have tried the biggest possible exponent: (%i31) float(rectform(1e154/(1e154+3/2*%i))); (%o31) 1.0  1.5e154 %i (%i32) float(rectform(10^154/(10^154+3/2*%i))); (%o32) 1.0  1.5e154 %i Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2974948&group_id=4933 
From: SourceForge.net <noreply@so...>  20110804 20:50:52

Bugs item #2998214, was opened at 20100507 17:18 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2998214&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: Pending >Resolution: Works For Me Priority: 5 Private: No Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: Cell can't end with a comment Initial Comment: When the last line in a cell is comment, maxima comlains about premature termination of input. To reproduce: Just enter a comment only in a input cell. /*Comment*/ and evaluate. Why do you need that? If you want to do somehting like this: a:99 /*assign a the magic value*/  >Comment By: Dieter Kaiser (crategus) Date: 20110804 22:50 Message: As written in one of the posts, it is save to use the syntax a:99 /* Comment */ ; in wxMaxima. The following syntax a:99; /* Comment */ might give an error in wxMaxima, because wxMaxima automatically appends an ;char and a comment followed by a ;char is not a valid command. Furthermore, I do not see a problem with a comment at the end of a file with SBCL. I have tried the Maxima versions 5.9, 5.10, 5.11, 5.12, 5.13, 5.14, 5.15, 5.16, 5.17, 5.18, 5.19, 5.20, and 5.25post and had no problems loading the test file of the last posting. I have tried several variations of the file and had a look at the code to search for a potential problem. Perhaps, we can close this bug report and open a new bug report, if we have a report with a known Lisp and Maxima version. Setting the status to pending and the resolution to "works for me". Dieter Kaiser  Comment By: Robert Dodier (robert_dodier) Date: 20100603 17:54 Message: This is actually a bug in the Maxima parser if I am not mistaken; it is not specific to wxMaxima so far as I know. At one time (maybe it is still true) the parser would barf on a batch script if it ended with a comment. I just tried that with Maxima 5.20.0 + GCL + Windows but the script was loaded successfully. Maybe the bug has been fixed? Need to test on other platforms and Lisp implementations. To the original poster: what does build_info(); report? Here is my test script: x:1; /* comment */  Comment By: Nobody/Anonymous (nobody) Date: 20100527 08:55 Message: Got it. A comment can't be the only text in a field. /*comment*/a:99 and a:99/*comment*/ both work ok. wxMaxima appears to append a semicolon to the line if it is absent. /*comment*/ gets a semicolon appended, then an error message "incorrect syntax: Premature termination of input at ;. /*comment*/; ^" With the caret indicating the closing slash of the comment. Using a $ terminator gets a similar message, no semicolon is appended, and the caret indicates the closing slash as before.  Comment By: Peter Cusack (pogo1) Date: 20100527 03:53 Message: I note that wmMaxima version 0.8.5 helpfully adds a semicolon and that fixes the problem sometimes. I'm still trying to make it reproducibly fail. If the line contains a $ terminator before the comment ( like this B0:i*Mu0/(2*a)/*Field at centre of single turn*/$) then adding the semicolon at the end makes Maxima hang.  Comment By: l_butler () Date: 20100520 10:16 Message: Hi, I think the problem is that Maxima input must be terminated with a semicolon: a:99; /* comment */ a:99 /* comment */; Both are valid inputs that assign 99 to a.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2998214&group_id=4933 
From: SourceForge.net <noreply@so...>  20110804 19:14:31

Bugs item #2995695, was opened at 20100503 04:37 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2995695&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: Documentation Group: None >Status: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: Documentation for resultant is confusing Initial Comment: The documentation for resultant starts:  Variable: resultant Computes the resultant of the two polynomials <p_1> and <p_2>, eliminating the variable <x>. The resultant is a determinant of the coefficients of <x> in <p_1> and <p_2>, which equals zero if and only if <p_1> and <p_2> have a nonconstant factor in common. The text talks about p_1 and p_2, but header says it's a variable. It's really talking about the function resultant. The variable resultant is discucssed later. This needs to be separated.  >Comment By: Dieter Kaiser (crategus) Date: 20110804 21:14 Message: Some work on the documentation has been don in Polynomials.texi revision 04.08.2011. Closing this bug report as fixed. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2995695&group_id=4933 
From: SourceForge.net <noreply@so...>  20110804 19:12:46

Bugs item #3386347, was opened at 20110804 21:12 Message generated for change (Tracker Item Submitted) made by stefano_ferri You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3386347&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: Stefano Ferri (stefano_ferri) Assigned to: Nobody/Anonymous (nobody) Summary: Rank leads to a "quotient by zero" error Initial Comment: I have a matrix, called K and listed below, and I need to compute its rank. Here a strange thing happens: (%i3) rank(K); `quotient' by `zero'  an error. To debug this try: debugmode(true); while (%i4) rank(ratsimp(K)); (%o4) 12 (%i5) as it should be, 12 is the correct answer. I would always expect a result from rank, not an error, therefore I think this is a bug in some matrixrelated subroutine, since the same error is present when I try to use K to solve a linear system (that do has a solution) with linsolve, on a system of equations generated starting from the same matrix. This error was generated with Maxima 5.23.2, compiled against gcl, on Slackware 13.1 (I reported this issue on the mailing list some time ago), but it is still present on Maxima 5.24. I'm sorry, the copy & paste for the following matrix (if somebody kindly wants to try) only works in wxMaxima, and not in a Maxima shell, due to the text format. For the note, in my code I assume l>0, but the error is always present Here is the matrix: K : matrix([1,0,0,0,0,0,0,0,0,0,0,0],[0,1,0,0,0,0,0,0,0,0,0,0], [0,0,1,0,0,0,0,0,0,0,0,0], [0,0,0, 6*l^2*E*I/((sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) +6*E*I/l^3 +A*E/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2))+A*E/(2*l), 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) 6*E*I/l^3 l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2))+A*E/(2*l), 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) +3*sqrt(2)*E*I/l^2, 6*l^2*E*I/((sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) A*E/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)), l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)), 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)),0,0,0], [0,0,0, 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) 6*E*I/l^3 l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2))+A*E/(2*l), 12*E*I/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) +6*E*I/l^3 +l^2*A*E/(2*(sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2))+A*E/(2*l), 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*sqrt(2)*E*I/l^2, l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)), 12*E*I/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) l^2*A*E/(2*(sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)),0,0,0], [0,0,0, 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) +3*sqrt(2)*E*I/l^2, 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*sqrt(2)*E*I/l^2, 4*E*I/sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)+4*E*I/l, 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 2*E*I/sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2),0,0,0], [0,0,0, 6*l^2*E*I/((sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) A*E/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)), l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)), 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 6*l^2*E*I/((sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) +6*l^2*E*I/((l/sqrt(2)sqrt(2)*l)^2 *(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)^(3/2)) +A*E/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) +A*E/((l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *sqrt((l/sqrt(2)sqrt(2)*l)^2+l^2/2)), 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) +3*2^(3/2)*l*E*I/((l/sqrt(2)sqrt(2)*l) *(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)^(3/2)) l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) l*A*E/(sqrt(2)*(l/sqrt(2)sqrt(2)*l) *(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *sqrt((l/sqrt(2)sqrt(2)*l)^2+l^2/2)), 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*sqrt(2)*l*E*I/((l/sqrt(2)sqrt(2)*l) *sqrt(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)),0,0,0], [0,0,0, l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)), 12*E*I/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) l^2*A*E/(2*(sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 3*2^(3/2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) +3*2^(3/2)*l*E*I/((l/sqrt(2)sqrt(2)*l) *(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)^(3/2)) l*A*E/(sqrt(2)*(sqrt(2)*ll/sqrt(2)) *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) l*A*E/(sqrt(2)*(l/sqrt(2)sqrt(2)*l) *(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *sqrt((l/sqrt(2)sqrt(2)*l)^2+l^2/2)), 12*E*I/((l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)^(3/2)) +12*E*I/((l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)^(3/2)) +l^2*A*E/(2*(sqrt(2)*ll/sqrt(2))^2 *(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2)) +l^2*A*E/(2*(l/sqrt(2)sqrt(2)*l)^2 *(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *sqrt((l/sqrt(2)sqrt(2)*l)^2+l^2/2)), 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 6*E*I/(sqrt(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)),0,0,0], [0,0,0, 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)), 2*E*I/sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2), 3*sqrt(2)*l*E*I/((sqrt(2)*ll/sqrt(2)) *sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 3*sqrt(2)*l*E*I/((l/sqrt(2)sqrt(2)*l) *sqrt(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)), 6*E*I/(sqrt(l^2/(2*(sqrt(2)*ll/sqrt(2))^2)+1) *((sqrt(2)*ll/sqrt(2))^2+l^2/2)) 6*E*I/(sqrt(l^2/(2*(l/sqrt(2)sqrt(2)*l)^2)+1) *((l/sqrt(2)sqrt(2)*l)^2+l^2/2)), 4*E*I/sqrt((sqrt(2)*ll/sqrt(2))^2+l^2/2) +4*E*I/sqrt((l/sqrt(2)sqrt(2)*l)^2+l^2/2),0,0,0], [0,0,0,0,0,0,0,0,0,1,0,0],[0,0,0,0,0,0,0,0,0,0,1,0], [0,0,0,0,0,0,0,0,0,0,0,1]);  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3386347&group_id=4933 
From: SourceForge.net <noreply@so...>  20110804 16:34:07

Bugs item #2983881, was opened at 20100408 15:25 Message generated for change (Settings changed) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2983881&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: Pending >Resolution: Works For Me Priority: 5 Private: No Submitted By: Stefano Ferri (stefano_ferri) Assigned to: Nobody/Anonymous (nobody) Summary: Triangularize does not preserve rank Initial Comment: I already reported this issue on the mailing list. Here there is an example where triangularize fails with a symbolic matrix. In a nutshell, doing triangularize over a symbolic matrix and then computing rank when the triangularized matrix is evaluated giving a value to a parameter, yields differents results than doing evaluation followed by triangularization. (%i1) display2d:false; (%o1) false (%i2) M:matrix([h1,h,h+2,h+1,1],[1h,2,h2,1,1],[h1,h,1,0,1],[0,h2,0,2,1]); (%o2) matrix([h1,h,h+2,h+1,1],[1h,2,h2,1,1],[h1,h,1,0,1],[0,h2,0,2,1]) (%i3) ev(M,h=1); (%o3) matrix([2,1,1,0,1],[2,2,1,1,1],[2,1,1,0,1],[0,1,0,2,1]) (%i4) rank(%); (%o4) 3 (%i5) triangularize(%o3); (%o5) matrix([2,1,1,0,1],[0,2,0,2,4],[0,0,0,2,6],[0,0,0,0,0]) (%i6) rank(%); (%o6) 3 that is correct, but: (%i7) Mt:triangularize(M); (%o7) matrix([h1,h,1,0,1],[0,h^2h+2,0,22*h,h1], [0,0,h^32*h^2+h+2,h^32*h^2+h+2,0], [0,0,0,h^42*h^3+h^2+2*h,3*h^36*h^2+3*h+6]) (%i8) ev(Mt,h=1); (%o8) matrix([2,1,1,0,1],[0,2,0,4,2],[0,0,0,0,0],[0,0,0,0,0]) (%i9) rank(%); (%o9) 2 %o9 is a wrong result, if h=1, rank of matrix M is 3  >Comment By: Dieter Kaiser (crategus) Date: 20110804 18:34 Message: I think the problem is not a bug, but the algorithm of the function triangularize, that is the Gaussian elimination, is not suited for the type of problem. The problem is that the matrix elements (h1), (h+1), or (h+2) are zero for the values h=1, h=1, or h=2. To get the triangular matrix the algorithm multiplies the rows with the factors (h1), (h+1), and (h+2), but e. g. for h=1 the factor (h+1) is zero. Therefore we are multiplying a row with a zero, which causes a wrong result when we insert h=1. The same happens when inserting h=1 or h=2. The following shows the factored result from the function triangularize: (%i1) m:matrix([h1,h,h+2,h+1,1],[1h,2,h2,1,1],[h1,h,1,0,1],[0,h2,0,2,1])$ (%i2) factor(triangularize(m)); (%o2) matrix([h1,h,1,0,1], [0,(h1)*(h+2),0,2*(h1),h1], [0,0,(h1)*(h+1)*(h+2),(h1)*(h+1)*(h+2),0], [0,0,0,(h1)*h*(h+1)*(h+2),3*(h1)*(h+1)*(h+2)]) The elements of the last two rows contain the factors (h1), (h+1), and (h+2) and will vanish when inserting one of the values h=1, h=1, or h=2. Therefore, I would argue that the triangular form of the matrix must not be correct for the special values h=1, h=1, and h=2. The algorithm of e.g. ptriangularize has not the problem for this example: (%i3) ptriangularize(m,h); (%o3) matrix([h1,h,1,0,1],[0,h2,0,2,1],[0,0,h1,1,3],[0,0,0,h,3]) At this point, I would suggest to close this bug report as "works for me". Dieter Kaiser  Comment By: Stefano Ferri (stefano_ferri) Date: 20100414 00:49 Message: Maybe I've found something. Triangularize dos not check if the matrix is already reduced, insisting to multiply for the pivot. Here is an example. (%i3) m:matrix([a,b],[c,d]); (%o3) matrix([a,b],[c,d]) (%i4) triangularize(m); (%o4) matrix([a,b],[0,a*db*c]) (%i5) triangularize(%); (%o5) matrix([a,b],[0,a^2*da*b*c]) (%i6) triangularize(%); (%o6) matrix([a,b],[0,a^3*da^2*b*c]) the matrix %o4 is already triangularized, but triangularize doesn't verify that. It continues to multiply for the pivot a in the first row. Maybe the problem is here. I don't understand why triangularize continues with such multiplications. In fact, in the above example, let a be the pivot: we want all the elements under it (only c here) vanish. To do that, we can sum to the 2nd row the first multiplied by c/a, that's what we see in %o4, that is correct. But if triangularization is repeated, triangularize multiplies for the pivot. Altough this behaviour should not change rank (it is a multiplication by a constant), maybe there are some issues related to it giving rise to errors? Surely, they give rise to high order polynomials (see ptriangularize for example, it checks if the matrix is already reduced and produces nicer expressions).  Comment By: Stefano Ferri (stefano_ferri) Date: 20100412 16:43 Message: Determinant and rank are different things: to reduce a matrix, one can multiply a row for a constant, or exchange the position of two rows: so, in general the determinant is changed, but this is normal. I don't think there is a documentation problem here, since this is a mathematical question. Maybe it could be mentioned in a note, but I think it is not necessary... But the rank must be preserved. This too is mathematics. What I'm seeing as a problem here is that triangularize has some problems handling symbolic matrices. As a consequence, evaluation after reduction yields differents results than reducing after evaluation. And with different results I don't mean a different appearence, but a different rank. Why in your example are you saying R2 < c*R1d*R2 is not rankpreserving?  Comment By: Barton Willis (willisbl) Date: 20100409 20:05 Message: Think about this: (%i7) triangularize(matrix([a,b],[c,d])); (%o7) matrix([a,b],[0,a*db*c]) There is no manifestly nonzero pivot, so triangularize uses a rank non preserving row operation R2 < c*R1d*R2. Even for matrices with rational entries, triangularize doesn't preserve the determinant: (%i13) determinant(triangularize(matrix([2,1],[6,7]))); (%o13) 16 (%i14) determinant(matrix([2,1],[6,7])); (%o14) 8 This too is a documentation problem, not a algorithm bug.  Comment By: Stefano Ferri (stefano_ferri) Date: 20100409 11:16 Message: I think it's not a documentationrelated bug or a lack of feature, but a real bug. What Maxima calls triangularize() is a matrix reduction by row. It doesn't matter how you are doing that, rank has to be preserved. There is a a very simple theorem that states that the rank of a reduced matrix is the same of the original matrix. As a consequence, the rank of a triangular matrix is equal the number of non null rows. Therefore, triangularize is doing something wrong in this case...  Comment By: Barton Willis (willisbl) Date: 20100409 02:47 Message: What do you expect for triangularize(matrix([a,b],[c,d])); A rank preserving triangularize would need to return a conditional statement. The user documentation doesn't say that triangularize preserves the rank. But it does say that triangularize does Gauss elimination, so it's reasonable to think that it is rankpreserving. This is more of a documentation bug than a algorithm bug, I think. You might like to try ptriangularize: (%i5) ptriangularize(matrix([8z, 7z],[1z, 3z]),z); (%o5) matrix([7,4],[0,17/7(3*z)/7]) (%i6) triangularize(matrix([8z, 7z],[1z, 3z]),z); (%o6) matrix([8z,7z],[0,173*z])  Comment By: Stefano Ferri (stefano_ferri) Date: 20100408 15:28 Message: Maxima version is 5.20.1. This problem is present both on Windows version, compiled with GCL, and Linux version, compiled with CLISP (Slackware package).  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2983881&group_id=4933 
From: SourceForge.net <noreply@so...>  20110804 08:11:32

Bugs item #2981631, was opened at 20100403 22:38 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2981631&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: Pending >Resolution: Works For Me Priority: 5 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Barton Willis (willisbl) Summary: abs_integrate bug Initial Comment: (%i3) load(abs_integrate)$ (%i4) integrate(min(x+1,1x,y+1,1y),y,1,1); first: empty argument. #0: abs_int_extra(q=x*signum(x)/4,x=y)(abs_integrate.mac line 215) #1: signum_int(q=x*signum(x)/4,x=y)(abs_integrate.mac line 136) #2: lambda([s],signum_int(s,x))(s=x*signum(x)/4)(abs_integrate.mac line 145)  an error. To debug this try: debugmode(true);  >Comment By: Dieter Kaiser (crategus) Date: 20110804 10:11 Message: The error of this bug report seems to be no longer present in the current Maxima version 5.25post. The result is a noun form. Maxima version: 5.25post Maxima build date: 21:58 8/3/2011 Host type: i686pclinuxgnu Lisp implementation type: SBCL Lisp implementation version: 1.0.45 (%i2) load(abs_integrate)$ (%i3) integrate(min(x+1,1x,y+1,1y),y,1,1); (%o3) 'integrate(min(x+1,1x,y+1,1y),y,1,1) Setting the resolution to "works for me" and the status to pending. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2981631&group_id=4933 