Thanks for the reply.
I tried both versions in a fresh maxima--same result.
in the first I tried
aa.aa$
%;
The error is triggered on %;
Similarly for "."(aa,aa)$
followed by ;$
I'm using clisp-2.39 and today's CVS for maxima.
David
On Wed, 2007-05-09 at 17:21 -0400, Stavros Macrakis wrote:
> On 5/9/07, David Ronis <ronis at ronispc.chem.mcgill.ca> wrote:
> The following gives an error with the cvs version of maxima:
>
> This appears to be a new bug. It works fine in Maxima 5.11.0 (GCL
> 2.6.8 / Windows) and in 5.11.0cvs (Clisp) built a few weeks ago.
>
> What version of Lisp are you using? I assume you don't have any
> special initializations. Have you tested this in a fresh Maxima (not
> after kill(all))?
>
> Please try the following things:
>
> Calculate aa.aa without displaying it using result: aa.aa$
> instead of aa.aa; . Then display it afterwards result; --
> this will check if it's a display problem.
>
> Input the product as "."(aa,aa) instead of using infix -- this will
> check if it's a parser problem.
>
> Thanks,
>
> -s
>
>
> s[0]:matrix([1,0],[0,1]);
> s[1]:matrix([1,0],[0,-1]);
> s[2]:matrix([0,1],[1,0]);
> s[3]:matrix([0,%i],[-%i,0]);
>
> aa:a0*s[0]+a1*s[1]+a2*s[2]+a3*s[3];
>
> aa.aa;
>
> gives:
> (%i1) s[0]:matrix([1,0],[0,1]);
> [ 1 0 ]
> (%o1) [ ]
> [ 0 1 ]
> (%i2) s[1]:matrix([1,0],[0,-1]);
> [ 1 0 ]
> (%o2) [ ]
> [ 0 - 1 ]
> (%i3) s[2]:matrix([0,1],[1,0]);
> [ 0 1 ]
> (%o3) [ ]
> [ 1 0 ]
> (%i4) s[3]:matrix([0,%i],[-%i,0]);
> [ 0 %i ]
> (%o4) [ ]
> [ - %i 0 ]
> (%i5) aa:a0*s[0]+a1*s[1]+a2*s[2]+a3*s[3];
> [ a1 + a0 %i a3 + a2 ]
> (%o5) [ ]
> [ a2 - %i a3 a0 - a1 ]
> (%i6) aa.aa;
> Maxima encountered a Lisp error:
>
>
> SYMBOL-NAME: 1 is not a symbol
>
> Automatically continuing.
> To reenable the Lisp debugger set *debugger-hook* to nil.
>
> There is a more general version of this problem that I'm
> trying to deal
> with. I have expressions that contain linear combinations of
> matrices
> like
>
> A[i]=a[i,1]*M[1] + a[i,2]*M[2] etc.
>
> where M[i] is a matrix and a[i,j] is a scalar. Maxima does
> seem to
> create the matrix as expected, but A_i.A_i behaves as if both
> the a's
> and the M's are some sort of array. I've tried scalar(a) as
> well as
> scalar(a[i,j]) constructs but this doesn't seem to help.
>
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima