| /phasta/shapeFunction/src/ |
| H A D | Bn.c | 21 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18; in Bn() local 134 t7 = s+r-1.0; in Bn() 135 t9 = t7*t7; in Bn() 136 f = (t1+6.0/11.0*s*t2+3.0/55.0*t4)*(t6+6.0/5.0*t*t7+3.0/10.0*t9); in Bn() 162 t7 = t5*t5; in Bn() 163 f = (r-2.0/13.0)*(s+2.0/9.0*r-2.0/9.0)*(t4+6.0/5.0*t*t5+3.0/10.0*t7); in Bn() 214 t7 = t4*t4; in Bn() 215 t11 = t7*t7; in Bn() 216 f = t2*t+30.0/11.0*t2*t4+30.0/11.0*t1*t*t7+40.0/33.0*t1*t7*t4+5.0/22.0* in Bn() 223 t7 = t5*t5; in Bn() [all …]
|
| H A D | EnDrv.c | 18 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18; in EnDrv() local 68 t7 = t5*t1*s; in EnDrv() 73 drv[0] = -21.0*t3+210.0*t4-525.0*t7+420.0*t9-105.0*t11+6.0*t12; in EnDrv() 74 drv[1] = 6.0*t3-105.0*t4+420.0*t7-525.0*t9+210.0*t11-21.0*t12; in EnDrv() 82 t7 = t6*t2; in EnDrv() 88 drv[0] = -28.0*t3+392.0*t5-1470.0*t7+1960.0*t10-980.0*t12+168.0*t14-7.0* in EnDrv() 90 drv[1] = 7.0*t3-168.0*t5+980.0*t7-1960.0*t10+1470.0*t12-392.0*t14+28.0* in EnDrv() 99 t7 = r*r; in EnDrv() 100 t9 = t7*t3*s; in EnDrv() 101 t10 = t7*r; in EnDrv() [all …]
|
| H A D | BnDrv.c | 20 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18; in BnDrv() local 92 t7 = t5*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0); in BnDrv() 93 drv[0] = 2.0/9.0*t1+4.0/15.0*t3+t4/15+t7; in BnDrv() 94 drv[1] = t1+6.0/5.0*t3+3.0/10.0*t4+t7; in BnDrv() 102 t7 = t5*t5; in BnDrv() 103 drv[0] = (2.0/3.0*s+r/6-1.0/6.0)*t2+2.0/3.0*t4+4.0/9.0*t6+t7/18; in BnDrv() 104 drv[1] = (2.0*s+2.0/3.0*r-2.0/3.0)*t2+2.0/3.0*t4+4.0/9.0*t6+t7/18; in BnDrv() 105 drv[2] = t4+2.0/3.0*t6+t7/12; in BnDrv() 121 t7 = t5*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0); in BnDrv() 122 drv[0] = t1+6.0/5.0*t*t2+3.0/10.0*t4+t7; in BnDrv() [all …]
|
| H A D | En.c | 19 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t18; in En() local 65 t7 = r*r; in En() 66 t10 = t7*r; in En() 67 t12 = t7*t7; in En() 68 f = t3*t2-28.0*r*t3*t1+196.0*t7*t3*s-490.0*t10*t3+490.0*t12*t2-196.0* in En() 69 t12*r*t1+28.0*t12*t7*s-t12*t10; in En() 76 t7 = r*r; in En() 77 t10 = t7*r; in En() 78 t13 = t7*t7; in En() 80 f = t3-36.0*r*t2*t4+336.0*t7*t2*t1-1176.0*t10*t2*s+1764.0*t13*t2-1176.0 in En() [all …]
|
| H A D | Fn.c | 20 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18; in Fn() local 202 t7 = t2*s; in Fn() 203 f = t1*t4-5.0/3.0*t1*t3+50.0/51.0*t1*t7-25.0/102.0*t1*t2+5.0/204.0*t1*s in Fn() 204 -t1/1428-r*t4/3+10.0/17.0*r*t3-25.0/68.0*r*t7+5.0/51.0*r*t2-5.0/476.0*r*s+r in Fn() 205 /3094+t4/51-5.0/136.0*t3+5.0/204.0*t7-5.0/714.0*t2+5.0/6188.0*s-1.0/37128.0; in Fn() 308 t7 = t3*t3; in Fn() 309 t12 = 1.0/352716.0-t2*t4+t1/1938-5.0/3876.0*t6+t3/1938-5.0/3876.0*t4+t7 in Fn() 310 /969-s/13566-2.0/57.0*r*t7-20.0/57.0*t6*t3+2.0/969.0*r*s+20.0/19.0*t6*t4; in Fn() 311 t24 = -5.0/323.0*r*t3+40.0/323.0*t1*t3-t6*t7+40.0/969.0*t6*s+6.0/19.0*t1* in Fn() 312 t7-20.0/57.0*t1*t4+t2*t7-5.0/323.0*t1*s+6.0/19.0*t2*t3-2.0/57.0*t2*s+40.0/969.0 in Fn() [all …]
|
| H A D | FnDrv.c | 19 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18; in FnDrv() local 115 t7 = r*r; in FnDrv() 118 drv[1] = 3.0*t7*t1-12.0/7.0*t7*s+18.0/91.0*t7-9.0/7.0*t4+72.0/91.0*t5-9.0 in FnDrv() 135 t7 = t1*t1; in FnDrv() 138 drv[1] = t7-10.0/7.0*t2+60.0/91.0*t1-10.0/91.0*r+5.0/1001.0; in FnDrv() 167 t7 = r*s; in FnDrv() 170 t7+3.0/91.0*r+t4/10-3.0/28.0*t1+3.0/91.0*s-1.0/364.0; in FnDrv() 171 drv[0] = 2.0*r*t2-5.0/2.0*t5+t6-t7/7+r/182-3.0/8.0*t2+t4/2-3.0/14.0*t1+ in FnDrv() 197 t7 = t1*s; in FnDrv() 199 t10 = 4.0*t2*t3-3.0/2.0*t5+t2/10-15.0/4.0*t1*t3+3.0/2.0*t7-3.0/28.0*t1+r* in FnDrv() [all …]
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 62 c$$$ t7 = y**2 63 c$$$ t8 = t7**2 64 c$$$ t9 = t8*t7 68 c$$$ t18 = t7*y 73 c$$$ t29 = t17*t7 76 c$$$ t46 = -40.D0*t6*t9-6.D0*t12*t7+92.D0*t15*t19+132.D0*t23*t24+80.D0* 92 c$$$ t131 = -8.D0*t57*t18-20.D0*t6*t7+56.D0*t77*t7-42.D0*t12*t8-24.D0*t 93 c$$$ #23*t17-48.D0*t50*t7+4.D0*t73*t7-16.D0*t73*t18+392.D0*t77*t8+14.D0* 94 c$$$ #t66*t7-84.D0*t66*t60+12.D0*t57*t7 140 c$$$ t7 = x**2 [all …]
|
| H A D | e3source.fRHSonly | 62 c$$$ t7 = y**2 63 c$$$ t8 = t7**2 64 c$$$ t9 = t8*t7 68 c$$$ t18 = t7*y 73 c$$$ t29 = t17*t7 76 c$$$ t46 = -40.D0*t6*t9-6.D0*t12*t7+92.D0*t15*t19+132.D0*t23*t24+80.D0* 92 c$$$ t131 = -8.D0*t57*t18-20.D0*t6*t7+56.D0*t77*t7-42.D0*t12*t8-24.D0*t 93 c$$$ #23*t17-48.D0*t50*t7+4.D0*t73*t7-16.D0*t73*t18+392.D0*t77*t8+14.D0* 94 c$$$ #t66*t7-84.D0*t66*t60+12.D0*t57*t7 140 c$$$ t7 = x**2 [all …]
|
| /phasta/phSolver/incompressible/ |
| H A D | e3source.f | 181 t7 = t2*x 182 t8 = 8*t7 189 t54 = -8/Re*(24.E0/5.E0*t4-12*t3+t8+2*(4*t7-6*t2+2*x)*(12 190 & *t13-2)+(24*x-12)*t21)-64*(t26/2-2*t3*t7+3*t30-2*t4+t3
|