Overture  Version 25
derivStencilCoeff.h
Go to the documentation of this file.
1 ! This file was generated by Overture/op/src/stencil.maple
2 ! Here are the coefficients of the ghost points (left side) in some stencils for derivatives.
3 ! Use -dr, -ds or -dt to get ghost points on the right side.
4 
5 
6 
7 ! *** orderOfAccuracy = 2
8 
9 #defineMacro xCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*aj ## rx/(dr) )
10 
11 #defineMacro xCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*aj ## sx/(ds) )
12 
13 #defineMacro yCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*aj ## ry/(dr) )
14 
15 #defineMacro yCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*aj ## sy/(ds) )
16 
17 #defineMacro lapCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(2.*aj ## rx**2-aj ## rxx*(dr)+2.*aj ## ry**2-aj ## ryy*(dr))/(dr)**2 )
18 
19 #defineMacro lapCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(-2.*aj ## sx**2+aj ## sxx*(ds)-2.*aj ## sy**2+aj ## syy*(ds))/(ds)**2 )
20 
21 #defineMacro xLapCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*(-2.*aj ## rx**3.*(ds)**2-6.*aj ## rx*aj ## sx**2.*(dr)**2-6.*aj ## rx*aj ## rxx*(dr)*(ds)**2+aj ## rxxx*(dr)**2.*(ds)**2-2.*aj ## ry**2.*aj ## rx*(ds)**2-4.*aj ## sy*(dr)**2.*aj ## ry*aj ## sx-2.*aj ## sy**2.*(dr)**2.*aj ## rx-2.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-4.*(dr)*(ds)**2.*aj ## ry*aj ## rxy+aj ## rxyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
22 
23 #defineMacro xLapCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## rx*(aj ## rx**2+aj ## ry**2)/(dr)**3 )
24 
25 #defineMacro xLapCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(-6.*aj ## rx**2.*aj ## sx*(ds)**2-2.*aj ## sx**3.*(dr)**2-6.*aj ## sxx*aj ## sx*(dr)**2.*(ds)+aj ## sxxx*(dr)**2.*(ds)**2-4.*aj ## ry*(ds)**2.*aj ## sy*aj ## rx-2.*aj ## ry**2.*(ds)**2.*aj ## sx-2.*aj ## sy**2.*aj ## sx*(dr)**2-2.*(dr)**2.*(ds)*aj ## syy*aj ## sx-4.*(dr)**2.*(ds)*aj ## sy*aj ## sxy+aj ## sxyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
26 
27 #defineMacro xLapCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sx*(aj ## sx**2+aj ## sy**2)/(ds)**3 )
28 
29 #defineMacro yLapCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*(-2.*aj ## ry*aj ## rx**2.*(ds)**2-4.*aj ## sx*(dr)**2.*aj ## sy*aj ## rx-2.*aj ## sx**2.*(dr)**2.*aj ## ry-4.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-2.*(dr)*(ds)**2.*aj ## ry*aj ## rxx+aj ## rxxy*(dr)**2.*(ds)**2-2.*aj ## ry**3.*(ds)**2-6.*aj ## ry*aj ## sy**2.*(dr)**2-6.*aj ## ry*aj ## ryy*(dr)*(ds)**2+aj ## ryyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
30 
31 #defineMacro yLapCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## ry*(aj ## rx**2+aj ## ry**2)/(dr)**3 )
32 
33 #defineMacro yLapCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(-4.*aj ## rx*(ds)**2.*aj ## ry*aj ## sx-2.*aj ## rx**2.*(ds)**2.*aj ## sy-2.*aj ## sy*aj ## sx**2.*(dr)**2-2.*(dr)**2.*(ds)*aj ## sy*aj ## sxx-4.*(dr)**2.*(ds)*aj ## sxy*aj ## sx+aj ## sxxy*(dr)**2.*(ds)**2-6.*aj ## ry**2.*aj ## sy*(ds)**2-2.*aj ## sy**3.*(dr)**2-6.*aj ## syy*aj ## sy*(dr)**2.*(ds)+aj ## syyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
34 
35 #defineMacro yLapCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sy*(aj ## sx**2+aj ## sy**2)/(ds)**3 )
36 
37 #defineMacro lapSqCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*(-6.*(dr)**2.*(ds)**2.*aj ## rxx**2-16.*(dr)**3.*aj ## sy*aj ## sx*aj ## rxy-12.*aj ## sx**2.*(dr)**3.*aj ## rxx-16.*(dr)*(ds)**2.*aj ## ry*aj ## rxy*aj ## rx-8.*(dr)**3.*aj ## sy*aj ## ry*aj ## sxx-4.*(dr)**3.*aj ## ryy*aj ## sx**2+24.*aj ## rx**2.*aj ## sx**2.*(dr)**2-12.*aj ## rx**2.*aj ## rxx*(dr)*(ds)**2+32.*(dr)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rx-12.*aj ## sy**2.*(dr)**3.*aj ## ryy+24.*aj ## ry**2.*aj ## sy**2.*(dr)**2-12.*aj ## ry**2.*aj ## ryy*(dr)*(ds)**2-4.*(dr)**3.*aj ## sy**2.*aj ## rxx+8.*(dr)**2.*aj ## sy**2.*aj ## rx**2+8.*(dr)**2.*aj ## ry**2.*aj ## sx**2+2.*aj ## rxxyy*(dr)**3.*(ds)**2+aj ## ryyyy*(dr)**3.*(ds)**2+8.*aj ## rx**4.*(ds)**2+8.*aj ## ry**4.*(ds)**2+aj ## rxxxx*(dr)**3.*(ds)**2-4.*(dr)*(ds)**2.*aj ## ry**2.*aj ## rxx-4.*(dr)*(ds)**2.*aj ## ryy*aj ## rx**2-8.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryyy-8.*(dr)**2.*(ds)**2.*aj ## rx*aj ## rxyy-4.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rxx-8.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxxy-8.*(dr)**2.*(ds)**2.*aj ## rx*aj ## rxxx-16.*(dr)**3.*aj ## sy*aj ## sxy*aj ## rx-16.*(dr)**3.*aj ## ry*aj ## sxy*aj ## sx-8.*(dr)**3.*aj ## syy*aj ## sx*aj ## rx-24.*aj ## sx*(dr)**3.*aj ## sxx*aj ## rx-24.*aj ## sy*(dr)**3.*aj ## syy*aj ## ry-6.*(dr)**2.*(ds)**2.*aj ## ryy**2-8.*(dr)**2.*(ds)**2.*aj ## rxy**2+16.*aj ## ry**2.*aj ## rx**2.*(ds)**2)/(dr)**4/(ds)**2 )
38 
39 #defineMacro lapSqCoeff2dOrder2Ghost20(aj,dr,ds) ( (-4.*(dr)*aj ## ry*aj ## rxy*aj ## rx-(dr)*aj ## ry**2.*aj ## rxx-(dr)*aj ## ryy*aj ## rx**2+aj ## rx**4+aj ## ry**4-3.*aj ## ry**2.*aj ## ryy*(dr)-3.*aj ## rx**2.*aj ## rxx*(dr)+2.*aj ## ry**2.*aj ## rx**2)/(dr)**4 )
40 
41 #defineMacro lapSqCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(16.*aj ## sy**2.*aj ## sx**2.*(dr)**2-4.*(dr)**2.*(ds)**2.*aj ## syy*aj ## sxx+24.*aj ## rx**2.*aj ## sx**2.*(ds)**2+24.*aj ## ry**2.*aj ## sy**2.*(ds)**2-12.*aj ## sy**2.*aj ## syy*(ds)*(dr)**2-12.*aj ## sx**2.*aj ## sxx*(ds)*(dr)**2+aj ## syyyy*(ds)**3.*(dr)**2-8.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxxy+32.*(ds)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rx-8.*(dr)**2.*(ds)**2.*aj ## sxxx*aj ## sx-6.*(dr)**2.*(ds)**2.*aj ## sxx**2-8.*(dr)**2.*(ds)**2.*aj ## sx*aj ## sxyy-4.*(dr)**2.*(ds)*aj ## syy*aj ## sx**2-4.*(dr)**2.*(ds)*aj ## sy**2.*aj ## sxx-4.*(ds)**3.*aj ## syy*aj ## rx**2+8.*aj ## sy**4.*(dr)**2+8.*aj ## sx**4.*(dr)**2+2.*aj ## sxxyy*(ds)**3.*(dr)**2+aj ## sxxxx*(ds)**3.*(dr)**2-8.*(dr)**2.*(ds)**2.*aj ## syyy*aj ## sy-6.*(dr)**2.*(ds)**2.*aj ## syy**2-8.*(dr)**2.*(ds)**2.*aj ## sxy**2+8.*(ds)**2.*aj ## sy**2.*aj ## rx**2+8.*(ds)**2.*aj ## ry**2.*aj ## sx**2-4.*(ds)**3.*aj ## ry**2.*aj ## sxx-12.*aj ## rx**2.*(ds)**3.*aj ## sxx-12.*aj ## ry**2.*(ds)**3.*aj ## syy-16.*(dr)**2.*(ds)*aj ## sy*aj ## sxy*aj ## sx-16.*(ds)**3.*aj ## sy*aj ## rxy*aj ## rx-8.*(ds)**3.*aj ## sy*aj ## ry*aj ## rxx-8.*(ds)**3.*aj ## ryy*aj ## sx*aj ## rx-16.*(ds)**3.*aj ## ry*aj ## sx*aj ## rxy-16.*(ds)**3.*aj ## ry*aj ## sxy*aj ## rx-24.*aj ## rx*(ds)**3.*aj ## sx*aj ## rxx-24.*aj ## ry*(ds)**3.*aj ## sy*aj ## ryy)/(ds)**4/(dr)**2 )
42 
43 #defineMacro lapSqCoeff2dOrder2Ghost02(aj,dr,ds) ( (-4.*(ds)*aj ## sy*aj ## sxy*aj ## sx-(ds)*aj ## sy**2.*aj ## sxx-(ds)*aj ## syy*aj ## sx**2+aj ## sy**4+aj ## sx**4-3.*aj ## sy**2.*aj ## syy*(ds)-3.*aj ## sx**2.*aj ## sxx*(ds)+2.*aj ## sy**2.*aj ## sx**2)/(ds)**4 )
44 
45 #defineMacro xLapSqCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(-96.*(dr)**3.*(ds)**2.*aj ## sy*aj ## rxy*aj ## sx*aj ## rx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## ryy*aj ## sy*aj ## sx-48.*(dr)**3.*(ds)**2.*aj ## syy*aj ## sy*aj ## ry*aj ## rx-48.*aj ## rx**2.*(dr)**2.*(ds)**2.*aj ## ry*aj ## sy*aj ## sx-96.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rx*aj ## sxy*aj ## sx-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rxx-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxx*aj ## rx-10.*aj ## ry**2.*aj ## rx**3.*(ds)**4-5.*aj ## ry**4.*aj ## rx*(ds)**4-80.*aj ## rxx*aj ## rx**3.*(dr)*(ds)**4+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rx*aj ## rxxy+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## rxy+8.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## ry*aj ## rx-5.*aj ## rx**5.*(ds)**4+2.*(dr)**3.*(ds)**4.*aj ## ryyyy*aj ## rx-48.*aj ## rx*(dr)*(ds)**4.*aj ## ry**2.*aj ## rxx-96.*aj ## rx**2.*(dr)*(ds)**4.*aj ## ry*aj ## rxy-48.*aj ## ry**2.*(dr)*(ds)**4.*aj ## ryy*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## rx*aj ## rxx+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxy*aj ## rxx+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ryy*aj ## sx+8.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxyyy+8.*(dr)**4.*(ds)**2.*aj ## ryyy*aj ## sy*aj ## sx+8.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxxy+4.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxxx+12.*(dr)**3.*(ds)**4.*aj ## rxx*aj ## rxyy+24.*(dr)**3.*(ds)**4.*aj ## rxy*aj ## rxxy+24.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sy*aj ## rxy+12.*(dr)**3.*(ds)**4.*aj ## rx*aj ## rxxyy+12.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxyy+8.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## rxy+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sxy+20.*(dr)**3.*(ds)**4.*aj ## rxx*aj ## rxxx+8.*(dr)**4.*(ds)**2.*aj ## ry*aj ## syyy*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ry*aj ## sxy+10.*(dr)**3.*(ds)**4.*aj ## rx*aj ## rxxxx+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxyy+8.*(dr)**4.*(ds)**2.*aj ## sy*aj ## syyy*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxyy+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sxx*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sx*aj ## rxx+4.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxxx+24.*(dr)**4.*(ds)**2.*aj ## sxyy*aj ## sx*aj ## rx+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxxy+12.*(dr)**2.*(ds)**4.*aj ## rxyy*aj ## rx**2+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rxy+6.*(dr)**2.*(ds)**4.*aj ## ryy**2.*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## sxx*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rxx+24.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxy*aj ## sxx+24.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxxy*aj ## sx+8.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxxx+40.*(dr)**4.*(ds)**2.*aj ## sxxx*aj ## sx*aj ## rx+30.*aj ## rx*(dr)**2.*(ds)**4.*aj ## rxx**2+60.*(dr)**4.*(ds)**2.*aj ## rxx*aj ## sx*aj ## sxx+20.*aj ## rx**2.*(dr)**2.*(ds)**4.*aj ## rxxx+48.*(dr)**4.*(ds)**2.*aj ## rxy*aj ## sxy*aj ## sx+24.*(dr)**2.*(ds)**4.*aj ## rxy**2.*aj ## rx+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxy*aj ## rx-24.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## sy**2.*aj ## rx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy**2.*aj ## rxy-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry**2.*aj ## sxy-24.*(dr)**3.*(ds)**2.*aj ## syy*aj ## ry**2.*aj ## sx+12.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxyy-24.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## sxx*aj ## sx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxy*aj ## sx**2-24.*(dr)**3.*(ds)**2.*aj ## sy**2.*aj ## rx*aj ## rxx-16.*aj ## rx**3.*(dr)*(ds)**4.*aj ## ryy-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rx**2-32.*aj ## ry**3.*(dr)*(ds)**4.*aj ## rxy-24.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## rx*aj ## sx**2-24.*(dr)**3.*(ds)**2.*aj ## syy*aj ## rx**2.*aj ## sx+12.*(dr)**4.*(ds)**2.*aj ## rxyy*aj ## sx**2-120.*aj ## sx*aj ## rx**2.*(dr)**3.*(ds)**2.*aj ## sxx-120.*aj ## sx**2.*aj ## rx*(dr)**3.*(ds)**2.*aj ## rxx+6.*(dr)**4.*(ds)**2.*aj ## syy**2.*aj ## rx+24.*(dr)**4.*(ds)**2.*aj ## sxy**2.*aj ## rx+4.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxxx-24.*aj ## sy**2.*aj ## ry**2.*(dr)**2.*(ds)**2.*aj ## rx-16.*aj ## sy*aj ## ry**3.*(dr)**2.*(ds)**2.*aj ## sx-24.*aj ## rx*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## sx**2+30.*(dr)**4.*(ds)**2.*aj ## rx*aj ## sxx**2-6.*aj ## sy**4.*(dr)**4.*aj ## rx+20.*(dr)**4.*(ds)**2.*aj ## rxxx*aj ## sx**2-8.*aj ## rx**3.*(dr)**2.*(ds)**2.*aj ## sy**2-24.*aj ## sy**3.*(dr)**4.*aj ## ry*aj ## sx-24.*aj ## sy*aj ## sx**3.*(dr)**4.*aj ## ry-36.*aj ## sy**2.*aj ## sx**2.*(dr)**4.*aj ## rx-40.*aj ## sx**2.*aj ## rx**3.*(dr)**2.*(ds)**2-30.*aj ## sx**4.*aj ## rx*(dr)**4)/(dr)**5/(ds)**4 )
46 
47 #defineMacro xLapSqCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*(-20.*aj ## sx**2.*aj ## rx**3.*(dr)**2-12.*aj ## rx*(dr)**2.*aj ## ry**2.*aj ## sx**2-20.*aj ## rxx*aj ## rx**3.*(dr)*(ds)**2-4.*aj ## ry**4.*aj ## rx*(ds)**2-8.*aj ## ry**2.*aj ## rx**3.*(ds)**2+6.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rx*aj ## rxx+12.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxy*aj ## rxx+12.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## rxxy+6.*(dr)**2.*(ds)**2.*aj ## rxyy*aj ## rx**2+12.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryy*aj ## rxy+4.*(dr)**2.*(ds)**2.*aj ## ryyy*aj ## ry*aj ## rx-24.*aj ## rx**2.*(dr)**2.*aj ## ry*aj ## sy*aj ## sx-12.*aj ## rx*(dr)*(ds)**2.*aj ## ry**2.*aj ## rxx-24.*aj ## rx**2.*(dr)*(ds)**2.*aj ## ry*aj ## rxy-12.*aj ## ry**2.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-4.*aj ## rx**3.*(dr)**2.*aj ## sy**2-4.*aj ## rx**5.*(ds)**2-12.*aj ## sy**2.*aj ## ry**2.*(dr)**2.*aj ## rx-8.*aj ## sy*aj ## ry**3.*(dr)**2.*aj ## sx+2.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxxx+3.*(dr)**2.*(ds)**2.*aj ## ryy**2.*aj ## rx+6.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxyy+12.*(dr)**2.*(ds)**2.*aj ## rxy**2.*aj ## rx+15.*aj ## rx*(dr)**2.*(ds)**2.*aj ## rxx**2+10.*aj ## rx**2.*(dr)**2.*(ds)**2.*aj ## rxxx-8.*aj ## ry**3.*(dr)*(ds)**2.*aj ## rxy-4.*aj ## rx**3.*(dr)*(ds)**2.*aj ## ryy)/(dr)**5/(ds)**2 )
48 
49 #defineMacro xLapSqCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## rx*(2.*aj ## ry**2.*aj ## rx**2+aj ## ry**4+aj ## rx**4)/(dr)**5 )
50 
51 #defineMacro xLapSqCoeff2dOrder2Ghost01(aj,dr,ds) ( 1/2.*(-96.*(dr)**2.*(ds)**3.*aj ## sy*aj ## rxy*aj ## sx*aj ## rx+24.*(ds)**3.*(dr)**4.*aj ## sxy*aj ## sxxy+12.*(ds)**3.*(dr)**4.*aj ## sxyy*aj ## sxx+8.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxyyy+4.*(ds)**3.*(dr)**4.*aj ## syy*aj ## sxxx+8.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxxxy+12.*(ds)**3.*(dr)**4.*aj ## sxxyy*aj ## sx+20.*(ds)**3.*(dr)**4.*aj ## sxx*aj ## sxxx+10.*(ds)**3.*(dr)**4.*aj ## sxxxx*aj ## sx+12.*(ds)**2.*(dr)**4.*aj ## sy**2.*aj ## sxyy+6.*(ds)**2.*(dr)**4.*aj ## syy**2.*aj ## sx+12.*(ds)**2.*(dr)**4.*aj ## sxyy*aj ## sx**2+4.*(ds)**2.*(dr)**4.*aj ## sy**2.*aj ## sxxx+12.*(ds)**2.*(dr)**4.*aj ## syy*aj ## sxx*aj ## sx+12.*(ds)**4.*(dr)**2.*aj ## ryy*aj ## sxx*aj ## rx+8.*(ds)**2.*(dr)**4.*aj ## syyy*aj ## sy*aj ## sx-16.*aj ## sx**3.*(ds)*(dr)**4.*aj ## syy+30.*(ds)**4.*(dr)**2.*aj ## sx*aj ## rxx**2-32.*aj ## sy**3.*(ds)*(dr)**4.*aj ## sxy-5.*aj ## sx**5.*(dr)**4-5.*aj ## sy**4.*aj ## sx*(dr)**4-10.*aj ## sy**2.*aj ## sx**3.*(dr)**4-30.*aj ## sx*aj ## rx**4.*(ds)**4-80.*aj ## sxx*aj ## sx**3.*(ds)*(dr)**4-40.*aj ## sx**3.*aj ## rx**2.*(ds)**2.*(dr)**2-48.*aj ## sy**2.*(ds)*(dr)**4.*aj ## syy*aj ## sx+40.*(ds)**4.*(dr)**2.*aj ## rxxx*aj ## sx*aj ## rx+60.*(ds)**4.*(dr)**2.*aj ## sxx*aj ## rx*aj ## rxx+4.*(ds)**4.*(dr)**2.*aj ## ry**2.*aj ## sxxx-96.*aj ## sx**2.*(ds)*(dr)**4.*aj ## sy*aj ## sxy-48.*aj ## sx*(ds)*(dr)**4.*aj ## sy**2.*aj ## sxx+20.*(ds)**4.*(dr)**2.*aj ## sxxx*aj ## rx**2+12.*(ds)**4.*(dr)**2.*aj ## syy*aj ## rx*aj ## rxx+24.*(ds)**4.*(dr)**2.*aj ## sy*aj ## rxy*aj ## rxx+24.*(ds)**4.*(dr)**2.*aj ## sy*aj ## rx*aj ## rxxy+12.*(ds)**4.*(dr)**2.*aj ## sxyy*aj ## rx**2+48.*(ds)**4.*(dr)**2.*aj ## sxy*aj ## rxy*aj ## rx+24.*(ds)**4.*(dr)**2.*aj ## rxyy*aj ## sx*aj ## rx+24.*(ds)**4.*(dr)**2.*aj ## sx*aj ## rxy**2+8.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sy*aj ## rxxx+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sxxy*aj ## rx+24.*(ds)**2.*(dr)**4.*aj ## sxy**2.*aj ## sx+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sx*aj ## rxxy+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sxy*aj ## rxx-24.*aj ## ry**3.*(ds)**4.*aj ## sy*aj ## rx+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sxx*aj ## rxy+12.*(ds)**4.*(dr)**2.*aj ## ryy*aj ## sx*aj ## rxx-36.*aj ## ry**2.*aj ## rx**2.*(ds)**4.*aj ## sx+8.*(ds)**4.*(dr)**2.*aj ## ryyy*aj ## sy*aj ## rx+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sy*aj ## rxyy+12.*(ds)**4.*(dr)**2.*aj ## ry**2.*aj ## sxyy+8.*(ds)**4.*(dr)**2.*aj ## ryyy*aj ## ry*aj ## sx+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## syy*aj ## rxy+24.*(ds)**4.*(dr)**2.*aj ## ryy*aj ## sy*aj ## rxy+8.*(ds)**4.*(dr)**2.*aj ## syyy*aj ## ry*aj ## rx-24.*aj ## ry*aj ## rx**3.*(ds)**4.*aj ## sy-6.*aj ## ry**4.*(ds)**4.*aj ## sx+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## ryy*aj ## sxy+12.*(ds)**4.*(dr)**2.*aj ## ryy*aj ## syy*aj ## rx+6.*(ds)**4.*(dr)**2.*aj ## ryy**2.*aj ## sx-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxy*aj ## rx**2-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry*aj ## sxx*aj ## rx-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry*aj ## sx*aj ## rxx-24.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## rx*aj ## rxx-96.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rx*aj ## sxy*aj ## sx-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rxy*aj ## sx**2-24.*(dr)**2.*(ds)**3.*aj ## ry**2.*aj ## sxx*aj ## sx-24.*(dr)**2.*(ds)**3.*aj ## syy*aj ## rx**2.*aj ## sx+30.*aj ## sx*(ds)**2.*(dr)**4.*aj ## sxx**2-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry**2.*aj ## sxy-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## rxy-8.*aj ## sx**3.*(ds)**2.*(dr)**2.*aj ## ry**2-48.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sy*aj ## ry*aj ## rx-24.*(dr)**2.*(ds)**3.*aj ## syy*aj ## ry**2.*aj ## sx-24.*(dr)**2.*(ds)**3.*aj ## ryy*aj ## rx*aj ## sx**2-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## ryy*aj ## sy*aj ## sx-24.*(dr)**2.*(ds)**3.*aj ## ryy*aj ## sy**2.*aj ## rx-120.*aj ## sx**2.*aj ## rx*(ds)**3.*(dr)**2.*aj ## rxx-24.*aj ## sx*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## rx**2+20.*aj ## sx**2.*(ds)**2.*(dr)**4.*aj ## sxxx-120.*aj ## sx*aj ## rx**2.*(ds)**3.*(dr)**2.*aj ## sxx-48.*aj ## sx**2.*(ds)**2.*(dr)**2.*aj ## ry*aj ## sy*aj ## rx-16.*aj ## ry*aj ## sy**3.*(ds)**2.*(dr)**2.*aj ## rx-24.*aj ## ry**2.*aj ## sy**2.*(ds)**2.*(dr)**2.*aj ## sx+24.*(ds)**2.*(dr)**4.*aj ## syy*aj ## sy*aj ## sxy+12.*(ds)**3.*(dr)**4.*aj ## syy*aj ## sxyy+8.*(ds)**3.*(dr)**4.*aj ## syyy*aj ## sxy+2.*(ds)**3.*(dr)**4.*aj ## syyyy*aj ## sx+24.*(ds)**2.*(dr)**4.*aj ## sy*aj ## sxxy*aj ## sx+24.*(ds)**2.*(dr)**4.*aj ## sy*aj ## sxy*aj ## sxx)/(ds)**5/(dr)**4 )
52 
53 #defineMacro xLapSqCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*(3.*(ds)**2.*(dr)**2.*aj ## syy**2.*aj ## sx-8.*aj ## sy**3.*(ds)*(dr)**2.*aj ## sxy-8.*aj ## sy**2.*aj ## sx**3.*(dr)**2-20.*aj ## sx**3.*aj ## rx**2.*(ds)**2-4.*aj ## sy**4.*aj ## sx*(dr)**2-4.*aj ## sx**5.*(dr)**2+4.*(ds)**2.*(dr)**2.*aj ## syyy*aj ## sy*aj ## sx+12.*(ds)**2.*(dr)**2.*aj ## syy*aj ## sy*aj ## sxy-12.*aj ## sy**2.*(ds)*(dr)**2.*aj ## syy*aj ## sx-20.*aj ## sxx*aj ## sx**3.*(ds)*(dr)**2+12.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sxxy*aj ## sx+12.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sxy*aj ## sxx+6.*(ds)**2.*(dr)**2.*aj ## syy*aj ## sxx*aj ## sx+12.*(ds)**2.*(dr)**2.*aj ## sxy**2.*aj ## sx+2.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## sxxx+6.*(ds)**2.*(dr)**2.*aj ## sxyy*aj ## sx**2-24.*aj ## sx**2.*(ds)**2.*aj ## ry*aj ## sy*aj ## rx-4.*aj ## sx**3.*(ds)**2.*aj ## ry**2-24.*aj ## sx**2.*(ds)*(dr)**2.*aj ## sy*aj ## sxy-12.*aj ## sx*(ds)*(dr)**2.*aj ## sy**2.*aj ## sxx-12.*aj ## sx*(ds)**2.*aj ## sy**2.*aj ## rx**2-4.*aj ## sx**3.*(ds)*(dr)**2.*aj ## syy+10.*aj ## sx**2.*(ds)**2.*(dr)**2.*aj ## sxxx+15.*aj ## sx*(ds)**2.*(dr)**2.*aj ## sxx**2+6.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## sxyy-8.*aj ## ry*aj ## sy**3.*(ds)**2.*aj ## rx-12.*aj ## ry**2.*aj ## sy**2.*(ds)**2.*aj ## sx)/(ds)**5/(dr)**2 )
54 
55 #defineMacro xLapSqCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sx*(aj ## sy**4+2.*aj ## sy**2.*aj ## sx**2+aj ## sx**4)/(ds)**5 )
56 
57 #defineMacro yLapSqCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(-5.*aj ## ry**5.*(ds)**4-16.*aj ## ry**3.*(dr)*(ds)**4.*aj ## rxx-30.*aj ## sy**4.*aj ## ry*(dr)**4-10.*aj ## ry**3.*aj ## rx**2.*(ds)**4-5.*aj ## ry*aj ## rx**4.*(ds)**4-40.*aj ## sy**2.*aj ## ry**3.*(dr)**2.*(ds)**2-80.*aj ## ry**3.*aj ## ryy*(dr)*(ds)**4-96.*aj ## ry**2.*(dr)*(ds)**4.*aj ## rxy*aj ## rx-48.*aj ## ry*(dr)*(ds)**4.*aj ## ryy*aj ## rx**2-48.*aj ## rx**2.*(dr)*(ds)**4.*aj ## ry*aj ## rxx+24.*(dr)**4.*(ds)**2.*aj ## rxx*aj ## sxy*aj ## sx+40.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sy*aj ## ry-32.*aj ## rx**3.*(dr)*(ds)**4.*aj ## rxy+12.*(dr)**4.*(ds)**2.*aj ## rxxy*aj ## sx**2+60.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## syy*aj ## sy+20.*(dr)**4.*(ds)**2.*aj ## ryyy*aj ## sy**2+30.*(dr)**4.*(ds)**2.*aj ## syy**2.*aj ## ry+8.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxxx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rxx+8.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxxx*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## sxxy*aj ## sx*aj ## rx+8.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxx*aj ## rx+24.*(dr)**4.*(ds)**2.*aj ## sxx*aj ## sxy*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxxy-120.*aj ## ry**2.*aj ## sy*(dr)**3.*(ds)**2.*aj ## syy-48.*aj ## ry**2.*(dr)**2.*(ds)**2.*aj ## sy*aj ## rx*aj ## sx-24.*aj ## ry*(dr)**2.*(ds)**2.*aj ## sy**2.*aj ## rx**2-24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## ryy*aj ## sx**2-120.*aj ## ry*aj ## sy**2.*(dr)**3.*(ds)**2.*aj ## ryy+24.*(dr)**4.*(ds)**2.*aj ## rxy*aj ## sx*aj ## sxx+6.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxx**2-8.*aj ## ry**3.*(dr)**2.*(ds)**2.*aj ## sx**2-24.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry**2.*aj ## sxx-24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy**2.*aj ## rxx-48.*(dr)**3.*(ds)**2.*aj ## sy**2.*aj ## rxy*aj ## rx-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sx*aj ## rx-96.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy*aj ## sx*aj ## rxy-48.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## sxy*aj ## sx-24.*aj ## sy**3.*aj ## sx*(dr)**4.*aj ## rx-24.*(dr)**3.*(ds)**2.*aj ## syy*aj ## sy*aj ## rx**2-96.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxy*aj ## rx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## syy*aj ## sx*aj ## rx-24.*aj ## sx**3.*(dr)**4.*aj ## sy*aj ## rx-48.*(dr)**3.*(ds)**2.*aj ## sxy*aj ## rx**2.*aj ## sx-48.*(dr)**3.*(ds)**2.*aj ## rxy*aj ## rx*aj ## sx**2-24.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rx**2-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sx*aj ## rx*aj ## rxx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rx*aj ## sxx*aj ## sx-24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxx*aj ## sx**2-36.*aj ## sy**2.*aj ## sx**2.*(dr)**4.*aj ## ry+48.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rxy+4.*(dr)**4.*(ds)**2.*aj ## ryyy*aj ## sx**2+20.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## ryy+12.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxxy+24.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxy**2+12.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxyy+24.*(dr)**3.*(ds)**4.*aj ## rxy*aj ## rxyy+8.*(dr)**3.*(ds)**4.*aj ## rx*aj ## rxyyy+12.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxxy+10.*(dr)**3.*(ds)**4.*aj ## ry*aj ## ryyyy+24.*(dr)**2.*(ds)**4.*aj ## rxy*aj ## rxx*aj ## rx+8.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rx*aj ## rxxx+4.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## rxx+12.*(dr)**2.*(ds)**4.*aj ## rxxy*aj ## rx**2+6.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxx**2+12.*(dr)**3.*(ds)**4.*aj ## rxxy*aj ## rxx+8.*(dr)**3.*(ds)**4.*aj ## rxxxy*aj ## rx+20.*aj ## ry**2.*(dr)**2.*(ds)**4.*aj ## ryyy+30.*aj ## ry*(dr)**2.*(ds)**4.*aj ## ryy**2+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxy**2+4.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## rx**2+2.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxxx+8.*(dr)**3.*(ds)**4.*aj ## rxxx*aj ## rxy+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rx*aj ## rxyy+24.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## rxy*aj ## rx-16.*aj ## rx**3.*aj ## sx*(dr)**2.*(ds)**2.*aj ## sy-24.*aj ## rx**2.*aj ## sx**2.*(dr)**2.*(ds)**2.*aj ## ry+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## rxx+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxyy*aj ## rx+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxyy+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sy*aj ## rxx+24.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxxy+24.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sx*aj ## rxy+24.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sxy*aj ## rx+24.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## sxy*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sx*aj ## sxyy+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sxx+8.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sx*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ry*aj ## sxx-6.*aj ## sx**4.*(dr)**4.*aj ## ry)/(dr)**5/(ds)**4 )
58 
59 #defineMacro yLapSqCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*(-20.*aj ## ry**3.*aj ## ryy*(dr)*(ds)**2-4.*aj ## ry**5.*(ds)**2+12.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxy**2+2.*(dr)**2.*(ds)**2.*aj ## ryyy*aj ## rx**2-8.*aj ## ry**3.*aj ## rx**2.*(ds)**2-4.*aj ## ry*aj ## rx**4.*(ds)**2-20.*aj ## sy**2.*aj ## ry**3.*(dr)**2+6.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxxy+12.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rxy*aj ## rx-8.*aj ## rx**3.*(dr)*(ds)**2.*aj ## rxy+12.*(dr)**2.*(ds)**2.*aj ## rxy*aj ## rxx*aj ## rx+4.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## rxxx-24.*aj ## ry**2.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-4.*aj ## ry**3.*(dr)*(ds)**2.*aj ## rxx-8.*aj ## rx**3.*aj ## sx*(dr)**2.*aj ## sy+6.*(dr)**2.*(ds)**2.*aj ## rxxy*aj ## rx**2+10.*aj ## ry**2.*(dr)**2.*(ds)**2.*aj ## ryyy+15.*aj ## ry*(dr)**2.*(ds)**2.*aj ## ryy**2+3.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxx**2+6.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryy*aj ## rxx+12.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## rxyy-12.*aj ## rx**2.*aj ## sx**2.*(dr)**2.*aj ## ry-24.*aj ## ry**2.*(dr)**2.*aj ## sy*aj ## rx*aj ## sx-4.*aj ## ry**3.*(dr)**2.*aj ## sx**2-12.*aj ## ry*(dr)*(ds)**2.*aj ## ryy*aj ## rx**2-12.*aj ## ry*(dr)**2.*aj ## sy**2.*aj ## rx**2-12.*aj ## rx**2.*(dr)*(ds)**2.*aj ## ry*aj ## rxx)/(dr)**5/(ds)**2 )
60 
61 #defineMacro yLapSqCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## ry*(2.*aj ## ry**2.*aj ## rx**2+aj ## ry**4+aj ## rx**4)/(dr)**5 )
62 
63 #defineMacro yLapSqCoeff2dOrder2Ghost01(aj,dr,ds) ( 1/2.*(8.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sx*aj ## rxxx+24.*(ds)**4.*(dr)**2.*aj ## sxx*aj ## rxy*aj ## rx-120.*aj ## ry*aj ## sy**2.*(ds)**3.*(dr)**2.*aj ## ryy-96.*aj ## sy**2.*(ds)*(dr)**4.*aj ## sxy*aj ## sx-6.*aj ## rx**4.*(ds)**4.*aj ## sy+6.*(ds)**4.*(dr)**2.*aj ## sy*aj ## rxx**2+12.*(ds)**4.*(dr)**2.*aj ## sxxy*aj ## rx**2-5.*aj ## sy**5.*(dr)**4+20.*(ds)**3.*(dr)**4.*aj ## syy*aj ## syyy-10.*aj ## sy**3.*aj ## sx**2.*(dr)**4-30.*aj ## sy*aj ## ry**4.*(ds)**4-80.*aj ## sy**3.*aj ## syy*(ds)*(dr)**4+4.*(ds)**4.*(dr)**2.*aj ## syyy*aj ## rx**2-40.*aj ## sy**3.*aj ## ry**2.*(ds)**2.*(dr)**2+24.*(ds)**4.*(dr)**2.*aj ## ryy*aj ## sx*aj ## rxy+12.*(ds)**4.*(dr)**2.*aj ## ry**2.*aj ## sxxy+24.*(ds)**4.*(dr)**2.*aj ## sy*aj ## rxy**2-5.*aj ## sy*aj ## sx**4.*(dr)**4+30.*(ds)**4.*(dr)**2.*aj ## sy*aj ## ryy**2+20.*(ds)**4.*(dr)**2.*aj ## syyy*aj ## ry**2+8.*(ds)**4.*(dr)**2.*aj ## ryyy*aj ## sx*aj ## rx+48.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sxy*aj ## rxy+12.*(ds)**4.*(dr)**2.*aj ## ry*aj ## ryy*aj ## sxx+24.*(ds)**4.*(dr)**2.*aj ## sy*aj ## rx*aj ## rxyy+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sx*aj ## rxyy+24.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sxyy*aj ## rx+24.*(ds)**4.*(dr)**2.*aj ## syy*aj ## rxy*aj ## rx+12.*(ds)**4.*(dr)**2.*aj ## sy*aj ## ryy*aj ## rxx+24.*(ds)**4.*(dr)**2.*aj ## sy*aj ## ry*aj ## rxxy+12.*(ds)**4.*(dr)**2.*aj ## syy*aj ## ry*aj ## rxx+24.*(ds)**4.*(dr)**2.*aj ## ryy*aj ## sxy*aj ## rx+60.*(ds)**4.*(dr)**2.*aj ## syy*aj ## ry*aj ## ryy+40.*(ds)**4.*(dr)**2.*aj ## ryyy*aj ## sy*aj ## ry+24.*(ds)**4.*(dr)**2.*aj ## rxxy*aj ## sx*aj ## rx+24.*(ds)**4.*(dr)**2.*aj ## sxy*aj ## rxx*aj ## rx+8.*(ds)**4.*(dr)**2.*aj ## sy*aj ## rx*aj ## rxxx+24.*(ds)**4.*(dr)**2.*aj ## rxy*aj ## sx*aj ## rxx+8.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sxxx*aj ## rx+2.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxxxx+8.*(ds)**3.*(dr)**4.*aj ## sx*aj ## sxyyy+4.*(ds)**3.*(dr)**4.*aj ## syyy*aj ## sxx+12.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxxyy+24.*(ds)**3.*(dr)**4.*aj ## sxyy*aj ## sxy-24.*aj ## ry**3.*aj ## rx*(ds)**4.*aj ## sx-24.*aj ## rx**3.*(ds)**4.*aj ## ry*aj ## sx+8.*(ds)**3.*(dr)**4.*aj ## sxxx*aj ## sxy+8.*(ds)**3.*(dr)**4.*aj ## sxxxy*aj ## sx+12.*(ds)**3.*(dr)**4.*aj ## sxx*aj ## sxxy+12.*(ds)**4.*(dr)**2.*aj ## ry*aj ## sxx*aj ## rxx-36.*aj ## ry**2.*aj ## rx**2.*(ds)**4.*aj ## sy+10.*(ds)**3.*(dr)**4.*aj ## sy*aj ## syyyy+12.*(ds)**3.*(dr)**4.*aj ## syy*aj ## sxxy-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sx*aj ## rx*aj ## rxx-24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxx*aj ## rx**2-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rx*aj ## sxx*aj ## sx-24.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rxx*aj ## sx**2-24.*(dr)**2.*(ds)**3.*aj ## ry*aj ## ryy*aj ## sx**2-120.*aj ## ry**2.*aj ## sy*(ds)**3.*(dr)**2.*aj ## syy-48.*(dr)**2.*(ds)**3.*aj ## sxy*aj ## rx**2.*aj ## sx-48.*(dr)**2.*(ds)**3.*aj ## rxy*aj ## rx*aj ## sx**2-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ryy*aj ## sx*aj ## rx-24.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sy*aj ## rx**2-96.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy*aj ## sxy*aj ## rx-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## syy*aj ## sx*aj ## rx-96.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy*aj ## sx*aj ## rxy-24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry**2.*aj ## sxx-24.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## rxx-48.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## rxy*aj ## rx+12.*(ds)**2.*(dr)**4.*aj ## sxxy*aj ## sx**2-48.*(dr)**2.*(ds)**3.*aj ## ry**2.*aj ## sxy*aj ## sx+24.*(ds)**2.*(dr)**4.*aj ## sxy*aj ## sx*aj ## sxx+8.*(ds)**2.*(dr)**4.*aj ## sy*aj ## sxxx*aj ## sx+6.*(ds)**2.*(dr)**4.*aj ## sy*aj ## sxx**2+12.*(ds)**2.*(dr)**4.*aj ## sy**2.*aj ## sxxy+24.*(ds)**2.*(dr)**4.*aj ## sy*aj ## sxy**2+24.*(ds)**2.*(dr)**4.*aj ## sy*aj ## sx*aj ## sxyy+12.*(ds)**2.*(dr)**4.*aj ## sy*aj ## syy*aj ## sxx+4.*(ds)**2.*(dr)**4.*aj ## syyy*aj ## sx**2+20.*aj ## sy**2.*(ds)**2.*(dr)**4.*aj ## syyy+24.*(ds)**2.*(dr)**4.*aj ## syy*aj ## sxy*aj ## sx+30.*aj ## sy*(ds)**2.*(dr)**4.*aj ## syy**2-8.*aj ## sy**3.*(ds)**2.*(dr)**2.*aj ## rx**2-32.*aj ## sx**3.*(ds)*(dr)**4.*aj ## sxy-16.*aj ## sy**3.*(ds)*(dr)**4.*aj ## sxx-24.*aj ## sy*(ds)**2.*(dr)**2.*aj ## ry**2.*aj ## sx**2-48.*aj ## sy**2.*(ds)**2.*(dr)**2.*aj ## ry*aj ## rx*aj ## sx-24.*aj ## rx**2.*aj ## sx**2.*(ds)**2.*(dr)**2.*aj ## sy-16.*aj ## rx*aj ## sx**3.*(ds)**2.*(dr)**2.*aj ## ry-48.*aj ## sx**2.*(ds)*(dr)**4.*aj ## sy*aj ## sxx-48.*aj ## sy*(ds)*(dr)**4.*aj ## syy*aj ## sx**2)/(ds)**5/(dr)**4 )
64 
65 #defineMacro yLapSqCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*(12.*(ds)**2.*(dr)**2.*aj ## sxy*aj ## sx*aj ## sxx-12.*aj ## sy*(ds)**2.*aj ## ry**2.*aj ## sx**2+10.*aj ## sy**2.*(ds)**2.*(dr)**2.*aj ## syyy+15.*aj ## sy*(ds)**2.*(dr)**2.*aj ## syy**2-4.*aj ## sy**5.*(dr)**2-4.*aj ## sy**3.*(ds)**2.*aj ## rx**2-8.*aj ## sy**3.*aj ## sx**2.*(dr)**2-20.*aj ## sy**3.*aj ## ry**2.*(ds)**2-20.*aj ## sy**3.*aj ## syy*(ds)*(dr)**2-4.*aj ## sy*aj ## sx**4.*(dr)**2-24.*aj ## sy**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## sx+12.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sx*aj ## sxyy+6.*(ds)**2.*(dr)**2.*aj ## sy*aj ## syy*aj ## sxx+12.*(ds)**2.*(dr)**2.*aj ## syy*aj ## sxy*aj ## sx-24.*aj ## sy**2.*(ds)*(dr)**2.*aj ## sxy*aj ## sx-12.*aj ## sy*(ds)*(dr)**2.*aj ## syy*aj ## sx**2-4.*aj ## sy**3.*(ds)*(dr)**2.*aj ## sxx-12.*aj ## sx**2.*(ds)*(dr)**2.*aj ## sy*aj ## sxx-12.*aj ## rx**2.*aj ## sx**2.*(ds)**2.*aj ## sy-8.*aj ## rx*aj ## sx**3.*(ds)**2.*aj ## ry+6.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## sxxy+12.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sxy**2+2.*(ds)**2.*(dr)**2.*aj ## syyy*aj ## sx**2+6.*(ds)**2.*(dr)**2.*aj ## sxxy*aj ## sx**2+4.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sxxx*aj ## sx+3.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sxx**2-8.*aj ## sx**3.*(ds)*(dr)**2.*aj ## sxy)/(ds)**5/(dr)**2 )
66 
67 #defineMacro yLapSqCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sy*(aj ## sy**4+2.*aj ## sy**2.*aj ## sx**2+aj ## sx**4)/(ds)**5 )
68 
69 #defineMacro lapCubedCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(180.*aj ## sy**4.*aj ## ry**2.*(dr)**4+30.*aj ## rx**6.*(ds)**4+90.*aj ## ry**4.*aj ## rx**2.*(ds)**4-75.*aj ## ryy*aj ## ry**4.*(dr)*(ds)**4-75.*aj ## rxx*aj ## rx**4.*(dr)*(ds)**4+180.*aj ## rx**2.*aj ## sx**4.*(dr)**4+30.*aj ## ry**6.*(ds)**4+90.*aj ## ry**2.*aj ## rx**4.*(ds)**4-18.*aj ## sy**4.*(dr)**5.*aj ## rxx+288.*aj ## sx**3.*(dr)**4.*aj ## ry*aj ## sy*aj ## rx+288.*aj ## sy**3.*(dr)**4.*aj ## ry*aj ## rx*aj ## sx-360.*aj ## sy**3.*(dr)**5.*aj ## syy*aj ## ry+240.*aj ## sy**2.*aj ## ry**4.*(dr)**2.*(ds)**2+30.*(dr)**4.*(ds)**4.*aj ## ryy*aj ## ryyyy+30.*(dr)**3.*(ds)**4.*aj ## rxxxx*aj ## rx**2+240.*aj ## rx**4.*aj ## sx**2.*(dr)**2.*(ds)**2+216.*aj ## sy**2.*(dr)**4.*aj ## ry**2.*aj ## sx**2+24.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rx*aj ## rxxx+144.*(dr)**3.*(ds)**4.*aj ## rxxy*aj ## rxy*aj ## rx+48.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxxy*aj ## rx+72.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxy*aj ## rxx+48.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxx*aj ## rxy+72.*(dr)**3.*(ds)**4.*aj ## rxyy*aj ## rxx*aj ## rx-432.*aj ## sy**2.*(dr)**5.*aj ## ry*aj ## sxy*aj ## sx-216.*aj ## sy**2.*(dr)**5.*aj ## syy*aj ## sx*aj ## rx-216.*aj ## sy*(dr)**5.*aj ## syy*aj ## ry*aj ## sx**2-288.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rx*aj ## rxyy-144.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## ry**2.*aj ## rxx-576.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## ry*aj ## rxy*aj ## rx-72.*(dr)**2.*(ds)**4.*aj ## ryy**2.*aj ## rx**2-96.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## ry*aj ## rx**2-216.*aj ## sx*(dr)**5.*aj ## sy**2.*aj ## rx*aj ## sxx-216.*aj ## sx**2.*(dr)**5.*aj ## ry*aj ## sy*aj ## sxx+48.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rx*aj ## rxyyy+48.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## rxy*aj ## rx-432.*aj ## sx**2.*(dr)**5.*aj ## sy*aj ## sxy*aj ## rx-120.*aj ## rx**3.*(dr)*(ds)**4.*aj ## ry*aj ## rxy+72.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rx*aj ## rxyy+72.*(dr)**3.*(ds)**4.*aj ## ry*aj ## ryy*aj ## rxxy+24.*(dr)**3.*(ds)**4.*aj ## ry*aj ## ryyy*aj ## rxx+144.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxy*aj ## rxyy-360.*aj ## rx**2.*aj ## sx**2.*(dr)**3.*(ds)**2.*aj ## rxx-120.*aj ## ry**3.*(dr)*(ds)**4.*aj ## rxy*aj ## rx-90.*aj ## rx**2.*(dr)*(ds)**4.*aj ## ry**2.*aj ## rxx-144.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sx*aj ## rx*aj ## rxx-36.*(dr)**4.*(ds)**2.*aj ## ry**2.*aj ## sxx**2-72.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sxx*aj ## rx**2-720.*(dr)**4.*(ds)**2.*aj ## sx*aj ## sxx*aj ## rx*aj ## rxx-144.*(dr)**4.*(ds)**2.*aj ## sx**2.*aj ## rxy**2-240.*(dr)**4.*(ds)**2.*aj ## rxxx*aj ## rx*aj ## sx**2-240.*(dr)**4.*(ds)**2.*aj ## sxxx*aj ## rx**2.*aj ## sx-240.*aj ## ry**3.*aj ## sy*(dr)**3.*(ds)**2.*aj ## syy-360.*aj ## ry**2.*aj ## sy**2.*(dr)**3.*(ds)**2.*aj ## ryy-240.*aj ## rx**3.*aj ## sx*(dr)**3.*(ds)**2.*aj ## sxx-96.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxxx*aj ## rx-576.*(dr)**4.*(ds)**2.*aj ## sx*aj ## sxy*aj ## rxy*aj ## rx-144.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxx*aj ## rxx-288.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxx*aj ## sxy*aj ## rx-288.*(dr)**4.*(ds)**2.*aj ## ry*aj ## rxy*aj ## sx*aj ## sxx-288.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxxy*aj ## sx*aj ## rx-288.*(dr)**4.*(ds)**2.*aj ## ry*aj ## rxx*aj ## sxy*aj ## sx-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## rxy*aj ## sx*aj ## rxx-288.*(dr)**2.*(ds)**4.*aj ## rxy**2.*aj ## rx**2-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rxy*aj ## rx-96.*(dr)**2.*(ds)**4.*aj ## rxyy*aj ## rx**3-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## rxxy*aj ## sx*aj ## rx-96.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sy*aj ## sx*aj ## rxxx+48.*aj ## ry**4.*(dr)**2.*(ds)**2.*aj ## sx**2-48.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rx*aj ## rxxx-144.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxy*aj ## rx**2-48.*(dr)**4.*(ds)**2.*aj ## ry**2.*aj ## sxxx*aj ## sx-144.*(dr)**4.*(ds)**2.*aj ## ry*aj ## rxxy*aj ## sx**2-144.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## rx*aj ## sxx*aj ## sx-36.*(dr)**4.*(ds)**2.*aj ## ryy**2.*aj ## sx**2-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rxx*aj ## rx-36.*(dr)**4.*(ds)**2.*aj ## syy**2.*aj ## rx**2-576.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxy*aj ## rxy+90.*(dr)**5.*(ds)**2.*aj ## rxx*aj ## sxx**2-96.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## ry*aj ## sx*aj ## rx-360.*aj ## sx**3.*(dr)**5.*aj ## sxx*aj ## rx-144.*(dr)**4.*(ds)**2.*aj ## rxyy*aj ## rx*aj ## sx**2-144.*(dr)**4.*(ds)**2.*aj ## sxyy*aj ## rx**2.*aj ## sx-72.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## rxx*aj ## sx**2-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sxy*aj ## rx+72.*(dr)**5.*(ds)**2.*aj ## ryy*aj ## sxy**2-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## syy*aj ## rxy*aj ## rx+18.*(dr)**5.*(ds)**2.*aj ## syy**2.*aj ## rxx-96.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ryyy*aj ## sx*aj ## rx+36.*(dr)**5.*(ds)**2.*aj ## sy**2.*aj ## rxxyy-144.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ry*aj ## ryy*aj ## sxx+6.*(dr)**5.*(ds)**2.*aj ## ryyyy*aj ## sx**2-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxyy*aj ## rx+36.*(dr)**5.*(ds)**2.*aj ## rxxyy*aj ## sx**2+120.*(dr)**3.*(ds)**4.*aj ## rxxx*aj ## rxx*aj ## rx-144.*(dr)**4.*(ds)**2.*aj ## sy*aj ## syy*aj ## ry*aj ## rxx+18.*(dr)**5.*(ds)**2.*aj ## ryy*aj ## sxx**2-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rxyy+6.*(dr)**5.*(ds)**2.*aj ## sy**2.*aj ## rxxxx-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sx*aj ## rxy+72.*(dr)**5.*(ds)**2.*aj ## sxy**2.*aj ## rxx-48.*(dr)**4.*(ds)**2.*aj ## ry*aj ## ryyy*aj ## sx**2-72.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## ryy*aj ## rxx-144.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rx*aj ## rxyy-144.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## ry*aj ## rxxy-144.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ry**2.*aj ## sxxy-48.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sy*aj ## rx**2-144.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ryy*aj ## sx*aj ## rx-144.*(dr)**4.*(ds)**2.*aj ## ry**2.*aj ## sxy**2-288.*(dr)**4.*(ds)**2.*aj ## ry*aj ## syy*aj ## sx*aj ## rxy-144.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxy**2-288.*(dr)**4.*(ds)**2.*aj ## ry*aj ## syy*aj ## sxy*aj ## rx-180.*(dr)**4.*(ds)**2.*aj ## ryy**2.*aj ## sy**2-288.*(dr)**4.*(ds)**2.*aj ## ry*aj ## ryy*aj ## sxy*aj ## sx+30.*(dr)**5.*(ds)**2.*aj ## ryyyy*aj ## sy**2+60.*(dr)**5.*(ds)**2.*aj ## syyyy*aj ## sy*aj ## ry+90.*(dr)**5.*(ds)**2.*aj ## ryy*aj ## syy**2-240.*(dr)**4.*(ds)**2.*aj ## ryyy*aj ## ry*aj ## sy**2-240.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## ry**2.*aj ## sy-720.*(dr)**4.*(ds)**2.*aj ## ry*aj ## ryy*aj ## syy*aj ## sy-180.*(dr)**4.*(ds)**2.*aj ## ry**2.*aj ## syy**2-72.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ry**2.*aj ## sxx-144.*(dr)**4.*(ds)**2.*aj ## ry**2.*aj ## sx*aj ## sxyy+12.*(dr)**5.*(ds)**2.*aj ## syyyy*aj ## sx*aj ## rx+120.*(dr)**5.*(ds)**2.*aj ## ryyy*aj ## sy*aj ## syy+120.*(dr)**5.*(ds)**2.*aj ## ry*aj ## syy*aj ## syyy+120.*(dr)**5.*(ds)**2.*aj ## ryy*aj ## syyy*aj ## sy+24.*(dr)**5.*(ds)**2.*aj ## sy*aj ## ryyy*aj ## sxx+72.*(dr)**5.*(ds)**2.*aj ## syy*aj ## sx*aj ## rxyy+72.*(dr)**5.*(ds)**2.*aj ## syy*aj ## sxyy*aj ## rx+144.*(dr)**5.*(ds)**2.*aj ## syy*aj ## sxy*aj ## rxy+24.*(dr)**5.*(ds)**2.*aj ## syyy*aj ## ry*aj ## sxx+48.*(dr)**5.*(ds)**2.*aj ## ryyy*aj ## sxy*aj ## sx+144.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rxyy+72.*(dr)**5.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxxyy+36.*(dr)**5.*(ds)**2.*aj ## syy*aj ## ryy*aj ## sxx+72.*(dr)**5.*(ds)**2.*aj ## ry*aj ## syy*aj ## sxxy+48.*(dr)**5.*(ds)**2.*aj ## ry*aj ## sx*aj ## sxyyy+72.*(dr)**5.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sxxy+144.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxyy*aj ## rxy+48.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxyyy+24.*(dr)**5.*(ds)**2.*aj ## sy*aj ## syyy*aj ## rxx+48.*(dr)**5.*(ds)**2.*aj ## syyy*aj ## sx*aj ## rxy+48.*(dr)**5.*(ds)**2.*aj ## syyy*aj ## sxy*aj ## rx+48.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxyyy*aj ## rx+72.*(dr)**5.*(ds)**2.*aj ## sy*aj ## syy*aj ## rxxy+72.*(dr)**5.*(ds)**2.*aj ## ryy*aj ## sx*aj ## sxyy+144.*(dr)**5.*(ds)**2.*aj ## ry*aj ## sxyy*aj ## sxy+72.*(dr)**5.*(ds)**2.*aj ## ry*aj ## sxx*aj ## sxxy+144.*(dr)**5.*(ds)**2.*aj ## sx*aj ## sxxy*aj ## rxy+144.*(dr)**5.*(ds)**2.*aj ## rxxy*aj ## sxy*aj ## sx+72.*(dr)**5.*(ds)**2.*aj ## sxxyy*aj ## sx*aj ## rx+24.*(dr)**5.*(ds)**2.*aj ## syy*aj ## sxxx*aj ## rx+24.*(dr)**5.*(ds)**2.*aj ## syy*aj ## sx*aj ## rxxx+36.*(dr)**5.*(ds)**2.*aj ## syy*aj ## sxx*aj ## rxx+48.*(dr)**5.*(ds)**2.*aj ## ry*aj ## sxxx*aj ## sxy+48.*(dr)**5.*(ds)**2.*aj ## ry*aj ## sxxxy*aj ## sx+12.*(dr)**5.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxxxx+144.*(dr)**5.*(ds)**2.*aj ## sxy*aj ## sxx*aj ## rxy+144.*(dr)**5.*(ds)**2.*aj ## sxy*aj ## sxxy*aj ## rx+72.*(dr)**5.*(ds)**2.*aj ## sxyy*aj ## sxx*aj ## rx+48.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxxx*aj ## rxy+48.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxxxy*aj ## rx+48.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rxxx-144.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## rx**2.*aj ## rxx-96.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rx*aj ## rxxx-576.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxy*aj ## rxx*aj ## rx-72.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxx**2-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxxy*aj ## rx**2+120.*(dr)**5.*(ds)**2.*aj ## rxx*aj ## sxxx*aj ## sx+120.*(dr)**5.*(ds)**2.*aj ## sxx*aj ## sxxx*aj ## rx+60.*(dr)**5.*(ds)**2.*aj ## sxxxx*aj ## sx*aj ## rx+120.*(dr)**5.*(ds)**2.*aj ## rxxx*aj ## sxx*aj ## sx+30.*(dr)**5.*(ds)**2.*aj ## rxxxx*aj ## sx**2+72.*(dr)**5.*(ds)**2.*aj ## sx*aj ## sxyy*aj ## rxx+72.*(dr)**5.*(ds)**2.*aj ## rxyy*aj ## sxx*aj ## sx+24.*(dr)**5.*(ds)**2.*aj ## ryy*aj ## sxxx*aj ## sx+48.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxxxy+72.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rxxy+72.*(dr)**5.*(ds)**2.*aj ## sy*aj ## sxxy*aj ## rxx-288.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry**2.*aj ## sx*aj ## rxy+30.*(dr)**3.*(ds)**4.*aj ## ryyyy*aj ## ry**2-288.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry**2.*aj ## sxy*aj ## rx-288.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy*aj ## ryy*aj ## sx*aj ## rx+384.*aj ## rx**3.*(dr)**2.*(ds)**2.*aj ## ry*aj ## sy*aj ## sx+48.*aj ## rx**4.*(dr)**2.*(ds)**2.*aj ## sy**2+288.*aj ## rx**2.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## sx**2+288.*aj ## ry**2.*(dr)**2.*(ds)**2.*aj ## sy**2.*aj ## rx**2+384.*aj ## ry**3.*(dr)**2.*(ds)**2.*aj ## sy*aj ## rx*aj ## sx-72.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## sy**2.*aj ## rxx-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry**3.*aj ## sxx-96.*(dr)**3.*(ds)**2.*aj ## ry**3.*aj ## sxy*aj ## sx-72.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## sy**2.*aj ## rx**2-72.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## ry**2.*aj ## sx**2-144.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy*aj ## syy*aj ## rx**2+24.*(dr)**4.*(ds)**4.*aj ## rxyy*aj ## rxxx-288.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy**2.*aj ## rxy*aj ## rx+36.*(dr)**4.*(ds)**4.*aj ## rxx*aj ## rxxyy-144.*(dr)**3.*(ds)**2.*aj ## syy*aj ## ry**2.*aj ## sx*aj ## rx+6.*(dr)**4.*(ds)**4.*aj ## ryy*aj ## rxxxx-144.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxx*aj ## rx**2+36.*(dr)**4.*(ds)**4.*aj ## ryy*aj ## rxxyy-288.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rx*aj ## rxx-96.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rx**3-288.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxy*aj ## rx**2+48.*(dr)**4.*(ds)**4.*aj ## rxy*aj ## rxxxy-72.*(dr)**3.*(ds)**2.*aj ## sy**2.*aj ## rx**2.*aj ## rxx-72.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## rxx*aj ## sx**2-144.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## rx*aj ## sxx*aj ## sx+6.*(dr)**4.*(ds)**4.*aj ## ryyyy*aj ## rxx-288.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxy*aj ## rx*aj ## sx**2+24.*(dr)**4.*(ds)**4.*aj ## ryyy*aj ## rxxy-288.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sxy*aj ## rx**2.*aj ## sx+48.*(dr)**4.*(ds)**4.*aj ## rxyyy*aj ## rxy-72.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## rx**2.*aj ## sx**2-48.*(dr)**3.*(ds)**2.*aj ## syy*aj ## rx**3.*aj ## sx+120.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## ryy*aj ## ry-90.*aj ## ry**2.*(dr)*(ds)**4.*aj ## ryy*aj ## rx**2+30.*(dr)**4.*(ds)**4.*aj ## rxx*aj ## rxxxx-108.*aj ## sy**2.*(dr)**5.*aj ## ryy*aj ## sx**2-96.*(dr)**2.*(ds)**4.*aj ## ry**3.*aj ## rxxy-288.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxy**2+216.*aj ## sx**2.*(dr)**4.*aj ## sy**2.*aj ## rx**2-144.*aj ## sy**3.*(dr)**5.*aj ## sx*aj ## rxy-72.*aj ## sy**3.*(dr)**5.*aj ## ry*aj ## sxx-144.*aj ## sy**3.*(dr)**5.*aj ## sxy*aj ## rx+36.*(dr)**3.*(ds)**4.*aj ## rxxyy*aj ## rx**2+6.*(dr)**3.*(ds)**4.*aj ## ry**2.*aj ## rxxxx-144.*aj ## sx**3.*(dr)**5.*aj ## ry*aj ## sxy+18.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxx**2+72.*(dr)**3.*(ds)**4.*aj ## rxy**2.*aj ## rxx-72.*aj ## sx**3.*(dr)**5.*aj ## syy*aj ## rx-144.*aj ## sx**3.*(dr)**5.*aj ## sy*aj ## rxy-108.*aj ## sx**2.*(dr)**5.*aj ## sy**2.*aj ## rxx+36.*(dr)**3.*(ds)**4.*aj ## ry**2.*aj ## rxxyy+72.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxy**2+18.*(dr)**3.*(ds)**4.*aj ## ryy**2.*aj ## rxx+6.*(dr)**3.*(ds)**4.*aj ## ryyyy*aj ## rx**2-360.*aj ## ry**2.*(dr)**2.*(ds)**4.*aj ## ryy**2-160.*aj ## ry**3.*(dr)**2.*(ds)**4.*aj ## ryyy-160.*aj ## rx**3.*(dr)**2.*(ds)**4.*aj ## rxxx-360.*aj ## rx**2.*(dr)**2.*(ds)**4.*aj ## rxx**2-15.*aj ## rx**4.*(dr)*(ds)**4.*aj ## ryy-15.*aj ## ry**4.*(dr)*(ds)**4.*aj ## rxx-36.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxx**2-180.*(dr)**4.*(ds)**2.*aj ## rx**2.*aj ## sxx**2-180.*(dr)**4.*(ds)**2.*aj ## sx**2.*aj ## rxx**2-144.*(dr)**4.*(ds)**2.*aj ## rx**2.*aj ## sxy**2-90.*aj ## sy**4.*(dr)**5.*aj ## ryy+20.*(dr)**4.*(ds)**4.*aj ## ryyy**2+20.*(dr)**4.*(ds)**4.*aj ## rxxx**2+36.*aj ## sy**4.*(dr)**4.*aj ## rx**2+30.*(dr)**3.*(ds)**4.*aj ## rxx**3+36.*aj ## sx**4.*(dr)**4.*aj ## ry**2-18.*aj ## sx**4.*(dr)**5.*aj ## ryy+36.*(dr)**4.*(ds)**4.*aj ## rxxy**2+36.*(dr)**4.*(ds)**4.*aj ## rxyy**2-90.*aj ## sx**4.*(dr)**5.*aj ## rxx+30.*(dr)**3.*(ds)**4.*aj ## ryy**3)/(dr)**6/(ds)**4 )
70 
71 #defineMacro lapCubedCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*(36.*(dr)**3.*(ds)**2.*aj ## rxy**2.*aj ## rxx-72.*(dr)**2.*(ds)**2.*aj ## rxy**2.*aj ## rx**2+18.*(dr)**3.*(ds)**2.*aj ## rxxyy*aj ## rx**2+3.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## rxxxx-24.*(dr)**2.*(ds)**2.*aj ## ry**3.*aj ## rxxy-72.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxy**2+3.*(dr)**3.*(ds)**2.*aj ## ryyyy*aj ## rx**2+18.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## rxxyy+12.*aj ## ry**6.*(ds)**2+12.*aj ## rx**6.*(ds)**2+36.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## rxy**2-60.*aj ## rxx*aj ## rx**4.*(dr)*(ds)**2-60.*aj ## ryy*aj ## ry**4.*(dr)*(ds)**2+60.*aj ## rx**4.*aj ## sx**2.*(dr)**2+9.*(dr)**3.*(ds)**2.*aj ## ryy**2.*aj ## rxx+36.*aj ## ry**4.*aj ## rx**2.*(ds)**2+36.*aj ## ry**2.*aj ## rx**4.*(ds)**2+12.*(dr)**3.*(ds)**2.*aj ## ry*aj ## ryyy*aj ## rxx+36.*(dr)**3.*(ds)**2.*aj ## ry*aj ## ryy*aj ## rxxy+36.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## rx*aj ## rxyy+24.*(dr)**3.*(ds)**2.*aj ## ryyy*aj ## rxy*aj ## rx+24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rx*aj ## rxyyy+72.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxy*aj ## rxyy+60.*aj ## sy**2.*aj ## ry**4.*(dr)**2+60.*(dr)**3.*(ds)**2.*aj ## ryyy*aj ## ryy*aj ## ry+15.*(dr)**3.*(ds)**2.*aj ## ryy**3+72.*(dr)**3.*(ds)**2.*aj ## rxxy*aj ## rxy*aj ## rx+12.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## rx*aj ## rxxx+15.*(dr)**3.*(ds)**2.*aj ## ryyyy*aj ## ry**2+24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxxx*aj ## rxy+36.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxxy*aj ## rxx+24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxxxy*aj ## rx+36.*(dr)**3.*(ds)**2.*aj ## rxyy*aj ## rxx*aj ## rx-96.*aj ## ry**3.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-96.*aj ## rx**3.*(dr)*(ds)**2.*aj ## ry*aj ## rxy-72.*aj ## rx**2.*(dr)*(ds)**2.*aj ## ry**2.*aj ## rxx-72.*aj ## ry**2.*(dr)*(ds)**2.*aj ## ryy*aj ## rx**2-12.*aj ## rx**4.*(dr)*(ds)**2.*aj ## ryy-12.*aj ## ry**4.*(dr)*(ds)**2.*aj ## rxx-40.*aj ## ry**3.*(dr)**2.*(ds)**2.*aj ## ryyy-90.*aj ## ry**2.*(dr)**2.*(ds)**2.*aj ## ryy**2+12.*aj ## ry**4.*(dr)**2.*aj ## sx**2-36.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rx**2.*aj ## rxx-18.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxx**2-24.*(dr)**2.*(ds)**2.*aj ## rxyy*aj ## rx**3-72.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxxy*aj ## rx**2-144.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxy*aj ## rxx*aj ## rx-24.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rx*aj ## rxxx-144.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## ry*aj ## rxy*aj ## rx-36.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## ry**2.*aj ## rxx-72.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rx*aj ## rxyy-24.*(dr)**2.*(ds)**2.*aj ## ryyy*aj ## ry*aj ## rx**2-90.*aj ## rx**2.*(dr)**2.*(ds)**2.*aj ## rxx**2+60.*(dr)**3.*(ds)**2.*aj ## rxxx*aj ## rxx*aj ## rx+15.*(dr)**3.*(ds)**2.*aj ## rxx**3-144.*(dr)**3.*aj ## sy*aj ## sx*aj ## rxy*aj ## rx**2-144.*(dr)**3.*aj ## sy*aj ## ry*aj ## sx*aj ## rx*aj ## rxx-72.*(dr)**3.*aj ## sy*aj ## ry*aj ## sxx*aj ## rx**2-36.*(dr)**3.*aj ## sy**2.*aj ## rx**2.*aj ## rxx-48.*(dr)**3.*aj ## sy*aj ## sxy*aj ## rx**3-36.*(dr)**3.*aj ## ry**2.*aj ## rxx*aj ## sx**2+12.*aj ## rx**4.*(dr)**2.*aj ## sy**2-72.*(dr)**3.*aj ## syy*aj ## ry**2.*aj ## sx*aj ## rx-144.*(dr)**3.*aj ## ry*aj ## sy**2.*aj ## rxy*aj ## rx-72.*(dr)**3.*aj ## ry*aj ## sy*aj ## syy*aj ## rx**2-144.*(dr)**3.*aj ## ry*aj ## sy*aj ## ryy*aj ## sx*aj ## rx-144.*(dr)**3.*aj ## sy*aj ## ry**2.*aj ## sxy*aj ## rx-144.*(dr)**3.*aj ## sy*aj ## ry**2.*aj ## sx*aj ## rxy+96.*aj ## rx**3.*(dr)**2.*aj ## ry*aj ## sy*aj ## sx-180.*aj ## rx**2.*aj ## sx**2.*(dr)**3.*aj ## rxx-120.*aj ## rx**3.*aj ## sx*(dr)**3.*aj ## sxx-24.*(dr)**3.*aj ## syy*aj ## rx**3.*aj ## sx-36.*(dr)**3.*aj ## ryy*aj ## rx**2.*aj ## sx**2-144.*(dr)**3.*aj ## ry*aj ## sxy*aj ## rx**2.*aj ## sx-144.*(dr)**3.*aj ## ry*aj ## rxy*aj ## rx*aj ## sx**2-72.*(dr)**3.*aj ## ry**2.*aj ## rx*aj ## sxx*aj ## sx-24.*(dr)**3.*aj ## sy*aj ## ry**3.*aj ## sxx-36.*(dr)**3.*aj ## ry**2.*aj ## sy**2.*aj ## rxx+72.*aj ## ry**2.*(dr)**2.*aj ## sy**2.*aj ## rx**2-36.*(dr)**3.*aj ## ryy*aj ## ry**2.*aj ## sx**2-36.*(dr)**3.*aj ## ryy*aj ## sy**2.*aj ## rx**2-48.*(dr)**3.*aj ## ry**3.*aj ## sxy*aj ## sx+96.*aj ## ry**3.*(dr)**2.*aj ## sy*aj ## rx*aj ## sx-18.*(dr)**2.*(ds)**2.*aj ## ryy**2.*aj ## rx**2+15.*(dr)**3.*(ds)**2.*aj ## rxxxx*aj ## rx**2+9.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## rxx**2-180.*aj ## ry**2.*aj ## sy**2.*(dr)**3.*aj ## ryy+72.*aj ## rx**2.*(dr)**2.*aj ## ry**2.*aj ## sx**2-120.*aj ## ry**3.*aj ## sy*(dr)**3.*aj ## syy-40.*aj ## rx**3.*(dr)**2.*(ds)**2.*aj ## rxxx)/(dr)**6/(ds)**2 )
72 
73 #defineMacro lapCubedCoeff2dOrder2Ghost30(aj,dr,ds) ( 1/2.*(2.*aj ## ry**6+2.*aj ## rx**6+6.*aj ## ry**2.*aj ## rx**4-18.*aj ## ry**2.*(dr)*aj ## ryy*aj ## rx**2-24.*aj ## ry**3.*(dr)*aj ## rxy*aj ## rx-3.*aj ## ry**4.*(dr)*aj ## rxx-3.*aj ## rx**4.*(dr)*aj ## ryy-24.*aj ## rx**3.*(dr)*aj ## ry*aj ## rxy-18.*aj ## rx**2.*(dr)*aj ## ry**2.*aj ## rxx-15.*aj ## ryy*aj ## ry**4.*(dr)+6.*aj ## ry**4.*aj ## rx**2-15.*aj ## rxx*aj ## rx**4.*(dr))/(dr)**6 )
74 
75 #defineMacro lapCubedCoeff2dOrder2Ghost01(aj,dr,ds) ( 1/2.*(-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxx*aj ## sxy*aj ## sx-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxxy*aj ## sx*aj ## rx-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxy*aj ## sx*aj ## sxx-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxx*aj ## sxy*aj ## rx-144.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sy*aj ## sxx*aj ## rxx+36.*aj ## ry**4.*(ds)**4.*aj ## sx**2-72.*(dr)**2.*(ds)**4.*aj ## syy*aj ## sxx*aj ## rx**2-36.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## sxx**2-36.*(dr)**2.*(ds)**4.*aj ## sy**2.*aj ## rxx**2-144.*(dr)**2.*(ds)**4.*aj ## rx**2.*aj ## sxy**2-180.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## syy**2+6.*(ds)**3.*(dr)**4.*aj ## syyyy*aj ## sx**2-576.*(dr)**2.*(ds)**4.*aj ## sx*aj ## sxy*aj ## rxy*aj ## rx-96.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sy*aj ## sxxx*aj ## rx-96.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sy*aj ## sx*aj ## rxxx-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rxxy*aj ## sx*aj ## rx-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## sxx*aj ## rxy*aj ## rx-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rxy*aj ## sx*aj ## rxx-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## sxy*aj ## rxx*aj ## rx-144.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## rx*aj ## sxx*aj ## sx-144.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxxy*aj ## sx**2-48.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## sxxx*aj ## sx-144.*(dr)**2.*(ds)**4.*aj ## sy*aj ## sxxy*aj ## rx**2-48.*(dr)**2.*(ds)**4.*aj ## sy**2.*aj ## rx*aj ## rxxx-72.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## rxx*aj ## sx**2-144.*(dr)**2.*(ds)**4.*aj ## sxyy*aj ## rx**2.*aj ## sx-144.*(dr)**2.*(ds)**4.*aj ## rxyy*aj ## rx*aj ## sx**2-240.*(dr)**2.*(ds)**4.*aj ## syyy*aj ## ry**2.*aj ## sy-240.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## ry*aj ## sy**2-180.*(dr)**2.*(ds)**4.*aj ## ryy**2.*aj ## sy**2+30.*(ds)**3.*(dr)**4.*aj ## sxxxx*aj ## sx**2+216.*aj ## ry**2.*(ds)**4.*aj ## sy**2.*aj ## rx**2+216.*aj ## rx**2.*(ds)**4.*aj ## ry**2.*aj ## sx**2+36.*(ds)**3.*(dr)**4.*aj ## sxxyy*aj ## sx**2+72.*(ds)**3.*(dr)**4.*aj ## sxx*aj ## sxy**2+6.*(ds)**3.*(dr)**4.*aj ## sy**2.*aj ## sxxxx+18.*(ds)**3.*(dr)**4.*aj ## syy*aj ## sxx**2-90.*aj ## rx**4.*(ds)**5.*aj ## sxx+48.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxxxy*aj ## sx+72.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxx*aj ## sxxy+24.*(ds)**3.*(dr)**4.*aj ## syy*aj ## sxxx*aj ## sx-288.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxy*aj ## sx**2-96.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## sxxx*aj ## sx-144.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sx**2.*aj ## sxx-288.*(dr)**4.*(ds)**2.*aj ## sxy**2.*aj ## sx**2+288.*aj ## rx**3.*(ds)**4.*aj ## ry*aj ## sy*aj ## sx+72.*(ds)**3.*(dr)**4.*aj ## sxyy*aj ## sx*aj ## sxx+144.*(ds)**3.*(dr)**4.*aj ## sxxy*aj ## sxy*aj ## sx+48.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxxx*aj ## sxy+288.*aj ## ry**3.*(ds)**4.*aj ## sy*aj ## rx*aj ## sx+36.*aj ## rx**4.*(ds)**4.*aj ## sy**2-96.*(dr)**4.*(ds)**2.*aj ## sxyy*aj ## sx**3-360.*aj ## sy**2.*(ds)**2.*(dr)**4.*aj ## syy**2-160.*aj ## sy**3.*(ds)**2.*(dr)**4.*aj ## syyy-160.*aj ## sx**3.*(ds)**2.*(dr)**4.*aj ## sxxx-360.*aj ## sx**2.*(ds)**2.*(dr)**4.*aj ## sxx**2-240.*aj ## rx*aj ## sx**3.*(ds)**3.*(dr)**2.*aj ## rxx-360.*aj ## rx**2.*aj ## sx**2.*(ds)**3.*(dr)**2.*aj ## sxx-240.*aj ## ry*aj ## sy**3.*(ds)**3.*(dr)**2.*aj ## ryy-360.*aj ## ry**2.*aj ## sy**2.*(ds)**3.*(dr)**2.*aj ## syy+6.*(ds)**4.*(dr)**4.*aj ## syy*aj ## sxxxx+36.*(ds)**4.*(dr)**4.*aj ## sxxyy*aj ## sxx+24.*(ds)**4.*(dr)**4.*aj ## sxyy*aj ## sxxx-288.*(dr)**2.*(ds)**3.*aj ## sy*aj ## rxy*aj ## rx*aj ## sx**2-288.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxy*aj ## rx**2.*aj ## sx+48.*aj ## sy**4.*(ds)**2.*(dr)**2.*aj ## rx**2-72.*(dr)**2.*(ds)**3.*aj ## ry**2.*aj ## sx**2.*aj ## sxx-48.*(dr)**2.*(ds)**3.*aj ## ryy*aj ## rx*aj ## sx**3-144.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## sx*aj ## rx*aj ## rxx-72.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## sxx*aj ## rx**2-288.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sxy*aj ## rx*aj ## sx**2-96.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rxy*aj ## sx**3-72.*(dr)**2.*(ds)**3.*aj ## syy*aj ## rx**2.*aj ## sx**2-288.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy*aj ## rx*aj ## sxx*aj ## sx-144.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy*aj ## rxx*aj ## sx**2-240.*(dr)**2.*(ds)**4.*aj ## sxxx*aj ## rx**2.*aj ## sx-240.*(dr)**2.*(ds)**4.*aj ## rxxx*aj ## rx*aj ## sx**2-720.*(dr)**2.*(ds)**4.*aj ## sx*aj ## sxx*aj ## rx*aj ## rxx-288.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## sx*aj ## rxy-288.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## sxy*aj ## rx-288.*(dr)**2.*(ds)**3.*aj ## ry*aj ## syy*aj ## sy*aj ## sx*aj ## rx-144.*(dr)**2.*(ds)**3.*aj ## ry*aj ## ryy*aj ## sy*aj ## sx**2-72.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## ry**2.*aj ## sxx-72.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## syy*aj ## rx**2+48.*aj ## sx**4.*(ds)**2.*(dr)**2.*aj ## ry**2+144.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxyy*aj ## rxy+48.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sx*aj ## rxyyy+24.*(ds)**5.*(dr)**2.*aj ## ry*aj ## syyy*aj ## rxx-96.*(dr)**2.*(ds)**3.*aj ## sy**3.*aj ## rxy*aj ## rx-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy**3.*aj ## rxx-144.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## ryy*aj ## sx*aj ## rx+384.*aj ## sx**3.*(ds)**2.*(dr)**2.*aj ## ry*aj ## sy*aj ## rx+288.*aj ## sx**2.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## rx**2+288.*aj ## sy**2.*(ds)**2.*(dr)**2.*aj ## ry**2.*aj ## sx**2+384.*aj ## sy**3.*(ds)**2.*(dr)**2.*aj ## ry*aj ## rx*aj ## sx-120.*aj ## sx**3.*(ds)*(dr)**4.*aj ## sy*aj ## sxy-90.*aj ## sx**2.*(ds)*(dr)**4.*aj ## sy**2.*aj ## sxx-15.*aj ## sx**4.*(ds)*(dr)**4.*aj ## syy-15.*aj ## sy**4.*(ds)*(dr)**4.*aj ## sxx-120.*aj ## sy**3.*(ds)*(dr)**4.*aj ## sxy*aj ## sx-90.*aj ## sy**2.*(ds)*(dr)**4.*aj ## syy*aj ## sx**2+120.*(ds)**3.*(dr)**4.*aj ## sxxx*aj ## sx*aj ## sxx-720.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## syy*aj ## sy+120.*(ds)**5.*(dr)**2.*aj ## syyy*aj ## ryy*aj ## ry-576.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ry*aj ## sxy*aj ## rxy-144.*(dr)**2.*(ds)**4.*aj ## sy**2.*aj ## rxy**2-72.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## sxx**2-288.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry**2.*aj ## sxy*aj ## sx+30.*(ds)**3.*(dr)**4.*aj ## sxx**3-72.*(dr)**2.*(ds)**3.*aj ## syy*aj ## ry**2.*aj ## sx**2-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ryy*aj ## sx*aj ## rxy-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ry*aj ## sx*aj ## rxyy-144.*(dr)**2.*(ds)**4.*aj ## sy*aj ## syy*aj ## ry*aj ## rxx-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ry*aj ## sxyy*aj ## rx-144.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ry*aj ## ryy*aj ## sxx+12.*(ds)**5.*(dr)**2.*aj ## ryyyy*aj ## sx*aj ## rx+72.*(ds)**5.*(dr)**2.*aj ## syy*aj ## rxy**2-96.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ryyy*aj ## sx*aj ## rx-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## syy*aj ## rxy*aj ## rx-288.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ryy*aj ## sxy*aj ## rx-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## sxy*aj ## sx-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## syy*aj ## sxy*aj ## rx-576.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy*aj ## sx*aj ## sxx-288.*(dr)**2.*(ds)**4.*aj ## ry*aj ## syy*aj ## sx*aj ## rxy-144.*(dr)**2.*(ds)**4.*aj ## syy*aj ## ryy*aj ## sx*aj ## rx-48.*(dr)**2.*(ds)**4.*aj ## syyy*aj ## sy*aj ## rx**2-144.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ry**2.*aj ## sxxy-144.*(dr)**2.*(ds)**4.*aj ## sy**2.*aj ## ry*aj ## rxxy-144.*(dr)**2.*(ds)**4.*aj ## sy**2.*aj ## rx*aj ## rxyy-96.*(dr)**2.*(ds)**4.*aj ## syyy*aj ## ry*aj ## sx*aj ## rx-72.*(dr)**2.*(ds)**4.*aj ## sy**2.*aj ## ryy*aj ## rxx-48.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryyy*aj ## sx**2-144.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## sx*aj ## sxyy-72.*(dr)**2.*(ds)**4.*aj ## syy*aj ## ry**2.*aj ## sxx-144.*(dr)**2.*(ds)**4.*aj ## syy*aj ## sx*aj ## rx*aj ## rxx-144.*(dr)**2.*(ds)**4.*aj ## sx**2.*aj ## rxy**2+180.*aj ## sy**2.*aj ## ry**4.*(ds)**4+180.*aj ## rx**4.*aj ## sx**2.*(ds)**4+30.*aj ## sy**6.*(dr)**4+30.*aj ## sx**6.*(dr)**4+18.*(ds)**5.*(dr)**2.*aj ## ryy**2.*aj ## sxx-288.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## sxy**2+6.*(ds)**5.*(dr)**2.*aj ## syyyy*aj ## rx**2+240.*aj ## sy**4.*aj ## ry**2.*(ds)**2.*(dr)**2+240.*aj ## rx**2.*aj ## sx**4.*(ds)**2.*(dr)**2+48.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxxxy*aj ## rx+48.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxxx*aj ## rxy-75.*aj ## syy*aj ## sy**4.*(ds)*(dr)**4-75.*aj ## sx**4.*aj ## sxx*(ds)*(dr)**4+18.*(ds)**5.*(dr)**2.*aj ## syy*aj ## rxx**2+36.*(ds)**5.*(dr)**2.*aj ## sxxyy*aj ## rx**2+72.*(ds)**5.*(dr)**2.*aj ## rxxyy*aj ## sx*aj ## rx+72.*(ds)**5.*(dr)**2.*aj ## sxyy*aj ## rxx*aj ## rx+30.*(ds)**3.*(dr)**4.*aj ## syy**3+144.*(ds)**5.*(dr)**2.*aj ## sxxy*aj ## rxy*aj ## rx+72.*(ds)**5.*(dr)**2.*aj ## sxx*aj ## rxy**2+72.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxx*aj ## rxxy+72.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxxy*aj ## rxx+48.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxy*aj ## rxxx+12.*(ds)**5.*(dr)**2.*aj ## sy*aj ## ry*aj ## rxxxx+48.*(ds)**5.*(dr)**2.*aj ## sy*aj ## rxxx*aj ## rxy+72.*(ds)**5.*(dr)**2.*aj ## sy*aj ## rxxy*aj ## rxx+36.*(ds)**4.*(dr)**4.*aj ## sxyy**2+144.*(ds)**5.*(dr)**2.*aj ## sy*aj ## rxy*aj ## rxyy+48.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sx*aj ## rxxxy+144.*(ds)**5.*(dr)**2.*aj ## rx*aj ## sxy*aj ## rxxy+72.*(ds)**5.*(dr)**2.*aj ## rx*aj ## sxx*aj ## rxyy+144.*(ds)**5.*(dr)**2.*aj ## rxxy*aj ## sx*aj ## rxy+72.*(ds)**5.*(dr)**2.*aj ## rxyy*aj ## sx*aj ## rxx+144.*(ds)**5.*(dr)**2.*aj ## rxx*aj ## sxy*aj ## rxy+24.*(ds)**5.*(dr)**2.*aj ## ryy*aj ## sxxx*aj ## rx+24.*(ds)**5.*(dr)**2.*aj ## ryy*aj ## sx*aj ## rxxx+36.*(ds)**5.*(dr)**2.*aj ## ryy*aj ## sxx*aj ## rxx+24.*(ds)**5.*(dr)**2.*aj ## syy*aj ## rx*aj ## rxxx+48.*(ds)**5.*(dr)**2.*aj ## syyy*aj ## rxy*aj ## rx+144.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxy*aj ## rxyy+36.*(ds)**4.*(dr)**4.*aj ## sxxy**2+72.*(ds)**5.*(dr)**2.*aj ## syy*aj ## rx*aj ## rxyy+36.*(ds)**5.*(dr)**2.*aj ## ry**2.*aj ## sxxyy+48.*(ds)**5.*(dr)**2.*aj ## sy*aj ## rx*aj ## rxyyy+24.*(ds)**5.*(dr)**2.*aj ## ry*aj ## ryyy*aj ## sxx+72.*(ds)**5.*(dr)**2.*aj ## ry*aj ## ryy*aj ## sxxy+72.*(ds)**5.*(dr)**2.*aj ## ryy*aj ## sx*aj ## rxyy+72.*(ds)**5.*(dr)**2.*aj ## ryy*aj ## sxyy*aj ## rx+144.*(ds)**5.*(dr)**2.*aj ## ryy*aj ## sxy*aj ## rxy+48.*(ds)**5.*(dr)**2.*aj ## ryyy*aj ## sx*aj ## rxy+48.*(ds)**5.*(dr)**2.*aj ## ryyy*aj ## sxy*aj ## rx+36.*(ds)**5.*(dr)**2.*aj ## syy*aj ## ryy*aj ## rxx+72.*(ds)**5.*(dr)**2.*aj ## syy*aj ## ry*aj ## rxxy+48.*(ds)**5.*(dr)**2.*aj ## ry*aj ## sxyyy*aj ## rx+72.*(ds)**5.*(dr)**2.*aj ## sy*aj ## ry*aj ## rxxyy+72.*(ds)**5.*(dr)**2.*aj ## sy*aj ## ryy*aj ## rxxy+24.*(ds)**5.*(dr)**2.*aj ## sy*aj ## ryyy*aj ## rxx+48.*(ds)**4.*(dr)**4.*aj ## sxy*aj ## sxyyy+36.*(ds)**4.*(dr)**4.*aj ## syy*aj ## sxxyy+24.*(ds)**4.*(dr)**4.*aj ## syyy*aj ## sxxy-96.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sy*aj ## sx**2-576.*(dr)**4.*(ds)**2.*aj ## sy*aj ## syy*aj ## sxy*aj ## sx-288.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## sx*aj ## sxyy-144.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sy**2.*aj ## sxx-72.*(dr)**4.*(ds)**2.*aj ## syy**2.*aj ## sx**2+48.*(ds)**4.*(dr)**4.*aj ## sxxxy*aj ## sxy-96.*(dr)**4.*(ds)**2.*aj ## sy**3.*aj ## sxxy+120.*(ds)**3.*(dr)**4.*aj ## syyy*aj ## sy*aj ## syy+60.*(ds)**5.*(dr)**2.*aj ## rxxxx*aj ## sx*aj ## rx+120.*(ds)**5.*(dr)**2.*aj ## sxxx*aj ## rxx*aj ## rx+120.*(ds)**5.*(dr)**2.*aj ## sxx*aj ## rx*aj ## rxxx+120.*(ds)**5.*(dr)**2.*aj ## sx*aj ## rxx*aj ## rxxx+30.*(ds)**4.*(dr)**4.*aj ## sxxxx*aj ## sxx+90.*(ds)**5.*(dr)**2.*aj ## sxx*aj ## rxx**2+30.*(ds)**5.*(dr)**2.*aj ## sxxxx*aj ## rx**2+30.*(ds)**5.*(dr)**2.*aj ## syyyy*aj ## ry**2+90.*(ds)**5.*(dr)**2.*aj ## syy*aj ## ryy**2-360.*aj ## rx**3.*(ds)**5.*aj ## sx*aj ## rxx-144.*aj ## rx**3.*(ds)**5.*aj ## sy*aj ## rxy+30.*(ds)**3.*(dr)**4.*aj ## syyyy*aj ## sy**2+120.*(ds)**5.*(dr)**2.*aj ## ry*aj ## syy*aj ## ryyy+6.*(ds)**4.*(dr)**4.*aj ## syyyy*aj ## sxx+120.*(ds)**5.*(dr)**2.*aj ## ryy*aj ## sy*aj ## ryyy+60.*(ds)**5.*(dr)**2.*aj ## ryyyy*aj ## sy*aj ## ry-18.*aj ## rx**4.*(ds)**5.*aj ## syy-216.*aj ## rx**2.*(ds)**5.*aj ## ry*aj ## sy*aj ## rxx-432.*aj ## rx**2.*(ds)**5.*aj ## ry*aj ## sx*aj ## rxy-144.*aj ## rx**3.*(ds)**5.*aj ## ry*aj ## sxy-108.*aj ## rx**2.*(ds)**5.*aj ## ry**2.*aj ## sxx-72.*aj ## rx**3.*(ds)**5.*aj ## ryy*aj ## sx-216.*aj ## rx*(ds)**5.*aj ## ry**2.*aj ## sx*aj ## rxx-432.*aj ## ry**2.*(ds)**5.*aj ## sy*aj ## rxy*aj ## rx-72.*aj ## ry**3.*(ds)**5.*aj ## sy*aj ## rxx-108.*aj ## ry**2.*(ds)**5.*aj ## syy*aj ## rx**2-144.*aj ## ry**3.*(ds)**5.*aj ## sx*aj ## rxy-18.*aj ## ry**4.*(ds)**5.*aj ## sxx+90.*aj ## sy**4.*aj ## sx**2.*(dr)**4+90.*aj ## sy**2.*aj ## sx**4.*(dr)**4+48.*(ds)**5.*(dr)**2.*aj ## sy*aj ## rxxxy*aj ## rx+6.*(ds)**5.*(dr)**2.*aj ## ry**2.*aj ## sxxxx+20.*(ds)**4.*(dr)**4.*aj ## sxxx**2-216.*aj ## ry*(ds)**5.*aj ## sy*aj ## ryy*aj ## rx**2-216.*aj ## ry**2.*(ds)**5.*aj ## ryy*aj ## sx*aj ## rx-144.*aj ## ry**3.*(ds)**5.*aj ## sxy*aj ## rx+18.*(ds)**3.*(dr)**4.*aj ## syy**2.*aj ## sxx+72.*(ds)**3.*(dr)**4.*aj ## syy*aj ## sxy**2+36.*(ds)**3.*(dr)**4.*aj ## sy**2.*aj ## sxxyy-360.*aj ## ry**3.*(ds)**5.*aj ## sy*aj ## ryy+72.*(ds)**3.*(dr)**4.*aj ## sy*aj ## syy*aj ## sxxy+48.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sx*aj ## sxyyy+24.*(ds)**3.*(dr)**4.*aj ## sy*aj ## syyy*aj ## sxx+144.*(ds)**3.*(dr)**4.*aj ## sy*aj ## sxyy*aj ## sxy+72.*(ds)**3.*(dr)**4.*aj ## syy*aj ## sx*aj ## sxyy+48.*(ds)**3.*(dr)**4.*aj ## syyy*aj ## sxy*aj ## sx+20.*(ds)**4.*(dr)**4.*aj ## syyy**2-90.*aj ## ry**4.*(ds)**5.*aj ## syy-180.*(dr)**2.*(ds)**4.*aj ## sx**2.*aj ## rxx**2-180.*(dr)**2.*(ds)**4.*aj ## rx**2.*aj ## sxx**2+30.*(ds)**4.*(dr)**4.*aj ## syy*aj ## syyyy-36.*(dr)**2.*(ds)**4.*aj ## ryy**2.*aj ## sx**2-36.*(dr)**2.*(ds)**4.*aj ## syy**2.*aj ## rx**2-144.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## sxy**2)/(ds)**6/(dr)**4 )
76 
77 #defineMacro lapCubedCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*(3.*(dr)**2.*(ds)**3.*aj ## syyyy*aj ## sx**2-144.*(ds)**3.*aj ## sy*aj ## rxy*aj ## rx*aj ## sx**2+60.*(dr)**2.*(ds)**3.*aj ## sxxx*aj ## sx*aj ## sxx+15.*(dr)**2.*(ds)**3.*aj ## sxxxx*aj ## sx**2+36.*(dr)**2.*(ds)**3.*aj ## sxyy*aj ## sx*aj ## sxx+18.*(dr)**2.*(ds)**3.*aj ## sxxyy*aj ## sx**2-96.*aj ## sx**3.*(ds)*(dr)**2.*aj ## sy*aj ## sxy+12.*aj ## sx**6.*(dr)**2+12.*aj ## sy**6.*(dr)**2+15.*(dr)**2.*(ds)**3.*aj ## syyyy*aj ## sy**2+60.*aj ## sy**4.*aj ## ry**2.*(ds)**2+60.*aj ## rx**2.*aj ## sx**4.*(ds)**2+36.*aj ## sy**2.*aj ## sx**4.*(dr)**2+15.*(dr)**2.*(ds)**3.*aj ## syy**3+36.*aj ## sy**4.*aj ## sx**2.*(dr)**2-60.*aj ## syy*aj ## sy**4.*(ds)*(dr)**2-60.*aj ## sx**4.*aj ## sxx*(ds)*(dr)**2+60.*(dr)**2.*(ds)**3.*aj ## syyy*aj ## sy*aj ## syy+36.*(dr)**2.*(ds)**3.*aj ## sy*aj ## syy*aj ## sxxy+24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sx*aj ## sxyyy+12.*(dr)**2.*(ds)**3.*aj ## sy*aj ## syyy*aj ## sxx+72.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxyy*aj ## sxy+36.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sx*aj ## sxyy+9.*(dr)**2.*(ds)**3.*aj ## syy**2.*aj ## sxx+24.*(dr)**2.*(ds)**3.*aj ## syyy*aj ## sxy*aj ## sx+36.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sxy**2+18.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## sxxyy+36.*(dr)**2.*(ds)**3.*aj ## sxx*aj ## sxy**2+9.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sxx**2-120.*aj ## rx*aj ## sx**3.*(ds)**3.*aj ## rxx-180.*aj ## rx**2.*aj ## sx**2.*(ds)**3.*aj ## sxx-120.*aj ## ry*aj ## sy**3.*(ds)**3.*aj ## ryy-144.*(ds)**3.*aj ## sy*aj ## sxy*aj ## rx**2.*aj ## sx-144.*(ds)**3.*aj ## ry*aj ## sxy*aj ## rx*aj ## sx**2-48.*(ds)**3.*aj ## ry*aj ## rxy*aj ## sx**3-36.*(ds)**3.*aj ## ry**2.*aj ## sx**2.*aj ## sxx-24.*(ds)**3.*aj ## ryy*aj ## rx*aj ## sx**3-36.*(ds)**3.*aj ## syy*aj ## rx**2.*aj ## sx**2+36.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxx*aj ## sxxy+12.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sxxx*aj ## sx-180.*aj ## ry**2.*aj ## sy**2.*(ds)**3.*aj ## syy+15.*(dr)**2.*(ds)**3.*aj ## sxx**3+72.*(dr)**2.*(ds)**3.*aj ## sxxy*aj ## sxy*aj ## sx+24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxxx*aj ## sxy+24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxxxy*aj ## sx+3.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## sxxxx-144.*(ds)**3.*aj ## ry*aj ## sy*aj ## rx*aj ## sxx*aj ## sx-72.*(ds)**3.*aj ## ry*aj ## sy*aj ## rxx*aj ## sx**2-72.*(ds)**3.*aj ## sy**2.*aj ## sx*aj ## rx*aj ## rxx-36.*(ds)**3.*aj ## sy**2.*aj ## sxx*aj ## rx**2-36.*(ds)**3.*aj ## sy**2.*aj ## ry**2.*aj ## sxx-36.*(ds)**3.*aj ## sy**2.*aj ## syy*aj ## rx**2-48.*(ds)**3.*aj ## sy**3.*aj ## rxy*aj ## rx-24.*(ds)**3.*aj ## ry*aj ## sy**3.*aj ## rxx-18.*(ds)**2.*(dr)**2.*aj ## syy**2.*aj ## sx**2-36.*(ds)**3.*aj ## syy*aj ## ry**2.*aj ## sx**2-24.*(ds)**2.*(dr)**2.*aj ## syyy*aj ## sy*aj ## sx**2-144.*(ds)**2.*(dr)**2.*aj ## sy*aj ## syy*aj ## sxy*aj ## sx-72.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## sx*aj ## sxyy-36.*(ds)**2.*(dr)**2.*aj ## syy*aj ## sy**2.*aj ## sxx-24.*(ds)**2.*(dr)**2.*aj ## sy**3.*aj ## sxxy-18.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## sxx**2-72.*(ds)**2.*(dr)**2.*aj ## sxy**2.*aj ## sx**2+12.*aj ## sx**4.*(ds)**2.*aj ## ry**2-144.*(ds)**3.*aj ## sy*aj ## ry**2.*aj ## sxy*aj ## sx-144.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## sx*aj ## rxy-144.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## sxy*aj ## rx-144.*(ds)**3.*aj ## ry*aj ## syy*aj ## sy*aj ## sx*aj ## rx-72.*(ds)**3.*aj ## ry*aj ## ryy*aj ## sy*aj ## sx**2-72.*(ds)**3.*aj ## sy**2.*aj ## ryy*aj ## sx*aj ## rx-72.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## sxy**2-144.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sxy*aj ## sx*aj ## sxx-72.*(ds)**2.*(dr)**2.*aj ## sy*aj ## sxxy*aj ## sx**2-24.*(ds)**2.*(dr)**2.*aj ## sy**2.*aj ## sxxx*aj ## sx-36.*(ds)**2.*(dr)**2.*aj ## syy*aj ## sx**2.*aj ## sxx-24.*(ds)**2.*(dr)**2.*aj ## sxyy*aj ## sx**3-90.*aj ## sy**2.*(ds)**2.*(dr)**2.*aj ## syy**2+96.*aj ## sx**3.*(ds)**2.*aj ## ry*aj ## sy*aj ## rx+12.*aj ## sy**4.*(ds)**2.*aj ## rx**2-40.*aj ## sy**3.*(ds)**2.*(dr)**2.*aj ## syyy+72.*aj ## sx**2.*(ds)**2.*aj ## sy**2.*aj ## rx**2+72.*aj ## sy**2.*(ds)**2.*aj ## ry**2.*aj ## sx**2-40.*aj ## sx**3.*(ds)**2.*(dr)**2.*aj ## sxxx-90.*aj ## sx**2.*(ds)**2.*(dr)**2.*aj ## sxx**2-12.*aj ## sx**4.*(ds)*(dr)**2.*aj ## syy-12.*aj ## sy**4.*(ds)*(dr)**2.*aj ## sxx+96.*aj ## sy**3.*(ds)**2.*aj ## ry*aj ## rx*aj ## sx-72.*aj ## sx**2.*(ds)*(dr)**2.*aj ## sy**2.*aj ## sxx-96.*aj ## sy**3.*(ds)*(dr)**2.*aj ## sxy*aj ## sx-72.*aj ## sy**2.*(ds)*(dr)**2.*aj ## syy*aj ## sx**2)/(ds)**6/(dr)**2 )
78 
79 #defineMacro lapCubedCoeff2dOrder2Ghost03(aj,dr,ds) ( 1/2.*(2.*aj ## sx**6-15.*aj ## syy*aj ## sy**4.*(ds)-24.*aj ## sy**3.*(ds)*aj ## sxy*aj ## sx-3.*aj ## sy**4.*(ds)*aj ## sxx-18.*aj ## sy**2.*(ds)*aj ## syy*aj ## sx**2-15.*aj ## sx**4.*aj ## sxx*(ds)+6.*aj ## sy**2.*aj ## sx**4+6.*aj ## sy**4.*aj ## sx**2-24.*aj ## sx**3.*(ds)*aj ## sy*aj ## sxy-18.*aj ## sx**2.*(ds)*aj ## sy**2.*aj ## sxx-3.*aj ## sx**4.*(ds)*aj ## syy+2.*aj ## sy**6)/(ds)**6 )
80 
81 #defineMacro xxxCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*(-2.*aj ## rx**3.*(ds)**2-6.*aj ## rx*aj ## sx**2.*(dr)**2-6.*aj ## rx*aj ## rxx*(dr)*(ds)**2+aj ## rxxx*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
82 
83 #defineMacro xxxCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## rx**3/(dr)**3 )
84 
85 #defineMacro xxxCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(-6.*aj ## rx**2.*aj ## sx*(ds)**2-2.*aj ## sx**3.*(dr)**2-6.*aj ## sxx*aj ## sx*(dr)**2.*(ds)+aj ## sxxx*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
86 
87 #defineMacro xxxCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sx**3/(ds)**3 )
88 
89 #defineMacro xxyCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*(-2.*aj ## ry*aj ## rx**2.*(ds)**2-2.*aj ## sx**2.*(dr)**2.*aj ## ry-4.*aj ## sx*(dr)**2.*aj ## sy*aj ## rx-4.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-2.*(dr)*(ds)**2.*aj ## ry*aj ## rxx+aj ## rxxy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
90 
91 #defineMacro xxyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## ry*aj ## rx**2/(dr)**3 )
92 
93 #defineMacro xxyCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(-2.*aj ## rx**2.*(ds)**2.*aj ## sy-4.*aj ## rx*(ds)**2.*aj ## ry*aj ## sx-2.*aj ## sy*aj ## sx**2.*(dr)**2-2.*(dr)**2.*(ds)*aj ## sy*aj ## sxx-4.*(dr)**2.*(ds)*aj ## sxy*aj ## sx+aj ## sxxy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
94 
95 #defineMacro xxyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sy*aj ## sx**2/(ds)**3 )
96 
97 #defineMacro xyyCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*(-2.*aj ## ry**2.*aj ## rx*(ds)**2-2.*aj ## sy**2.*(dr)**2.*aj ## rx-4.*aj ## sy*(dr)**2.*aj ## ry*aj ## sx-2.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-4.*(dr)*(ds)**2.*aj ## ry*aj ## rxy+aj ## rxyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
98 
99 #defineMacro xyyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## ry**2.*aj ## rx/(dr)**3 )
100 
101 #defineMacro xyyCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(-2.*aj ## ry**2.*(ds)**2.*aj ## sx-4.*aj ## ry*(ds)**2.*aj ## sy*aj ## rx-2.*aj ## sy**2.*aj ## sx*(dr)**2-2.*(dr)**2.*(ds)*aj ## syy*aj ## sx-4.*(dr)**2.*(ds)*aj ## sy*aj ## sxy+aj ## sxyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
102 
103 #defineMacro xyyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sy**2.*aj ## sx/(ds)**3 )
104 
105 #defineMacro yyyCoeff2dOrder2Ghost10(aj,dr,ds) ( -1/2.*(-2.*aj ## ry**3.*(ds)**2-6.*aj ## ry*aj ## sy**2.*(dr)**2-6.*aj ## ry*aj ## ryy*(dr)*(ds)**2+aj ## ryyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
106 
107 #defineMacro yyyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## ry**3/(dr)**3 )
108 
109 #defineMacro yyyCoeff2dOrder2Ghost01(aj,dr,ds) ( -1/2.*(-6.*aj ## ry**2.*aj ## sy*(ds)**2-2.*aj ## sy**3.*(dr)**2-6.*aj ## syy*aj ## sy*(dr)**2.*(ds)+aj ## syyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
110 
111 #defineMacro yyyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sy**3/(ds)**3 )
112 
113 #defineMacro xxxxxCoeff2dOrder2Ghost10(aj,dr,ds) ( 5/2.*(-aj ## rx**5.*(ds)**4-8.*aj ## sx**2.*aj ## rx**3.*(dr)**2.*(ds)**2-6.*aj ## sx**4.*aj ## rx*(dr)**4-16.*aj ## rxx*aj ## rx**3.*(dr)*(ds)**4-24.*aj ## rx*aj ## sx**2.*(dr)**3.*(ds)**2.*aj ## rxx-24.*aj ## rx**2.*aj ## sx*(dr)**3.*(ds)**2.*aj ## sxx+4.*aj ## rx**2.*(dr)**2.*(ds)**4.*aj ## rxxx+6.*aj ## rx*(dr)**2.*(ds)**4.*aj ## rxx**2+12.*(dr)**4.*(ds)**2.*aj ## rxx*aj ## sx*aj ## sxx+8.*(dr)**4.*(ds)**2.*aj ## sxxx*aj ## sx*aj ## rx+6.*(dr)**4.*(ds)**2.*aj ## rx*aj ## sxx**2+4.*(dr)**4.*(ds)**2.*aj ## rxxx*aj ## sx**2+4.*(dr)**3.*(ds)**4.*aj ## rxx*aj ## rxxx+2.*(dr)**3.*(ds)**4.*aj ## rx*aj ## rxxxx)/(dr)**5/(ds)**4 )
114 
115 #defineMacro xxxxxCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## rx*(-4.*aj ## rx**4.*(ds)**2-20.*aj ## sx**2.*aj ## rx**2.*(dr)**2-20.*aj ## rxx*aj ## rx**2.*(dr)*(ds)**2+10.*(dr)**2.*(ds)**2.*aj ## rx*aj ## rxxx+15.*(dr)**2.*(ds)**2.*aj ## rxx**2)/(dr)**5/(ds)**2 )
116 
117 #defineMacro xxxxxCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## rx**5/(dr)**5 )
118 
119 #defineMacro xxxxxCoeff2dOrder2Ghost01(aj,dr,ds) ( 5/2.*(-6.*aj ## sx*aj ## rx**4.*(ds)**4-8.*aj ## sx**3.*aj ## rx**2.*(dr)**2.*(ds)**2-aj ## sx**5.*(dr)**4-24.*aj ## rx*aj ## sx**2.*(dr)**2.*(ds)**3.*aj ## rxx-24.*aj ## rx**2.*aj ## sx*(dr)**2.*(ds)**3.*aj ## sxx-16.*aj ## sxx*aj ## sx**3.*(dr)**4.*(ds)+4.*(dr)**2.*(ds)**4.*aj ## sxxx*aj ## rx**2+8.*(dr)**2.*(ds)**4.*aj ## rxxx*aj ## sx*aj ## rx+6.*(dr)**2.*(ds)**4.*aj ## sx*aj ## rxx**2+12.*(dr)**2.*(ds)**4.*aj ## sxx*aj ## rx*aj ## rxx+4.*aj ## sx**2.*(dr)**4.*(ds)**2.*aj ## sxxx+6.*aj ## sx*(dr)**4.*(ds)**2.*aj ## sxx**2+2.*(dr)**4.*(ds)**3.*aj ## sxxxx*aj ## sx+4.*(dr)**4.*(ds)**3.*aj ## sxx*aj ## sxxx)/(dr)**4/(ds)**5 )
120 
121 #defineMacro xxxxxCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sx*(-20.*aj ## sx**2.*aj ## rx**2.*(ds)**2-4.*aj ## sx**4.*(dr)**2-20.*aj ## sxx*aj ## sx**2.*(dr)**2.*(ds)+10.*(dr)**2.*(ds)**2.*aj ## sxxx*aj ## sx+15.*(dr)**2.*(ds)**2.*aj ## sxx**2)/(dr)**2/(ds)**5 )
122 
123 #defineMacro xxxxxCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sx**5/(ds)**5 )
124 
125 #defineMacro xxxxyCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(-5.*aj ## ry*aj ## rx**4.*(ds)**4-6.*aj ## sx**4.*(dr)**4.*aj ## ry-48.*aj ## rx**2.*(dr)*(ds)**4.*aj ## ry*aj ## rxx-48.*(dr)**3.*(ds)**2.*aj ## sxy*aj ## rx**2.*aj ## sx-48.*(dr)**3.*(ds)**2.*aj ## rxy*aj ## rx*aj ## sx**2-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sx*aj ## rx*aj ## rxx-24.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rx**2-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rx*aj ## sxx*aj ## sx-16.*aj ## rx**3.*aj ## sx*(dr)**2.*(ds)**2.*aj ## sy-24.*aj ## rx**2.*aj ## sx**2.*(dr)**2.*(ds)**2.*aj ## ry-24.*aj ## sx**3.*(dr)**4.*aj ## sy*aj ## rx-32.*aj ## rx**3.*(dr)*(ds)**4.*aj ## rxy+24.*(dr)**4.*(ds)**2.*aj ## rxx*aj ## sxy*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## sxxy*aj ## sx*aj ## rx+24.*(dr)**4.*(ds)**2.*aj ## rxy*aj ## sx*aj ## sxx+24.*(dr)**4.*(ds)**2.*aj ## sxx*aj ## sxy*aj ## rx+8.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxx*aj ## rx+12.*(dr)**3.*(ds)**4.*aj ## rxxy*aj ## rxx+24.*(dr)**2.*(ds)**4.*aj ## rxy*aj ## rxx*aj ## rx+8.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rx*aj ## rxxx-24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxx*aj ## sx**2+12.*(dr)**2.*(ds)**4.*aj ## rxxy*aj ## rx**2+6.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxx**2+6.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxx**2+12.*(dr)**4.*(ds)**2.*aj ## rxxy*aj ## sx**2+2.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxxx+8.*(dr)**3.*(ds)**4.*aj ## rxxx*aj ## rxy+8.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxxx+8.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxxx*aj ## sx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rxx+8.*(dr)**3.*(ds)**4.*aj ## rxxxy*aj ## rx)/(dr)**5/(ds)**4 )
126 
127 #defineMacro xxxxyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*(-4.*aj ## ry*aj ## rx**4.*(ds)**2-8.*aj ## rx**3.*aj ## sx*(dr)**2.*aj ## sy-12.*aj ## rx**2.*aj ## sx**2.*(dr)**2.*aj ## ry-8.*aj ## rx**3.*(dr)*(ds)**2.*aj ## rxy-12.*aj ## rx**2.*(dr)*(ds)**2.*aj ## ry*aj ## rxx+6.*(dr)**2.*(ds)**2.*aj ## rxxy*aj ## rx**2+12.*(dr)**2.*(ds)**2.*aj ## rxy*aj ## rxx*aj ## rx+4.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## rxxx+3.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxx**2)/(dr)**5/(ds)**2 )
128 
129 #defineMacro xxxxyCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## ry*aj ## rx**4/(dr)**5 )
130 
131 #defineMacro xxxxyCoeff2dOrder2Ghost01(aj,dr,ds) ( 1/2.*(-24.*aj ## rx**3.*(ds)**4.*aj ## ry*aj ## sx-32.*aj ## sx**3.*(dr)**4.*(ds)*aj ## sxy+12.*(dr)**2.*(ds)**4.*aj ## sxxy*aj ## rx**2+6.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rxx**2-6.*aj ## rx**4.*(ds)**4.*aj ## sy-24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxx*aj ## rx**2-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rx*aj ## sxx*aj ## sx-24.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rxx*aj ## sx**2-48.*aj ## sx**2.*(dr)**4.*(ds)*aj ## sy*aj ## sxx+24.*(dr)**2.*(ds)**4.*aj ## sxx*aj ## rxy*aj ## rx+24.*(dr)**2.*(ds)**4.*aj ## rxxy*aj ## sx*aj ## rx+8.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rx*aj ## rxxx-16.*aj ## rx*aj ## sx**3.*(dr)**2.*(ds)**2.*aj ## ry-24.*aj ## rx**2.*aj ## sx**2.*(dr)**2.*(ds)**2.*aj ## sy-5.*aj ## sy*aj ## sx**4.*(dr)**4+24.*(dr)**2.*(ds)**4.*aj ## rxy*aj ## sx*aj ## rxx+8.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxxx*aj ## rx+8.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sx*aj ## rxxx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxx*aj ## rxx+24.*(dr)**2.*(ds)**4.*aj ## sxy*aj ## rxx*aj ## rx+8.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxx*aj ## sx+6.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxx**2+8.*(dr)**4.*(ds)**3.*aj ## sxxx*aj ## sxy+8.*(dr)**4.*(ds)**3.*aj ## sxxxy*aj ## sx+12.*(dr)**4.*(ds)**3.*aj ## sxx*aj ## sxxy+2.*(dr)**4.*(ds)**3.*aj ## sy*aj ## sxxxx+24.*(dr)**4.*(ds)**2.*aj ## sxy*aj ## sx*aj ## sxx+12.*(dr)**4.*(ds)**2.*aj ## sxxy*aj ## sx**2-48.*(dr)**2.*(ds)**3.*aj ## sxy*aj ## rx**2.*aj ## sx-48.*(dr)**2.*(ds)**3.*aj ## rxy*aj ## rx*aj ## sx**2-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sx*aj ## rx*aj ## rxx)/(dr)**4/(ds)**5 )
132 
133 #defineMacro xxxxyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*(-8.*aj ## rx*aj ## sx**3.*(ds)**2.*aj ## ry-12.*aj ## rx**2.*aj ## sx**2.*(ds)**2.*aj ## sy-4.*aj ## sy*aj ## sx**4.*(dr)**2-8.*aj ## sx**3.*(dr)**2.*(ds)*aj ## sxy-12.*aj ## sx**2.*(dr)**2.*(ds)*aj ## sy*aj ## sxx+12.*(dr)**2.*(ds)**2.*aj ## sxy*aj ## sx*aj ## sxx+6.*(dr)**2.*(ds)**2.*aj ## sxxy*aj ## sx**2+4.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxxx*aj ## sx+3.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxx**2)/(dr)**2/(ds)**5 )
134 
135 #defineMacro xxxxyCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sy*aj ## sx**4/(ds)**5 )
136 
137 #defineMacro xxxyyCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(-12.*aj ## sy*aj ## sx**3.*(dr)**4.*aj ## ry-18.*aj ## sy**2.*aj ## sx**2.*(dr)**4.*aj ## rx-8.*aj ## rx**3.*(dr)*(ds)**4.*aj ## ryy+6.*(dr)**3.*(ds)**4.*aj ## rxx*aj ## rxyy+12.*(dr)**3.*(ds)**4.*aj ## rxy*aj ## rxxy+6.*(dr)**3.*(ds)**4.*aj ## rx*aj ## rxxyy+12.*(dr)**4.*(ds)**2.*aj ## sxy**2.*aj ## rx+2.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxxx+4.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxxy-24.*aj ## rx*(dr)*(ds)**4.*aj ## ry**2.*aj ## rxx-48.*aj ## rx**2.*(dr)*(ds)**4.*aj ## ry*aj ## rxy-24.*aj ## rx**2.*(dr)**2.*(ds)**2.*aj ## ry*aj ## sy*aj ## sx-4.*aj ## rx**3.*(dr)**2.*(ds)**2.*aj ## sy**2-12.*aj ## rx*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## sx**2-5.*aj ## ry**2.*aj ## rx**3.*(ds)**4-24.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxx*aj ## rx+6.*(dr)**4.*(ds)**2.*aj ## rxyy*aj ## sx**2-24.*(dr)**3.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rx**2-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## rxy*aj ## sx*aj ## rx+2.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxxx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rx*aj ## sxy*aj ## sx+2.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxxx-12.*(dr)**3.*(ds)**2.*aj ## sy**2.*aj ## rx*aj ## rxx-24.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rxx+12.*(dr)**2.*(ds)**4.*aj ## rxy**2.*aj ## rx-12.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## sxx*aj ## sx-24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## rxy*aj ## sx**2-12.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## rx*aj ## sx**2-12.*(dr)**3.*(ds)**2.*aj ## syy*aj ## rx**2.*aj ## sx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxxy+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxx*aj ## rxy+12.*(dr)**4.*(ds)**2.*aj ## sxyy*aj ## sx*aj ## rx+6.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## sxx*aj ## sx+6.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## rx*aj ## rxx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxy*aj ## rxx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rx*aj ## rxxy+6.*(dr)**2.*(ds)**4.*aj ## rxyy*aj ## rx**2+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rxx+12.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxy*aj ## sxx+6.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sxx*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxxy*aj ## sx+4.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxxx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxy*aj ## rx+6.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sx*aj ## rxx+24.*(dr)**4.*(ds)**2.*aj ## rxy*aj ## sxy*aj ## sx)/(dr)**5/(ds)**4 )
138 
139 #defineMacro xxxyyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*(-4.*aj ## ry**2.*aj ## rx**3.*(ds)**2-6.*aj ## rx*(dr)**2.*aj ## ry**2.*aj ## sx**2-12.*aj ## rx**2.*(dr)**2.*aj ## ry*aj ## sy*aj ## sx-2.*aj ## rx**3.*(dr)**2.*aj ## sy**2-12.*aj ## rx**2.*(dr)*(ds)**2.*aj ## ry*aj ## rxy-6.*aj ## rx*(dr)*(ds)**2.*aj ## ry**2.*aj ## rxx-2.*aj ## rx**3.*(dr)*(ds)**2.*aj ## ryy+6.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## rxxy+6.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxy*aj ## rxx+3.*(dr)**2.*(ds)**2.*aj ## rxyy*aj ## rx**2+(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxxx+6.*(dr)**2.*(ds)**2.*aj ## rxy**2.*aj ## rx+3.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rx*aj ## rxx)/(dr)**5/(ds)**2 )
140 
141 #defineMacro xxxyyCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## ry**2.*aj ## rx**3/(dr)**5 )
142 
143 #defineMacro xxxyyCoeff2dOrder2Ghost01(aj,dr,ds) ( 1/2.*(-12.*aj ## sx*(dr)**2.*(ds)**2.*aj ## sy**2.*aj ## rx**2-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## rxy*aj ## sx*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## sx*aj ## rxy**2-24.*aj ## sx**2.*(dr)**2.*(ds)**2.*aj ## ry*aj ## sy*aj ## rx-5.*aj ## sy**2.*aj ## sx**3.*(dr)**4-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rx*aj ## sxy*aj ## sx-24.*(dr)**2.*(ds)**3.*aj ## ry*aj ## rxy*aj ## sx**2-24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry*aj ## sx*aj ## rxx-12.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## rx*aj ## rxx-8.*aj ## sx**3.*(dr)**4.*(ds)*aj ## syy-24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## sxy*aj ## rx**2-24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry*aj ## sxx*aj ## rx+6.*(dr)**2.*(ds)**4.*aj ## sxyy*aj ## rx**2-24.*aj ## sx*(dr)**4.*(ds)*aj ## sy**2.*aj ## sxx+6.*(dr)**2.*(ds)**4.*aj ## syy*aj ## rx*aj ## rxx-12.*(dr)**2.*(ds)**3.*aj ## ryy*aj ## rx*aj ## sx**2-48.*aj ## sx**2.*(dr)**4.*(ds)*aj ## sy*aj ## sxy-12.*(dr)**2.*(ds)**3.*aj ## ry**2.*aj ## sxx*aj ## sx-12.*(dr)**2.*(ds)**3.*aj ## syy*aj ## rx**2.*aj ## sx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxy*aj ## rxx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxx*aj ## rxy+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxxy*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sx*aj ## rxxy+24.*(dr)**2.*(ds)**4.*aj ## sxy*aj ## rxy*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## rxyy*aj ## sx*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rx*aj ## rxxy+12.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rxy*aj ## rxx+6.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sxx*aj ## sx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxxy*aj ## sx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy*aj ## sxx+4.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sy*aj ## rxxx+12.*(dr)**4.*(ds)**2.*aj ## sxy**2.*aj ## sx+6.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## sx*aj ## rxx+6.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## sxx*aj ## rx-18.*aj ## ry**2.*aj ## rx**2.*(ds)**4.*aj ## sx-12.*aj ## ry*aj ## rx**3.*(ds)**4.*aj ## sy-4.*aj ## sx**3.*(dr)**2.*(ds)**2.*aj ## ry**2+2.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## sxxx+2.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## sxxx+6.*(dr)**4.*(ds)**3.*aj ## sxyy*aj ## sxx+6.*(dr)**4.*(ds)**3.*aj ## sxxyy*aj ## sx+6.*(dr)**4.*(ds)**2.*aj ## sxyy*aj ## sx**2+4.*(dr)**4.*(ds)**3.*aj ## sy*aj ## sxxxy+12.*(dr)**4.*(ds)**3.*aj ## sxy*aj ## sxxy+2.*(dr)**4.*(ds)**3.*aj ## syy*aj ## sxxx)/(dr)**4/(ds)**5 )
144 
145 #defineMacro xxxyyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*(-12.*aj ## sx**2.*(ds)**2.*aj ## ry*aj ## sy*aj ## rx-2.*aj ## sx**3.*(ds)**2.*aj ## ry**2-6.*aj ## sx*(ds)**2.*aj ## sy**2.*aj ## rx**2-4.*aj ## sy**2.*aj ## sx**3.*(dr)**2-12.*aj ## sx**2.*(dr)**2.*(ds)*aj ## sy*aj ## sxy-2.*aj ## sx**3.*(dr)**2.*(ds)*aj ## syy-6.*aj ## sx*(dr)**2.*(ds)*aj ## sy**2.*aj ## sxx+6.*(dr)**2.*(ds)**2.*aj ## sxy**2.*aj ## sx+6.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxxy*aj ## sx+(dr)**2.*(ds)**2.*aj ## sy**2.*aj ## sxxx+6.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxy*aj ## sxx+3.*(dr)**2.*(ds)**2.*aj ## sxyy*aj ## sx**2+3.*(dr)**2.*(ds)**2.*aj ## syy*aj ## sxx*aj ## sx)/(dr)**2/(ds)**5 )
146 
147 #defineMacro xxxyyCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sy**2.*aj ## sx**3/(ds)**5 )
148 
149 #defineMacro xxyyyCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(-24.*(dr)**3.*(ds)**2.*aj ## sy**2.*aj ## rxy*aj ## rx-24.*(dr)**3.*(ds)**2.*aj ## ry*aj ## syy*aj ## sx*aj ## rx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxy*aj ## rx-12.*(dr)**3.*(ds)**2.*aj ## syy*aj ## sy*aj ## rx**2-12.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy**2.*aj ## rxx-12.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry**2.*aj ## sxx-12.*(dr)**3.*(ds)**2.*aj ## ry*aj ## ryy*aj ## sx**2-48.*aj ## ry**2.*(dr)*(ds)**4.*aj ## rxy*aj ## rx-24.*aj ## ry*(dr)*(ds)**4.*aj ## ryy*aj ## rx**2-24.*(dr)**3.*(ds)**2.*aj ## ry**2.*aj ## sxy*aj ## sx-24.*aj ## ry**2.*(dr)**2.*(ds)**2.*aj ## sy*aj ## rx*aj ## sx-12.*aj ## ry*(dr)**2.*(ds)**2.*aj ## sy**2.*aj ## rx**2-5.*aj ## ry**3.*aj ## rx**2.*(ds)**4+12.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## rxy*aj ## rx+6.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## rxx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rx*aj ## rxyy-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy*aj ## sx*aj ## rxy-24.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sx*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sx*aj ## rxyy+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy*aj ## rxy+6.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ry*aj ## sxx+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sxy*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sx*aj ## sxyy+6.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sxx+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxyy*aj ## rx+4.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sx*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## sxy*aj ## sx+6.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sy*aj ## rxx+12.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sy*aj ## sxxy+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sx*aj ## rxy-4.*aj ## ry**3.*(dr)**2.*(ds)**2.*aj ## sx**2+6.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxxy-8.*aj ## ry**3.*(dr)*(ds)**4.*aj ## rxx-18.*aj ## sy**2.*aj ## sx**2.*(dr)**4.*aj ## ry-12.*aj ## sy**3.*aj ## sx*(dr)**4.*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## rxy**2+2.*(dr)**4.*(ds)**2.*aj ## ryyy*aj ## sx**2+6.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxxy+2.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## rx**2+12.*(dr)**4.*(ds)**2.*aj ## ry*aj ## sxy**2+12.*(dr)**3.*(ds)**4.*aj ## rxy*aj ## rxyy+6.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxxyy+2.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## rxx+6.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxxy+4.*(dr)**3.*(ds)**4.*aj ## rx*aj ## rxyyy)/(dr)**5/(ds)**4 )
150 
151 #defineMacro xxyyyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*(-4.*aj ## ry**3.*aj ## rx**2.*(ds)**2-12.*aj ## ry**2.*(dr)**2.*aj ## sy*aj ## rx*aj ## sx-2.*aj ## ry**3.*(dr)**2.*aj ## sx**2-6.*aj ## ry*(dr)**2.*aj ## sy**2.*aj ## rx**2-12.*aj ## ry**2.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-2.*aj ## ry**3.*(dr)*(ds)**2.*aj ## rxx-6.*aj ## ry*(dr)*(ds)**2.*aj ## ryy*aj ## rx**2+6.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rxy*aj ## rx+3.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryy*aj ## rxx+6.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## rxyy+3.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxxy+6.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxy**2+(dr)**2.*(ds)**2.*aj ## ryyy*aj ## rx**2)/(dr)**5/(ds)**2 )
152 
153 #defineMacro xxyyyCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## ry**3.*aj ## rx**2/(dr)**5 )
154 
155 #defineMacro xxyyyCoeff2dOrder2Ghost01(aj,dr,ds) ( 1/2.*(-5.*aj ## sy**3.*aj ## sx**2.*(dr)**4-12.*aj ## sy*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## sx**2-24.*aj ## sy**2.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## sx-4.*aj ## sy**3.*(dr)**2.*(ds)**2.*aj ## rx**2-12.*aj ## ry**3.*aj ## rx*(ds)**4.*aj ## sx-18.*aj ## ry**2.*aj ## rx**2.*(ds)**4.*aj ## sy-24.*(dr)**2.*(ds)**3.*aj ## ry**2.*aj ## sxy*aj ## sx-12.*(dr)**2.*(ds)**3.*aj ## ry*aj ## ryy*aj ## sx**2-12.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sy*aj ## rx**2-12.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## rxx-12.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry**2.*aj ## sxx-24.*(dr)**2.*(ds)**3.*aj ## sy**2.*aj ## rxy*aj ## rx-24.*(dr)**2.*(ds)**3.*aj ## ry*aj ## syy*aj ## sx*aj ## rx-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy*aj ## sxy*aj ## rx-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy*aj ## sx*aj ## rxy-24.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ryy*aj ## sx*aj ## rx-8.*aj ## sy**3.*(dr)**4.*(ds)*aj ## sxx-48.*aj ## sy**2.*(dr)**4.*(ds)*aj ## sxy*aj ## sx-24.*aj ## sy*(dr)**4.*(ds)*aj ## syy*aj ## sx**2+2.*(dr)**2.*(ds)**4.*aj ## syyy*aj ## rx**2+6.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## sxxy+12.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rxy**2+6.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## sxxy+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sxy**2+12.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ry*aj ## rxxy+6.*(dr)**2.*(ds)**4.*aj ## syy*aj ## ry*aj ## rxx+12.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## sx*aj ## rxy+12.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## sxy*aj ## rx+6.*(dr)**4.*(ds)**3.*aj ## syy*aj ## sxxy+12.*(dr)**2.*(ds)**4.*aj ## syy*aj ## rxy*aj ## rx+4.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## sx*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## sy*aj ## rx*aj ## rxyy+6.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ryy*aj ## rxx+4.*(dr)**4.*(ds)**3.*aj ## sx*aj ## sxyyy+12.*(dr)**4.*(ds)**2.*aj ## sy*aj ## sx*aj ## sxyy+6.*(dr)**4.*(ds)**2.*aj ## sy*aj ## syy*aj ## sxx+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sxy*aj ## sx+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxy*aj ## rxy+6.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## sxx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sxyy*aj ## rx+12.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sx*aj ## rxyy+2.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sx**2+2.*(dr)**4.*(ds)**3.*aj ## syyy*aj ## sxx+12.*(dr)**4.*(ds)**3.*aj ## sxyy*aj ## sxy+6.*(dr)**4.*(ds)**3.*aj ## sy*aj ## sxxyy)/(dr)**4/(ds)**5 )
156 
157 #defineMacro xxyyyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*(-6.*aj ## sy*(ds)**2.*aj ## ry**2.*aj ## sx**2-12.*aj ## sy**2.*(ds)**2.*aj ## ry*aj ## rx*aj ## sx-2.*aj ## sy**3.*(ds)**2.*aj ## rx**2-4.*aj ## sy**3.*aj ## sx**2.*(dr)**2-12.*aj ## sy**2.*(dr)**2.*(ds)*aj ## sxy*aj ## sx-2.*aj ## sy**3.*(dr)**2.*(ds)*aj ## sxx-6.*aj ## sy*(dr)**2.*(ds)*aj ## syy*aj ## sx**2+3.*(dr)**2.*(ds)**2.*aj ## sy**2.*aj ## sxxy+6.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sx*aj ## sxyy+3.*(dr)**2.*(ds)**2.*aj ## sy*aj ## syy*aj ## sxx+6.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxy**2+(dr)**2.*(ds)**2.*aj ## syyy*aj ## sx**2+6.*(dr)**2.*(ds)**2.*aj ## syy*aj ## sxy*aj ## sx)/(dr)**2/(ds)**5 )
158 
159 #defineMacro xxyyyCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sy**3.*aj ## sx**2/(ds)**5 )
160 
161 #defineMacro xyyyyCoeff2dOrder2Ghost10(aj,dr,ds) ( 1/2.*(-24.*aj ## sy**3.*(dr)**4.*aj ## ry*aj ## sx-6.*aj ## sy**4.*(dr)**4.*aj ## rx-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## ryy*aj ## sy*aj ## sx-24.*(dr)**3.*(ds)**2.*aj ## ryy*aj ## sy**2.*aj ## rx+8.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## ry*aj ## rx+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## rxy-48.*(dr)**3.*(ds)**2.*aj ## syy*aj ## sy*aj ## ry*aj ## rx-24.*(dr)**3.*(ds)**2.*aj ## syy*aj ## ry**2.*aj ## sx-48.*(dr)**3.*(ds)**2.*aj ## sy*aj ## ry**2.*aj ## sxy-48.*(dr)**3.*(ds)**2.*aj ## ry*aj ## sy**2.*aj ## rxy-48.*aj ## ry**2.*(dr)*(ds)**4.*aj ## ryy*aj ## rx-32.*aj ## ry**3.*(dr)*(ds)**4.*aj ## rxy+12.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## rxyy+6.*(dr)**2.*(ds)**4.*aj ## ryy**2.*aj ## rx+6.*(dr)**4.*(ds)**2.*aj ## syy**2.*aj ## rx+12.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## rxyy-16.*aj ## sy*aj ## ry**3.*(dr)**2.*(ds)**2.*aj ## sx-24.*aj ## sy**2.*aj ## ry**2.*(dr)**2.*(ds)**2.*aj ## rx-5.*aj ## ry**4.*aj ## rx*(ds)**4+8.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## rxy+12.*(dr)**3.*(ds)**4.*aj ## ryy*aj ## rxyy+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ryy*aj ## sxy+24.*(dr)**4.*(ds)**2.*aj ## sy*aj ## ry*aj ## sxyy+8.*(dr)**4.*(ds)**2.*aj ## sy*aj ## syyy*aj ## rx+8.*(dr)**4.*(ds)**2.*aj ## ry*aj ## syyy*aj ## sx+8.*(dr)**3.*(ds)**4.*aj ## ry*aj ## rxyyy+2.*(dr)**3.*(ds)**4.*aj ## ryyyy*aj ## rx+8.*(dr)**4.*(ds)**2.*aj ## ryyy*aj ## sy*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sy*aj ## rxy+12.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ryy*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## syy*aj ## ry*aj ## sxy)/(dr)**5/(ds)**4 )
162 
163 #defineMacro xyyyyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*(-4.*aj ## ry**4.*aj ## rx*(ds)**2-8.*aj ## sy*aj ## ry**3.*(dr)**2.*aj ## sx-12.*aj ## sy**2.*aj ## ry**2.*(dr)**2.*aj ## rx-12.*aj ## ry**2.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-8.*aj ## ry**3.*(dr)*(ds)**2.*aj ## rxy+4.*(dr)**2.*(ds)**2.*aj ## ryyy*aj ## ry*aj ## rx+6.*(dr)**2.*(ds)**2.*aj ## ry**2.*aj ## rxyy+12.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryy*aj ## rxy+3.*(dr)**2.*(ds)**2.*aj ## ryy**2.*aj ## rx)/(dr)**5/(ds)**2 )
164 
165 #defineMacro xyyyyCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## ry**4.*aj ## rx/(dr)**5 )
166 
167 #defineMacro xyyyyCoeff2dOrder2Ghost01(aj,dr,ds) ( 1/2.*(-5.*aj ## sy**4.*aj ## sx*(dr)**4-6.*aj ## ry**4.*(ds)**4.*aj ## sx-24.*aj ## ry**3.*(ds)**4.*aj ## sy*aj ## rx-24.*(dr)**2.*(ds)**3.*aj ## syy*aj ## ry**2.*aj ## sx-48.*(dr)**2.*(ds)**3.*aj ## sy*aj ## ry**2.*aj ## sxy-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## sy**2.*aj ## rxy-48.*(dr)**2.*(ds)**3.*aj ## ry*aj ## ryy*aj ## sy*aj ## sx-24.*(dr)**2.*(ds)**3.*aj ## ryy*aj ## sy**2.*aj ## rx-32.*aj ## sy**3.*(dr)**4.*(ds)*aj ## sxy-16.*aj ## ry*aj ## sy**3.*(dr)**2.*(ds)**2.*aj ## rx-24.*aj ## ry**2.*aj ## sy**2.*(dr)**2.*(ds)**2.*aj ## sx-48.*(dr)**2.*(ds)**3.*aj ## syy*aj ## sy*aj ## ry*aj ## rx-48.*aj ## sy**2.*(dr)**4.*(ds)*aj ## syy*aj ## sx+8.*(dr)**2.*(ds)**4.*aj ## syyy*aj ## ry*aj ## rx+8.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## sy*aj ## rx+8.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## ry*aj ## sx+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## syy*aj ## rxy+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## sy*aj ## rxyy+24.*(dr)**2.*(ds)**4.*aj ## ry*aj ## ryy*aj ## sxy+24.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## sy*aj ## rxy+12.*(dr)**2.*(ds)**4.*aj ## ryy*aj ## syy*aj ## rx+6.*(dr)**2.*(ds)**4.*aj ## ryy**2.*aj ## sx+12.*(dr)**2.*(ds)**4.*aj ## ry**2.*aj ## sxyy+6.*(dr)**4.*(ds)**2.*aj ## syy**2.*aj ## sx+8.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sy*aj ## sx+24.*(dr)**4.*(ds)**2.*aj ## syy*aj ## sy*aj ## sxy+12.*(dr)**4.*(ds)**2.*aj ## sy**2.*aj ## sxyy+8.*(dr)**4.*(ds)**3.*aj ## sy*aj ## sxyyy+2.*(dr)**4.*(ds)**3.*aj ## syyyy*aj ## sx+8.*(dr)**4.*(ds)**3.*aj ## syyy*aj ## sxy+12.*(dr)**4.*(ds)**3.*aj ## syy*aj ## sxyy)/(dr)**4/(ds)**5 )
168 
169 #defineMacro xyyyyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*(-8.*aj ## ry*aj ## sy**3.*(ds)**2.*aj ## rx-12.*aj ## ry**2.*aj ## sy**2.*(ds)**2.*aj ## sx-4.*aj ## sy**4.*aj ## sx*(dr)**2-12.*aj ## sy**2.*(dr)**2.*(ds)*aj ## syy*aj ## sx-8.*aj ## sy**3.*(dr)**2.*(ds)*aj ## sxy+4.*(dr)**2.*(ds)**2.*aj ## syyy*aj ## sy*aj ## sx+12.*(dr)**2.*(ds)**2.*aj ## syy*aj ## sy*aj ## sxy+3.*(dr)**2.*(ds)**2.*aj ## syy**2.*aj ## sx+6.*(dr)**2.*(ds)**2.*aj ## sy**2.*aj ## sxyy)/(dr)**2/(ds)**5 )
170 
171 #defineMacro xyyyyCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sy**4.*aj ## sx/(ds)**5 )
172 
173 #defineMacro yyyyyCoeff2dOrder2Ghost10(aj,dr,ds) ( 5/2.*(-aj ## ry**5.*(ds)**4-8.*aj ## sy**2.*aj ## ry**3.*(dr)**2.*(ds)**2-6.*aj ## sy**4.*aj ## ry*(dr)**4-16.*aj ## ry**3.*aj ## ryy*(dr)*(ds)**4-24.*aj ## ry*aj ## sy**2.*(dr)**3.*(ds)**2.*aj ## ryy-24.*aj ## ry**2.*aj ## sy*(dr)**3.*(ds)**2.*aj ## syy+4.*aj ## ry**2.*(dr)**2.*(ds)**4.*aj ## ryyy+6.*aj ## ry*(dr)**2.*(ds)**4.*aj ## ryy**2+12.*(dr)**4.*(ds)**2.*aj ## ryy*aj ## syy*aj ## sy+4.*(dr)**4.*(ds)**2.*aj ## ryyy*aj ## sy**2+8.*(dr)**4.*(ds)**2.*aj ## syyy*aj ## sy*aj ## ry+6.*(dr)**4.*(ds)**2.*aj ## syy**2.*aj ## ry+2.*(dr)**3.*(ds)**4.*aj ## ry*aj ## ryyyy+4.*(dr)**3.*(ds)**4.*aj ## ryyy*aj ## ryy)/(dr)**5/(ds)**4 )
174 
175 #defineMacro yyyyyCoeff2dOrder2Ghost20(aj,dr,ds) ( -1/2.*aj ## ry*(-4.*aj ## ry**4.*(ds)**2-20.*aj ## sy**2.*aj ## ry**2.*(dr)**2-20.*aj ## ry**2.*aj ## ryy*(dr)*(ds)**2+10.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryyy+15.*(dr)**2.*(ds)**2.*aj ## ryy**2)/(dr)**5/(ds)**2 )
176 
177 #defineMacro yyyyyCoeff2dOrder2Ghost30(aj,dr,ds) ( -1/2.*aj ## ry**5/(dr)**5 )
178 
179 #defineMacro yyyyyCoeff2dOrder2Ghost01(aj,dr,ds) ( 5/2.*(-6.*aj ## sy*aj ## ry**4.*(ds)**4-8.*aj ## sy**3.*aj ## ry**2.*(dr)**2.*(ds)**2-aj ## sy**5.*(dr)**4-24.*aj ## ry*aj ## sy**2.*(dr)**2.*(ds)**3.*aj ## ryy-24.*aj ## ry**2.*aj ## sy*(dr)**2.*(ds)**3.*aj ## syy-16.*aj ## sy**3.*aj ## syy*(dr)**4.*(ds)+12.*(dr)**2.*(ds)**4.*aj ## syy*aj ## ry*aj ## ryy+8.*(dr)**2.*(ds)**4.*aj ## ryyy*aj ## sy*aj ## ry+6.*(dr)**2.*(ds)**4.*aj ## sy*aj ## ryy**2+4.*(dr)**2.*(ds)**4.*aj ## syyy*aj ## ry**2+4.*aj ## sy**2.*(dr)**4.*(ds)**2.*aj ## syyy+6.*aj ## sy*(dr)**4.*(ds)**2.*aj ## syy**2+4.*(dr)**4.*(ds)**3.*aj ## syy*aj ## syyy+2.*(dr)**4.*(ds)**3.*aj ## sy*aj ## syyyy)/(dr)**4/(ds)**5 )
180 
181 #defineMacro yyyyyCoeff2dOrder2Ghost02(aj,dr,ds) ( -1/2.*aj ## sy*(-20.*aj ## sy**2.*aj ## ry**2.*(ds)**2-4.*aj ## sy**4.*(dr)**2-20.*aj ## sy**2.*aj ## syy*(dr)**2.*(ds)+10.*(dr)**2.*(ds)**2.*aj ## syyy*aj ## sy+15.*(dr)**2.*(ds)**2.*aj ## syy**2)/(dr)**2/(ds)**5 )
182 
183 #defineMacro yyyyyCoeff2dOrder2Ghost03(aj,dr,ds) ( -1/2.*aj ## sy**5/(ds)**5 )
184 
185 
186 ! *** orderOfAccuracy = 4
187 
188 #defineMacro xCoeff2dOrder4Ghost10(aj,dr,ds) ( -2/3.*aj ## rx/(dr) )
189 
190 #defineMacro xCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/12.*aj ## rx/(dr) )
191 
192 #defineMacro xCoeff2dOrder4Ghost01(aj,dr,ds) ( -2/3.*aj ## sx/(ds) )
193 
194 #defineMacro xCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/12.*aj ## sx/(ds) )
195 
196 #defineMacro yCoeff2dOrder4Ghost10(aj,dr,ds) ( -2/3.*aj ## ry/(dr) )
197 
198 #defineMacro yCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/12.*aj ## ry/(dr) )
199 
200 #defineMacro yCoeff2dOrder4Ghost01(aj,dr,ds) ( -2/3.*aj ## sy/(ds) )
201 
202 #defineMacro yCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/12.*aj ## sy/(ds) )
203 
204 #defineMacro lapCoeff2dOrder4Ghost10(aj,dr,ds) ( 2/3.*(2.*aj ## rx**2-aj ## rxx*(dr)+2.*aj ## ry**2-aj ## ryy*(dr))/(dr)**2 )
205 
206 #defineMacro lapCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/12.*(-aj ## rx**2+aj ## rxx*(dr)-aj ## ry**2+aj ## ryy*(dr))/(dr)**2 )
207 
208 #defineMacro lapCoeff2dOrder4Ghost01(aj,dr,ds) ( -2/3.*(-2.*aj ## sx**2+aj ## sxx*(ds)-2.*aj ## sy**2+aj ## syy*(ds))/(ds)**2 )
209 
210 #defineMacro lapCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/12.*(-aj ## sx**2+aj ## sxx*(ds)-aj ## sy**2+aj ## syy*(ds))/(ds)**2 )
211 
212 #defineMacro xLapCoeff2dOrder4Ghost10(aj,dr,ds) ( -1/24.*(-39.*aj ## rx**3.*(ds)**2-120.*aj ## rx*aj ## sx**2.*(dr)**2-96.*aj ## rx*aj ## rxx*(dr)*(ds)**2+16.*aj ## rxxx*(dr)**2.*(ds)**2-39.*aj ## ry**2.*aj ## rx*(ds)**2-80.*aj ## sy*(dr)**2.*aj ## ry*aj ## sx-40.*aj ## sy**2.*(dr)**2.*aj ## rx-32.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-64.*(dr)*(ds)**2.*aj ## ry*aj ## rxy+16.*aj ## rxyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
213 
214 #defineMacro xLapCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/24.*(-24.*aj ## rx**3.*(ds)**2-15.*aj ## rx*aj ## sx**2.*(dr)**2-6.*aj ## rx*aj ## rxx*(dr)*(ds)**2+2.*aj ## rxxx*(dr)**2.*(ds)**2-24.*aj ## ry**2.*aj ## rx*(ds)**2-10.*aj ## sy*(dr)**2.*aj ## ry*aj ## sx-5.*aj ## sy**2.*(dr)**2.*aj ## rx-2.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-4.*(dr)*(ds)**2.*aj ## ry*aj ## rxy+2.*aj ## rxyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
215 
216 #defineMacro xLapCoeff2dOrder4Ghost30(aj,dr,ds) ( 1/8.*aj ## rx*(aj ## rx**2+aj ## ry**2)/(dr)**3 )
217 
218 #defineMacro xLapCoeff2dOrder4Ghost01(aj,dr,ds) ( -1/24.*(-120.*aj ## rx**2.*aj ## sx*(ds)**2-39.*aj ## sx**3.*(dr)**2-96.*aj ## sxx*aj ## sx*(dr)**2.*(ds)+16.*aj ## sxxx*(dr)**2.*(ds)**2-80.*aj ## ry*(ds)**2.*aj ## sy*aj ## rx-40.*aj ## ry**2.*(ds)**2.*aj ## sx-39.*aj ## sy**2.*aj ## sx*(dr)**2-32.*(dr)**2.*(ds)*aj ## syy*aj ## sx-64.*(dr)**2.*(ds)*aj ## sy*aj ## sxy+16.*aj ## sxyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
219 
220 #defineMacro xLapCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/24.*(-15.*aj ## rx**2.*aj ## sx*(ds)**2-24.*aj ## sx**3.*(dr)**2-6.*aj ## sxx*aj ## sx*(dr)**2.*(ds)+2.*aj ## sxxx*(dr)**2.*(ds)**2-10.*aj ## ry*(ds)**2.*aj ## sy*aj ## rx-5.*aj ## ry**2.*(ds)**2.*aj ## sx-24.*aj ## sy**2.*aj ## sx*(dr)**2-2.*(dr)**2.*(ds)*aj ## syy*aj ## sx-4.*(dr)**2.*(ds)*aj ## sy*aj ## sxy+2.*aj ## sxyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
221 
222 #defineMacro xLapCoeff2dOrder4Ghost03(aj,dr,ds) ( 1/8.*aj ## sx*(aj ## sx**2+aj ## sy**2)/(ds)**3 )
223 
224 #defineMacro yLapCoeff2dOrder4Ghost10(aj,dr,ds) ( -1/24.*(-39.*aj ## ry*aj ## rx**2.*(ds)**2-80.*aj ## sx*(dr)**2.*aj ## sy*aj ## rx-40.*aj ## sx**2.*(dr)**2.*aj ## ry-64.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-32.*(dr)*(ds)**2.*aj ## ry*aj ## rxx+16.*aj ## rxxy*(dr)**2.*(ds)**2-39.*aj ## ry**3.*(ds)**2-120.*aj ## ry*aj ## sy**2.*(dr)**2-96.*aj ## ry*aj ## ryy*(dr)*(ds)**2+16.*aj ## ryyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
225 
226 #defineMacro yLapCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/24.*(-24.*aj ## ry*aj ## rx**2.*(ds)**2-10.*aj ## sx*(dr)**2.*aj ## sy*aj ## rx-5.*aj ## sx**2.*(dr)**2.*aj ## ry-4.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-2.*(dr)*(ds)**2.*aj ## ry*aj ## rxx+2.*aj ## rxxy*(dr)**2.*(ds)**2-24.*aj ## ry**3.*(ds)**2-15.*aj ## ry*aj ## sy**2.*(dr)**2-6.*aj ## ry*aj ## ryy*(dr)*(ds)**2+2.*aj ## ryyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
227 
228 #defineMacro yLapCoeff2dOrder4Ghost30(aj,dr,ds) ( 1/8.*aj ## ry*(aj ## rx**2+aj ## ry**2)/(dr)**3 )
229 
230 #defineMacro yLapCoeff2dOrder4Ghost01(aj,dr,ds) ( -1/24.*(-80.*aj ## rx*(ds)**2.*aj ## ry*aj ## sx-40.*aj ## rx**2.*(ds)**2.*aj ## sy-39.*aj ## sy*aj ## sx**2.*(dr)**2-32.*(dr)**2.*(ds)*aj ## sy*aj ## sxx-64.*(dr)**2.*(ds)*aj ## sxy*aj ## sx+16.*aj ## sxxy*(dr)**2.*(ds)**2-120.*aj ## ry**2.*aj ## sy*(ds)**2-39.*aj ## sy**3.*(dr)**2-96.*aj ## syy*aj ## sy*(dr)**2.*(ds)+16.*aj ## syyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
231 
232 #defineMacro yLapCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/24.*(-10.*aj ## rx*(ds)**2.*aj ## ry*aj ## sx-5.*aj ## rx**2.*(ds)**2.*aj ## sy-24.*aj ## sy*aj ## sx**2.*(dr)**2-2.*(dr)**2.*(ds)*aj ## sy*aj ## sxx-4.*(dr)**2.*(ds)*aj ## sxy*aj ## sx+2.*aj ## sxxy*(dr)**2.*(ds)**2-15.*aj ## ry**2.*aj ## sy*(ds)**2-24.*aj ## sy**3.*(dr)**2-6.*aj ## syy*aj ## sy*(dr)**2.*(ds)+2.*aj ## syyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
233 
234 #defineMacro yLapCoeff2dOrder4Ghost03(aj,dr,ds) ( 1/8.*aj ## sy*(aj ## sx**2+aj ## sy**2)/(ds)**3 )
235 
236 #defineMacro lapSqCoeff2dOrder4Ghost10(aj,dr,ds) ( -1/12.*(8.*aj ## ryyyy*(dr)**3.*(ds)**2+16.*aj ## rxxyy*(dr)**3.*(ds)**2+8.*aj ## rxxxx*(dr)**3.*(ds)**2+78.*aj ## rx**4.*(ds)**2+78.*aj ## ry**4.*(ds)**2-160.*(dr)**3.*aj ## sy*aj ## sx*aj ## rxy-80.*(dr)**3.*aj ## sy*aj ## ry*aj ## sxx-160.*(dr)**3.*aj ## ry*aj ## sxy*aj ## sx-160.*(dr)**3.*aj ## sy*aj ## sxy*aj ## rx-240.*aj ## sy*(dr)**3.*aj ## syy*aj ## ry-80.*(dr)**3.*aj ## syy*aj ## sx*aj ## rx-64.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxxy-32.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rxx-64.*(dr)**2.*(ds)**2.*aj ## rx*aj ## rxyy-240.*aj ## sx*(dr)**3.*aj ## sxx*aj ## rx-64.*(dr)**2.*(ds)**2.*aj ## rx*aj ## rxxx-39.*(dr)*(ds)**2.*aj ## ry**2.*aj ## rxx-48.*(dr)**2.*(ds)**2.*aj ## rxx**2-64.*(dr)**2.*(ds)**2.*aj ## rxy**2-64.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryyy-120.*aj ## sx**2.*(dr)**3.*aj ## rxx-120.*aj ## sy**2.*(dr)**3.*aj ## ryy+80.*(dr)**2.*aj ## ry**2.*aj ## sx**2+80.*(dr)**2.*aj ## sy**2.*aj ## rx**2-40.*(dr)**3.*aj ## ryy*aj ## sx**2-40.*(dr)**3.*aj ## sy**2.*aj ## rxx-39.*(dr)*(ds)**2.*aj ## ryy*aj ## rx**2-48.*(dr)**2.*(ds)**2.*aj ## ryy**2+240.*aj ## ry**2.*aj ## sy**2.*(dr)**2+240.*aj ## rx**2.*aj ## sx**2.*(dr)**2-117.*aj ## rx**2.*aj ## rxx*(dr)*(ds)**2-117.*aj ## ry**2.*aj ## ryy*(dr)*(ds)**2+156.*aj ## ry**2.*aj ## rx**2.*(ds)**2+320.*(dr)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rx-156.*(dr)*(ds)**2.*aj ## ry*aj ## rxy*aj ## rx)/(dr)**4/(ds)**2 )
237 
238 #defineMacro lapSqCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/12.*(-5.*(dr)**3.*aj ## ryy*aj ## sx**2-5.*(dr)**3.*aj ## sy**2.*aj ## rxx-15.*aj ## sx**2.*(dr)**3.*aj ## rxx-15.*aj ## sy**2.*(dr)**3.*aj ## ryy+5.*(dr)**2.*aj ## sy**2.*aj ## rx**2-4.*(dr)**2.*(ds)**2.*aj ## rxy**2+5.*(dr)**2.*aj ## ry**2.*aj ## sx**2-3.*(dr)**2.*(ds)**2.*aj ## ryy**2-3.*(dr)**2.*(ds)**2.*aj ## rxx**2+aj ## rxxxx*(dr)**3.*(ds)**2+aj ## ryyyy*(dr)**3.*(ds)**2+2.*aj ## rxxyy*(dr)**3.*(ds)**2-30.*aj ## sy*(dr)**3.*aj ## syy*aj ## ry+24.*aj ## rx**4.*(ds)**2+24.*aj ## ry**4.*(ds)**2-10.*(dr)**3.*aj ## sy*aj ## ry*aj ## sxx-30.*aj ## sx*(dr)**3.*aj ## sxx*aj ## rx-20.*(dr)**3.*aj ## sy*aj ## sx*aj ## rxy-20.*(dr)**3.*aj ## ry*aj ## sxy*aj ## sx-20.*(dr)**3.*aj ## sy*aj ## sxy*aj ## rx-10.*(dr)**3.*aj ## syy*aj ## sx*aj ## rx-2.*(dr)**2.*(ds)**2.*aj ## ryy*aj ## rxx-4.*(dr)**2.*(ds)**2.*aj ## rx*aj ## rxyy-4.*(dr)**2.*(ds)**2.*aj ## ry*aj ## rxxy-4.*(dr)**2.*(ds)**2.*aj ## ry*aj ## ryyy-4.*(dr)**2.*(ds)**2.*aj ## rx*aj ## rxxx-24.*(dr)*(ds)**2.*aj ## ryy*aj ## rx**2-24.*(dr)*(ds)**2.*aj ## ry**2.*aj ## rxx+48.*aj ## ry**2.*aj ## rx**2.*(ds)**2-72.*aj ## ry**2.*aj ## ryy*(dr)*(ds)**2-72.*aj ## rx**2.*aj ## rxx*(dr)*(ds)**2+15.*aj ## rx**2.*aj ## sx**2.*(dr)**2-96.*(dr)*(ds)**2.*aj ## ry*aj ## rxy*aj ## rx+20.*(dr)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rx+15.*aj ## ry**2.*aj ## sy**2.*(dr)**2)/(dr)**4/(ds)**2 )
239 
240 #defineMacro lapSqCoeff2dOrder4Ghost30(aj,dr,ds) ( -1/12.*(-9.*aj ## rx**2.*aj ## rxx*(dr)-9.*aj ## ry**2.*aj ## ryy*(dr)+4.*aj ## ry**2.*aj ## rx**2+2.*aj ## ry**4+2.*aj ## rx**4-12.*(dr)*aj ## ry*aj ## rxy*aj ## rx-3.*(dr)*aj ## ry**2.*aj ## rxx-3.*(dr)*aj ## ryy*aj ## rx**2)/(dr)**4 )
241 
242 #defineMacro lapSqCoeff2dOrder4Ghost01(aj,dr,ds) ( -1/12.*(-48.*(dr)**2.*(ds)**2.*aj ## sxx**2+80.*(ds)**2.*aj ## ry**2.*aj ## sx**2+80.*(ds)**2.*aj ## sy**2.*aj ## rx**2-64.*(dr)**2.*(ds)**2.*aj ## sxy**2-48.*(dr)**2.*(ds)**2.*aj ## syy**2-40.*(ds)**3.*aj ## ry**2.*aj ## sxx-40.*(ds)**3.*aj ## syy*aj ## rx**2-120.*aj ## rx**2.*(ds)**3.*aj ## sxx-120.*aj ## ry**2.*(ds)**3.*aj ## syy+8.*aj ## sxxxx*(ds)**3.*(dr)**2+16.*aj ## sxxyy*(ds)**3.*(dr)**2+8.*aj ## syyyy*(ds)**3.*(dr)**2+240.*aj ## rx**2.*aj ## sx**2.*(ds)**2+156.*aj ## sy**2.*aj ## sx**2.*(dr)**2-156.*(dr)**2.*(ds)*aj ## sy*aj ## sxy*aj ## sx-117.*aj ## sy**2.*aj ## syy*(ds)*(dr)**2-117.*aj ## sx**2.*aj ## sxx*(ds)*(dr)**2+240.*aj ## ry**2.*aj ## sy**2.*(ds)**2+78.*aj ## sy**4.*(dr)**2+320.*(ds)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rx-64.*(dr)**2.*(ds)**2.*aj ## syyy*aj ## sy+78.*aj ## sx**4.*(dr)**2-39.*(dr)**2.*(ds)*aj ## sy**2.*aj ## sxx-64.*(dr)**2.*(ds)**2.*aj ## sx*aj ## sxyy-32.*(dr)**2.*(ds)**2.*aj ## syy*aj ## sxx-64.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxxy-240.*aj ## ry*(ds)**3.*aj ## sy*aj ## ryy-64.*(dr)**2.*(ds)**2.*aj ## sxxx*aj ## sx-39.*(dr)**2.*(ds)*aj ## syy*aj ## sx**2-160.*(ds)**3.*aj ## ry*aj ## sx*aj ## rxy-80.*(ds)**3.*aj ## ryy*aj ## sx*aj ## rx-80.*(ds)**3.*aj ## sy*aj ## ry*aj ## rxx-160.*(ds)**3.*aj ## sy*aj ## rxy*aj ## rx-240.*aj ## rx*(ds)**3.*aj ## sx*aj ## rxx-160.*(ds)**3.*aj ## ry*aj ## sxy*aj ## rx)/(ds)**4/(dr)**2 )
243 
244 #defineMacro lapSqCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/12.*(-72.*aj ## sx**2.*aj ## sxx*(ds)*(dr)**2+48.*aj ## sy**2.*aj ## sx**2.*(dr)**2-72.*aj ## sy**2.*aj ## syy*(ds)*(dr)**2-4.*(dr)**2.*(ds)**2.*aj ## sy*aj ## sxxy-2.*(dr)**2.*(ds)**2.*aj ## syy*aj ## sxx-4.*(dr)**2.*(ds)**2.*aj ## sx*aj ## sxyy-96.*(dr)**2.*(ds)*aj ## sy*aj ## sxy*aj ## sx-24.*(dr)**2.*(ds)*aj ## sy**2.*aj ## sxx-24.*(dr)**2.*(ds)*aj ## syy*aj ## sx**2-10.*(ds)**3.*aj ## sy*aj ## ry*aj ## rxx-10.*(ds)**3.*aj ## ryy*aj ## sx*aj ## rx-20.*(ds)**3.*aj ## sy*aj ## rxy*aj ## rx+20.*(ds)**2.*aj ## sy*aj ## ry*aj ## sx*aj ## rx-4.*(dr)**2.*(ds)**2.*aj ## sxxx*aj ## sx-30.*aj ## rx*(ds)**3.*aj ## sx*aj ## rxx-20.*(ds)**3.*aj ## ry*aj ## sx*aj ## rxy-20.*(ds)**3.*aj ## ry*aj ## sxy*aj ## rx-30.*aj ## ry*(ds)**3.*aj ## sy*aj ## ryy+24.*aj ## sy**4.*(dr)**2+aj ## sxxxx*(ds)**3.*(dr)**2+2.*aj ## sxxyy*(ds)**3.*(dr)**2+aj ## syyyy*(ds)**3.*(dr)**2-3.*(dr)**2.*(ds)**2.*aj ## syy**2-5.*(ds)**3.*aj ## syy*aj ## rx**2-15.*aj ## rx**2.*(ds)**3.*aj ## sxx-3.*(dr)**2.*(ds)**2.*aj ## sxx**2+5.*(ds)**2.*aj ## ry**2.*aj ## sx**2+5.*(ds)**2.*aj ## sy**2.*aj ## rx**2-4.*(dr)**2.*(ds)**2.*aj ## sxy**2-15.*aj ## ry**2.*(ds)**3.*aj ## syy-5.*(ds)**3.*aj ## ry**2.*aj ## sxx+24.*aj ## sx**4.*(dr)**2-4.*(dr)**2.*(ds)**2.*aj ## syyy*aj ## sy+15.*aj ## rx**2.*aj ## sx**2.*(ds)**2+15.*aj ## ry**2.*aj ## sy**2.*(ds)**2)/(ds)**4/(dr)**2 )
245 
246 #defineMacro lapSqCoeff2dOrder4Ghost03(aj,dr,ds) ( -1/12.*(-12.*(ds)*aj ## sy*aj ## sxy*aj ## sx-3.*(ds)*aj ## sy**2.*aj ## sxx-3.*(ds)*aj ## syy*aj ## sx**2+2.*aj ## sy**4+2.*aj ## sx**4-9.*aj ## sy**2.*aj ## syy*(ds)-9.*aj ## sx**2.*aj ## sxx*(ds)+4.*aj ## sy**2.*aj ## sx**2)/(ds)**4 )
247 
248 #defineMacro xxxCoeff2dOrder4Ghost10(aj,dr,ds) ( -1/24.*(-39.*aj ## rx**3.*(ds)**2-120.*aj ## rx*aj ## sx**2.*(dr)**2-96.*aj ## rx*aj ## rxx*(dr)*(ds)**2+16.*aj ## rxxx*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
249 
250 #defineMacro xxxCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/24.*(-24.*aj ## rx**3.*(ds)**2-15.*aj ## rx*aj ## sx**2.*(dr)**2-6.*aj ## rx*aj ## rxx*(dr)*(ds)**2+2.*aj ## rxxx*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
251 
252 #defineMacro xxxCoeff2dOrder4Ghost30(aj,dr,ds) ( 1/8.*aj ## rx**3/(dr)**3 )
253 
254 #defineMacro xxxCoeff2dOrder4Ghost01(aj,dr,ds) ( -1/24.*(-120.*aj ## rx**2.*aj ## sx*(ds)**2-39.*aj ## sx**3.*(dr)**2-96.*aj ## sxx*aj ## sx*(dr)**2.*(ds)+16.*aj ## sxxx*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
255 
256 #defineMacro xxxCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/24.*(-15.*aj ## rx**2.*aj ## sx*(ds)**2-24.*aj ## sx**3.*(dr)**2-6.*aj ## sxx*aj ## sx*(dr)**2.*(ds)+2.*aj ## sxxx*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
257 
258 #defineMacro xxxCoeff2dOrder4Ghost03(aj,dr,ds) ( 1/8.*aj ## sx**3/(ds)**3 )
259 
260 #defineMacro xxyCoeff2dOrder4Ghost10(aj,dr,ds) ( -1/24.*(-39.*aj ## ry*aj ## rx**2.*(ds)**2-40.*aj ## sx**2.*(dr)**2.*aj ## ry-80.*aj ## sx*(dr)**2.*aj ## sy*aj ## rx-64.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-32.*(dr)*(ds)**2.*aj ## ry*aj ## rxx+16.*aj ## rxxy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
261 
262 #defineMacro xxyCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/24.*(-24.*aj ## ry*aj ## rx**2.*(ds)**2-5.*aj ## sx**2.*(dr)**2.*aj ## ry-10.*aj ## sx*(dr)**2.*aj ## sy*aj ## rx-4.*(dr)*(ds)**2.*aj ## rxy*aj ## rx-2.*(dr)*(ds)**2.*aj ## ry*aj ## rxx+2.*aj ## rxxy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
263 
264 #defineMacro xxyCoeff2dOrder4Ghost30(aj,dr,ds) ( 1/8.*aj ## ry*aj ## rx**2/(dr)**3 )
265 
266 #defineMacro xxyCoeff2dOrder4Ghost01(aj,dr,ds) ( -1/24.*(-40.*aj ## rx**2.*(ds)**2.*aj ## sy-80.*aj ## rx*(ds)**2.*aj ## ry*aj ## sx-39.*aj ## sy*aj ## sx**2.*(dr)**2-32.*(dr)**2.*(ds)*aj ## sy*aj ## sxx-64.*(dr)**2.*(ds)*aj ## sxy*aj ## sx+16.*aj ## sxxy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
267 
268 #defineMacro xxyCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/24.*(-5.*aj ## rx**2.*(ds)**2.*aj ## sy-10.*aj ## rx*(ds)**2.*aj ## ry*aj ## sx-24.*aj ## sy*aj ## sx**2.*(dr)**2-2.*(dr)**2.*(ds)*aj ## sy*aj ## sxx-4.*(dr)**2.*(ds)*aj ## sxy*aj ## sx+2.*aj ## sxxy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
269 
270 #defineMacro xxyCoeff2dOrder4Ghost03(aj,dr,ds) ( 1/8.*aj ## sy*aj ## sx**2/(ds)**3 )
271 
272 #defineMacro xyyCoeff2dOrder4Ghost10(aj,dr,ds) ( -1/24.*(-39.*aj ## ry**2.*aj ## rx*(ds)**2-40.*aj ## sy**2.*(dr)**2.*aj ## rx-80.*aj ## sy*(dr)**2.*aj ## ry*aj ## sx-32.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-64.*(dr)*(ds)**2.*aj ## ry*aj ## rxy+16.*aj ## rxyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
273 
274 #defineMacro xyyCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/24.*(-24.*aj ## ry**2.*aj ## rx*(ds)**2-5.*aj ## sy**2.*(dr)**2.*aj ## rx-10.*aj ## sy*(dr)**2.*aj ## ry*aj ## sx-2.*(dr)*(ds)**2.*aj ## ryy*aj ## rx-4.*(dr)*(ds)**2.*aj ## ry*aj ## rxy+2.*aj ## rxyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
275 
276 #defineMacro xyyCoeff2dOrder4Ghost30(aj,dr,ds) ( 1/8.*aj ## ry**2.*aj ## rx/(dr)**3 )
277 
278 #defineMacro xyyCoeff2dOrder4Ghost01(aj,dr,ds) ( -1/24.*(-40.*aj ## ry**2.*(ds)**2.*aj ## sx-80.*aj ## ry*(ds)**2.*aj ## sy*aj ## rx-39.*aj ## sy**2.*aj ## sx*(dr)**2-32.*(dr)**2.*(ds)*aj ## syy*aj ## sx-64.*(dr)**2.*(ds)*aj ## sy*aj ## sxy+16.*aj ## sxyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
279 
280 #defineMacro xyyCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/24.*(-5.*aj ## ry**2.*(ds)**2.*aj ## sx-10.*aj ## ry*(ds)**2.*aj ## sy*aj ## rx-24.*aj ## sy**2.*aj ## sx*(dr)**2-2.*(dr)**2.*(ds)*aj ## syy*aj ## sx-4.*(dr)**2.*(ds)*aj ## sy*aj ## sxy+2.*aj ## sxyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
281 
282 #defineMacro xyyCoeff2dOrder4Ghost03(aj,dr,ds) ( 1/8.*aj ## sy**2.*aj ## sx/(ds)**3 )
283 
284 #defineMacro yyyCoeff2dOrder4Ghost10(aj,dr,ds) ( -1/24.*(-39.*aj ## ry**3.*(ds)**2-120.*aj ## ry*aj ## sy**2.*(dr)**2-96.*aj ## ry*aj ## ryy*(dr)*(ds)**2+16.*aj ## ryyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
285 
286 #defineMacro yyyCoeff2dOrder4Ghost20(aj,dr,ds) ( 1/24.*(-24.*aj ## ry**3.*(ds)**2-15.*aj ## ry*aj ## sy**2.*(dr)**2-6.*aj ## ry*aj ## ryy*(dr)*(ds)**2+2.*aj ## ryyy*(dr)**2.*(ds)**2)/(dr)**3/(ds)**2 )
287 
288 #defineMacro yyyCoeff2dOrder4Ghost30(aj,dr,ds) ( 1/8.*aj ## ry**3/(dr)**3 )
289 
290 #defineMacro yyyCoeff2dOrder4Ghost01(aj,dr,ds) ( -1/24.*(-120.*aj ## ry**2.*aj ## sy*(ds)**2-39.*aj ## sy**3.*(dr)**2-96.*aj ## syy*aj ## sy*(dr)**2.*(ds)+16.*aj ## syyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
291 
292 #defineMacro yyyCoeff2dOrder4Ghost02(aj,dr,ds) ( 1/24.*(-15.*aj ## ry**2.*aj ## sy*(ds)**2-24.*aj ## sy**3.*(dr)**2-6.*aj ## syy*aj ## sy*(dr)**2.*(ds)+2.*aj ## syyy*(dr)**2.*(ds)**2)/(dr)**2/(ds)**3 )
293 
294 #defineMacro yyyCoeff2dOrder4Ghost03(aj,dr,ds) ( 1/8.*aj ## sy**3/(ds)**3 )
295 
296 
297 ! *** orderOfAccuracy = 6
298 
299 #defineMacro xCoeff2dOrder6Ghost10(aj,dr,ds) ( -3/4.*aj ## rx/(dr) )
300 
301 #defineMacro xCoeff2dOrder6Ghost20(aj,dr,ds) ( 3/20.*aj ## rx/(dr) )
302 
303 #defineMacro xCoeff2dOrder6Ghost30(aj,dr,ds) ( -1/60.*aj ## rx/(dr) )
304 
305 #defineMacro xCoeff2dOrder6Ghost01(aj,dr,ds) ( -3/4.*aj ## sx/(ds) )
306 
307 #defineMacro xCoeff2dOrder6Ghost02(aj,dr,ds) ( 3/20.*aj ## sx/(ds) )
308 
309 #defineMacro xCoeff2dOrder6Ghost03(aj,dr,ds) ( -1/60.*aj ## sx/(ds) )
310 
311 #defineMacro yCoeff2dOrder6Ghost10(aj,dr,ds) ( -3/4.*aj ## ry/(dr) )
312 
313 #defineMacro yCoeff2dOrder6Ghost20(aj,dr,ds) ( 3/20.*aj ## ry/(dr) )
314 
315 #defineMacro yCoeff2dOrder6Ghost30(aj,dr,ds) ( -1/60.*aj ## ry/(dr) )
316 
317 #defineMacro yCoeff2dOrder6Ghost01(aj,dr,ds) ( -3/4.*aj ## sy/(ds) )
318 
319 #defineMacro yCoeff2dOrder6Ghost02(aj,dr,ds) ( 3/20.*aj ## sy/(ds) )
320 
321 #defineMacro yCoeff2dOrder6Ghost03(aj,dr,ds) ( -1/60.*aj ## sy/(ds) )
322 
323 #defineMacro lapCoeff2dOrder6Ghost10(aj,dr,ds) ( 3/4.*(2.*aj ## rx**2-aj ## rxx*(dr)+2.*aj ## ry**2-aj ## ryy*(dr))/(dr)**2 )
324 
325 #defineMacro lapCoeff2dOrder6Ghost20(aj,dr,ds) ( 3/20.*(-aj ## rx**2+aj ## rxx*(dr)-aj ## ry**2+aj ## ryy*(dr))/(dr)**2 )
326 
327 #defineMacro lapCoeff2dOrder6Ghost30(aj,dr,ds) ( -1/180.*(-2.*aj ## rx**2+3.*aj ## rxx*(dr)-2.*aj ## ry**2+3.*aj ## ryy*(dr))/(dr)**2 )
328 
329 #defineMacro lapCoeff2dOrder6Ghost01(aj,dr,ds) ( -3/4.*(-2.*aj ## sx**2+aj ## sxx*(ds)-2.*aj ## sy**2+aj ## syy*(ds))/(ds)**2 )
330 
331 #defineMacro lapCoeff2dOrder6Ghost02(aj,dr,ds) ( 3/20.*(-aj ## sx**2+aj ## sxx*(ds)-aj ## sy**2+aj ## syy*(ds))/(ds)**2 )
332 
333 #defineMacro lapCoeff2dOrder6Ghost03(aj,dr,ds) ( 1/180.*(2.*aj ## sx**2-3.*aj ## sxx*(ds)+2.*aj ## sy**2-3.*aj ## syy*(ds))/(ds)**2 )
334 
Generated on Fri Jan 4 2013 10:17:50 for Overture by   doxygen 1.8.3