# Ticket #141: shimorin@kth.se_18Jun2009.txt

File shimorin@kth.se_18Jun2009.txt, 2.3 KB (added by seelisch, 14 years ago)

Singular code producing error

Line
1monitor(":w ResultGHB.txt","o");
2ring r=5,(u223,u222,u12,u23,u1,u2,delta2v,deltav1,deltav2,deltav3,v11,v12,v13,v22,v23,v33,v1,v2,v3,v),wp(1024,1024,256,256,64,64,4,3,3,3,2,2,2,2,2,2,1,1,1,0);
3
4//i to je chto j//
5
6//Simplifying assumptions://
7number u111=1;
8number u112=0;
9number u113=0;
10number u123=0;
11
12//extra//
13//number u222=1;
14number u122=-1;
15number u11=1;
16number u13=0;
17number u22=-1;
18//number u23=-1;
19//number u1=1;
20number u3=1;
21
22
23poly Rv=v-u1^2-u2^2-u3^2;
24poly Rv1=v1-2*(u11*u1+u12*u2+u13*u3);
25poly Rv2=v2-2*(u12*u1+u22*u2+u23*u3);
26poly Rv3=v3-2*(u13*u1+u23*u2+(-u11-u22)*u3);
27poly Rv11=v11-2*(u111*u1+u112*u2+u113*u3+u11*u11+u12*u12+u13*u13);
28poly Rv12=v12-2*(u112*u1+u122*u2+u123*u3+u11*u12+u12*u22+u13*u23);
29poly Rv13=v13-2*(u113*u1+u123*u2+(-u111-u122)*u3+u11*u13+u12*u23+u13*(-u11-u22));
30poly Rv22=v22-2*(u122*u1+u222*u2+u223*u3+u12*u12+u22*u22+u23*u23);
31poly Rv23=v23-2*(u123*u1+u223*u2+(-u112-u222)*u3+u12*u13+u22*u23+u23*(-u11-u22));
32poly Rv33=v33-2*((-u111-u122)*u1+(-u112-u222)*u2+(-u113-u223)*u3+u13*u13+u23*u23+(-u11-u22)*(-u11-u22));
33
34poly Rdeltav1=deltav1-4*(u11*u111+2*u12*u112+2*u13*u113+u22*u122+2*u23*u123+(-u11-u22)*(-u111-u122));
35poly Rdeltav2=deltav2-4*(u11*u112+2*u12*u122+2*u13*u123+u22*u222+2*u23*u223+(-u11-u22)*(-u112-u222));
36poly Rdeltav3=deltav3-4*(u11*u113+2*u12*u123+2*u13*(-u111-u122)+u22*u223+2*u23*(-u112-u222)+(-u11-u22)*(-u113-u223));
37
39//poly Rdeltav=deltav-(v11+v22+v33);
40
42
44
45//poly RdetD2v=detD2v-(v11*v22*v33+2*v12*v23*v13-v13*v22*v13-v12*v12*v13-v11*v23*v23);
46
47//poly RD2vHS=D2vHS-(v11^2+2*v12^2+2*v13^2+v22^2+2*v23^2+v33^2);
48