Thank you Frederic. It does work charmingly. Thanks a lot. subsop works as well. Pari works to the extent that I checked (I tried to initialize a complicated elliptic curve with ellinit and that worked). Thanks to you, now i would not need to import pari from another library.
I just discovered a problem with this version that I did not have with version 0.4.4. A simple example that will illustrate this problem is the following code that tries to symbolically compute the nullspace of a complicated matrix:
Code : Tout sélectionner
from giacpy import matrix, ker
m=matrix("[(-2*l1+2)/2,(2*l1+2)/2,0,0,(l1*a2+a2)/2,(l1*a2-a2)/2,(l1*d2-d2)/2,(-l1*d2-d2)/2,0,0,0,0,0,0,0,0],[0,0,(-2*l1+2)/2,(2*l1+2)/2,(-l1*d2+d2)/2,(l1*d2+d2)/2,(-l1*a2-a2)/2,(-l1*a2+a2)/2,(-2*l1+2)/2,(2*l1+2)/2,0,0,(l1*a2+a2)/2,(l1*a2-a2)/2,(l1*d2-d2)/2,(-l1*d2-d2)/2],[0,0,0,0,0,0,0,0,0,0,(-2*l1+2)/2,(2*l1+2)/2,(-l1*d2+d2)/2,(l1*d2+d2)/2,(-l1*a2-a2)/2,(-l1*a2+a2)/2],[0,0,0,0,0,0,(l1+1)/2,(l1-1)/2,0,0,0,0,0,0,0,0],[0,0,0,0,(l1+1)/2,(l1-1)/2,0,0,0,0,0,0,0,0,(l1+1)/2,(l1-1)/2],[0,0,0,0,0,0,0,0,0,0,0,0,(l1+1)/2,(l1-1)/2,0,0],[0,0,(2*l1+2)/2,(-2*l1+2)/2,(l1*d2+d2)/2,(-l1*d2+d2)/2,(l1*a2-a2)/2,(l1*a2+a2)/2,0,0,0,0,0,0,0,0],[(-2*l1-2)/2,(2*l1-2)/2,0,0,(l1*a2-a2)/2,(l1*a2+a2)/2,(l1*d2+d2)/2,(-l1*d2+d2)/2,0,0,(2*l1+2)/2,(-2*l1+2)/2,(l1*d2+d2)/2,(-l1*d2+d2)/2,(l1*a2-a2)/2,(l1*a2+a2)/2],[0,0,0,0,0,0,0,0,(-2*l1-2)/2,(2*l1-2)/2,0,0,(l1*a2-a2)/2,(l1*a2+a2)/2,(l1*d2+d2)/2,(-l1*d2+d2)/2],[0,0,0,0,(-l1+1)/2,(-l1-1)/2,0,0,0,0,0,0,0,0,0,0],[0,0,0,0,0,0,(l1-1)/2,(l1+1)/2,0,0,0,0,(-l1+1)/2,(-l1-1)/2,0,0],[0,0,0,0,0,0,0,0,0,0,0,0,0,0,(l1-1)/2,(l1+1)/2]")
ker(m)
This could be more of a giac problem than a giacpy problem. Could you reproduce this error as well?
With giacpy version 0.4.4 I get the nullspace (I haven't verify its correctness but I would assume it is correct).