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}
(21) 
_{Jul}
(19) 
_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

From: Robert Dodier <robert_dodier@us...>  20170304 06:13:45

I don't see that behavior. I tried it with both Clisp (Maxima 5.39.0) and GCL (post5.39) on Ubuntu 14.04. Perhaps you can retry it with a new version? ~~~~ $ maxima Maxima 5.39.0 http://maxima.sourceforge.net using Lisp CLISP 2.49 (20100707) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) N:20; (%o1) 20 (%i2) F:taylor((1+z)^(1/3),z,0,N); 2 3 4 5 6 7 8 z z 5 z 10 z 22 z 154 z 374 z 935 z (%o2)/T/ 1 +    +    +    +    3 9 81 243 729 6561 19683 59049 9 10 11 12 13 14 21505 z 55913 z 147407 z 1179256 z 3174920 z 8617640 z +    +    +    1594323 4782969 14348907 129140163 387420489 1162261467 15 16 17 18 70664648 z 194327782 z 537259162 z 13431479050 z +    +    10460353203 31381059609 94143178827 2541865828329 19 20 37466757350 z 104906920580 z +    + . . . 7625597484987 22876792454961 (%i3) tmp:coeff(F,z,1); 1 (%o3)/R/  3 (%i4) F/tmp; 2 3 4 5 6 7 8 z 5 z 10 z 22 z 154 z 374 z 935 z (%o4)/T/ 3 + 1 z   +    +    +    3 27 81 243 2187 6561 19683 9 10 11 12 13 14 21505 z 55913 z 147407 z 1179256 z 3174920 z 8617640 z +    +    +    531441 1594323 4782969 43046721 129140163 387420489 15 16 17 18 70664648 z 194327782 z 537259162 z 13431479050 z +    +    3486784401 10460353203 31381059609 847288609443 19 20 37466757350 z 104906920580 z +    + . . . 2541865828329 7625597484987~~~~ ~~~~ $ sh maximalocal l gcl Maxima branch_5_39_base_82_g59007f6_dirty http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.12 Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) N:20; (%o1) 20 (%i2) F:taylor((1+z)^(1/3),z,0,N); 2 3 4 5 6 7 8 z z 5 z 10 z 22 z 154 z 374 z 935 z (%o2)/T/ 1 +    +    +    +    3 9 81 243 729 6561 19683 59049 9 10 11 12 13 14 21505 z 55913 z 147407 z 1179256 z 3174920 z 8617640 z +    +    +    1594323 4782969 14348907 129140163 387420489 1162261467 15 16 17 18 70664648 z 194327782 z 537259162 z 13431479050 z +    +    10460353203 31381059609 94143178827 2541865828329 19 20 37466757350 z 104906920580 z +    + . . . 7625597484987 22876792454961 (%i3) tmp:coeff(F,z,1); 1 (%o3)/R/  3 (%i4) F/tmp; 2 3 4 5 6 7 8 z 5 z 10 z 22 z 154 z 374 z 935 z (%o4)/T/ 3 + 1 z   +    +    +    3 27 81 243 2187 6561 19683 9 10 11 12 13 14 21505 z 55913 z 147407 z 1179256 z 3174920 z 8617640 z +    +    +    531441 1594323 4782969 43046721 129140163 387420489 15 16 17 18 70664648 z 194327782 z 537259162 z 13431479050 z +    +    3486784401 10460353203 31381059609 847288609443 19 20 37466757350 z 104906920580 z +    + . . . 2541865828329 7625597484987 ~~~~  ** [bugs:#3290] taylor expansion issue** **Status:** open **Group:** None **Created:** Fri Mar 03, 2017 06:58 PM UTC by igorpesando **Last Updated:** Fri Mar 03, 2017 06:58 PM UTC **Owner:** nobody Maxima version 5.37.2 using Lisp GNU Common Lisp (GCL) GCL 2.6.12 in ubuntu 16.04.2 Example: N:20; F:taylor((1+z)^(1/3),z,0,N); tmp:coeff(F,z,1); F/tmp; the last command gives an error (with segfault if one tries to repeat it) Condition in MACSYMATOPLEVEL [or a callee]: INTERNALSIMPLEERROR: Caught fatal error [memory may be damaged]  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: igorpesando <igorpesando@us...>  20170303 18:58:41

 ** [bugs:#3290] taylor expansion issue** **Status:** open **Group:** None **Created:** Fri Mar 03, 2017 06:58 PM UTC by igorpesando **Last Updated:** Fri Mar 03, 2017 06:58 PM UTC **Owner:** nobody Maxima version 5.37.2 using Lisp GNU Common Lisp (GCL) GCL 2.6.12 in ubuntu 16.04.2 Example: N:20; F:taylor((1+z)^(1/3),z,0,N); tmp:coeff(F,z,1); F/tmp; the last command gives an error (with segfault if one tries to repeat it) Condition in MACSYMATOPLEVEL [or a callee]: INTERNALSIMPLEERROR: Caught fatal error [memory may be damaged]  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Van <lvwarren@us...>  20170302 23:54:41

Command Line Maxima works. 2+2; works running more than a four line program crashes, see attached text file which can be transferred to wxMaxima bug tracker.  ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** wontfix **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Mon Feb 27, 2017 01:48 AM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170227 02:04:28

 ** [bugs:#3289] creating, removing rule causes symbol to stay on props list even after kill(all) and kill(props)** **Status:** open **Group:** None **Labels:** kill tellsimp tellsimpafter props **Created:** Mon Feb 27, 2017 02:04 AM UTC by Robert Dodier **Last Updated:** Mon Feb 27, 2017 02:04 AM UTC **Owner:** nobody `props` is the list of symbols which have userdefined properties. `kill(all)` and `kill(props)` is supposed to remove such properties and leave `props` empty. However it is possible to create a symbol which can't be removed from `props`. ~~~~ (%i1) tellsimp(f(1), 0); (%o1) [frule1, false] (%i2) put(f, a, b); (%o2) a (%i3) props; (%o3) [f] (%i4) kill(all); (%o0) done (%i1) props; (%o1) [f] ~~~~ I believe the problem is that killing the rule leaves the OPERATORS property behind, and this is enough to prevent `kill(props)` from removing `f` from `props`. ~~~~ (%i2) :lisp (symbolplist '$f) (OPERATORS NIL) ~~~~ It seems like `kill` could just `(remprop '$f 'operators)` instead of setting it to NIL. I think that would have the same effect, and also avoid confusion about whether or not there is really a property. Am I understanding this correctly? Other notes: Running `tests/rtest_rules.mac` creates two such symbols, `foo` and `f`. There is another test script (not sure which one) which puts un`kill`able MTIMES on `props`.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170227 01:48:43

 **status**: open > wontfix  **Comment**: I'm closing this report as it appears to be a wxMaxima problem on the basis of the evidence so far. I have opened a new issue on the wxMaxima bug tracker: https://github.com/andrejv/wxmaxima/issues/858 This item can be reopened if it turns out to be a Maxima problem after all. In the meantime please direct further discussion to the wxMaxima issue.  ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** wontfix **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Thu Feb 23, 2017 08:33 PM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Daniel Volinski <danielvolinski@us...>  20170225 08:58:27

Hi Robert, Please take a look at this problem:Based on Leon Magiera and Josef Böhm, CAS for Physics Examples, Problem I.4 Page 9http://rfdz.phnoe.ac.at/fileadmin/Mathematik_Uploads/ACDCA/materialien/CAS_Physics_1.pdf load("vect")$ e:ident(3)$ a_:[a_1,a_2,a_3]$ b_:[b_1,b_2,b_3]$ I'm trying to define the cross product using index notation First I define the LeviCivita Symbol ε(i,j,k):=e[i].express(e[j]~e[k])$ The result of ε(i,j,k) should be a scalar Some examples: ε(1,2,3); => 1ε(1,3,2); => 1ε(1,2,2); => [0,0,0]express(ε(1,2,3)*a_[2]*b_[3]); => a_2*b_3express(ε(1,3,2)*a_[2]*b_[3]); => a_2*b_3express(ε(1,2,2)*a_[2]*b_[2]); => [0,0,0] This should give the cross product makelist(sum(sum(ε(i,j,k)*a_[j]*b_[k],j,1,3),k,1,3),i,3); but the result is a list of three lists, this is wrong! I have to define my own cross product: mycross(va,vb):=[va[2]*vb[3]va[3]*vb[2], va[3]*vb[1]va[1]*vb[3], va[1]*vb[2]va[2]*vb[1]]$ Then the LeviCivita Symbol again: ε(i,j,k):=e[i].mycross(e[j],e[k])$ And now I get the correct result (the cross product): makelist(sum(sum(ε(i,j,k)*a_[j]*b_[k],j,1,3),k,1,3),i,3); [a_2*b_3a_3*b_2,a_3*b_1a_1*b_3,a_1*b_2a_2*b_1] Daniel El Sábado 25 de febrero de 2017 8:00, Robert Dodier <robert_dodier@...> escribió: I'm not sure this is a bug. 0 has the appropriate property that x + 0 is x even when x is a list. If you change your interpretation of 0 to mean "the additive identity", does anything work incorrectly? Are there algebraic identities which do not hold? [bugs:#3287] The cross product is wrong when I use the same vectorStatus: open Group: None Created: Mon Feb 20, 2017 05:11 PM UTC by Daniel Volinski Last Updated: Mon Feb 20, 2017 05:11 PM UTC Owner: nobodyThe cross product is wrong when I use the same vector:load(vect)$ express([0,1,1]~[0,1,2]);=> [1,0,0] / Correct! / express([0,1,1]~[0,1,1]);=> 0 / Incorrect, it should be a vector! /DanielSent from sourceforge.net because you indicated interest in https://sourceforge.net/p/maxima/bugs/3287/To unsubscribe from further messages, please visit https://sourceforge.net/auth/subscriptions/  ** [bugs:#3287] The cross product is wrong when I use the same vector** **Status:** open **Group:** None **Created:** Mon Feb 20, 2017 05:11 PM UTC by Daniel Volinski **Last Updated:** Sat Feb 25, 2017 06:00 AM UTC **Owner:** nobody The cross product is wrong when I use the same vector: load(vect)$ express([0,1,1]~[0,1,2]);=> [1,0,0] /* Correct! */ express([0,1,1]~[0,1,1]);=> 0 /* Incorrect, it should be a vector! */ Daniel  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170225 06:00:45

I'm not sure this is a bug. 0 has the appropriate property that `x + 0` is `x` even when `x` is a list. If you change your interpretation of 0 to mean "the additive identity", does anything work incorrectly? Are there algebraic identities which do not hold?  ** [bugs:#3287] The cross product is wrong when I use the same vector** **Status:** open **Group:** None **Created:** Mon Feb 20, 2017 05:11 PM UTC by Daniel Volinski **Last Updated:** Mon Feb 20, 2017 05:11 PM UTC **Owner:** nobody The cross product is wrong when I use the same vector: load(vect)$ express([0,1,1]~[0,1,2]);=> [1,0,0] /* Correct! */ express([0,1,1]~[0,1,1]);=> 0 /* Incorrect, it should be a vector! */ Daniel  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170223 20:33:56

It appears from the log file you posted that wxMaxima has crashed (with a segfault). Does command line Maxima work OK?  ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** open **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Tue Feb 21, 2017 06:45 AM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Van <lvwarren@us...>  20170221 06:46:00

After each Maxima crash, a common lisp process 'sbcl' is left running on the machine and this process has very high CPU utilization. If not detected several of these processes taken together can overheat some machines, Stuxnet style. I have to keep on an icepack to keep the CPU temperature below 210 degrees. Bringing up the activity manager or doing a ps au shows the zombie common lisps which can be killed one at a time.  ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** open **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Tue Feb 21, 2017 06:32 AM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Van <lvwarren@us...>  20170221 06:32:32

 ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** open **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Tue Feb 21, 2017 06:32 AM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Daniel Volinski <danielvolinski@us...>  20170220 17:11:09

 ** [bugs:#3287] The cross product is wrong when I use the same vector** **Status:** open **Group:** None **Created:** Mon Feb 20, 2017 05:11 PM UTC by Daniel Volinski **Last Updated:** Mon Feb 20, 2017 05:11 PM UTC **Owner:** nobody The cross product is wrong when I use the same vector: load(vect)$ express([0,1,1]~[0,1,2]);=> [1,0,0] /* Correct! */ express([0,1,1]~[0,1,1]);=> 0 /* Incorrect, it should be a vector! */ Daniel  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Gilles Schintgen <gschint@us...>  20170205 18:33:35

 ** [bugs:#3286] implicit_plot treats some discontinuities as sign changes** **Status:** open **Group:** None **Created:** Sun Feb 05, 2017 06:33 PM UTC by Gilles Schintgen **Last Updated:** Sun Feb 05, 2017 06:33 PM UTC **Owner:** nobody It seems that implicit_plot treats a change from positive to complex as a sign change. Here's a simple example: load(implicit_plot)$ implicit_plot(1/sqrt(x^2+y^21),[x,2,2],[y,2,2])$ Result: the unit circle is drawn even though the expression to be plotted is either positive or nonreal. It's not even approaching zero when (x,y) is close to the unit circle. A cursory look at the code seems to confirm that implicit_plot tries to detect sign changes but is actually looking for changes from larger than ip_epsilon to "everything else" (in function sampledata).  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: David Scherfgen <tomasriker@us...>  20170203 08:05:15

 Description has changed: Diff: ~~~~  old +++ new @@ 1,4 +1,4 @@ `radcan` fails to simplify an expression to zero that contains no variables, only logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: +`radcan` fails to simplify an expression to zero that contains no variables, only exponentials and logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: ~~~ (%i1) X : 2^(log(3)/(log(2)log(3)))3^(log(3)/(log(2)log(3))+1)$ (%i2) radcan(X); ~~~~  ** [bugs:#3285] radcan: failure to simplify to zero depending on slight change and/or ?expsumsplit** **Status:** open **Group:** None **Labels:** radcan **Created:** Fri Feb 03, 2017 08:04 AM UTC by David Scherfgen **Last Updated:** Fri Feb 03, 2017 08:04 AM UTC **Owner:** nobody `radcan` fails to simplify an expression to zero that contains no variables, only exponentials and logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: ~~~ (%i1) X : 2^(log(3)/(log(2)log(3)))3^(log(3)/(log(2)log(3))+1)$ (%i2) radcan(X); (%o2) (%e^(log(3)^2/(log(3)log(2)))*(%e^(log(3)^2/(log(3)log(2)))3*2^(log(3)/(log(3)log(2)))))/2^(log(3)/(log(3)log(2))) (%i3) radcan(logcontract(X)); (%o3) 0 (%i4) Y : 2^(log(3)/(log(2)log(3)))1/3^(log(2)/(log(3)log(2)))$ (%i5) radcan(Y); (%o5) 0 (%i6) radcan(X), ?expsumsplit:false; (%o6) 0 ~~~ Thanks to Richard Fateman for investigating.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: David Scherfgen <tomasriker@us...>  20170203 08:04:37

 ** [bugs:#3285] radcan: failure to simplify to zero depending on slight change and/or ?expsumsplit** **Status:** open **Group:** None **Labels:** radcan **Created:** Fri Feb 03, 2017 08:04 AM UTC by David Scherfgen **Last Updated:** Fri Feb 03, 2017 08:04 AM UTC **Owner:** nobody `radcan` fails to simplify an expression to zero that contains no variables, only logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: ~~~ (%i1) X : 2^(log(3)/(log(2)log(3)))3^(log(3)/(log(2)log(3))+1)$ (%i2) radcan(X); (%o2) (%e^(log(3)^2/(log(3)log(2)))*(%e^(log(3)^2/(log(3)log(2)))3*2^(log(3)/(log(3)log(2)))))/2^(log(3)/(log(3)log(2))) (%i3) radcan(logcontract(X)); (%o3) 0 (%i4) Y : 2^(log(3)/(log(2)log(3)))1/3^(log(2)/(log(3)log(2)))$ (%i5) radcan(Y); (%o5) 0 (%i6) radcan(X), ?expsumsplit:false; (%o6) 0 ~~~ Thanks to Richard Fateman for investigating.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Jinsong Zhao <jszhao@us...>  20170131 09:59:05

To the issue 4. I have wrote a Lisp script to replace the python script used during compiling the document. I attach the extract_categories1.awk and extract_categories1.sed here, FYI. The only issue with this approach is I don't know how to set correct Lisp program name in extract_categories.sh. ~~~  > tmpmakecategories.py + > tmpmakecategories.lisp ${PYTHONBIN:python} tmpmakecategories.py +#echo ${SBCL_NAME} +/d/Program/Steel\ Bank\ Common\ Lisp/1.3.10/sbcl.exe script tmpmakecategories.lisp ~~~ I meant here, now, it's hard coded in the script, it's not the right way. Best regards, Jinsong Attachments:  [extract_categories1.awk](https://sourceforge.net/p/maxima/bugs/_discuss/thread/9d22061f/886c/attachment/extract_categories1.awk) (1.3 kB; application/octetstream)  [extract_categories1.sed](https://sourceforge.net/p/maxima/bugs/_discuss/thread/9d22061f/886c/attachment/extract_categories1.sed) (617 Bytes; application/octetstream)  ** [bugs:#3284] maxima build system** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 03:04 PM UTC by Jinsong Zhao **Last Updated:** Sun Jan 29, 2017 03:04 PM UTC **Owner:** nobody There are several issues that related with the maxima building system on Windows platform. I don't know if there are bugs, however, it makes trouble or confusion when I try to make a maxima installer for Windows. I posted them in the maximadiscuss mailing list. 1) After './configure...' and 'make', 'make install' gives the following error message: /usr/bin/install: cannot create regular file 'C:/maxima5.39.0sbcl/bin/xmaxima': File exists 'xmaxima' in interface/xmaxima seems to be a script, and not need by XMaxima on Windows platform. On Windows, 'xmaxima.exe' is the correct executable file. This error prevents the running of make install. I have to changed the Makefile.in in interface/xmaxima from $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  exit $$?; \ to $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  echo "xxx" $$?; \ according to https://sourceforge.net/projects/maxima/files/MaximaWindows/5.XX.XWindows%20Build%20System/README.md Then make install could be run successfully. If possible, could the 'xmaxima' script be removed from makefile for Windows platform? 2) Could the 'xmaxima' script be removed from the 'maxima.iss' on Windows platform? 3) There appears a directory named ".cirtcl" outside the maxima toplevel directory during the making process. In the directory, there are two files, i.e., mgnuplot.bat and xmaixma.bat. I don't know if they are needed, and why there are be put outside the toplevel directory. It seems they are related with xmaxima interface. 4) Could we not compile the document again? In the doc/, there are correct html and info files, why we generate them again? I encountered the most difficult thing is the python script 'tmpmakecategories.py does not work in MSYS2 system. I don't know why. So, if possible, could configure provide a option that disable doc generation. Best regards, Jinsong  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170130 01:48:23

 **labels**: gamma, gamma_greek, gamma_incomp > gamma, gamma_greek, gamma_incomplete  ** [bugs:#3277] no numerical evaluation for gamma_greek** **Status:** open **Group:** None **Labels:** gamma gamma_greek gamma_incomplete **Created:** Thu Jan 19, 2017 12:50 AM UTC by kofoid **Last Updated:** Mon Jan 30, 2017 01:47 AM UTC **Owner:** nobody I am running maxima 5.38.0 under OS X 10.10.5. Of the various gammaassociated functions, only "gamma" itself produces output. As an example, attempting to run "gamma_greek(10,1)" produces itself as output.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170130 01:47:21

 **labels**: > gamma, gamma_greek, gamma_incomp  **summary**: gammarelated functions under OS X > no numerical evaluation for gamma_greek  **Comment**: Adjusted title to focus on gamma_greek and omit OS.  ** [bugs:#3277] no numerical evaluation for gamma_greek** **Status:** open **Group:** None **Labels:** gamma gamma_greek gamma_incomp **Created:** Thu Jan 19, 2017 12:50 AM UTC by kofoid **Last Updated:** Sun Jan 22, 2017 10:23 PM UTC **Owner:** nobody I am running maxima 5.38.0 under OS X 10.10.5. Of the various gammaassociated functions, only "gamma" itself produces output. As an example, attempting to run "gamma_greek(10,1)" produces itself as output.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170129 20:04:39

 **labels**: > lsquares, share  **status**: open > closed  **Comment**: Fixed by commit b58a547. Thanks for the pointing it out. Closing this report as fixed.  ** [bugs:#3282] lsquares and lists of list of data to be fitted on** **Status:** closed **Group:** None **Labels:** lsquares share **Created:** Tue Jan 24, 2017 06:58 PM UTC by Gunter Königsmann **Last Updated:** Thu Jan 26, 2017 11:46 PM UTC **Owner:** nobody I have a list of data sets that are huge => running lsquares_estimates() on them to fit them on an equation would need too big amounts of time. Normally lsquares_estimates_approximate() helps in this case as it skips the step that tries to find an exact solution. But as I have not a single data set, but a list of data sets that somehow backfires: ~~~ (%i34) load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; lsquares_mse(data[1],[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o30) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) [matrix( [0, 0], [1, 1], [2, 2] ),matrix( [2, 2], [1, 1], [2, 2] )] (eqtn) y=a*x^2+b*x+c (%o33) sum((m6[i,2]a*m6[i,1]^2b*m6[i,1]c)^2,i,1,1) Maxima encountered a Lisp error: The value ((MTIMES SIMP) 2.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2) ((MPLUS SIMP) 1.0 ((MTIMES SIMP) 1.0 ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1)) ((MTIMES SIMP) 1.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2)) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 2))) is not of type DOUBLEFLOAT. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. ~~~ Only running lsquares_estimates on a single entry of the list works: ~~~ (%i29) load("lsquares"); data:matrix([0,0],[1,1],[2,2]); eqtn:y=a*x^2+b*x+c; lsquares_mse(data,[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o25) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) matrix( [0, 0], [1, 1], [2, 2] ) (eqtn) y=a*x^2+b*x+c (%o28) sum((data[i,2]a*data[i,1]^2b*data[i,1]c)^2,i,1,3)/3 ************************************************* N= 3 NUMBER OF CORRECTIONS=25 INITIAL VALUES F= 1.000000000000000D+01 GNORM= 1.753726191728787D+01 ************************************************* I NFN FUNC GNORM STEPLENGTH 1 2 2.796937667198332D01 1.961076084276657D+00 5.702144409522791D02 2 3 1.464269617620298D01 4.983306664572200D01 1.000000000000000D+00 ~~~ snip ~~~ 13 14 3.397773998012277D06 4.742450181229817D03 1.000000000000000D+00 14 15 2.643557109421608D08 2.597101638571121D04 1.000000000000000D+00 THE MINIMIZATION TERMINATED WITHOUT DETECTING ERRORS. IFLAG = 0 (%o29) [[a=3.405898581949557*10^4,b=1.000676542954337,c=1.342533507359797*10^4]] ~~~ I have been able to ship around this problem by putting the data set I am working on into a temporary variable before running lsquares: ~~~ load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; makelist(block([tmp], tmp:data[i], lsquares_estimates_approximate(lsquares_mse(tmp,[x,y],eqtn),[a,b,c]) ), i,1,length(data) ); ...but I am convinced that the fact that I needed this is a bug. ~~~  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Stavros Macrakis <macrakis@us...>  20170129 19:48:03

More, simpler, examples: taylor(rat(1/z+1),z,0,1) => unfamiliar singularity error taylor( 1/z+1,z,0,1) => 1/z+1 taylor(rat(z+1),z,inf,0) => unfamiliar singularity error taylor( z+1 ,z,inf,0) => z+1 Maxima branch_5_37_base_528_g802ae7f using Lisp SBCL 1.2.11  ** [bugs:#3068] taylor of CRE fails** **Status:** open **Group:** None **Created:** Tue Jan 12, 2016 07:50 PM UTC by Stavros Macrakis **Last Updated:** Tue Jan 12, 2016 07:50 PM UTC **Owner:** nobody taylor(rat(x/(1+x),x,inf,2) => taylor: encountered an unfamiliar singularity though taylor((x/(1+x),x,inf,2) => 1  1/x + 1/x^2 + ...  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Jinsong Zhao <jszhao@us...>  20170129 15:04:46

 ** [bugs:#3284] maxima build system** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 03:04 PM UTC by Jinsong Zhao **Last Updated:** Sun Jan 29, 2017 03:04 PM UTC **Owner:** nobody There are several issues that related with the maxima building system on Windows platform. I don't know if there are bugs, however, it makes trouble or confusion when I try to make a maxima installer for Windows. I posted them in the maximadiscuss mailing list. 1) After './configure...' and 'make', 'make install' gives the following error message: /usr/bin/install: cannot create regular file 'C:/maxima5.39.0sbcl/bin/xmaxima': File exists 'xmaxima' in interface/xmaxima seems to be a script, and not need by XMaxima on Windows platform. On Windows, 'xmaxima.exe' is the correct executable file. This error prevents the running of make install. I have to changed the Makefile.in in interface/xmaxima from $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  exit $$?; \ to $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  echo "xxx" $$?; \ according to https://sourceforge.net/projects/maxima/files/MaximaWindows/5.XX.XWindows%20Build%20System/README.md Then make install could be run successfully. If possible, could the 'xmaxima' script be removed from makefile for Windows platform? 2) Could the 'xmaxima' script be removed from the 'maxima.iss' on Windows platform? 3) There appears a directory named ".cirtcl" outside the maxima toplevel directory during the making process. In the directory, there are two files, i.e., mgnuplot.bat and xmaixma.bat. I don't know if they are needed, and why there are be put outside the toplevel directory. It seems they are related with xmaxima interface. 4) Could we not compile the document again? In the doc/, there are correct html and info files, why we generate them again? I encountered the most difficult thing is the python script 'tmpmakecategories.py does not work in MSYS2 system. I don't know why. So, if possible, could configure provide a option that disable doc generation. Best regards, Jinsong  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Dima Pasechnik <dimpase@us...>  20170129 11:49:27

Here is how to get this with Maxima only (5.39.0 with SBCL or ECL): ~~~ s1: 2*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)*(cos(z)*sin(y) + cos(y)*sin(z))*((cos(1/2*x)*cos(z)  cos(y)*cos(z) + sin(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)*cos(z)*sin(y) + cos(y)^2*sin(z)  (cos(z)*sin(y) + cos(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))/((cos(1/2*x)*cos(y)  1)*((cos(y)*cos(z)  sin(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)*cos(y)  1) + (cos(z)*sin(y) + cos(y)*sin(z))*((cos(1/2*x)*cos(z)  cos(y)*cos(z) + sin(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)*cos(z)*sin(y) + cos(y)^2*sin(z)  (cos(z)*sin(y) + cos(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))*sin(1/2*x)/((cos(1/2*x)*cos(y)  1)*((cos(z)*sin(y)  cos(1/2*x)*sin(z) + cos(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)^2*cos(z)  cos(y)*sin(y)*sin(z)  (cos(y)*cos(z)  sin(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))))*((cos(z)*sin(y)  cos(1/2*x)*sin(z) + cos(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)^2*cos(z)  cos(y)*sin(y)*sin(z)  (cos(y)*cos(z)  sin(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))^2*sin(y))  2*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/(((cos(z)*sin(y)  cos(1/2*x)*sin(z) + cos(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)^2*cos(z)  cos(y)*sin(y)*sin(z)  (cos(y)*cos(z)  sin(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))*sin(1/2*x)*sin(y))$ trigsimp(s1); s2: 2*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)*(cos(b)*sin(a) + cos(a)*sin(b))*((cos(1/2*x)*cos(b)  cos(a)*cos(b) + sin(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)*cos(b)*sin(a) + cos(a)^2*sin(b)  (cos(b)*sin(a) + cos(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))/((cos(1/2*x)*cos(a)  1)*((cos(a)*cos(b)  sin(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)*cos(a)  1) + (cos(b)*sin(a) + cos(a)*sin(b))*((cos(1/2*x)*cos(b)  cos(a)*cos(b) + sin(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)*cos(b)*sin(a) + cos(a)^2*sin(b)  (cos(b)*sin(a) + cos(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))*sin(1/2*x)/((cos(1/2*x)*cos(a)  1)*((cos(b)*sin(a)  cos(1/2*x)*sin(b) + cos(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)^2*cos(b)  cos(a)*sin(a)*sin(b)  (cos(a)*cos(b)  sin(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))))*((cos(b)*sin(a)  cos(1/2*x)*sin(b) + cos(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)^2*cos(b)  cos(a)*sin(a)*sin(b)  (cos(a)*cos(b)  sin(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))^2*sin(a))  2*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/(((cos(b)*sin(a)  cos(1/2*x)*sin(b) + cos(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)^2*cos(b)  cos(a)*sin(a)*sin(b)  (cos(a)*cos(b)  sin(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))*sin(1/2*x)*sin(a))$ trigsimp(s2); ~~~ gives ~~~ Maxima 5.39.0 http://maxima.sourceforge.net using Lisp SBCL 1.3.12 Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) (%i2) x (%o2) (2 cos(y)  2 cos()) sin(y) sin(z) 2 x 2 + (2 cos() cos(y)  2 cos (y)) cos(z) 2 (%i3) (%i4) Quotient by a polynomial of higher degree  an error. To debug this try: debugmode(true); ~~~ Expression s2 is obtained from s1 by renaming y>a, z>b.  ** [bugs:#3283] trigsimp fails for some order of the variables** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 10:26 AM UTC by Eric Gourgoulhon **Last Updated:** Sun Jan 29, 2017 10:26 AM UTC **Owner:** nobody Hi, trigsimp may fail depending on the (alphabetical?) order of the variables; please see this SageMath bug report: https://groups.google.com/forum/#!topic/sagedevel/bJBgJhvRX4  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Eric Gourgoulhon <egourgoulhon@us...>  20170129 10:26:26

 ** [bugs:#3283] trigsimp fails for some order of the variables** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 10:26 AM UTC by Eric Gourgoulhon **Last Updated:** Sun Jan 29, 2017 10:26 AM UTC **Owner:** nobody Hi, trigsimp may fail depending on the (alphabetical?) order of the variables; please see this SageMath bug report: https://groups.google.com/forum/#!topic/sagedevel/bJBgJhvRX4  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170126 23:46:18

Yes, this is a bug in `lsquares_mse`. The problem is that it introduces a temporary variable and tries to assign the data to the temporary variable, but it is incorrect. I think I see how to fix it.  ** [bugs:#3282] lsquares and lists of list of data to be fitted on** **Status:** open **Group:** None **Created:** Tue Jan 24, 2017 06:58 PM UTC by Gunter Königsmann **Last Updated:** Tue Jan 24, 2017 06:58 PM UTC **Owner:** nobody I have a list of data sets that are huge => running lsquares_estimates() on them to fit them on an equation would need too big amounts of time. Normally lsquares_estimates_approximate() helps in this case as it skips the step that tries to find an exact solution. But as I have not a single data set, but a list of data sets that somehow backfires: ~~~ (%i34) load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; lsquares_mse(data[1],[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o30) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) [matrix( [0, 0], [1, 1], [2, 2] ),matrix( [2, 2], [1, 1], [2, 2] )] (eqtn) y=a*x^2+b*x+c (%o33) sum((m6[i,2]a*m6[i,1]^2b*m6[i,1]c)^2,i,1,1) Maxima encountered a Lisp error: The value ((MTIMES SIMP) 2.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2) ((MPLUS SIMP) 1.0 ((MTIMES SIMP) 1.0 ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1)) ((MTIMES SIMP) 1.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2)) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 2))) is not of type DOUBLEFLOAT. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. ~~~ Only running lsquares_estimates on a single entry of the list works: ~~~ (%i29) load("lsquares"); data:matrix([0,0],[1,1],[2,2]); eqtn:y=a*x^2+b*x+c; lsquares_mse(data,[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o25) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) matrix( [0, 0], [1, 1], [2, 2] ) (eqtn) y=a*x^2+b*x+c (%o28) sum((data[i,2]a*data[i,1]^2b*data[i,1]c)^2,i,1,3)/3 ************************************************* N= 3 NUMBER OF CORRECTIONS=25 INITIAL VALUES F= 1.000000000000000D+01 GNORM= 1.753726191728787D+01 ************************************************* I NFN FUNC GNORM STEPLENGTH 1 2 2.796937667198332D01 1.961076084276657D+00 5.702144409522791D02 2 3 1.464269617620298D01 4.983306664572200D01 1.000000000000000D+00 ~~~ snip ~~~ 13 14 3.397773998012277D06 4.742450181229817D03 1.000000000000000D+00 14 15 2.643557109421608D08 2.597101638571121D04 1.000000000000000D+00 THE MINIMIZATION TERMINATED WITHOUT DETECTING ERRORS. IFLAG = 0 (%o29) [[a=3.405898581949557*10^4,b=1.000676542954337,c=1.342533507359797*10^4]] ~~~ I have been able to ship around this problem by putting the data set I am working on into a temporary variable before running lsquares: ~~~ load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; makelist(block([tmp], tmp:data[i], lsquares_estimates_approximate(lsquares_mse(tmp,[x,y],eqtn),[a,b,c]) ), i,1,length(data) ); ...but I am convinced that the fact that I needed this is a bug. ~~~  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Peter Gasparovic <buteo@us...>  20170126 10:03:36

I confirm, that your example is correctly solved in Maxima 5.21.1. Problem with solve inconsistency is present with newest Maxima 5.38.1. I have the same problem with: eq: y=(x^2/2+1)^k; sol_a1: solve([eq], [x]); (I.e. correctly solved in 5.21.1, problem in 5.38.1.) In new version you can get the correct result by repeating solve sommand. In my case the first solve command works when I assign as a first command: simp:false; However that doesn't work on your example. I described it on Stackexchange: [Why “simp:false” helps Maxima to solve equation?](http://math.stackexchange.com/questions/2106603/whysimpfalsehelpsmaximatosolveequation)  ** [bugs:#3070] solve() inconsistency** **Status:** open **Group:** None **Created:** Wed Jan 13, 2016 11:35 AM UTC by Jaime E. Villate **Last Updated:** Thu Jan 14, 2016 08:19 PM UTC **Owner:** nobody The result of the following two commands should be the same: ~~~~ (%i2) solve(sqrt((x^2+1)/2) = 1); (%o2) [x = 1,x = 1] (%i3) solve(sqrt(x^2/2+1/2) = 1); (%o3) [sqrt(x^2+1)/sqrt(2) = 1] ~~~~ %o3 used to be the same as %o2 in older versions (at least in 5.21).  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Gunter Königsmann <peterpall@us...>  20170124 18:58:05

 ** [bugs:#3282] lsquares and lists of list of data to be fitted on** **Status:** open **Group:** None **Created:** Tue Jan 24, 2017 06:58 PM UTC by Gunter Königsmann **Last Updated:** Tue Jan 24, 2017 06:58 PM UTC **Owner:** nobody I have a list of data sets that are huge => running lsquares_estimates() on them to fit them on an equation would need too big amounts of time. Normally lsquares_estimates_approximate() helps in this case as it skips the step that tries to find an exact solution. But as I have not a single data set, but a list of data sets that somehow backfires: ~~~ (%i34) load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; lsquares_mse(data[1],[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o30) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) [matrix( [0, 0], [1, 1], [2, 2] ),matrix( [2, 2], [1, 1], [2, 2] )] (eqtn) y=a*x^2+b*x+c (%o33) sum((m6[i,2]a*m6[i,1]^2b*m6[i,1]c)^2,i,1,1) Maxima encountered a Lisp error: The value ((MTIMES SIMP) 2.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2) ((MPLUS SIMP) 1.0 ((MTIMES SIMP) 1.0 ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1)) ((MTIMES SIMP) 1.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2)) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 2))) is not of type DOUBLEFLOAT. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. ~~~ Only running lsquares_estimates on a single entry of the list works: ~~~ (%i29) load("lsquares"); data:matrix([0,0],[1,1],[2,2]); eqtn:y=a*x^2+b*x+c; lsquares_mse(data,[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o25) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) matrix( [0, 0], [1, 1], [2, 2] ) (eqtn) y=a*x^2+b*x+c (%o28) sum((data[i,2]a*data[i,1]^2b*data[i,1]c)^2,i,1,3)/3 ************************************************* N= 3 NUMBER OF CORRECTIONS=25 INITIAL VALUES F= 1.000000000000000D+01 GNORM= 1.753726191728787D+01 ************************************************* I NFN FUNC GNORM STEPLENGTH 1 2 2.796937667198332D01 1.961076084276657D+00 5.702144409522791D02 2 3 1.464269617620298D01 4.983306664572200D01 1.000000000000000D+00 ~~~ snip ~~~ 13 14 3.397773998012277D06 4.742450181229817D03 1.000000000000000D+00 14 15 2.643557109421608D08 2.597101638571121D04 1.000000000000000D+00 THE MINIMIZATION TERMINATED WITHOUT DETECTING ERRORS. IFLAG = 0 (%o29) [[a=3.405898581949557*10^4,b=1.000676542954337,c=1.342533507359797*10^4]] ~~~ I have been able to ship around this problem by putting the data set I am working on into a temporary variable before running lsquares: ~~~ load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; makelist(block([tmp], tmp:data[i], lsquares_estimates_approximate(lsquares_mse(tmp,[x,y],eqtn),[a,b,c]) ), i,1,length(data) ); ...but I am convinced that the fact that I needed this is a bug. ~~~  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 