| /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 233 t8 = s+r-1.0; in Bn() 234 t10 = t8*t8; in Bn() 235 f = (t1+6.0/13.0*t2*s+t4/26)*(t6*t+12.0/7.0*t6*t8+6.0/7.0*t*t10+4.0/ in Bn() 236 35.0*t10*t8); in Bn() 278 t8 = t6*t6; in Bn() 280 *t8+4.0/35.0*t8*t6); in Bn() 286 t8 = t*t; in Bn() 289 f = (r-2.0/15.0)*(t2+6.0/11.0*s*t3+3.0/55.0*t5)*(t8+6.0/5.0*t*t9+3.0/ in Bn() 320 t8 = t6*t6; in Bn() [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 113 t8 = r*r; in En() 114 t11 = t8*r; in En() 115 t13 = t8*t8; in En() 117 f = t2*t4-66.0*r*t4*t1+1210.0*t8*t4*s-9075.0*t11*t4+32670.0*t13*t3*t2 in En() 118 -60984.0*t13*r*t3*t1+60984.0*t13*t8*t3*s-32670.0*t13*t11*t3+9075.0*t24*t2 in En() 119 -1210.0*t24*r*t1+66.0*t24*t8*s-t24*t11; in En() 126 t8 = r*r; in En() 127 t11 = t8*r; in En() 128 t14 = t8*t8; in En() [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 173 t8 = 20.0/9.0*t2+10.0/3.0*t4+10.0/7.0*t6+10.0/63.0*t7; in BnDrv() 174 drv[0] = t8; in BnDrv() 175 drv[1] = t8; in BnDrv() 186 t8 = s+2.0/11.0*r-2.0/11.0; in BnDrv() 188 t11 = t8*(12.0/7.0*t1+12.0/7.0*t9+12.0/35.0*t5); in BnDrv() 191 drv[2] = t8*(3.0*t1+24.0/7.0*t9+6.0/7.0*t5); in BnDrv() 198 t8 = s*s; in BnDrv() 201 t12 = t8+6.0/11.0*s*t9+3.0/55.0*t11; in BnDrv() 213 t8 = t1*s; in BnDrv() [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 126 t8 = r*t1; in FnDrv() 129 drv[1] = 2.0*t8*s-3.0/7.0*t8-12.0/7.0*t4+36.0/91.0*t1+36.0/91.0*t6-9.0/ in FnDrv() 181 t8 = r*t2; in FnDrv() 184 t8-9.0/35.0*t9+6.0/455.0*r+3.0/20.0*t3-9.0/70.0*t2+27.0/910.0*s-3.0/1820.0; in FnDrv() 187 70.0*t1+9.0/20.0*t8-9.0/35.0*t9+27.0/910.0*r-3.0/140.0*t2+6.0/455.0*s-3.0/ in FnDrv() 247 t8 = r*t1; in FnDrv() 249 t10 = 2.0*r*t3-10.0/3.0*t5+100.0/51.0*t7-25.0/51.0*t8+5.0/102.0*t9-r/714- in FnDrv() 253 204.0*t11-5.0/3.0*t5+40.0/17.0*t7-75.0/68.0*t8+10.0/51.0*t9-5.0/476.0*r+5.0/ in FnDrv() 265 t8 = t1*s; in FnDrv() [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 57 t8 = t5*r*s; in EnDrv() 59 drv[0] = -15.0*t2+100.0*t4-150.0*t6+60.0*t8-5.0*t9; in EnDrv() 60 drv[1] = 5.0*t2-60.0*t4+150.0*t6-100.0*t8+15.0*t9; in EnDrv() 165 t8 = t7*t3; in EnDrv() 178 t28 = -66.0*t4+2420.0*t6-27225.0*t8+130680.0*t12-304920.0*t15+365904.0* in EnDrv() 180 t29 = 11.0*t4-660.0*t6+10890.0*t8-72600.0*t12+228690.0*t15-365904.0*t18+ in EnDrv() 192 t8 = r*r; in EnDrv() 193 t10 = t8*t4*s; in EnDrv() 194 t11 = t8*r; in EnDrv() [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 192 t8 = t1*s; in Fn() 193 f = r*t3-7.0/3.0*r*t5+35.0/17.0*r*t2-175.0/204.0*r*t8+35.0/204.0*r*t1-r in Fn() 194 *s/68+r/2652-t3/9+14.0/51.0*t5-35.0/136.0*t2+35.0/306.0*t8-5.0/204.0*t1+s/442 in Fn() 222 t8 = t1*r; in Fn() 223 f = t2*t4-2.0/3.0*t2*t3+2.0/17.0*t2*s-t2/204-10.0/9.0*t8*t4+40.0/51.0* in Fn() 224 t8*t3-5.0/34.0*t8*s+t8/153+20.0/51.0*t1*t4-5.0/17.0*t1*t3+t1*s/17-t1/357-5.0/ in Fn() 242 t8 = t1*r; in Fn() 244 -175.0/204.0*t8*s+35.0/306.0*t8+35.0/204.0*t1*s-5.0/204.0*t1-r*s/68+r/442+s in Fn() 268 t8 = t3*s; in Fn() [all …]
|
| H A D | Lagrange.c | 24 double r=par[0],s=par[1],t=par[2],t8,t4; in TetQuadLag() local 25 t8 = 1.0-r-s-t; in TetQuadLag() 29 N[3] = (1.0-2.0*r-2.0*s-2.0*t)*t8; in TetQuadLag() 33 N[7] = 4.0*r*t8; in TetQuadLag() 34 N[8] = 4.0*s*t8; in TetQuadLag() 35 N[9] = 4.0*t*t8; in TetQuadLag()
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fRHSonly | 63 c$$$ t8 = t7**2 64 c$$$ t9 = t8*t7 81 c$$$ t60 = t8*y 87 c$$$ #D0*t57*y+36.D0*t12*t60+192.D0*t50*t18+98.D0*t66*t8+24.D0*t12*t18+2 89 c$$$ t106 = -336.D0*t50*t8-12.D0*t12*t9+8.D0*t73*t9-140.D0*t6*t8+28.D0* 90 c$$$ #t73*t8+12.D0*t15*t17+12.D0*t43*t17+288.D0*t50*t60+112.D0*t77*t9-22 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* 182 t8 = 8*t7 189 t54 = -8/Re*(24.E0/5.E0*t4-12*t3+t8+2*(4*t7-6*t2+2*x)*(12 [all …]
|
| H A D | e3source.fimplicit | 63 c$$$ t8 = t7**2 64 c$$$ t9 = t8*t7 81 c$$$ t60 = t8*y 87 c$$$ #D0*t57*y+36.D0*t12*t60+192.D0*t50*t18+98.D0*t66*t8+24.D0*t12*t18+2 89 c$$$ t106 = -336.D0*t50*t8-12.D0*t12*t9+8.D0*t73*t9-140.D0*t6*t8+28.D0* 90 c$$$ #t73*t8+12.D0*t15*t17+12.D0*t43*t17+288.D0*t50*t60+112.D0*t77*t9-22 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* 182 t8 = 8*t7 189 t54 = -8/Re*(24.E0/5.E0*t4-12*t3+t8+2*(4*t7-6*t2+2*x)*(12 [all …]
|
| /phasta/phSolver/incompressible/ |
| H A D | e3source.f | 182 t8 = 8*t7 189 t54 = -8/Re*(24.E0/5.E0*t4-12*t3+t8+2*(4*t7-6*t2+2*x)*(12 192 & *t4-14*t3+t8-2*t2)
|