Screenshot instructions:
Windows
Mac
Red Hat Linux
Ubuntu
Click URL instructions:
Rightclick on ad, choose "Copy Link", then paste here →
(This may not be possible with some types of ads)
You can subscribe to this list here.
2014 
_{Jan}

_{Feb}
(232) 
_{Mar}
(323) 
_{Apr}
(383) 
_{May}
(359) 
_{Jun}
(435) 
_{Jul}
(252) 
_{Aug}
(172) 
_{Sep}
(265) 
_{Oct}
(263) 
_{Nov}
(350) 
_{Dec}
(359) 

2015 
_{Jan}
(267) 
_{Feb}
(220) 
_{Mar}
(311) 
_{Apr}
(269) 
_{May}
(388) 
_{Jun}
(403) 
_{Jul}
(172) 
_{Aug}
(399) 
_{Sep}
(364) 
_{Oct}
(269) 
_{Nov}
(357) 
_{Dec}
(468) 
2016 
_{Jan}
(618) 
_{Feb}
(592) 
_{Mar}
(625) 
_{Apr}
(516) 
_{May}
(375) 
_{Jun}
(155) 
_{Jul}
(346) 
_{Aug}
(262) 
_{Sep}
(346) 
_{Oct}
(291) 
_{Nov}
(333) 
_{Dec}
(335) 
2017 
_{Jan}
(436) 
_{Feb}
(460) 
_{Mar}
(370) 
_{Apr}
(189) 
_{May}
(252) 
_{Jun}
(272) 
_{Jul}
(286) 
_{Aug}
(293) 
_{Sep}
(303) 
_{Oct}
(331) 
_{Nov}
(346) 
_{Dec}
(273) 
2018 
_{Jan}
(295) 
_{Feb}
(217) 
_{Mar}

_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 



1
(3) 
2
(6) 
3
(32) 
4
(23) 
5
(16) 
6
(16) 
7
(17) 
8
(11) 
9
(9) 
10
(14) 
11
(8) 
12
(4) 
13
(3) 
14
(10) 
15
(4) 
16
(6) 
17
(7) 
18
(10) 
19
(7) 
20
(14) 
21
(12) 
22
(26) 
23
(31) 
24
(28) 
25
(23) 
26
(12) 
27
(12) 
28
(11) 
29
(4) 
30
(4) 



From: Panagiotis Papasotiriou <p.j.papasot@gm...>  20140423 22:43:04

Ted, for some reason I have some issues while loading your webpage (it used to work perfectly in the past.) It takes quite some time to load, and sometimes it even fails for no apparent reason (my internet connection works as expected otherwise, so it can't be a problem on my side.) Also, are you sure cp2.pdf is up? I can read cp1, project1, but not cp2: I get a 404 error, "The requested URL /~woollett/cp2.pdf was not found on this server." 20140423 23:04 GMT+03:00 Edwin Woollett <woollett@...>: > > Chapter 2, Initial Value Problems (Computational Physics Using > R and Maxima) has been posted to my web page. > > From that web page: >  > Chapter 2, Initial Value Problems > > A variety of initial value problems are solved using both > R and Maxima. Homemade codes for the Euler method and > the classical fourth order RungeKutta method are > demonstrated with numerous examples. The standard > Maxima methods rk and rkf45 are discussed, and the > R deSolve package ode method is extensively discussed > with many examples. Finally, methods for solving initial > value problems in which the derivatives are discontinuous > and in which the dependent variables are discontinous are discussed. > > cp2.pdf : Ch. 2, April 23, 2014, Maxima 5.31.0, R 3.0.2, 71 pages > cp2.tex : Ch. 2 Latex code file > k2util.mac : Maxima utility functions for Ch. 2, April 23, 2014 > myode.mac : Maxima code for Ch. 2 , April 23, 2014, Maxima 5.31 > myode.R : R code for Ch. 2, April 23, 2014, R 3.0.2 > myrkf45.mac : Maxima code for Ch. 2 , April 23, 2014, Maxima 5.31 > > Chapter 2 Topics > > 1 The Euler Method, Truncation Errors, Roundoff Errors, Instability > 2 4th Order RungeKutta Code for R and Maxima > 3 The Maxima Methods rk and rkf45 > 4 The R package deSolve and Use of ode > 5 Solving O.D.E.'s with External Forcing Data > 6 Integrating O.D.E.'s with Discontinuous Derivatives > 7 Integrating O.D.E.'s with Discontinuous Dependent Variables Using R > 1 Events Specified with a data.frame > 2 Intravenous Drug Injection Model > 3 Using an Event Function at Specified Times > 4 Using an Event Function when a Root Condition is Satisfied > 5 Use of a Switching Parameter as a Dependent Variable >  > > Ted Woollett > > http://www.csulb.edu/~woollett/ > > > > >  > Start Your Social Network Today  Download eXo Platform > Build your Enterprise Intranet with eXo Platform Software > Java Based Open Source Intranet  Social, Extensible, Cloud Ready > Get Started Now And Turn Your Intranet Into A Collaboration Platform > http://p.sf.net/sfu/ExoPlatform > _______________________________________________ > Maximadiscuss mailing list > Maximadiscuss@... > https://lists.sourceforge.net/lists/listinfo/maximadiscuss > 
From: Stavros Macrakis <macrakis@al...>  20140423 21:41:50

Well, let's list the ways they get into trouble and see what we can do about them. Here are some off the top of my head: * not all of Maxima allows x[1] where it allows a variable ** for x[1] thru 3 => syntax error ** function, block, and lambda formals (do we really want x[i] and x[2] as formals?) ** ... I'm sure there are more * the relationship between x[i] and i, x[1] and x[i], or x[i] and x[j] ** diff(x[i],x[j]) => 0 instead of kron_delta(i,j) ** diff(x[y],y) =>0 ... should be a nounform (or an error) ** diff(sum(x[i],i,0,n),x[j]) => 0 ... should be (if j>=0 and j<=n then 1 else 0) ** etc. ... these should all be fixable * you need to be sure to use mapatom instead of atom in many cases * what else? On Wed, Apr 23, 2014 at 5:16 PM, Barton Willis <willisb@...> wrote: > > > > I > agree with Robert that we shouldn't be warning users against using > features we think are working correctly. If we think they aren't working > correctly, > > we should fix or remove them. > And I'm not sure what "uncritical use" is anyway. > > Users sometimes use x[1] only because they like the way it prints. Then > they get into trouble. > > bw > 
From: Barton Willis <willisb@un...>  20140423 21:17:05

?? > I ? agree with Robert that we shouldn't be warning users against using features we think are working correctly. If we think they aren't working correctly, > we should fix or remove them. And I'm not sure what "uncritical use" is anyway. Users sometimes use x[1] only because they like the way it prints. Then they get into trouble. bw 
From: Stavros Macrakis <macrakis@al...>  20140423 20:55:02

On Tue, Apr 22, 2014 at 2:59 PM, Robert Dodier <robert.dodier@...>wrote: > On 20140421, Stavros Macrakis <macrakis@...> wrote: > > > See Technical Issues of Separation in Function Cells and Value > > Cells<http://www.nhplace.com/kent/Papers/TechnicalIssues.html>; (often > > called "Lisp 1 vs. Lisp 2") by R.P. Gabriel for a discussion of the > issues > > in Lisp. Maxima, with its "variables without values represent themselves" > > policy, adds another twist to the story. > > > > In short, this is not a "bug", but a design choice, and whether you are > > surprised by it depends on what you're familiar with.... > > It doesn't really have much to do with Lisp1 vs Lisp2; more > precisely, there was no clear choice of one which causes heartburn > for partisans of the other camp. Rather, the problem is that Maxima > doesn't clearly prefer one or the other, which makes it difficult > to predict what's going to happen, hence the > > surprise. > Aside from allowing functions in both value and function cells (in Lisp > terms), Maxima resolves the ambiguity of a symbol which has both in > favor of the function cell, which has the unfortunate effect of favoring > a global declaration to a local one; such action at a distance is, > I guess, what makes the suprise unpleasant. > This is the standard behavior in Lisp2 systems. > > Personally, I would support the syntax (f)(x) to indicate that the 'f' is > > to be evaluated. > > Some hacking suggests that to ensure that f is evaluated, it's enough to > construct an MPROGN expression containing f. At present the parser > elides MPROGN when it has just one argument (look for "Optimize out > MPROGN" in src/nparse.lisp). To implement the policy you suggest (not > convinced yet), I'm guessing it's enough to cut out that optimization. > Yes, I think so. best > > Robert Dodier > > > >  > Start Your Social Network Today  Download eXo Platform > Build your Enterprise Intranet with eXo Platform Software > Java Based Open Source Intranet  Social, Extensible, Cloud Ready > Get Started Now And Turn Your Intranet Into A Collaboration Platform > http://p.sf.net/sfu/ExoPlatform > _______________________________________________ > Maximadiscuss mailing list > Maximadiscuss@... > https://lists.sourceforge.net/lists/listinfo/maximadiscuss > 
From: Stavros Macrakis <macrakis@al...>  20140423 20:52:47

Re Our user documentation should have a stern warning against uncritical use of undeclared arrays. I agree with Robert that we shouldn't be warning users against using features we think are working correctly. If we think they aren't working correctly, we should fix or remove them. And I'm not sure what "uncritical use" is anyway. And I think undeclared arrays, a.k.a. hasharrays, are a powerful and valuable feature of Maxima, very consistent with Maxima's general style of assigning generic meaning to things without preliminaries such as declarations. I would not like to ask users to do something like array(x,hash) before using x[i]. I agree that it is strange that the x in x[i] is interpreted as a formal array name if there are any other x[j]'s in existence, and as a variable if there are not. This is some sort of weird Lisp3 (where the 3rd value of a variable is its array value). I can think of a couple of consistent ways of handling this: 1) x is *always* evaluated as a variable before being tested for array reference (a Lisp1 sort of approach); 2) x is *never* evaluated as a variable before being tested for array reference (a Lisp2 sort of approach). Both seem reasonable, but (2) is more consistent with Maxima's general Lisp2 style. If you wanted it to be evaluated, you'd use the same solution I proposed for functions, namely (x)[i]. s On Sun, Apr 20, 2014 at 6:22 PM, Barton Willis <willisb@...> wrote: > A simple example of the somewhat unfortunate behavior: > > (%i1) f(x,i) := x[i]$ > > OK: > > (%i2) f(z,8); > (%o2) z[8] > > (%i3) x[%pi] : 42$ > > unexpected (?): > > (%i4) f(z,8); > (%o4) x[8] > > A cure is to use localit " removes any properties before evaluating > other expressions, and > restores any saved properties on exit from the block or other ..." > > (%i5) g(x,i) := (local(x), x[i])$ > (%i6) g(z,8); > (%o6) z[8] > > Our user documentation should have a stern warning against uncritical use > of undeclared arrays. > If you use them because you like the they display, consider changing your > preferences. One > not so great thing about undeclared arrays is that some (all?) > declarations don't work: > > (%i8) declare(xx[1],integer); > declare: improper argument: xx[1] > > _______________________________________ > > >  > Learn Graph Databases  Download FREE O'Reilly Book > "Graph Databases" is the definitive new guide to graph databases and their > applications. Written by three acclaimed leaders in the field, > this first edition is now available. Download your free book today! > http://p.sf.net/sfu/NeoTech > _______________________________________________ > Maximadiscuss mailing list > Maximadiscuss@... > https://lists.sourceforge.net/lists/listinfo/maximadiscuss > 
From: Robert Dodier <robert.dodier@gm...>  20140423 20:25:06

On 20140423, Helmut Jarausch <jarausch@...> wrote: > thanks again for this precious blex. > Unfortunately, it disables source debugging. > Can this be fixed or is it some inevitable collateral damage? Looks like Maxima's substitution functions, which are called from blex, clobber the debugging data. I'm pretty sure that's fixable, but I will have to look at it more carefully. Sorry for the bother, and thanks for the encouragement. Robert Dodier 
From: Edwin Woollett <woollett@ch...>  20140423 20:03:56

Chapter 2, Initial Value Problems (Computational Physics Using R and Maxima) has been posted to my web page. From that web page:  Chapter 2, Initial Value Problems A variety of initial value problems are solved using both R and Maxima. Homemade codes for the Euler method and the classical fourth order RungeKutta method are demonstrated with numerous examples. The standard Maxima methods rk and rkf45 are discussed, and the R deSolve package ode method is extensively discussed with many examples. Finally, methods for solving initial value problems in which the derivatives are discontinuous and in which the dependent variables are discontinous are discussed. cp2.pdf : Ch. 2, April 23, 2014, Maxima 5.31.0, R 3.0.2, 71 pages cp2.tex : Ch. 2 Latex code file k2util.mac : Maxima utility functions for Ch. 2, April 23, 2014 myode.mac : Maxima code for Ch. 2 , April 23, 2014, Maxima 5.31 myode.R : R code for Ch. 2, April 23, 2014, R 3.0.2 myrkf45.mac : Maxima code for Ch. 2 , April 23, 2014, Maxima 5.31 Chapter 2 Topics 1 The Euler Method, Truncation Errors, Roundoff Errors, Instability 2 4th Order RungeKutta Code for R and Maxima 3 The Maxima Methods rk and rkf45 4 The R package deSolve and Use of ode 5 Solving O.D.E.'s with External Forcing Data 6 Integrating O.D.E.'s with Discontinuous Derivatives 7 Integrating O.D.E.'s with Discontinuous Dependent Variables Using R 1 Events Specified with a data.frame 2 Intravenous Drug Injection Model 3 Using an Event Function at Specified Times 4 Using an Event Function when a Root Condition is Satisfied 5 Use of a Switching Parameter as a Dependent Variable  Ted Woollett http://www.csulb.edu/~woollett/ 
From: Stavros Macrakis <macrakis@al...>  20140423 19:48:07

You can calculate the correct rounded result using bfloats. The easiest way to do that is a:50000000.0b0$ b:60000000.0b0$ round ((6^a2^b3^b)/(6^(a7)4^b5^b)); This gives 0. You can tell Maxima how many digits to carry in the calculation using fpprec:1000, for example, to calculate to 1000 digits of precision. This will get slow if you want, say, 10000 digits of precision, but that's probably not necessary. The result of 'round' may also turn out to be a very large integer, which may be too big to manipulate (e.g. you might get an out of memory error) and for that matter most of its digits may be meaningless. Why do you want to do these calculations? Another way to put it: what are you going to do with the results of these calculations? If you're looking at the behavior of these sorts of expressions as a and b get large, the 'taylor' function may help. For example, if you assume that b=6/5*a, then as a gets large, the expression has the value (3/5)^(6/5*a). s On Wed, Apr 23, 2014 at 3:09 PM, Scott Hughes <hughesscott@...>wrote: > An example... > > > > a: 50000000; > > b: 60000000; > > c: 40000000; > > round ((6^a2^b3^b)/(6^(a7)4^b5^b)); > > > > *From:* Gunter Königsmann [mailto:gunter@...] > *Sent:* April2314 2:58 PM > *To:* fateman@...; Scott Hughes; < > maximadiscuss@...> > *Subject:* Re: [Maximadiscuss] Maxima (lisp on Mac)  heap exhausted > > > > Another question: > How do you calculate the powers? I get a different error message, if I Try. > > > > On 23. April 2014 20:11:26 MESZ, Richard Fateman <fateman@...> > wrote: > > what do you plan to do with a number raised to the power 60,000,000? > > If you are computing it in order to then reduce it modulo some other > number, > then there is a better way to do it. > > > > On 4/23/2014 6:34 AM, Scott Hughes wrote: > > Firstly, I’m not a programmer and only recently started to use Maxima > (command line). I don’t know anything about lisp. I’m running Maxima on a > Mac but also on a PC through Windows. Both are the command line version. > > > > I only need to use Maxima for a few basic operations, primarily small > numbers raised to large (at least 8) exponents. Aside from that, I’m using > the round and asksign commands. > > > > I’m encountering “heap exhausted” errors once my exponents get beyond > about 60 000 000 and have looked in manual and forums for help. There is > reference to changing heap size but I didn’t understand the command or even > how to enter it. Can anyone help me with very explicit instructions ? I’m > not concerned with the length of time the calculations take, just don’t > want them to crash before they complete. > > > > > > > > > >  > > > Start Your Social Network Today  Download eXo Platform > Build your Enterprise Intranet with eXo Platform Software > Java Based Open Source Intranet  Social, Extensible, Cloud Ready > Get Started Now And Turn Your Intranet Into A Collaboration Platform > http://p.sf.net/sfu/ExoPlatform > >  > > > Maximadiscuss mailing list > Maximadiscuss@... > https://lists.sourceforge.net/lists/listinfo/maximadiscuss > > >  > Diese Nachricht wurde von meinem Mobiltelefon mit Kaiten Mail gesendet. > > >  > Start Your Social Network Today  Download eXo Platform > Build your Enterprise Intranet with eXo Platform Software > Java Based Open Source Intranet  Social, Extensible, Cloud Ready > Get Started Now And Turn Your Intranet Into A Collaboration Platform > http://p.sf.net/sfu/ExoPlatform > _______________________________________________ > Maximadiscuss mailing list > Maximadiscuss@... > https://lists.sourceforge.net/lists/listinfo/maximadiscuss > > 
From: Robert Dodier <robert.dodier@gm...>  20140423 19:38:50

On 20140423, Scott Hughes <hughesscott@...> wrote: > I'm encountering "heap exhausted" errors once my exponents get beyond about > 60000000 and have looked in manual and forums for help. "Heap exhausted" means that the underlying Lisp implementation has run out of memory trying to accomplish some operation. In this case, it may be trying to construct extremely large integers (just guessing). Different Lisp implementations manage memory differently. My experience is that GCL runs into outofmemory errors more often than other implementations. Unfortunately the Windows installers are built with GCL so that is the usual implementation. What Lisp implementation is your Mac using? build_info(); tells that. By the way, can you give an example that provokes the error so we can try it? best Robert Dodier 
From: Richard Fateman <fateman@be...>  20140423 19:26:49

On 4/23/2014 12:09 PM, Scott Hughes wrote: > > An example... > > a: 50000000; > > b: 60000000; > > c: 40000000; > > round ((6^a2^b3^b)/(6^(a7)4^b5^b)); > The result of the calculation, before rounding, is appx 1.815603675821045b3030638 this is easily computed using 2.0b0 etc instead of 2 etc. That is, it has over 3 million digits . Before computing and displaying them, it may make sense to ask.... what do you plan to do with them? > *From:*Gunter Königsmann [mailto:gunter@...] > *Sent:* April2314 2:58 PM > *To:* fateman@...; Scott Hughes; > <maximadiscuss@...> > *Subject:* Re: [Maximadiscuss] Maxima (lisp on Mac)  heap exhausted > > Another question: > How do you calculate the powers? I get a different error message, if I > Try. > > On 23. April 2014 20:11:26 MESZ, Richard Fateman <fateman@... > <mailto:fateman@...>> wrote: > > what do you plan to do with a number raised to the power 60,000,000? > > If you are computing it in order to then reduce it modulo some other > number, > then there is a better way to do it. > > > > On 4/23/2014 6:34 AM, Scott Hughes wrote: > > Firstly, I’m not a programmer and only recently started to use > Maxima (command line). I don’t know anything about lisp. I’m > running Maxima on a Mac but also on a PC through Windows. Both are > the command line version. > > I only need to use Maxima for a few basic operations, primarily > small numbers raised to large (at least 8) exponents. Aside from > that, I’m using the round and asksign commands. > > I’m encountering “heap exhausted” errors once my exponents get > beyond about 60 000 000 and have looked in manual and forums for > help. There is reference to changing heap size but I didn’t > understand the command or even how to enter it. Can anyone help me > with very explicit instructions ? I’m not concerned with the > length of time the calculations take, just don’t want them to > crash before they complete. > >  > Start Your Social Network Today  Download eXo Platform > Build your Enterprise Intranet with eXo Platform Software > Java Based Open Source Intranet  Social, Extensible, Cloud Ready > Get Started Now And Turn Your Intranet Into A Collaboration Platform > http://p.sf.net/sfu/ExoPlatform >  > Maximadiscuss mailing list > Maximadiscuss@... <mailto:Maximadiscuss@...> > https://lists.sourceforge.net/lists/listinfo/maximadiscuss > > >  > Diese Nachricht wurde von meinem Mobiltelefon mit Kaiten Mail gesendet. > 
From: Richard Fateman <fateman@be...>  20140423 19:14:38

On 4/23/2014 11:31 AM, Scott Hughes wrote: > > Thanks for answering. I am comparing the results to other similar > calculations. Unfortunately I need to be able to review both rounded > and full results. eg. (7^12345678+8^2345678) / (9^5678901 + 3^2334555) > You can get an approximate result by using bigfloats. If you write (7.0b0^12345678+8.0b0^2345678) / (9.0b0^5678901 + 3.0b0^2334555); The answer comes out 3.462179964245316b5014259 That is, 3.46 ... X 10^5014259 If you compute the numerator and denominator separately, you learn they have 10433308 and 5014259 decimal digits, respectively. So an exact answer will probably display with some 15 million characters. How do you plan to compare these characters with other "similar calculations"? You can easily calculate the last few digits, say: modulus:100 rat(7)^12345678 gives 11 so the last 2 digits are 89. Good luck with your numerology. RJF > *From:*Richard Fateman [mailto:fateman@...] > *Sent:* April2314 2:11 PM > *To:* Scott Hughes; <maximadiscuss@...> > *Subject:* Re: [Maximadiscuss] Maxima (lisp on Mac)  heap exhausted > > what do you plan to do with a number raised to the power 60,000,000? > > If you are computing it in order to then reduce it modulo some other > number, > then there is a better way to do it. > > > > On 4/23/2014 6:34 AM, Scott Hughes wrote: > > Firstly, I'm not a programmer and only recently started to use > Maxima (command line). I don't know anything about lisp. I'm > running Maxima on a Mac but also on a PC through Windows. Both are > the command line version. > > I only need to use Maxima for a few basic operations, primarily > small numbers raised to large (at least 8) exponents. Aside from > that, I'm using the round and asksign commands. > > I'm encountering "heap exhausted" errors once my exponents get > beyond about 60 000 000 and have looked in manual and forums for > help. There is reference to changing heap size but I didn't > understand the command or even how to enter it. Can anyone help me > with very explicit instructions ? I'm not concerned with the > length of time the calculations take, just don't want them to > crash before they complete. > 
From: Leo Butler <l_butler@us...>  20140423 19:13:33

Helmut Jarausch <jarausch@...> writes: > On 04/23/2014 08:17:38 PM, Leo Butler wrote: >> Leo Butler <l_butler@...> writes: >> >> >> As for your message: note that I simply load the lisp file attached to >> my previous message. I do not modify src/nparse.lisp and I do not >> recompile all of Maxima  there is no need to do either to test the >> change. >> >> If you want to try the patched version out for a longer period and >> have it automatically loaded, you can add the line >> >> load("nparsepatch.lisp"); > > I've understood now, still it doesn't do what you made me expect. > >> Maxima 5.33.0 http://maxima.sourceforge.net > using Lisp SBCL 1.1.16 > 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) load("nparsepatch.lisp")$ > STYLEWARNING: redefining MAXIMA::MREADSYNERR in DEFUN > > (%i2) load("Eulix_T5.mac"); > > /Src/Src/Math/Maxima/HJ/Eulix_T5.mac:221:incorrect syntax: !! is not a > prefix operator > ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ > this file name is an improvement, but 221 is still the byte offset, > not the line. You need to turn on line reporting, as my examples posted earlier showed. See below. > > Furthermore, if, .e.g, Eulix_T5.mac has a > load("Eulix.mac"); > which contains a syntax error, I still get "stdin" as filename and the > byte offset > as it was without your nparsepatch. Please post a simple example; I do not get this behaviour. Here is what I get when I turn on both flags that control syntax error reporting: (%i2) report_synerr_info : report_synerr_line : true$ (%i3) load("/tmp/err.mac"); /tmp/err2.mac:line 2:incorrect syntax: Y is not an infix operator \nx:1 y; ^ Here are the contents of the two files: err.mac: x:2* y; load("/tmp/err2.mac"); err2.mac: w : sin(1); x:1 y; > Thanks, > Helmut > >  > Start Your Social Network Today  Download eXo Platform > Build your Enterprise Intranet with eXo Platform Software > Java Based Open Source Intranet  Social, Extensible, Cloud Ready > Get Started Now And Turn Your Intranet Into A Collaboration Platform > http://p.sf.net/sfu/ExoPlatform  Leo Butler <l_butler@...> SDF Public Access UNIX System  http://sdf.lonestar.org 
From: Scott Hughes <hughesscott@ho...>  20140423 19:09:42

An example... a: 50000000; b: 60000000; c: 40000000; round ((6^a2^b3^b)/(6^(a7)4^b5^b)); From: Gunter Königsmann [mailto:gunter@...] Sent: April2314 2:58 PM To: fateman@...; Scott Hughes; <maximadiscuss@...> Subject: Re: [Maximadiscuss] Maxima (lisp on Mac)  heap exhausted Another question: How do you calculate the powers? I get a different error message, if I Try. On 23. April 2014 20:11:26 MESZ, Richard Fateman <fateman@...> wrote: what do you plan to do with a number raised to the power 60,000,000? If you are computing it in order to then reduce it modulo some other number, then there is a better way to do it. On 4/23/2014 6:34 AM, Scott Hughes wrote: Firstly, I’m not a programmer and only recently started to use Maxima (command line). I don’t know anything about lisp. I’m running Maxima on a Mac but also on a PC through Windows. Both are the command line version. I only need to use Maxima for a few basic operations, primarily small numbers raised to large (at least 8) exponents. Aside from that, I’m using the round and asksign commands. I’m encountering “heap exhausted” errors once my exponents get beyond about 60 000 000 and have looked in manual and forums for help. There is reference to changing heap size but I didn’t understand the command or even how to enter it. Can anyone help me with very explicit instructions ? I’m not concerned with the length of time the calculations take, just don’t want them to crash before they complete. _____ Start Your Social Network Today  Download eXo Platform Build your Enterprise Intranet with eXo Platform Software Java Based Open Source Intranet  Social, Extensible, Cloud Ready Get Started Now And Turn Your Intranet Into A Collaboration Platform http://p.sf.net/sfu/ExoPlatform _____ Maximadiscuss mailing list Maximadiscuss@... https://lists.sourceforge.net/lists/listinfo/maximadiscuss  Diese Nachricht wurde von meinem Mobiltelefon mit Kaiten Mail gesendet. 
From: Gunter Königsmann <gunter@pe...>  20140423 18:58:04

Another question: How do you calculate the powers? I get a different error message, if I Try. On 23. April 2014 20:11:26 MESZ, Richard Fateman <fateman@...> wrote: >what do you plan to do with a number raised to the power 60,000,000? > >If you are computing it in order to then reduce it modulo some other >number, >then there is a better way to do it. > > > >On 4/23/2014 6:34 AM, Scott Hughes wrote: >> >> Firstly, I'm not a programmer and only recently started to use Maxima > >> (command line). I don't know anything about lisp. I'm running Maxima >> on a Mac but also on a PC through Windows. Both are the command line >> version. >> >> I only need to use Maxima for a few basic operations, primarily small > >> numbers raised to large (at least 8) exponents. Aside from that, I'm >> using the round and asksign commands. >> >> I'm encountering "heap exhausted" errors once my exponents get beyond > >> about 60 000 000 and have looked in manual and forums for help. There > >> is reference to changing heap size but I didn't understand the >command >> or even how to enter it. Can anyone help me with very explicit >> instructions ? I'm not concerned with the length of time the >> calculations take, just don't want them to crash before they >complete. >> >> > > > > > > >Start Your Social Network Today  Download eXo Platform >Build your Enterprise Intranet with eXo Platform Software >Java Based Open Source Intranet  Social, Extensible, Cloud Ready >Get Started Now And Turn Your Intranet Into A Collaboration Platform >http://p.sf.net/sfu/ExoPlatform > > > >_______________________________________________ >Maximadiscuss mailing list >Maximadiscuss@... >https://lists.sourceforge.net/lists/listinfo/maximadiscuss  Diese Nachricht wurde von meinem Mobiltelefon mit Kaiten Mail gesendet. 
From: Gunter Königsmann <gunter@pe...>  20140423 18:57:02

On windows I seem to remember that you can use a whole gigabyte before maxima gives up. This could mean you are dealing with numbers that might exceed any amount of RAM your computer might have installed. If you just want to know the length of the result of the calculation you might want to calculate a logarithm of the result: the logarithm atva base of ten tells you how many digits a number needs. If your problem is that you want to keep too many but numbers in memory you might want to look at the nolabels flag, too. On 23. April 2014 20:11:26 MESZ, Richard Fateman <fateman@...> wrote: >what do you plan to do with a number raised to the power 60,000,000? > >If you are computing it in order to then reduce it modulo some other >number, >then there is a better way to do it. > > > >On 4/23/2014 6:34 AM, Scott Hughes wrote: >> >> Firstly, I'm not a programmer and only recently started to use Maxima > >> (command line). I don't know anything about lisp. I'm running Maxima >> on a Mac but also on a PC through Windows. Both are the command line >> version. >> >> I only need to use Maxima for a few basic operations, primarily small > >> numbers raised to large (at least 8) exponents. Aside from that, I'm >> using the round and asksign commands. >> >> I'm encountering "heap exhausted" errors once my exponents get beyond > >> about 60 000 000 and have looked in manual and forums for help. There > >> is reference to changing heap size but I didn't understand the >command >> or even how to enter it. Can anyone help me with very explicit >> instructions ? I'm not concerned with the length of time the >> calculations take, just don't want them to crash before they >complete. >> >> > > > > > > >Start Your Social Network Today  Download eXo Platform >Build your Enterprise Intranet with eXo Platform Software >Java Based Open Source Intranet  Social, Extensible, Cloud Ready >Get Started Now And Turn Your Intranet Into A Collaboration Platform >http://p.sf.net/sfu/ExoPlatform > > > >_______________________________________________ >Maximadiscuss mailing list >Maximadiscuss@... >https://lists.sourceforge.net/lists/listinfo/maximadiscuss  Diese Nachricht wurde von meinem Mobiltelefon mit Kaiten Mail gesendet. 
From: Helmut Jarausch <jarausch@ig...>  20140423 18:54:50

On 04/23/2014 08:17:38 PM, Leo Butler wrote: > Leo Butler <l_butler@...> writes: > > > As for your message: note that I simply load the lisp file attached to > my previous message. I do not modify src/nparse.lisp and I do not > recompile all of Maxima  there is no need to do either to test the > change. > > If you want to try the patched version out for a longer period and > have it automatically loaded, you can add the line > > load("nparsepatch.lisp"); I've understood now, still it doesn't do what you made me expect. > Maxima 5.33.0 http://maxima.sourceforge.net using Lisp SBCL 1.1.16 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) load("nparsepatch.lisp")$ STYLEWARNING: redefining MAXIMA::MREADSYNERR in DEFUN (%i2) load("Eulix_T5.mac"); /Src/Src/Math/Maxima/HJ/Eulix_T5.mac:221:incorrect syntax: !! is not a prefix operator ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ this file name is an improvement, but 221 is still the byte offset, not the line. Furthermore, if, .e.g, Eulix_T5.mac has a load("Eulix.mac"); which contains a syntax error, I still get "stdin" as filename and the byte offset as it was without your nparsepatch. Thanks, Helmut 
From: Raymond Toy <toy.raymond@gm...>  20140423 18:34:21

>>>>> "Pankaj" == Pankaj Sejwal <pankajsejwal@...> writes: Pankaj> >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Pankaj> Don't know what ctrl+k does, but that 2 pages of documentation is probably all that is available Pankaj> in maxima. The huge long list of lisp files are the files needed to implement the functions Pankaj> supported by Pankaj> maxima. There is a good chance that many more functions could be exposed to maxima, but no one Pankaj> has expressed a need for any additional lapack functions, so the interface hasn't been written Not Pankaj> 100% sure, but I think that all of the lapack files in maxima have Pankaj> been succesfully translated. >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> . Pankaj> If most of lapack has been translated than I think 99% Pankaj> work should be over and it will be great if someone who Pankaj> knows how to write this interface could devote some of his Pankaj> time. We will not have to write any doc for lapack Or if I would start by looking at, say, share/lapack/dgeqrf.lisp or eigensys.lisp to see how to implement an interface to the underlying lapack routines. This requires some knowledge of lisp, some knowledge of how Fortran works, some knowledge of lapack, and some knowledge of how Maxima stores matrices. None of these are hard, but does take a little time. If you are interested in a particular lapack function that doesn't have an interface for maxima, start with implementing that. I can help with that. Actually, there's a huge amount of boilerplate stuff and perhaps a lot of this could be hidden behind some macros to make it more uniform and consistent. This will also help with writing new interfaces. I can try implementing that. Ray 
From: Leo Butler <l_butler@us...>  20140423 18:17:58

From: Robert Dodier <robert.dodier@gm...>  20140423 18:16:03

On 20140423, 本田康晃 <yasuaki.honda@...> wrote: > (%i2) qe([[E,x],[E,y]], x^2+2*y^2=3 %and y=x/2+m); > > (%o2) (2*m3<=0) %and (2*m+3>=0) Please help me understand this result. By eliminating x and y, we get a result in terms of m alone. But what problem does that solve? If we are interested in intersections of an ellipse with a line, surely we need to know about x and y as well as m. I don't know anything about quantifier elimination in general so I am probably missing the point. Thanks for any information. best, Robert Dodier 
From: Richard Fateman <fateman@be...>  20140423 18:11:42

what do you plan to do with a number raised to the power 60,000,000? If you are computing it in order to then reduce it modulo some other number, then there is a better way to do it. On 4/23/2014 6:34 AM, Scott Hughes wrote: > > Firstly, I'm not a programmer and only recently started to use Maxima > (command line). I don't know anything about lisp. I'm running Maxima > on a Mac but also on a PC through Windows. Both are the command line > version. > > I only need to use Maxima for a few basic operations, primarily small > numbers raised to large (at least 8) exponents. Aside from that, I'm > using the round and asksign commands. > > I'm encountering "heap exhausted" errors once my exponents get beyond > about 60 000 000 and have looked in manual and forums for help. There > is reference to changing heap size but I didn't understand the command > or even how to enter it. Can anyone help me with very explicit > instructions ? I'm not concerned with the length of time the > calculations take, just don't want them to crash before they complete. > > 
From: Stavros Macrakis <macrakis@al...>  20140423 17:58:35

On Wed, Apr 23, 2014 at 1:14 PM, Pankaj Sejwal <pankajsejwal@...>wrote: > Hi, > 1) `member` function works only on first > > level of list, ex, > ... > In my opinion this is not completely correct. > That is the correct behavior  member is not membership in a tree, but in a list. This is similar to the similar settheoretic concept of membership: a is a member of {a,b,c}, but not a member of {{a},b,c}. You are of course free to write your own tree_member function if you find it useful. Similarly with your atlevel function. By the way, for the function to select nonlists, I think you'd be better off using 'not listp(x)' rather than trying to enumerate all the things that aren't lists. Perhaps you could explain why you find nestedlist membership and levelextraction so useful  maybe there is a better way to solve your original problem? s 2) > I found that there is no function to reach to some specific level in > nested lists and collect the same level lists from various sublists, so I > tried, >  > atlevel(lis,n):=block([li:lis,elem:[]],load(basic),for i:1 thru n do > (li:apply(append,li), > push(sublist(li,lambda([x],(symbolp(x) or > integerp(x) or operatorp(x) or stringp(x) or floatnump(x)))),elem), > li:sublist(li,listp)),[reverse(elem),li])$ > > 
From: Pankaj Sejwal <pankajsejwal@gm...>  20140423 17:15:43

