| /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 66 t5 = t3*t3; in Bn() 67 f = t1*t+12.0/7.0*t3*t1+6.0/7.0*t*t5+4.0/35.0*t3*t5; in Bn() 72 t5 = t3*t3; in Bn() 73 f = (s+2.0/9.0*r-2.0/9.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5); in Bn() 84 t5 = t3*t3; in Bn() 85 f = t1*s+4.0/3.0*t1*t3+s*t5/2+t5*t3/21; in Bn() 90 t5 = t3*t3; in Bn() 91 f = (r-2.0/11.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5); in Bn() 99 t5 = t3*t3; in Bn() [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 145 t5 = t2*s; in Fn() 146 f = t1*t3-5.0/4.0*t1*t5+t1*t2/2-t1*s/14+t1/364-3.0/8.0*r*t3+r*t5/2-3.0/ in Fn() 147 14.0*r*t2+3.0/91.0*r*s-r/728+t3/40-t5/28+3.0/182.0*t2-s/364+1.0/8008.0; in Fn() 191 t5 = t2*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() 241 t5 = r*t2; in Fn() 243 f = t3*s-t3/9-7.0/3.0*t5*s+14.0/51.0*t5+35.0/17.0*t2*s-35.0/136.0*t2 in Fn() 292 t5 = t4*s; in Fn() [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 81 t5 = 12.0/7.0*t1+12.0/7.0*t3+12.0/35.0*t4; in BnDrv() 82 drv[0] = t5; in BnDrv() 83 drv[1] = t5; in BnDrv() 91 t5 = s+2.0/9.0*r-2.0/9.0; in BnDrv() 92 t7 = t5*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0); in BnDrv() 95 drv[2] = t5*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0); in BnDrv() 100 t5 = r-1.0; in BnDrv() 101 t6 = s*t5; in BnDrv() 102 t7 = t5*t5; 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 82 t5 = r*r; in FnDrv() 84 drv[1] = 2.0*t5*s-t5/2-t3+3.0/11.0*r+s/11-3.0/110.0; in FnDrv() 114 t5 = r*s; in FnDrv() 116 drv[0] = 2.0*r*t2-12.0/7.0*t4+36.0/91.0*t5-2.0/91.0*r-3.0/7.0*t2+36.0/ 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() 165 t5 = r*t4; in FnDrv() 169 t13 = 4.0*t9*t4-15.0/4.0*t9*t1+t9*s-t9/14-3.0/2.0*t5+3.0/2.0*t6-3.0/7.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() 179 t5 = t1*t2; 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 46 t5 = t4*s; in EnDrv() 48 drv[0] = -10.0*t2+40.0*t3-30.0*t5+4.0*t6; in EnDrv() 49 drv[1] = 4.0*t2-30.0*t3+40.0*t5-10.0*t6; in EnDrv() 55 t5 = r*r; in EnDrv() 56 t6 = t5*t1; in EnDrv() 57 t8 = t5*r*s; in EnDrv() 58 t9 = t5*t5; in EnDrv() 67 t5 = r*r; in EnDrv() 68 t7 = t5*t1*s; in EnDrv() [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 42 t5 = r*r; in En() 43 t9 = t5*t5; in En() 44 f = t2-10.0*r*t1*s+20.0*t5*t1-10.0*t5*r*s+t9; in En() 49 t5 = r*r; in En() 50 t10 = t5*t5; in En() 51 f = t2*s-15.0*t2*r+50.0*t5*t1*s-50.0*t5*r*t1+15.0*t10*s-t10*r; in En() 125 t5 = t1*s; in En() 130 f = t3*t2-78.0*r*t3*t5+1716.0*t8*t3*t1-15730.0*t11*t3*s+70785.0*t14*t3 in En() 131 -169884.0*t14*r*t2*t5+226512.0*t14*t8*t2*t1-169884.0*t14*t11*t2*s+70785.0*t25* in En() [all …]
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fRHSonly | 60 c$$$ t5 = dexp(14.D0*x) 61 c$$$ t6 = t2*x*t5 65 c$$$ t12 = t2*t5 79 c$$$ t50 = t2*t22*t5 83 c$$$ t66 = t65*t5 84 c$$$ t73 = t22*t5 85 c$$$ t77 = t2*t1*t5 101 c$$$ t5 = dexp(7.D0*x) 102 c$$$ t6 = t5*y 116 c$$$ t50 = -84.D0*t3*t6+343.D0*t11*t5*t2+66.D0*t11*t5*x-98.D0*t11*t5*t1 [all …]
|
| H A D | e3source.fimplicit | 60 c$$$ t5 = dexp(14.D0*x) 61 c$$$ t6 = t2*x*t5 65 c$$$ t12 = t2*t5 79 c$$$ t50 = t2*t22*t5 83 c$$$ t66 = t65*t5 84 c$$$ t73 = t22*t5 85 c$$$ t77 = t2*t1*t5 101 c$$$ t5 = dexp(7.D0*x) 102 c$$$ t6 = t5*y 116 c$$$ t50 = -84.D0*t3*t6+343.D0*t11*t5*t2+66.D0*t11*t5*x-98.D0*t11*t5*t1 [all …]
|
| /phasta/phSolver/incompressible/ |
| H A D | stats.f | 217 real*8 u1, u2, u3, r0, r1, r2, r3, r4, r5, t3, t4, t5 local 286 t5 = CInv(6) * r3 + CInv(5) * r4 + CInv(3) * r5 289 & + r0 - DInv(1) * (B(1) * t3 + B(3) * t5) 293 & + r2 - DInv(3) * (B(3) * t5 + B(2) * t4) 297 stsVelSqr(i,6) = stsVelSqr(i,6) + t5 316 t5 = CInv(6) * r3 + CInv(5) * r4 + CInv(3) * r5 319 & + r0 - DInv(1) * (B(1) * t3 + B(3) * t5) 323 & + r2 - DInv(3) * (B(3) * t5 + B(2) * t4) 326 stsStress(i,6) = stsStress(i,6) + t5
|