| /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 103 t3 = t2*s; in Fn() 104 f = t1*t3-6.0/7.0*t1*t2+18.0/91.0*t1*s-t1/91-3.0/7.0*r*t3+36.0/91.0*r* in Fn() 105 t2-9.0/91.0*r*s+6.0/1001.0*r+3.0/91.0*t3-3.0/91.0*t2+9.0/1001.0*s-3.0/5005.0; in Fn() 110 t3 = s*s; in Fn() 111 f = t2*t3-3.0/7.0*t2*s+3.0/91.0*t2-6.0/7.0*t1*t3+36.0/91.0*t1*s-3.0/ in Fn() 112 91.0*t1+18.0/91.0*r*t3-9.0/91.0*r*s+9.0/1001.0*r-t3/91+6.0/1001.0*s-3.0/5005.0; in Fn() 136 t3 = t2*s; in Fn() 138 f = r*t3-15.0/8.0*r*t2+5.0/4.0*r*t6-5.0/14.0*r*t1+15.0/364.0*r*s-r/728- in Fn() 139 t3/8+t2/4-5.0/28.0*t6+5.0/91.0*t1-5.0/728.0*s+1.0/4004.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 60 t3 = r*r; in FnDrv() 62 drv[1] = t3-3.0/5.0*r+1.0/15.0; in FnDrv() 81 t3 = r*s; in FnDrv() 83 drv[0] = 2.0*r*t1-t3+r/11-t1/2+3.0/11.0*s-3.0/110.0; 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() 105 t3 = t1*s; in FnDrv() 106 drv[0] = t2-10.0/7.0*t3+60.0/91.0*t1-10.0/91.0*s+5.0/1001.0; in FnDrv() 107 drv[1] = 4.0*r*t3-30.0/7.0*r*t1+120.0/91.0*r*s-10.0/91.0*r-4.0/7.0*t3+ in FnDrv() 178 t3 = t2*s; in FnDrv() [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 31 t3 = r*r; in En() 32 f = t1-3.0*r*s+t3; in En() 64 t3 = t1*t1; 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() 74 t3 = t2*t2; 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() 86 t3 = t2*t2; in En() 92 f = t3*s-45.0*t3*r+540.0*t6*t2*t7-2520.0*t10*t2*t1+5292.0*t13*t2*s in En() 98 t3 = t2*t2; in En() [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 65 t3 = s+r-1.0; in Bn() 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() 71 t3 = s+r-1.0; 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() 83 t3 = r-1.0; 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() [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 37 t3 = r*r; in EnDrv() 38 drv[0] = -6.0*t1+12.0*t2-3.0*t3; in EnDrv() 39 drv[1] = 3.0*t1-12.0*t2+6.0*t3; in EnDrv() 44 t3 = r*t1; 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() 65 t3 = t2*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() [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 79 t3 = t*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() 89 t3 = t*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() 110 t3 = s*t2; in BnDrv() 112 drv[0] = 4.0/3.0*t1+t3+t4/7; in BnDrv() 113 drv[1] = 3.0*t1+8.0/3.0*t3+t4/2; in BnDrv() [all …]
|
| /phasta/phSolver/incompressible/ |
| H A D | e3lhs.f | 45 dimension t1(npro,3), t2(npro,3), t3(npro,3), 124 t3(:,1) = t2(:,1) + tmp1 * u1 + tmp2 * rLui(:,1) 125 t3(:,2) = t2(:,2) + tmp1 * u2 + tmp2 * rLui(:,2) 126 t3(:,3) = t2(:,3) + tmp1 * u3 + tmp2 * rLui(:,3) 137 tmp = t3(:,1) * shg(:,aa,1) 138 & + t3(:,2) * shg(:,aa,2) 139 & + t3(:,3) * shg(:,aa,3) 170 tmp = t3(:,1) * shg(:,aa,1) 171 & + t3(:,2) * shg(:,aa,2) 172 & + t3(:,3) * shg(:,aa,3)
|
| H A D | e3source.f | 179 t3 = t2**2 180 t4 = t3*x 186 t26 = t3**2 187 t30 = t3*t2 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 | stats.f | 217 real*8 u1, u2, u3, r0, r1, r2, r3, r4, r5, t3, t4, t5 local 284 t3 = CInv(1) * r3 + CInv(4) * r4 + CInv(6) * r5 289 & + r0 - DInv(1) * (B(1) * t3 + B(3) * t5) 291 & + r1 - DInv(2) * (B(2) * t4 + B(1) * t3) 295 stsVelSqr(i,4) = stsVelSqr(i,4) + t3 314 t3 = CInv(1) * r3 + CInv(4) * r4 + CInv(6) * r5 319 & + r0 - DInv(1) * (B(1) * t3 + B(3) * t5) 321 & + r1 - DInv(2) * (B(2) * t4 + B(1) * t3) 324 stsStress(i,4) = stsStress(i,4) + t3
|
| /phasta/phSolver/common/ |
| H A D | blowerControl.f90 | 101 real*8 :: t0, t1, t2, t3, t4, vmax, vmin local 129 t3 = t2 + blower(iBlower)%t_fallTime 140 elseif(t2 <= t .and. t < t3) then !falling edge 141 vmagref = (vmax - vmin)*(t - t3)/(t2 - t3) + vmin 142 elseif(t3 <= t .and. t < t4) then !low
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 100 c$$$ t3 = nu*t2 116 c$$$ t50 = -84.D0*t3*t6+343.D0*t11*t5*t2+66.D0*t11*t5*x-98.D0*t11*t5*t1 118 c$$$ #*t3*t5-56.D0*t39*t40-28.D0*t43*t40+8.D0*t43*t47 131 c$$$ #2*t92+427.D0*t3*t89+8.D0*t39*t46+4.D0*t43*t46 133 c$$$ #8.D0*t52*t40-24.D0*t52*t46-686.D0*t3*t92+2.D0*t86*t5-40.D0*t39*t10 138 c$$$ t3 = om3*x 143 c$$$ fxa = 2.D0*om2*w+2.D0*t3*t5*(2.D0+x-10.D0*t7+7.D0*t7*x)*t12*t15+om 144 c$$$ #2*(om1*y-om2*x)-om3*(t3-om1*z) 154 c$$$ t3 = dexp(7.D0*x) 159 c$$$ fza = -2.D0*om1*x*t3*(2.D0+x-10.D0*t5+7.D0*t5*x)*t10*t13-4.D0*om2* [all …]
|
| H A D | e3source.fRHSonly | 100 c$$$ t3 = nu*t2 116 c$$$ t50 = -84.D0*t3*t6+343.D0*t11*t5*t2+66.D0*t11*t5*x-98.D0*t11*t5*t1 118 c$$$ #*t3*t5-56.D0*t39*t40-28.D0*t43*t40+8.D0*t43*t47 131 c$$$ #2*t92+427.D0*t3*t89+8.D0*t39*t46+4.D0*t43*t46 133 c$$$ #8.D0*t52*t40-24.D0*t52*t46-686.D0*t3*t92+2.D0*t86*t5-40.D0*t39*t10 138 c$$$ t3 = om3*x 143 c$$$ fxa = 2.D0*om2*w+2.D0*t3*t5*(2.D0+x-10.D0*t7+7.D0*t7*x)*t12*t15+om 144 c$$$ #2*(om1*y-om2*x)-om3*(t3-om1*z) 154 c$$$ t3 = dexp(7.D0*x) 159 c$$$ fza = -2.D0*om1*x*t3*(2.D0+x-10.D0*t5+7.D0*t5*x)*t10*t13-4.D0*om2* [all …]
|