Hi all, I used the following code in notebook but it fails to give results.
Any suggestions?
var('a11,a12,a13,a14,a15,a16,a22,a23,a24,a25,a26,a33,a34,a35,a36,a44,a45,a46,a55,a56,a66,p,q,r',
domain='real')
g11=a11*p*p+a66*q*q+a55*r*r+2.0*a16*p*q+2.0*a15*p*r+2.0*a56*q*r;
g22=a66*p*p+a22*q*q+a44*r*r+2.0*a26*p*q+2.0*a46*p*r+2.0*a24*q*r;
g33=a55*p*p+a44*q*q+a33*r*r+2.0*a45*p*q+2.0*a35*p*r+2.0*a34*q*r;
g12=a16*p*p+a26*q*q+a45*r*r+(a12+a66)*p*q+(a14+a56)*p*r+(a46+a25)*q*r;
g13=a15*p*p+a46*q*q+a35*r*r+(a14+a56)*p*q+(a13+a55)*p*r+(a36+a45)*q*r;
g23=a56*p*p+a24*q*q+a34*r*r+(a46+a25)*p*q+(a36+a45)*p*r+(a23+a44)*q*r;
g0=matrix(3,3,[g11,g12,g13,g12,g22,g23,g13,g12,g33]).eigenvectors_right()
It shows error "TypeError: ECL says: Memory limit reached. Please jump to
an outer pointer, quit program and enlarge the memory limits before
executing the program again." I have 9GB ram, and it fails before
exhausting the ram.
The reason I insist on calculating symbolic expressions first instead of
numerical solution is that (1) aij and p,q,r are varying, (2) I want the
eigenvalue associated with the same eigenvector or eigensystem for
different aij, p,q,r. If I use numerical solution I lose the track of the
eigensystem. Thanks very much.
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To post to this group, send email to sage-support@googlegroups.com.
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.
