| /phasta/shapeFunction/src/ |
| 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 96 t4 = t1*s; in Fn() 97 f = r*t2-10.0/7.0*r*t4+60.0/91.0*r*t1-10.0/91.0*r*s+5.0/1001.0*r-t2/7+ in Fn() 98 20.0/91.0*t4-10.0/91.0*t1+20.0/1001.0*s-1.0/1001.0; in Fn() 117 t4 = r*t1; in Fn() 118 f = t2*s-t2/7-10.0/7.0*t4*s+20.0/91.0*t4+60.0/91.0*t1*s-10.0/91.0*t1 in Fn() 153 t4 = t3*s; in Fn() 154 f = t2*t4-3.0/4.0*t2*t3+3.0/20.0*t2*s-t2/140-3.0/4.0*t1*t4+3.0/5.0*t1* in Fn() 155 t3-9.0/70.0*t1*s+3.0/455.0*t1+3.0/20.0*r*t4-9.0/70.0*r*t3+27.0/910.0*r*s-3.0/ in Fn() 156 1820.0*r-t4/140+3.0/455.0*t3-3.0/1820.0*s+1.0/10010.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 113 t4 = r*t1; 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() 124 t4 = t1*s; in FnDrv() 127 drv[0] = 3.0*t1*t2-9.0/7.0*t4+9.0/91.0*t1-12.0/7.0*r*t2+72.0/91.0*t6-6.0/ 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() 156 t4 = t1*s; in FnDrv() 157 drv[0] = t2*s-15.0/8.0*t2+5.0/4.0*t4-5.0/14.0*t1+15.0/364.0*s-1.0/728.0; in FnDrv() 158 drv[1] = 5.0*r*t2-15.0/2.0*r*t4+15.0/4.0*r*t1-5.0/7.0*r*s+15.0/364.0*r in FnDrv() [all …]
|
| 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 41 t4 = t2*t2; in Bn() 42 f = t1+6.0/5.0*t*t2+3.0/10.0*t4; in Bn() 50 t4 = t2*t2; in Bn() 51 f = t1+6.0/7.0*t2*s+t4/7; in Bn() 78 t4 = t2*t2; in Bn() 79 f = (t1+2.0/3.0*s*t2+t4/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); in Bn() 117 t4 = s+r-1.0; in Bn() 118 t6 = t4*t4; in Bn() 120 f = t2+20.0/9.0*t1*t*t4+5.0/3.0*t1*t6+10.0/21.0*t*t6*t4+5.0/126.0*t10; in Bn() [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 80 t4 = t2*t2; in BnDrv() 81 t5 = 12.0/7.0*t1+12.0/7.0*t3+12.0/35.0*t4; in BnDrv() 84 drv[2] = 3.0*t1+24.0/7.0*t3+6.0/7.0*t4; in BnDrv() 90 t4 = t2*t2; 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() 99 t4 = s*s; 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() [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 45 t4 = r*r; in EnDrv() 46 t5 = t4*s; in EnDrv() 47 t6 = t4*r; in EnDrv() 54 t4 = r*t1*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() 66 t4 = r*t2; 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() [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 36 t4 = r*r; in En() 37 f = t1*s-6.0*r*t1+6.0*t4*s-t4*r; in En() 75 t4 = t1*s; 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() 81 *t13*r*t4+336.0*t13*t7*t1-36.0*t13*t10*s+t21; in En() 112 t4 = t3*t3; 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() 138 t4 = t2*t2; in En() 145 f = t4*t3-91.0*r*t2*t4+2366.0*t8*t4*t9-26026.0*t12*t4*t1+143143.0*t15* in En() [all …]
|
| H A D | Lagrange.c | 24 double r=par[0],s=par[1],t=par[2],t8,t4; in TetQuadLag() local 37 t4 = -3.0+4.0*r+4.0*s+4.0*t; in TetQuadLag() 47 dN[3][0] = t4; in TetQuadLag() 48 dN[3][1] = t4; in TetQuadLag() 49 dN[3][2] = t4; in TetQuadLag()
|
| /phasta/phSolver/incompressible/ |
| H A D | stats.f | 217 real*8 u1, u2, u3, r0, r1, r2, r3, r4, r5, t3, t4, t5 local 285 t4 = CInv(4) * r3 + CInv(2) * r4 + CInv(5) * r5 291 & + r1 - DInv(2) * (B(2) * t4 + B(1) * t3) 293 & + r2 - DInv(3) * (B(3) * t5 + B(2) * t4) 296 stsVelSqr(i,5) = stsVelSqr(i,5) + t4 315 t4 = CInv(4) * r3 + CInv(2) * r4 + CInv(5) * r5 321 & + r1 - DInv(2) * (B(2) * t4 + B(1) * t3) 323 & + r2 - DInv(3) * (B(3) * t5 + B(2) * t4) 325 stsStress(i,5) = stsStress(i,5) + t4
|
| H A D | e3source.f | 180 t4 = t3*x 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 192 & *t4-14*t3+t8-2*t2)
|
| /phasta/phSolver/common/ |
| H A D | blowerControl.f90 | 101 real*8 :: t0, t1, t2, t3, t4, vmax, vmin local 130 t4 = blower(iBlower)%t_cycle 134 t = mod(BC_t - blower(iBlower)%t0, t4) 142 elseif(t3 <= t .and. t < t4) then !low
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 147 c$$$ t4 = (-1.D0+x)**2 150 c$$$ fya = 4.D0*om3*t1*t4*t7*y*(1.D0-3.D0*y+2.D0*t10)-2.D0*om1*w+om3*(o 180 t4 = t3*x 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 192 & *t4-14*t3+t8-2*t2)
|
| H A D | e3source.fRHSonly | 147 c$$$ t4 = (-1.D0+x)**2 150 c$$$ fya = 4.D0*om3*t1*t4*t7*y*(1.D0-3.D0*y+2.D0*t10)-2.D0*om1*w+om3*(o 180 t4 = t3*x 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 192 & *t4-14*t3+t8-2*t2)
|