Commit f240b67c authored by Philippe Helluy's avatar Philippe Helluy

kernel flux and volume written

parent a745cc3e
......@@ -75,33 +75,35 @@ Italic" -1 43 "Times" 1 12 0 0 0 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "Maple
0 0 0 1 }{CSTYLE "2D Math Italic Small208" -1 201 "Times" 1 1 0 0 0 1
1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic Small207" -1 202 "Times"
1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Comment" -1 21 "Courie
r" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Caption Reference" -1
203 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Input Pla
ceholder" -1 204 "Courier New" 1 12 200 0 200 1 2 1 2 2 1 2 0 0 0 1 }
{CSTYLE "Code" -1 205 "Courier New" 1 12 255 0 0 1 2 2 2 2 2 2 0 0 0 1
}{CSTYLE "2D Math Italic Small206" -1 206 "Times" 1 1 0 0 0 1 1 2 2 2
2 2 0 0 0 1 }{CSTYLE "2D Math Italic Small205" -1 207 "Times" 1 1 0 0
0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Inert Output" -1 208 "Times" 1
12 144 144 144 1 2 2 2 2 1 2 0 0 0 1 }{CSTYLE "2D Math Italic Small"
-1 209 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Comment"
-1 18 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underlin
ed Bold" -1 41 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 0 0 0 1 }{CSTYLE "Copy
right" -1 34 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Hyperl
ink" -1 17 "Times" 1 12 0 128 128 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "Help
Underlined" -1 44 "Times" 1 12 0 0 0 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "
Prompt" -1 1 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D M
ath" -1 2 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Input"
-1 19 "Times" 1 12 0 0 0 1 2 2 2 2 1 2 0 0 0 1 }{CSTYLE "Header and F
ooter" -1 210 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Text"
-1 211 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Notes"
-1 37 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Bold"
-1 39 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Equation Labe
l" -1 212 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "LaTeX" -1
32 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Popup" -1 31 "T
imes" 1 12 0 128 128 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "Dictionary Hyperl
ink" -1 45 "Times" 1 12 147 0 15 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "Capti
on Text" -1 213 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Hel
p Fixed" -1 23 "Courier" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }}
r" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic Small209"
-1 203 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Caption Refe
rence" -1 204 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Maple
Input Placeholder" -1 205 "Courier New" 1 12 200 0 200 1 2 1 2 2 1 2
0 0 0 1 }{CSTYLE "Code" -1 206 "Courier New" 1 12 255 0 0 1 2 2 2 2 2
2 0 0 0 1 }{CSTYLE "2D Math Italic Small206" -1 207 "Times" 1 1 0 0 0
1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic Small205" -1 208 "Times
" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Inert Output" -1 209 "T
imes" 1 12 144 144 144 1 2 2 2 2 1 2 0 0 0 1 }{CSTYLE "2D Math Italic \+
Small" -1 210 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Com
ment" -1 18 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Un
derlined Bold" -1 41 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 0 0 0 1 }{CSTYLE
"Copyright" -1 34 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "
Hyperlink" -1 17 "Times" 1 12 0 128 128 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE
"Help Underlined" -1 44 "Times" 1 12 0 0 0 1 2 2 1 2 2 2 0 0 0 1 }
{CSTYLE "Prompt" -1 1 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }
{CSTYLE "2D Math" -1 2 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }
{CSTYLE "2D Input" -1 19 "Times" 1 12 0 0 0 1 2 2 2 2 1 2 0 0 0 1 }
{CSTYLE "Header and Footer" -1 211 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0
0 0 1 }{CSTYLE "Text" -1 212 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1
}{CSTYLE "Help Notes" -1 37 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }
{CSTYLE "Help Bold" -1 39 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }
{CSTYLE "Equation Label" -1 213 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0
1 }{CSTYLE "LaTeX" -1 32 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }
{CSTYLE "Popup" -1 31 "Times" 1 12 0 128 128 1 1 2 1 2 2 2 0 0 0 1 }
{CSTYLE "Dictionary Hyperlink" -1 45 "Times" 1 12 147 0 15 1 2 2 1 2 2
2 0 0 0 1 }{CSTYLE "Caption Text" -1 214 "Times" 1 12 0 0 0 1 2 1 2 2
2 2 0 0 0 1 }{CSTYLE "Help Fixed" -1 23 "Courier" 1 10 0 0 0 1 2 2 2
2 2 2 0 0 0 1 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "restart:with(linalg
):with(ListTools);with(CodeGeneration):" }}{PARA 11 "" 1 "" {XPPMATH
20 "7>I,BinaryPlaceG6\"I-BinarySearchGF$I+CategorizeGF$I+DotProductGF$
......@@ -124,10 +126,10 @@ F$7%F$F$F&7%F&F$F&7%F&F&F&7%F$F&F&" }}}{EXCHG {PARA 0 "> " 0 ""
e2[j],node),j=1..8)];C(permut);" }}{PARA 11 "" 1 "" {XPPMATH 20 "7*7%
\"\"!F$F$7%F$F$\"\"\"7%F&F$F&7%F&F$F$7%F$F&F$7%F$F&F&7%F&F&F&7%F&F&F$"
}}{PARA 11 "" 1 "" {XPPMATH 20 "7*\"\"\"\"\"&\"\"'\"\"#\"\"%\"\")\"\"
(\"\"$" }}{PARA 6 "" 1 "" {TEXT 214 10 "cg[0] = 1;" }{TEXT 214 11 "\nc
g[1] = 5;" }{TEXT 214 11 "\ncg[2] = 6;" }{TEXT 214 11 "\ncg[3] = 2;" }
{TEXT 214 11 "\ncg[4] = 4;" }{TEXT 214 11 "\ncg[5] = 8;" }{TEXT 214
11 "\ncg[6] = 7;" }{TEXT 214 11 "\ncg[7] = 3;" }}}{EXCHG {PARA 0 "> "
(\"\"$" }}{PARA 6 "" 1 "" {TEXT 215 10 "cg[0] = 1;" }{TEXT 215 11 "\nc
g[1] = 5;" }{TEXT 215 11 "\ncg[2] = 6;" }{TEXT 215 11 "\ncg[3] = 2;" }
{TEXT 215 11 "\ncg[4] = 4;" }{TEXT 215 11 "\ncg[5] = 8;" }{TEXT 215
11 "\ncg[6] = 7;" }{TEXT 215 11 "\ncg[7] = 3;" }}}{EXCHG {PARA 0 "> "
0 "" {MPLTEXT 1 0 56 "P:=vector(8,[1,x,y,z,x*y,x*z,y*z,x*y*z]);phi:=ve
ctor(8);" }}{PARA 11 "" 1 "" {XPPMATH 20 "=6\"6#;\"\"\"\"\")E\\[l)F&F&
\"\"#I\"xGF#\"\"$I\"yGF#\"\"%I\"zGF#\"\"&*&F*F&F,F&\"\"'*&F*F&F.F&\"\"
......@@ -153,83 +155,89 @@ F%F&,&!\"\"F&I\"yGF%F&F&F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "*(I\"xG6\"
{MPLTEXT 1 0 23 "#Dd:=matrix(8,8,(i,j)->" }{MPLTEXT 1 0 19 "subs(x=nod
e[j,1],y=" }{MPLTEXT 1 0 21 "node[j,2],z=node[j,3]" }{MPLTEXT 1 0 18 "
,phi[i])-dk(i,j));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "with(
CodeGeneration):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "C(phi);"
}}{PARA 6 "" 1 "" {TEXT 214 41 "phi[0] = -(-1 + z) * (-1 + y) * (-1 +
x);" }}{PARA 6 "" 1 "" {TEXT 214 33 "phi[1] = x * (-1 + z) * (-1 + y)
;" }{TEXT 214 28 "\nphi[2] = -x * y * (-1 + z);" }{TEXT 214 34 "\nphi[
3] = y * (-1 + z) * (-1 + x);" }{TEXT 214 34 "\nphi[4] = z * (-1 + y) \+
* (-1 + x);" }{TEXT 214 28 "\nphi[5] = -x * z * (-1 + y);" }{TEXT 214
20 "\nphi[6] = x * y * z;" }{TEXT 214 28 "\nphi[7] = -y * z * (-1 + x)
;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "tau:=vector(3,[0,0,0])
;p:=array(1..3*8);" }}{PARA 11 "" 1 "" {XPPMATH 20 "=6\"6#;\"\"\"\"\"$
E\\[l$F&\"\"!\"\"#F)F'F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "=6\"6#;\"\"
\"\"#CE\\[l!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "#?submatrix
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "i:='i';for k to 3 do ta
u[k]:=sum(phi[i]*p[3*(i-1)+k],i=1..8);od;" }}{PARA 11 "" 1 ""
{XPPMATH 20 "I\"iG6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 ",2**,&!\"\"\"\"
\"I\"zG6\"F&F&,&F%F&I\"yGF(F&F&,&F%F&I\"xGF(F&F&&I\"pGF(6#F&F&F%**F,F&
F$F&F)F&&F.6#\"\"%F&F&**F,F&F*F&F$F&&F.6#\"\"(F&F%**F*F&F$F&F+F&&F.6#
\"#5F&F&**F'F&F)F&F+F&&F.6#\"#8F&F&**F,F&F'F&F)F&&F.6#\"#;F&F%**F,F&F*
F&F'F&&F.6#\"#>F&F&**F*F&F'F&F+F&&F.6#\"#AF&F%" }}{PARA 11 "" 1 ""
{XPPMATH 20 ",2**,&!\"\"\"\"\"I\"zG6\"F&F&,&F%F&I\"yGF(F&F&,&F%F&I\"xG
F(F&F&&I\"pGF(6#\"\"#F&F%**F,F&F$F&F)F&&F.6#\"\"&F&F&**F,F&F*F&F$F&&F.
6#\"\")F&F%**F*F&F$F&F+F&&F.6#\"#6F&F&**F'F&F)F&F+F&&F.6#\"#9F&F&**F,F
&F'F&F)F&&F.6#\"#<F&F%**F,F&F*F&F'F&&F.6#\"#?F&F&**F*F&F'F&F+F&&F.6#\"
#BF&F%" }}{PARA 11 "" 1 "" {XPPMATH 20 ",2**,&!\"\"\"\"\"I\"zG6\"F&F&,
&F%F&I\"yGF(F&F&,&F%F&I\"xGF(F&F&&I\"pGF(6#\"\"$F&F%**F,F&F$F&F)F&&F.6
#\"\"'F&F&**F,F&F*F&F$F&&F.6#\"\"*F&F%**F*F&F$F&F+F&&F.6#\"#7F&F&**F'F
&F)F&F+F&&F.6#\"#:F&F&**F,F&F'F&F)F&&F.6#\"#=F&F%**F,F&F*F&F'F&&F.6#\"
#@F&F&**F*F&F'F&F+F&&F.6#\"#CF&F%" }}}{EXCHG {PARA 0 "> " 0 ""
{MPLTEXT 1 0 28 "dtau:=jacobian(tau,[x,y,z]);" }}{PARA 11 "" 1 ""
{XPPMATH 20 "=6\"6$;\"\"\"\"\"$F%E\\[l*6$F'F',2*(,&!\"\"F&I\"yGF#F&F&,
&F-F&I\"xGF#F&F&&I\"pGF#6#F'F&F-*(F0F&F,F&&F26#\"\"'F&F&*(F0F&F.F&&F26
#\"\"*F&F-*(F.F&F/F&&F26#\"#7F&F&*(F,F&F/F&&F26#\"#:F&F&*(F0F&F,F&&F26
#\"#=F&F-*(F0F&F.F&&F26#\"#@F&F&*(F.F&F/F&&F26#\"#CF&F-6$F'F&,2*(,&F-F
&I\"zGF#F&F&F,F&F1F&F-*(FSF&F,F&F5F&F&*(F.F&FSF&F9F&F-*(F.F&FSF&F=F&F&
*(FTF&F,F&FAF&F&*(FTF&F,F&FEF&F-*(F.F&FTF&FIF&F&*(F.F&FTF&FMF&F-6$F&F&
,2*(FSF&F,F&&F26#F&F&F-*(FSF&F,F&&F26#\"\"%F&F&*(F.F&FSF&&F26#\"\"(F&F
-*(F.F&FSF&&F26#\"#5F&F&*(FTF&F,F&&F26#\"#8F&F&*(FTF&F,F&&F26#\"#;F&F-
*(F.F&FTF&&F26#\"#>F&F&*(F.F&FTF&&F26#\"#AF&F-6$\"\"#Fhp,2*(FSF&F/F&&F
26#FhpF&F-*(F0F&FSF&&F26#\"\"&F&F&*(F0F&FSF&&F26#\"\")F&F-*(FSF&F/F&&F
26#\"#6F&F&*(FTF&F/F&&F26#\"#9F&F&*(F0F&FTF&&F26#\"#<F&F-*(F0F&FTF&&F2
6#\"#?F&F&*(FTF&F/F&&F26#\"#BF&F-6$F'Fhp,2*(FSF&F/F&F1F&F-*(F0F&FSF&F5
F&F&*(F0F&FSF&F9F&F-*(FSF&F/F&F=F&F&*(FTF&F/F&FAF&F&*(F0F&FTF&FEF&F-*(
F0F&FTF&FIF&F&*(FTF&F/F&FMF&F-6$F&F',2*(F,F&F/F&FinF&F-*(F0F&F,F&F\\oF
&F&*(F0F&F.F&F`oF&F-*(F.F&F/F&FdoF&F&*(F,F&F/F&FhoF&F&*(F0F&F,F&F\\pF&
F-*(F0F&F.F&F`pF&F&*(F.F&F/F&FdpF&F-6$F&Fhp,2*(FSF&F/F&FinF&F-*(F0F&FS
F&F\\oF&F&*(F0F&FSF&F`oF&F-*(FSF&F/F&FdoF&F&*(FTF&F/F&FhoF&F&*(F0F&FTF
&F\\pF&F-*(F0F&FTF&F`pF&F&*(FTF&F/F&FdpF&F-6$FhpF',2*(F,F&F/F&F[qF&F-*
(F0F&F,F&F^qF&F&*(F0F&F.F&FbqF&F-*(F.F&F/F&FfqF&F&*(F,F&F/F&FjqF&F&*(F
0F&F,F&F^rF&F-*(F0F&F.F&FbrF&F&*(F.F&F/F&FfrF&F-6$FhpF&,2*(FSF&F,F&F[q
F&F-*(FSF&F,F&F^qF&F&*(F.F&FSF&FbqF&F-*(F.F&FSF&FfqF&F&*(FTF&F,F&FjqF&
F&*(FTF&F,F&F^rF&F-*(F.F&FTF&FbrF&F&*(F.F&FTF&FfrF&F-" }}}{EXCHG
{PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "C(dtau,optimize);" }}{PARA 6 "" 1 "
" {TEXT 214 12 "t1 = -1 + z;" }}{PARA 6 "" 1 "" {TEXT 214 12 "t2 = -1 \+
+ y;" }{TEXT 214 14 "\nt3 = t1 * t2;" }{TEXT 214 13 "\nt6 = y * t1;" }
{TEXT 214 13 "\nt9 = z * t2;" }{TEXT 214 13 "\nt12 = y * z;" }{TEXT
214 14 "\nt16 = -1 + x;" }{TEXT 214 16 "\nt17 = t1 * t16;" }{TEXT 214
14 "\nt19 = x * t1;" }{TEXT 214 15 "\nt23 = z * t16;" }{TEXT 214 13 "
\nt25 = x * z;" }{TEXT 214 16 "\nt30 = t2 * t16;" }{TEXT 214 14 "\nt32
= x * t2;" }{TEXT 214 13 "\nt34 = x * y;" }{TEXT 214 15 "\nt36 = y * \+
t16;" }{TEXT 214 115 "\ndtau[0][0] = -t3 * p[0] + t3 * p[3] - t6 * p[6
] + t6 * p[9] + t9 * p[12] - t9 * p[15] + t12 * p[18] - t12 * p[21];"
}{TEXT 214 121 "\ndtau[0][1] = -t17 * p[0] + t19 * p[3] - t19 * p[6] +
t17 * p[9] + t23 * p[12] - t25 * p[15] + t25 * p[18] - t23 * p[21];"
}{TEXT 214 121 "\ndtau[0][2] = -t30 * p[0] + t32 * p[3] - t34 * p[6] +
t36 * p[9] + t30 * p[12] - t32 * p[15] + t34 * p[18] - t36 * p[21];"
}{TEXT 214 116 "\ndtau[1][0] = -t3 * p[1] + t3 * p[4] - t6 * p[7] + t6
* p[10] + t9 * p[13] - t9 * p[16] + t12 * p[19] - t12 * p[22];" }
{TEXT 214 122 "\ndtau[1][1] = -t17 * p[1] + t19 * p[4] - t19 * p[7] + \+
t17 * p[10] + t23 * p[13] - t25 * p[16] + t25 * p[19] - t23 * p[22];"
}{TEXT 214 122 "\ndtau[1][2] = -t30 * p[1] + t32 * p[4] - t34 * p[7] +
t36 * p[10] + t30 * p[13] - t32 * p[16] + t34 * p[19] - t36 * p[22];"
}{TEXT 214 116 "\ndtau[2][0] = -t3 * p[2] + t3 * p[5] - t6 * p[8] + t
6 * p[11] + t9 * p[14] - t9 * p[17] + t12 * p[20] - t12 * p[23];" }
{TEXT 214 122 "\ndtau[2][1] = -t17 * p[2] + t19 * p[5] - t19 * p[8] + \+
t17 * p[11] + t23 * p[14] - t25 * p[17] + t25 * p[20] - t23 * p[23];"
}{TEXT 214 122 "\ndtau[2][2] = -t30 * p[2] + t32 * p[5] - t34 * p[8] +
t36 * p[11] + t30 * p[14] - t32 * p[17] + t34 * p[20] - t36 * p[23];"
}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 2 "?C" }}}{EXCHG {PARA 0 ">
" 0 "" {MPLTEXT 1 0 0 "" }}}}
CodeGeneration):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "C(phi, \+
optimize);" }}{PARA 6 "" 1 "" {TEXT 215 12 "t1 = -1 + z;" }}{PARA 6 ""
1 "" {TEXT 215 12 "t2 = -1 + y;" }{TEXT 215 13 "\nt4 = -1 + x;" }
{TEXT 215 12 "\nt8 = x * y;" }{TEXT 215 24 "\nphi[0] = -t1 * t2 * t4;"
}{TEXT 215 22 "\nphi[1] = x * t1 * t2;" }{TEXT 215 19 "\nphi[2] = -t8
* t1;" }{TEXT 215 22 "\nphi[3] = y * t1 * t4;" }{TEXT 215 22 "\nphi[4
] = z * t2 * t4;" }{TEXT 215 22 "\nphi[5] = -x * z * t2;" }{TEXT 215
17 "\nphi[6] = t8 * z;" }{TEXT 215 22 "\nphi[7] = -y * z * t4;" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "tau:=vector(3,[0,0,0]);p:=ar
ray(1..3*8);" }}{PARA 11 "" 1 "" {XPPMATH 20 "=6\"6#;\"\"\"\"\"$E\\[l$
F&\"\"!\"\"#F)F'F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "=6\"6#;\"\"\"\"#CE
\\[l!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "#?submatrix" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "i:='i';for k to 3 do tau[k]:
=sum(phi[i]*p[3*(i-1)+k],i=1..8);od;" }}{PARA 11 "" 1 "" {XPPMATH 20 "
I\"iG6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 ",2**,&\"\"\"!\"\"I\"zG6\"F%F
%,&F%F&I\"yGF(F%F%,&F%F&I\"xGF(F%F%&I\"pGF(6#F%F%F&**F,F%F$F%F)F%&F.6#
\"\"%F%F%**F,F%F*F%F$F%&F.6#\"\"(F%F&**F*F%F$F%F+F%&F.6#\"#5F%F%**F'F%
F)F%F+F%&F.6#\"#8F%F%**F,F%F'F%F)F%&F.6#\"#;F%F&**F,F%F*F%F'F%&F.6#\"#
>F%F%**F*F%F'F%F+F%&F.6#\"#AF%F&" }}{PARA 11 "" 1 "" {XPPMATH 20 ",2**
,&\"\"\"!\"\"I\"zG6\"F%F%,&F%F&I\"yGF(F%F%,&F%F&I\"xGF(F%F%&I\"pGF(6#
\"\"#F%F&**F,F%F$F%F)F%&F.6#\"\"&F%F%**F,F%F*F%F$F%&F.6#\"\")F%F&**F*F
%F$F%F+F%&F.6#\"#6F%F%**F'F%F)F%F+F%&F.6#\"#9F%F%**F,F%F'F%F)F%&F.6#\"
#<F%F&**F,F%F*F%F'F%&F.6#\"#?F%F%**F*F%F'F%F+F%&F.6#\"#BF%F&" }}{PARA
11 "" 1 "" {XPPMATH 20 ",2**,&\"\"\"!\"\"I\"zG6\"F%F%,&F%F&I\"yGF(F%F%
,&F%F&I\"xGF(F%F%&I\"pGF(6#\"\"$F%F&**F,F%F$F%F)F%&F.6#\"\"'F%F%**F,F%
F*F%F$F%&F.6#\"\"*F%F&**F*F%F$F%F+F%&F.6#\"#7F%F%**F'F%F)F%F+F%&F.6#\"
#:F%F%**F,F%F'F%F)F%&F.6#\"#=F%F&**F,F%F*F%F'F%&F.6#\"#@F%F%**F*F%F'F%
F+F%&F.6#\"#CF%F&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "dtau:=
jacobian(tau,[x,y,z]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "=6\"6$;\"\"\"
\"\"$F%E\\[l*6$F'F',2*(,&F&!\"\"I\"yGF#F&F&,&F&F-I\"xGF#F&F&&I\"pGF#6#
F'F&F-*(F0F&F,F&&F26#\"\"'F&F&*(F0F&F.F&&F26#\"\"*F&F-*(F.F&F/F&&F26#
\"#7F&F&*(F,F&F/F&&F26#\"#:F&F&*(F0F&F,F&&F26#\"#=F&F-*(F0F&F.F&&F26#
\"#@F&F&*(F.F&F/F&&F26#\"#CF&F-6$F'F&,2*(,&F&F-I\"zGF#F&F&F,F&F1F&F-*(
FSF&F,F&F5F&F&*(F.F&FSF&F9F&F-*(F.F&FSF&F=F&F&*(FTF&F,F&FAF&F&*(FTF&F,
F&FEF&F-*(F.F&FTF&FIF&F&*(F.F&FTF&FMF&F-6$F&F&,2*(FSF&F,F&&F26#F&F&F-*
(FSF&F,F&&F26#\"\"%F&F&*(F.F&FSF&&F26#\"\"(F&F-*(F.F&FSF&&F26#\"#5F&F&
*(FTF&F,F&&F26#\"#8F&F&*(FTF&F,F&&F26#\"#;F&F-*(F.F&FTF&&F26#\"#>F&F&*
(F.F&FTF&&F26#\"#AF&F-6$\"\"#Fhp,2*(FSF&F/F&&F26#FhpF&F-*(F0F&FSF&&F26
#\"\"&F&F&*(F0F&FSF&&F26#\"\")F&F-*(FSF&F/F&&F26#\"#6F&F&*(FTF&F/F&&F2
6#\"#9F&F&*(F0F&FTF&&F26#\"#<F&F-*(F0F&FTF&&F26#\"#?F&F&*(FTF&F/F&&F26
#\"#BF&F-6$F'Fhp,2*(FSF&F/F&F1F&F-*(F0F&FSF&F5F&F&*(F0F&FSF&F9F&F-*(FS
F&F/F&F=F&F&*(FTF&F/F&FAF&F&*(F0F&FTF&FEF&F-*(F0F&FTF&FIF&F&*(FTF&F/F&
FMF&F-6$F&F',2*(F,F&F/F&FinF&F-*(F0F&F,F&F\\oF&F&*(F0F&F.F&F`oF&F-*(F.
F&F/F&FdoF&F&*(F,F&F/F&FhoF&F&*(F0F&F,F&F\\pF&F-*(F0F&F.F&F`pF&F&*(F.F
&F/F&FdpF&F-6$F&Fhp,2*(FSF&F/F&FinF&F-*(F0F&FSF&F\\oF&F&*(F0F&FSF&F`oF
&F-*(FSF&F/F&FdoF&F&*(FTF&F/F&FhoF&F&*(F0F&FTF&F\\pF&F-*(F0F&FTF&F`pF&
F&*(FTF&F/F&FdpF&F-6$FhpF',2*(F,F&F/F&F[qF&F-*(F0F&F,F&F^qF&F&*(F0F&F.
F&FbqF&F-*(F.F&F/F&FfqF&F&*(F,F&F/F&FjqF&F&*(F0F&F,F&F^rF&F-*(F0F&F.F&
FbrF&F&*(F.F&F/F&FfrF&F-6$FhpF&,2*(FSF&F,F&F[qF&F-*(FSF&F,F&F^qF&F&*(F
.F&FSF&FbqF&F-*(F.F&FSF&FfqF&F&*(FTF&F,F&FjqF&F&*(FTF&F,F&F^rF&F-*(F.F
&FTF&FbrF&F&*(F.F&FTF&FfrF&F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1
0 17 "C(dtau,optimize);" }}{PARA 6 "" 1 "" {TEXT 215 12 "t1 = -1 + z;"
}}{PARA 6 "" 1 "" {TEXT 215 12 "t2 = -1 + y;" }{TEXT 215 1 "\n" }
{TEXT 215 13 "t3 = t1 * t2;" }{TEXT 215 1 "\n" }{TEXT 215 12 "t6 = y *
t1;" }{TEXT 215 1 "\n" }{TEXT 215 12 "t9 = z * t2;" }{TEXT 215 1 "\n"
}{TEXT 215 12 "t12 = y * z;" }{TEXT 215 1 "\n" }{TEXT 215 13 "t16 = -
1 + x;" }{TEXT 215 1 "\n" }{TEXT 215 15 "t17 = t1 * t16;" }{TEXT 215
1 "\n" }{TEXT 215 13 "t19 = x * t1;" }{TEXT 215 1 "\n" }{TEXT 215 14 "
t23 = z * t16;" }{TEXT 215 1 "\n" }{TEXT 215 12 "t25 = x * z;" }{TEXT
215 1 "\n" }{TEXT 215 15 "t30 = t2 * t16;" }{TEXT 215 1 "\n" }{TEXT
215 13 "t32 = x * t2;" }{TEXT 215 1 "\n" }{TEXT 215 12 "t34 = x * y;"
}{TEXT 215 1 "\n" }{TEXT 215 14 "t36 = y * t16;" }{TEXT 215 1 "\n" }
{TEXT 215 114 "dtau[0][0] = -t3 * p[0] + t3 * p[3] - t6 * p[6] + t6 * \+
p[9] + t9 * p[12] - t9 * p[15] + t12 * p[18] - t12 * p[21];" }{TEXT
215 1 "\n" }{TEXT 215 120 "dtau[0][1] = -t17 * p[0] + t19 * p[3] - t19
* p[6] + t17 * p[9] + t23 * p[12] - t25 * p[15] + t25 * p[18] - t23 *
p[21];" }{TEXT 215 1 "\n" }{TEXT 215 120 "dtau[0][2] = -t30 * p[0] + \+
t32 * p[3] - t34 * p[6] + t36 * p[9] + t30 * p[12] - t32 * p[15] + t34
* p[18] - t36 * p[21];" }{TEXT 215 1 "\n" }{TEXT 215 115 "dtau[1][0] \+
= -t3 * p[1] + t3 * p[4] - t6 * p[7] + t6 * p[10] + t9 * p[13] - t9 * \+
p[16] + t12 * p[19] - t12 * p[22];" }{TEXT 215 1 "\n" }{TEXT 215 121 "
dtau[1][1] = -t17 * p[1] + t19 * p[4] - t19 * p[7] + t17 * p[10] + t23
* p[13] - t25 * p[16] + t25 * p[19] - t23 * p[22];" }{TEXT 215 1 "\n"
}{TEXT 215 121 "dtau[1][2] = -t30 * p[1] + t32 * p[4] - t34 * p[7] + \+
t36 * p[10] + t30 * p[13] - t32 * p[16] + t34 * p[19] - t36 * p[22];"
}{TEXT 215 1 "\n" }{TEXT 215 115 "dtau[2][0] = -t3 * p[2] + t3 * p[5] \+
- t6 * p[8] + t6 * p[11] + t9 * p[14] - t9 * p[17] + t12 * p[20] - t12
* p[23];" }{TEXT 215 1 "\n" }{TEXT 215 121 "dtau[2][1] = -t17 * p[2] \+
+ t19 * p[5] - t19 * p[8] + t17 * p[11] + t23 * p[14] - t25 * p[17] + \+
t25 * p[20] - t23 * p[23];" }{TEXT 215 1 "\n" }{TEXT 215 121 "dtau[2][
2] = -t30 * p[2] + t32 * p[5] - t34 * p[8] + t36 * p[11] + t30 * p[14]
- t32 * p[17] + t34 * p[20] - t36 * p[23];" }}}{EXCHG {PARA 0 "> " 0
"" {MPLTEXT 1 0 2 "?C" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }
}}}
{MARK "0 0 0" 0 }{VIEWOPTS 1 1 0 125 100 1803 1 1 1 1 }{PAGENUMBERS 0
1 2 33 1 1 }
\ No newline at end of file
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment