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
(30) |
Aug
(48) |
Sep
(39) |
Oct
(30) |
Nov
(75) |
Dec
(28) |
2018 |
Jan
(70) |
Feb
|
Mar
|
Apr
(1) |
May
|
Jun
|
Jul
|
Aug
|
Sep
|
Oct
|
Nov
|
Dec
(1) |
2019 |
Jan
(19) |
Feb
(61) |
Mar
(14) |
Apr
(7) |
May
(5) |
Jun
(17) |
Jul
(5) |
Aug
(7) |
Sep
(11) |
Oct
(2) |
Nov
(17) |
Dec
(9) |
2020 |
Jan
(8) |
Feb
(8) |
Mar
(12) |
Apr
(17) |
May
(2) |
Jun
(10) |
Jul
(24) |
Aug
(6) |
Sep
(16) |
Oct
(3) |
Nov
(10) |
Dec
(40) |
2021 |
Jan
(53) |
Feb
(18) |
Mar
(20) |
Apr
(11) |
May
(23) |
Jun
(37) |
Jul
(28) |
Aug
(32) |
Sep
(105) |
Oct
(81) |
Nov
(109) |
Dec
(41) |
2022 |
Jan
(139) |
Feb
(82) |
Mar
(96) |
Apr
(51) |
May
(58) |
Jun
(104) |
Jul
(32) |
Aug
(61) |
Sep
(37) |
Oct
(25) |
Nov
(94) |
Dec
(81) |
2023 |
Jan
(113) |
Feb
(77) |
Mar
(98) |
Apr
(43) |
May
(48) |
Jun
(28) |
Jul
(72) |
Aug
(40) |
Sep
(44) |
Oct
(61) |
Nov
(70) |
Dec
(94) |
2024 |
Jan
(101) |
Feb
(21) |
Mar
(66) |
Apr
(88) |
May
(55) |
Jun
(109) |
Jul
(57) |
Aug
(79) |
Sep
|
Oct
|
Nov
|
Dec
|
From: Raymond T. <rt...@us...> - 2024-08-25 00:11:18
|
That's a pretty old bug, going back to 2010. The solution there was wrong. I have a fix coming soon. Fairly easy, but needs some more testing. --- **[bugs:#4352] elliptic_e(1,1.23) --> lisp numeric error (introduces complex %i)** **Status:** open **Group:** None **Labels:** elliptic functions **Created:** Sat Aug 24, 2024 04:48 PM UTC by Richard Fateman **Last Updated:** Sun Aug 25, 2024 12:05 AM UTC **Owner:** Raymond Toy (%i39) elliptic_e(1,1.23); Maxima encountered a Lisp error: The value ((MPLUS SIMP) 0.822727348982166 ((MTIMES SIMP) 0.16698253315507405 $%I)) is not of type NUMBER Automatically continuing. ............ build_info(version="5.47.0",timestamp="2023-06-02 20:31:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.3.2", --- Sent from sourceforge.net because max...@li... 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: Raymond T. <rt...@us...> - 2024-08-25 00:05:42
|
- **labels**: --> elliptic functions - **assigned_to**: Raymond Toy --- **[bugs:#4352] elliptic_e(1,1.23) --> lisp numeric error (introduces complex %i)** **Status:** open **Group:** None **Labels:** elliptic functions **Created:** Sat Aug 24, 2024 04:48 PM UTC by Richard Fateman **Last Updated:** Sat Aug 24, 2024 04:48 PM UTC **Owner:** Raymond Toy (%i39) elliptic_e(1,1.23); Maxima encountered a Lisp error: The value ((MPLUS SIMP) 0.822727348982166 ((MTIMES SIMP) 0.16698253315507405 $%I)) is not of type NUMBER Automatically continuing. ............ build_info(version="5.47.0",timestamp="2023-06-02 20:31:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.3.2", --- Sent from sourceforge.net because max...@li... 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: Richard F. <fa...@us...> - 2024-08-24 16:48:25
|
--- **[bugs:#4352] elliptic_e(1,1.23) --> lisp numeric error (introduces complex %i)** **Status:** open **Group:** None **Created:** Sat Aug 24, 2024 04:48 PM UTC by Richard Fateman **Last Updated:** Sat Aug 24, 2024 04:48 PM UTC **Owner:** nobody (%i39) elliptic_e(1,1.23); Maxima encountered a Lisp error: The value ((MPLUS SIMP) 0.822727348982166 ((MTIMES SIMP) 0.16698253315507405 $%I)) is not of type NUMBER Automatically continuing. ............ build_info(version="5.47.0",timestamp="2023-06-02 20:31:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.3.2", --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-22 22:00:44
|
- **status**: open --> closed - **Comment**: Fixed by commit [ 9a2271d ]. --- **[bugs:#4350] arrayinfo complains "not an array" when supplied a Lisp array or hash table** **Status:** closed **Group:** None **Labels:** arrayinfo **Created:** Mon Aug 19, 2024 08:59 PM UTC by Robert Dodier **Last Updated:** Mon Aug 19, 2024 08:59 PM UTC **Owner:** nobody `arrayinfo` complains "not an array" when its argument is a Lisp array or Lisp hash table. However, when its argument is a symbol which has the Lisp array or hash table bound to it, then `arrayinfo` knows what to do. ``` (%i33) aa: make_array (any, 4); (%o33) {Lisp Array: #(NIL NIL NIL NIL)} (%i34) arrayinfo (aa); (%o34) [declared, 1, [3]] (%i35) apply (arrayinfo, [aa]); arrayinfo: {Lisp Array: #(NIL NIL NIL NIL)} is not an array. -- an error. To debug this try: debugmode(true); (%i36) hh: make_array (hashed, 1); (%o36) {Lisp Array: #<HASH-TABLE :TEST EQUAL :COUNT 1 {1004C322F3}>} (%i37) hh["foo"]: "bar"; (%o37) bar (%i38) arrayinfo (hh); (%o38) [hash_table, 1, foo] (%i39) apply (arrayinfo, [hh]); arrayinfo: {Lisp Array: #<HASH-TABLE :TEST EQUAL :COUNT 2 {1004C322F3}>} is not an array. -- an error. To debug this try: debugmode(true); ``` I bumped into this problem in the following context. I have a list of arrays, and I'm trying to call `arrayinfo` on each one. %o41 fails because getting the value of `L[1]` would require calling MEVAL, which, I'm happy to report, `arrayinfo` doesn't call. OK so far. The standard workaround is to call `apply(foo, [...])` for an argument-quoting function `foo`. But that doesn't work either, for `arrayinfo`: ``` (%i40) L: [aa]; (%o40) [{Lisp Array: #(NIL NIL NIL NIL)}] (%i41) arrayinfo (L[1]); arrayinfo: L is not an array. 1 -- an error. To debug this try: debugmode(true); (%i42) apply (arrayinfo, [L[1]]); arrayinfo: {Lisp Array: #(NIL NIL NIL NIL)} is not an array. -- an error. To debug this try: debugmode(true); ``` The fix is just to handle Lisp array and hash table values in addition to symbols. --- Sent from sourceforge.net because max...@li... 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: Raymond T. <rt...@us...> - 2024-08-21 13:19:27
|
- Description has changed: Diff: ~~~~ --- old +++ new @@ -2,8 +2,8 @@ Not positively sure this is a bug but anyway: ~~~ -~~~ + [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; -~~~ + ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? ~~~~ - **Comment**: Consider using the preview button to see what your message looks like before saving/sending the message. There's also a "Formatting Help" button with useful information. I have to use that fairly often because I can't remember. --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Wed Aug 21, 2024 12:41 PM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: Barton W. <wil...@us...> - 2024-08-21 12:41:21
|
As far as I know, regardless of your operating system there is no requirement to uninstall Maxima before installing a new Maxima version. You can have multiple versions of Maxima installed and you can even run them simultaneously. If disk space is short, you can either delete the old version or install a new version on a flash drive (at least the last I tried, this works). Of course if you have tweaked the Maxima source code or have appended your own files to Maxima, you'll need to manually preserve these before deleting Maxima. I'm not sure about Maxima initialization files--if you have something important in these files, be safe and copy them before deleting. If for some reason, the new version isn't satisfactory, you can download and install your previous version. Thus there isn't much risk in deleting an old version and installing a new. --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Wed Aug 21, 2024 11:43 AM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: dan h. <zm...@us...> - 2024-08-21 11:43:16
|
to Barton Williams. ok i posted that last before i saw your latest post. Anyway with regard to that should i uninstall my present version of maxima first ? or just what is best way to do that and not waste a lot of diskspace. --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Wed Aug 21, 2024 11:33 AM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: dan h. <zm...@us...> - 2024-08-21 11:33:20
|
But strangely enough the following give correct answer. **``** [rombergmin:4,romberg(abs(sin(x)),x,-1,1),romberg(sqrt(sin(x)^2),x,-1,1) ,romberg(signum(sin(x))*sin(x),x,-1,1)]; while replacing sin(x) by x or x^3 all gave wrong answer of 0. But plain integrate(...) or absint(...) all gave wrong answer but ~~~ [load(fourie),absint(signum(sin(x))*sin(x),x,-1,1) ,"ck"=ev(2*integrate(sin(x),x,0,1),numer) ] ; ~~~ Did not give answer 0 but gave a positive wrong answer and can one explain this; ~~~ [ev(absint(sqrt(sin(x)^2),x,-1,0),numer),ev(absint(sqrt(sin(x)^2),x,0,1) ,numer)]; ~~~ both gave different positive answers with of course the 2nd one correct. --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Wed Aug 21, 2024 11:24 AM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: Barton W. <wil...@us...> - 2024-08-21 11:24:57
|
Using Maxima 5.47.0, these bugs seem to be fixed: ~~~ (%i1) load(fourie)$ (%i2) integrate(abs(x),x,-1,1); (%o2) 1 (%i3) absint(abs(x),x,-1,1); (%o3) 1 (%i4) romberg(abs(x),x,-11,11); (%o4) 121.0 (%i5) romberg(sqrt(x^2),x,-11,11); (%o5) 121.0 (%i6) romberg(signum(x)*x,x,-11,11); (%o6) 121.0 ~~~ Could you try installing the latest Maxima and trying again? --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Wed Aug 21, 2024 03:31 AM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: dan h. <zm...@us...> - 2024-08-21 03:31:20
|
But now this is pretty bad. All the following also give 0.0 ~~~ ~~~ [romberg(abs(x),x,-11,11),romberg(sqrt(x^2),x,-11,11),romberg(signum(x)*x,x,-11,11)]; ~~~ ~~~ --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Wed Aug 21, 2024 03:06 AM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: dan h. <zm...@us...> - 2024-08-21 03:06:57
|
~~~ [load(fourie),integrate(sqrt(x^2),x,-1,1),absint(sqrt(x^2),x,-1,1),integrate(signum(x)*x,x,-1,1),absint(signum(x)*x,x,-1,1) ] ; ~~~ All give 0,no comments nor anything etc. --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Tue Aug 20, 2024 08:42 PM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-21 01:57:59
|
- **status**: open --> closed - **Comment**: Fixed by commit d1f6d5. --- **[bugs:#4349] user-defined rules apply to expressions with square brackets as well as parentheses** **Status:** closed **Group:** None **Labels:** rules tellsimp tellsimpafter defrule defmatch **Created:** Sun Aug 18, 2024 02:06 AM UTC by Robert Dodier **Last Updated:** Sun Aug 18, 2024 03:02 AM UTC **Owner:** nobody Rules defined via `tellsimp`, `tellsimpafter`, `defrule`, and `defmatch` all ignore whether the operator is an ordinary function (with parentheses) or a subscripted operator (with square brackets). Same thing happens with all of them. In every case the generated code is probably just checking whether the expression caar is equal to required one, and not looking at the presence or absence of the ARRAY flag in the expression car. The rule should require the ARRAY flag to be the same as in the pattern for the rule (either ARRAY is present in both pattern and exemplar or absent in both). Here's an example with `tellsimpafter`. Same bug present in `tellsimp`, `defrule`, and `defmatch` as well. ``` (%i4) matchdeclare(aa, integerp); (%o4) done (%i5) tellsimpafter(foo(aa), 1 - aa); (%o5) [foorule1, false] (%i6) foo(10); (%o6) - 9 (%i7) foo[10]; (%o7) - 9 ``` --- Sent from sourceforge.net because max...@li... 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 M. <mac...@us...> - 2024-08-20 20:42:12
|
- **labels**: --> integrate - **Comment**: Agreed, **integrate** should return at the very least return a noun-form in this case. --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Labels:** integrate **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Tue Aug 20, 2024 07:51 PM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: dan h. <zm...@us...> - 2024-08-20 19:51:19
|
ok after i did load(fourie) the absint(...) Did do the function but again wrong answer of 0 and no comments or anything) --- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Tue Aug 20, 2024 07:38 PM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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: dan h. <zm...@us...> - 2024-08-20 19:38:18
|
--- **[bugs:#4351] absolute value, integral** **Status:** open **Group:** None **Created:** Tue Aug 20, 2024 07:38 PM UTC by dan hayes **Last Updated:** Tue Aug 20, 2024 07:38 PM UTC **Owner:** nobody "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") Not positively sure this is a bug but anyway: ~~~ ~~~ [integrate(abs(x),x,-1,1),absint(abs(x),x,-1,1)]; ~~~ ~~~ For the first one it seems to ignore the abs, absolute value prefix and gives wrong answer of 0, no warnings,no comments etc. Second one absint in my help file is supposed to be a function . Doesn't say i have to load fourier or anything like that but it just gives a repeat of the exact words of the function absint(...) without doing any evaluation. Am i missing something or do later versions of maxima do the same thing ? --- Sent from sourceforge.net because max...@li... 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 M. <mac...@us...> - 2024-08-19 22:23:29
|
- **Comment**: The original sin seems to be with **arrayinfo**, which has weird semantics (because array evaluation has weird semantics) ~~~ ar:make_array(fixnum,2) => {Lisp Array: #(0 0)} arl:[ar] => [{Lisp Array: #(0 0)}] arrayinfo(ar) => [declared, 1, [1]] arrayinfo(arl[1]) => ERROR arrayinfo: arl[1] is not an array. ev(arrayinfo(i),i=arl[1]) => [declared, 1, [1]] block([i:arl[1]], arrayinfo(i)) => [declared, 1, [1]] arrayinfo(''ar) => ERROR arrayinfo: {Lisp Array: #(0 0)} is not an array. ~~~ That said, there is *also* something wrong with **translate** of **makelist**: ~~~ test():= block([ar,arl,i], [ ar:make_array('fixnum,2), arl:[ar], errcatch(arrayinfo(ar)), errcatch(arrayinfo(arl[1])), /* ERROR */ errcatch(ev(arrayinfo(i),i=arl[1])), errcatch(apply('arrayinfo,[ar])), /* ERROR */ errcatch(block([i:ar], arrayinfo(i))), errcatch(block([i:arl[1]], arrayinfo(i))), errcatch(makelist(arrayinfo(i),i,arl)) /* ERROR when translated only */ ] )$ test(); arrayinfo: arl[1] is not an array. arrayinfo: "{Lisp Array: #(0 0)}" is not an array. => ["{Lisp Array: #(0 0)}",["{Lisp Array: #(0 0)}"],[[declared,1,[1]]],[],[[declared,1,[1]]],[],[[declared,1,[1]]],[[declared,1,[1]]],[[[declared,1,[1]]]]] translate(test)$ test(); arrayinfo: arl[1] is not an array. arrayinfo: "{Lisp Array: #(0 0)}" is not an array. arrayinfo: i is not an array. => ["{Lisp Array: #(0 0)}",["{Lisp Array: #(0 0)}"],[[declared,1,[1]]],[],[[declared,1,[1]]],[],[[declared,1,[1]]],[[declared,1,[1]]],[]] ~~~ Maxima 5.47.0 SBCL 2.4.7 MacOS --- **[bugs:#3296] makelist differs slightly between compiler and interpreter.** **Status:** open **Group:** None **Labels:** compiler arrayinfo makelist **Created:** Mon Mar 20, 2017 06:56 AM UTC by Gunter Königsmann **Last Updated:** Mon Aug 19, 2024 08:59 PM UTC **Owner:** nobody Another response to the call for differences between compiler and interpreter in maxima. First define a list of arrays: ~~~ (%i208) arr:makelist(make_array(flonum,10),i,1,10); (arr) [Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\]] ~~~ Then try to map a function on each individual array element using makelist - which fails if this is done by compiled code: ~~~ (%i211) test(x):=makelist(arrayinfo(i),i,x); compile(test); test(arr); (%o209) test(x):=makelist(arrayinfo(i),i,x) (%o210) [test] arrayinfo: i is not an array. -- an error. To debug this try: debugmode(true); ~~~ Using interpreted code this won't fail, though: ~~~ (%i213) test(x):=makelist(arrayinfo(i),i,x); test(arr); (%o212) test(x):=makelist(arrayinfo(i),i,x) (%o213) [[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]]] (%i207) build_info(); (%o207) build_info(version="branch_5_39_base_137_g3987a2f",timestamp="2017-03-10 08:59:49",host="x86_64-pc-linux-gnu",lisp_name="CLISP",lisp_version="2.49 (2010-07-07) (built 3605577779) (memory 3698121592)") ~~~ --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-19 21:00:05
|
See also SF bug #4350 for a possibly-related bug. --- **[bugs:#3296] makelist differs slightly between compiler and interpreter.** **Status:** open **Group:** None **Labels:** compiler arrayinfo makelist **Created:** Mon Mar 20, 2017 06:56 AM UTC by Gunter Königsmann **Last Updated:** Mon Aug 19, 2024 08:48 PM UTC **Owner:** nobody Another response to the call for differences between compiler and interpreter in maxima. First define a list of arrays: ~~~ (%i208) arr:makelist(make_array(flonum,10),i,1,10); (arr) [Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\]] ~~~ Then try to map a function on each individual array element using makelist - which fails if this is done by compiled code: ~~~ (%i211) test(x):=makelist(arrayinfo(i),i,x); compile(test); test(arr); (%o209) test(x):=makelist(arrayinfo(i),i,x) (%o210) [test] arrayinfo: i is not an array. -- an error. To debug this try: debugmode(true); ~~~ Using interpreted code this won't fail, though: ~~~ (%i213) test(x):=makelist(arrayinfo(i),i,x); test(arr); (%o212) test(x):=makelist(arrayinfo(i),i,x) (%o213) [[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]]] (%i207) build_info(); (%o207) build_info(version="branch_5_39_base_137_g3987a2f",timestamp="2017-03-10 08:59:49",host="x86_64-pc-linux-gnu",lisp_name="CLISP",lisp_version="2.49 (2010-07-07) (built 3605577779) (memory 3698121592)") ~~~ --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-19 20:59:16
|
--- **[bugs:#4350] arrayinfo complains "not an array" when supplied a Lisp array or hash table** **Status:** open **Group:** None **Labels:** arrayinfo **Created:** Mon Aug 19, 2024 08:59 PM UTC by Robert Dodier **Last Updated:** Mon Aug 19, 2024 08:59 PM UTC **Owner:** nobody `arrayinfo` complains "not an array" when its argument is a Lisp array or Lisp hash table. However, when its argument is a symbol which has the Lisp array or hash table bound to it, then `arrayinfo` knows what to do. ``` (%i33) aa: make_array (any, 4); (%o33) {Lisp Array: #(NIL NIL NIL NIL)} (%i34) arrayinfo (aa); (%o34) [declared, 1, [3]] (%i35) apply (arrayinfo, [aa]); arrayinfo: {Lisp Array: #(NIL NIL NIL NIL)} is not an array. -- an error. To debug this try: debugmode(true); (%i36) hh: make_array (hashed, 1); (%o36) {Lisp Array: #<HASH-TABLE :TEST EQUAL :COUNT 1 {1004C322F3}>} (%i37) hh["foo"]: "bar"; (%o37) bar (%i38) arrayinfo (hh); (%o38) [hash_table, 1, foo] (%i39) apply (arrayinfo, [hh]); arrayinfo: {Lisp Array: #<HASH-TABLE :TEST EQUAL :COUNT 2 {1004C322F3}>} is not an array. -- an error. To debug this try: debugmode(true); ``` I bumped into this problem in the following context. I have a list of arrays, and I'm trying to call `arrayinfo` on each one. %o41 fails because getting the value of `L[1]` would require calling MEVAL, which, I'm happy to report, `arrayinfo` doesn't call. OK so far. The standard workaround is to call `apply(foo, [...])` for an argument-quoting function `foo`. But that doesn't work either, for `arrayinfo`: ``` (%i40) L: [aa]; (%o40) [{Lisp Array: #(NIL NIL NIL NIL)}] (%i41) arrayinfo (L[1]); arrayinfo: L is not an array. 1 -- an error. To debug this try: debugmode(true); (%i42) apply (arrayinfo, [L[1]]); arrayinfo: {Lisp Array: #(NIL NIL NIL NIL)} is not an array. -- an error. To debug this try: debugmode(true); ``` The fix is just to handle Lisp array and hash table values in addition to symbols. --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-19 20:48:34
|
- **labels**: compiler --> compiler, arrayinfo, makelist - **Comment**: `arrayinfo` seems to require its argument to be a symbol; the symbol's value causes an error about "not an array". Not sure whether that is the whole problem here, or if there is something about `makelist` in translated code, or what. --- **[bugs:#3296] makelist differs slightly between compiler and interpreter.** **Status:** open **Group:** None **Labels:** compiler arrayinfo makelist **Created:** Mon Mar 20, 2017 06:56 AM UTC by Gunter Königsmann **Last Updated:** Mon Mar 20, 2017 06:56 AM UTC **Owner:** nobody Another response to the call for differences between compiler and interpreter in maxima. First define a list of arrays: ~~~ (%i208) arr:makelist(make_array(flonum,10),i,1,10); (arr) [Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\],Lisp array \[10\]] ~~~ Then try to map a function on each individual array element using makelist - which fails if this is done by compiled code: ~~~ (%i211) test(x):=makelist(arrayinfo(i),i,x); compile(test); test(arr); (%o209) test(x):=makelist(arrayinfo(i),i,x) (%o210) [test] arrayinfo: i is not an array. -- an error. To debug this try: debugmode(true); ~~~ Using interpreted code this won't fail, though: ~~~ (%i213) test(x):=makelist(arrayinfo(i),i,x); test(arr); (%o212) test(x):=makelist(arrayinfo(i),i,x) (%o213) [[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]],[declared,1,[9]]] (%i207) build_info(); (%o207) build_info(version="branch_5_39_base_137_g3987a2f",timestamp="2017-03-10 08:59:49",host="x86_64-pc-linux-gnu",lisp_name="CLISP",lisp_version="2.49 (2010-07-07) (built 3605577779) (memory 3698121592)") ~~~ --- Sent from sourceforge.net because max...@li... 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: dan h. <zm...@us...> - 2024-08-18 18:26:34
|
yea apparently it did work but one would never know from looking in the help file for rootsconmode as that option 'super' is not listed for that at least in my help file. i wonder if it is listed in anyone elses more recent versions of maxima help ? On Saturday, August 17, 2024 at 11:51:11 AM CDT, Stavros Macrakis via Maxima-bugs <max...@li...> wrote: rootsconmode does not support super. I was suggesting that it should use that to handle this case. [bugs:#4347] maxima would not simplify nor rationalize the simplest expression Status: open Group: None Created: Thu Aug 15, 2024 02:40 PM UTC by dan hayes Last Updated: Sat Aug 17, 2024 02:09 PM UTC Owner: nobody build_info() or bug_report() "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") (assume_pos:true, file_output_append:true, ratprint:false, showtime:true,load(simplify_sum) , simpsum:true, load("lrats"),letrat:true,ratfac:true,algebraic:true,fpprintprec:4); (tt:solve([x+y-u,x*y-v],[x,y]),ta:rationalize(factor(lratsubst(part(tt,1),(1-x)^a(1+x)^b(1-y)^a(1+y)^b))) ,tb:factor(ratsimp( ((2-u)^2-u^2+4v)^a((2+u)^2-u^2+4v)^b)*(2^(-2))^(b+a)),disp([ta,tb]) ); ok i know from help it says rationalize is for floating numbers but i tried everything else so i tried that and maxima would not simplify the expression ta so i had to manually do it myself to get rid of irrational expression. Clearly maxima should have been able to do that itself and end up with my expression tb without me having to do it myself. Why did it not? and if not why is there no way to implement that. There certainly should be. Sent from sourceforge.net because max...@li... 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. _______________________________________________ Maxima-bugs mailing list Max...@li... https://lists.sourceforge.net/lists/listinfo/maxima-bugs --- **[bugs:#4347] maxima would not simplify nor rationalize the simplest expression** **Status:** open **Group:** None **Created:** Thu Aug 15, 2024 02:40 PM UTC by dan hayes **Last Updated:** Sun Aug 18, 2024 12:24 AM UTC **Owner:** nobody build_info() or bug_report() "branch_5_44_base_231_g5c411f69f",timestamp="2021-01-12 23:51:42",host="x86_64-w64-mingw32",lisp_name="SBCL",lisp_version="2.0.0",maxima_userdir="C:/Users/zmth1/maxima",maxima_tempdir="C:/Users/zmth1/AppData/Local/Temp",maxima_objdir="C:/Users/zmth1/maxima/binary/branch_5_44_base_231_g5c411f69f/sbcl/2_0_0",maxima_frontend="wxMaxima",maxima_frontend_version="20.12.2-DevelopmentSnapshot_MSW_OpenMP201511+Locks") (assume_pos:true, file_output_append:true, ratprint:false, showtime:true,load(simplify_sum) , simpsum:true, load("lrats"),letrat:true,ratfac:true,algebraic:true,fpprintprec:4); (tt:solve([x+y-u,x*y-v],[x,y]),ta:rationalize(factor(lratsubst(part(tt,1),(1-x)^a*(1+x)^b*(1-y)^a*(1+y)^b))) ,tb:factor(ratsimp( ((2-u)^2-u^2+4*v)^a*((2+u)^2-u^2+4*v)^b)*(2^(-2))^(b+a)),disp([ta,tb]) ); ok i know from help it says rationalize is for floating numbers but i tried everything else so i tried that and maxima would not simplify the expression ta so i had to manually do it myself to get rid of irrational expression. Clearly maxima should have been able to do that itself and end up with my expression tb without me having to do it myself. Why did it not? and if not why is there no way to implement that. There certainly should be. --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-18 03:02:26
|
- **summary**: user-defined rules apply to expressions with square brackets as parentheses --> user-defined rules apply to expressions with square brackets as well as parentheses --- **[bugs:#4349] user-defined rules apply to expressions with square brackets as well as parentheses** **Status:** open **Group:** None **Labels:** rules tellsimp tellsimpafter defrule defmatch **Created:** Sun Aug 18, 2024 02:06 AM UTC by Robert Dodier **Last Updated:** Sun Aug 18, 2024 02:06 AM UTC **Owner:** nobody Rules defined via `tellsimp`, `tellsimpafter`, `defrule`, and `defmatch` all ignore whether the operator is an ordinary function (with parentheses) or a subscripted operator (with square brackets). Same thing happens with all of them. In every case the generated code is probably just checking whether the expression caar is equal to required one, and not looking at the presence or absence of the ARRAY flag in the expression car. The rule should require the ARRAY flag to be the same as in the pattern for the rule (either ARRAY is present in both pattern and exemplar or absent in both). Here's an example with `tellsimpafter`. Same bug present in `tellsimp`, `defrule`, and `defmatch` as well. ``` (%i4) matchdeclare(aa, integerp); (%o4) done (%i5) tellsimpafter(foo(aa), 1 - aa); (%o5) [foorule1, false] (%i6) foo(10); (%o6) - 9 (%i7) foo[10]; (%o7) - 9 ``` --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-18 02:06:15
|
--- **[bugs:#4349] user-defined rules apply to expressions with square brackets as parentheses** **Status:** open **Group:** None **Labels:** rules tellsimp tellsimpafter defrule defmatch **Created:** Sun Aug 18, 2024 02:06 AM UTC by Robert Dodier **Last Updated:** Sun Aug 18, 2024 02:06 AM UTC **Owner:** nobody Rules defined via `tellsimp`, `tellsimpafter`, `defrule`, and `defmatch` all ignore whether the operator is an ordinary function (with parentheses) or a subscripted operator (with square brackets). Same thing happens with all of them. In every case the generated code is probably just checking whether the expression caar is equal to required one, and not looking at the presence or absence of the ARRAY flag in the expression car. The rule should require the ARRAY flag to be the same as in the pattern for the rule (either ARRAY is present in both pattern and exemplar or absent in both). Here's an example with `tellsimpafter`. Same bug present in `tellsimp`, `defrule`, and `defmatch` as well. ``` (%i4) matchdeclare(aa, integerp); (%o4) done (%i5) tellsimpafter(foo(aa), 1 - aa); (%o5) [foorule1, false] (%i6) foo(10); (%o6) - 9 (%i7) foo[10]; (%o7) - 9 ``` --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-18 01:47:21
|
- **labels**: testsuite, rtest16 --> testsuite, rtest16, integrate, defint - **summary**: rtest16 #538 --> rtest16 #538 (definite integral) --- **[bugs:#4339] rtest16 #538 (definite integral)** **Status:** open **Group:** None **Labels:** testsuite rtest16 integrate defint **Created:** Mon Aug 05, 2024 12:44 PM UTC by Barton Willis **Last Updated:** Sun Aug 18, 2024 01:45 AM UTC **Owner:** nobody Test 538 in rtest16 is ``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ (declare(m2726, even), block([tmp: integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4)], sign (subst(m2726 = 4, tmp)))); pos$ ``` First, I don't think that `m2767` should need to be declared to be even; and second, the integral diverges when `m2726 <= -1`. I think the test should be modified to ``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ block([tmp, m2726], assume(m2726+1 > 0), tmp : integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4), tmp : sign(subst(m2726 = 4, tmp)), forget(m2726+1 > 0), tmp); ``` Now Maxima 5.47 passes this test OK. But if `domain` is complex, the test will ask if `m2726` is an integer. Actually, with Maxima 5.47, we have ``` (%i9) integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4); (%o9) (-(2^(m2726/2-1/2)*(-1)^m2726*gamma_incomplete((m2726+1)/2,1/32)) -2^(m2726/2-1/2)*gamma_incomplete((m2726+1)/2,1/32) +2^((m2726-1)/2)*(-1)^m2726*gamma_incomplete((m2726+1)/2,0) +gamma_incomplete(m2726/2+1/2,0)*2^(m2726/2-1/2)) /(sqrt(2)*sqrt(%pi)) ``` The definite integral is OK, I think, but it's undefined when `m2726 <= -1`. --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-18 01:45:49
|
I think the idea with the `declare(..., even)` is that the integral is zero when the exponent is an odd positive integer. It wouldn't hurt to test both odd and even, I guess. The test for the result being positive seems like a pretty weak test -- seems much more informative to test for a specific result. Aside from that, we could test with a variable declared even or odd, and also with specific even or odd values -- those should come out the same, of course. I see from looking at the original bug report (https://sourceforge.net/p/maxima/bugs/2726/) that the declaration is said to be a workaround -- maybe we want to test with and, especially, without the declaration. --- **[bugs:#4339] rtest16 #538** **Status:** open **Group:** None **Labels:** testsuite rtest16 **Created:** Mon Aug 05, 2024 12:44 PM UTC by Barton Willis **Last Updated:** Sun Aug 18, 2024 01:38 AM UTC **Owner:** nobody Test 538 in rtest16 is ``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ (declare(m2726, even), block([tmp: integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4)], sign (subst(m2726 = 4, tmp)))); pos$ ``` First, I don't think that `m2767` should need to be declared to be even; and second, the integral diverges when `m2726 <= -1`. I think the test should be modified to ``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ block([tmp, m2726], assume(m2726+1 > 0), tmp : integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4), tmp : sign(subst(m2726 = 4, tmp)), forget(m2726+1 > 0), tmp); ``` Now Maxima 5.47 passes this test OK. But if `domain` is complex, the test will ask if `m2726` is an integer. Actually, with Maxima 5.47, we have ``` (%i9) integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4); (%o9) (-(2^(m2726/2-1/2)*(-1)^m2726*gamma_incomplete((m2726+1)/2,1/32)) -2^(m2726/2-1/2)*gamma_incomplete((m2726+1)/2,1/32) +2^((m2726-1)/2)*(-1)^m2726*gamma_incomplete((m2726+1)/2,0) +gamma_incomplete(m2726/2+1/2,0)*2^(m2726/2-1/2)) /(sqrt(2)*sqrt(%pi)) ``` The definite integral is OK, I think, but it's undefined when `m2726 <= -1`. --- Sent from sourceforge.net because max...@li... 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 D. <rob...@us...> - 2024-08-18 01:38:39
|
- Description has changed: Diff: ~~~~ --- old +++ new @@ -1,16 +1,15 @@ Test 538 in rtest16 is - -~~~ +``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ (declare(m2726, even), block([tmp: integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4)], sign (subst(m2726 = 4, tmp)))); pos$ -~~~ +``` First, I don't think that `m2767` should need to be declared to be even; and second, the integral diverges when `m2726 <= -1`. I think the test should be modified to -~~~ +``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ block([tmp, m2726], assume(m2726+1 > 0), @@ -18,13 +17,13 @@ tmp : sign(subst(m2726 = 4, tmp)), forget(m2726+1 > 0), tmp); -~~~ +``` Now Maxima 5.47 passes this test OK. But if `domain` is complex, the test will ask if `m2726` is an integer. Actually, with Maxima 5.47, we have -~~~ +``` (%i9) integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4); (%o9) (-(2^(m2726/2-1/2)*(-1)^m2726*gamma_incomplete((m2726+1)/2,1/32)) @@ -35,6 +34,6 @@ +gamma_incomplete(m2726/2+1/2,0)*2^(m2726/2-1/2)) /(sqrt(2)*sqrt(%pi)) - ~~~ +``` The definite integral is OK, I think, but it's undefined when `m2726 <= -1`. ~~~~ --- **[bugs:#4339] rtest16 #538** **Status:** open **Group:** None **Labels:** testsuite rtest16 **Created:** Mon Aug 05, 2024 12:44 PM UTC by Barton Willis **Last Updated:** Mon Aug 05, 2024 12:44 PM UTC **Owner:** nobody Test 538 in rtest16 is ``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ (declare(m2726, even), block([tmp: integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4)], sign (subst(m2726 = 4, tmp)))); pos$ ``` First, I don't think that `m2767` should need to be declared to be even; and second, the integral diverges when `m2726 <= -1`. I think the test should be modified to ``` /* #2726: Integrate produces wrong answer for Gaussian Moments */ block([tmp, m2726], assume(m2726+1 > 0), tmp : integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4), tmp : sign(subst(m2726 = 4, tmp)), forget(m2726+1 > 0), tmp); ``` Now Maxima 5.47 passes this test OK. But if `domain` is complex, the test will ask if `m2726` is an integer. Actually, with Maxima 5.47, we have ``` (%i9) integrate(exp(-x^2/2)/sqrt(2*%pi) * x^m2726, x, -1/4, 1/4); (%o9) (-(2^(m2726/2-1/2)*(-1)^m2726*gamma_incomplete((m2726+1)/2,1/32)) -2^(m2726/2-1/2)*gamma_incomplete((m2726+1)/2,1/32) +2^((m2726-1)/2)*(-1)^m2726*gamma_incomplete((m2726+1)/2,0) +gamma_incomplete(m2726/2+1/2,0)*2^(m2726/2-1/2)) /(sqrt(2)*sqrt(%pi)) ``` The definite integral is OK, I think, but it's undefined when `m2726 <= -1`. --- Sent from sourceforge.net because max...@li... 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. |