Re: [Maxima-discuss] ratsimp not working as expected

[andre m. <and...@gm...>]

On 03/08/2014 04:32 AM, Thomas D. Dean wrote:
> On 03/07/14 10:20, andre maute wrote:
>> @Thomas: what sort of linux distribution do you have?
>>

*snip*
Thomas' detailed description.
*snip*

>   > maxima
> Maxima 5.32.1 http://maxima.sourceforge.net
> using Lisp GNU Common Lisp (GCL) GCL 2.6.10 (a.k.a. GCL)
> 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) display2d : false;
>
> (%o1) false
> (%i2) d : 4$
> (%i3) A : genmatrix(lambda([k,l],v[[k,l]]-v[[0,l]]),d,d)$
> (%i4) detA : factor(determinant(A))$
> (%i5) AinvT : transpose(invert(A))$
> (%i6) AinvT : ratsimp(AinvT)$ /* error trigger */
> (%i7) zz : [0,-1,0,0]$
> (%i8) h1 : transpose(matrix(zz))$
> (%i9) h2 : AinvT.h1$
> (%i10) h3 : sum(h2[m,1]^2,m,1,d)*detA^2$
> (%i11) res : ratsimp(h3)$ /* chokes if trigger is commented */
> (%i12) res:ratsimp(h3)$
>
> Is this the correct sequence of input to generate the problem?
No, it is not. The first ratsimp line has to be commented.
But for GCL this code will even segfault for the invert() in (%i5).

Andre