Hi, 1) `member` function works only on first level of list, ex,  %IN:: li : [[a, b, c],[2,3,4,[3,4,5,[00,99,88,[77,66,55,"+"]]]],[ff,gg,tt], [[t, y, u,[jj,kk,ll]]]] $ %IN:: member([a,b,c],li) /*True*/ %IN:: member(b,li); /*False*/  In my opinion this is not completely correct. So, I tried writing this,  member1(lis,l):=block([li:lis,p:false],load(basic),if member(l,li) then p:true else (while(not li=[]) do (li:apply(append,li), if(member(l,li)) then p:true,li:sublist(li,listp))),p)$ Usage:: member1(b,li); /*True*/ member1([a,b,c],li)/*True*/ member1 ("+",li)/*True*/  Should it not be included in maxima ? 2) I found that there is no function to reach to some specific level in nested lists and collect the same level lists from various sublists, so I tried,  atlevel(lis,n):=block([li:lis,elem:[]],load(basic),for i:1 thru n do (li:apply(append,li), push(sublist(li,lambda([x],(symbolp(x) or integerp(x) or operatorp(x) or stringp(x) or floatnump(x)))),elem), li:sublist(li,listp)),[reverse(elem),li])$ Usage:: atlevel(li,2); /* [[[a,b,c,2,3,4,ff,gg,tt],[3,4,5,t,y,u]],[[0,99,88,[77,66,55,"a"]],[jj,kk,ll]]] */ %IN:: atlevel(li,4) /* [[[a,b,c,2,3,4,ff,gg,tt],[3,4,5,t,y,u],[0,99,88,jj,kk,ll],[77,66,55,"a"]],[]]  It collects two lists, one of them collects elements in each level before the level specified and next collects the sublists at specified level. These are very tiny functions but can be helpful in list manipulation. Please try these and suggest optimizations.  Regards, Pankaj Sejwal _______________________________________________ "The more I read, the more I acquire, the more certain I am that I know nothing.”  Voltaire <http://www.goodreads.com/author/show/5754446.Voltaire>; 
From: Pankaj Sejwal <pankajsejwal@gm...>  20140423 16:30:03

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Don't know what ctrl+k does, but that 2 pages of documentation is probably all that is available in maxima. The huge long list of lisp files are the files needed to implement the functions supported by maxima. There is a good chance that many more functions could be exposed to maxima, but no one has expressed a need for any additional lapack functions, so the interface hasn't been written Not 100% sure, but I think that all of the lapack files in maxima have been succesfully translated. >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>. If most of lapack has been translated than I think 99% work should be over and it will be great if someone who knows how to write this interface could devote some of his time. We will not have to write any doc for lapack Or if there is some sample guidance doc I would love to do it. Ctrl+k is function suggestion keyword.  Regards, Pankaj Sejwal _______________________________________________ "The more I read, the more I acquire, the more certain I am that I know nothing.”  Voltaire 
From: 本田康晃 <yasuaki.honda@gm...>  20140423 16:20:16

Dear Maixma users, Berns san, You may already noticed that the new version (v2.2) of Maxima on Android was pushed to Google Play (a week ago). This one includes the latest Maxima 5.33.0. Also, it includes Qepcad B and Qepmax for quantifier elimination. (%i1) load(qepmax)$ (%i2) qe([[E,x],[E,y]], x^2+2*y^2=3 %and y=x/2+m); (%o2) (2*m3<=0) %and (2*m+3>=0) on your mobile phone and tablet. Thanks and best regards, Yasuaki Honda, Chiba, Japan 20140406 21:50 GMT+09:00 Berns Buenaobra <berns.buenaobra@...>: > Honda San looking forward for the Android version on Google Play soon! My > students will be awed! > > Thanks, > Berns B. > > Bernardino (Berns) J. Buenaobra > Faculty at the Dept. of Physics  University of San Carlos > Contact information: > For immediate reply on issues on grades,ISMIS, exams, > attendance, emergencies: > Mobile +639054732772 > For online academic consulting in assisted chat, voice, video: > Google+ berns.buenaobra@... > Skype ID: poormanphysics > Professional reference Research and Education link: > http://ph. linkedin.com/pub/bernardinobuenaobra/16/42/637 > > > On Apr 5, 2014 1:13 PM, "本田康晃" <yasuaki.honda@...> wrote: > >> Dear Robert san, >> >> Thanks for the new release. >> >> I just would like to report that the following combinations work: >> Mac OS X 10.9.2, ECL12.12.1, maxima 5.33.0 >> imaxima works fine, too, on this setting. >> Android 4.4.2, ECL 12.12.1, maxima 5.33.0 >> I will create Maxima on Android package within couple of weeks. >> >> Thanks and best regards, >> Yasuaki Honda, Chiba, Japan >> >> >> 20140403 4:34 GMT+09:00 Robert Dodier <robert.dodier@...>: >> >>> Hey, I've tagged 5.33.0 in Git (on branch5_33) and uploaded rpms and >>> tar.gz to the Sourceforge file manager: >>> http://sourceforge.net/p/maxima/files >>> >>> Please give it a try. If someone can upload a Windows installer, that >>> would be terrific. >>> >>> best, >>> >>> Robert Dodier >>> >>> >>> >>>  >>> _______________________________________________ >>> Maximadiscuss mailing list >>> Maximadiscuss@... >>> https://lists.sourceforge.net/lists/listinfo/maximadiscuss >>> >> >> >> >>  >> >> _______________________________________________ >> Maximadiscuss mailing list >> Maximadiscuss@... >> https://lists.sourceforge.net/lists/listinfo/maximadiscuss >> >> 
From: Stavros Macrakis <macrakis@al...>  20140423 14:03:37

The multiplication operator "." denotes noncommutative multiplication; a^^n == a . ... . a <n times>. There is a large number of global variables controlling simplification of ".", all starting with "dot". The command [?? dot] will list them. If you set dotscrules:true, then <scalar>.A and A.<scalar> simplify to <scalar>*A. Numbers and constants are always scalar, and <number>.A always simplifies to <number>*A (regardless of dotscrules). You can explicitly declare variables scalar and nonscalar. s On Tue, Apr 22, 2014 at 9:55 AM, Henry Baker <hbaker1@...> wrote: > I'd like to fiddle around with some _matrix_ equations, with > _matrix_ unknowns. > > Is there a way in Maxima to declare a noncommutative variable? > > > >  > Start Your Social Network Today  Download eXo Platform > Build your Enterprise Intranet with eXo Platform Software > Java Based Open Source Intranet  Social, Extensible, Cloud Ready > Get Started Now And Turn Your Intranet Into A Collaboration Platform > http://p.sf.net/sfu/ExoPlatform > _______________________________________________ > Maximadiscuss mailing list > Maximadiscuss@... > https://lists.sourceforge.net/lists/listinfo/maximadiscuss > 