| /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
63 t2 = t1*t1; in Fn()
64 f = t2-5.0/3.0*t1*s+10.0/11.0*t1-2.0/11.0*s+1.0/99.0; in Fn()
68 t2 = t1*s; in Fn()
69 f = r*t2-r*t1+3.0/11.0*r*s-r/55-t2/6+2.0/11.0*t1-3.0/55.0*s+2.0/495.0; in Fn()
73 t2 = s*s; in Fn()
74 f = t1*t2-t1*s/2+t1/22-r*t2/2+3.0/11.0*r*s-3.0/110.0*r+t2/22-3.0/110.0*s in Fn()
79 t2 = r*t1; in Fn()
80 f = t2*s-t2/6-t1*s+2.0/11.0*t1+3.0/11.0*r*s-3.0/55.0*r-s/55+2.0/495.0; in Fn()
84 t2 = t1*t1; 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
98 t2 = t1*t1; in FnDrv()
100 drv[1] = 5.0*t2-60.0/7.0*t1*s+450.0/91.0*t1-100.0/91.0*s+75.0/1001.0; in FnDrv()
104 t2 = t1*t1; 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()
112 t2 = t1*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()
123 t2 = s*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()
128 91.0*r+18.0/91.0*t2-9.0/91.0*s+9.0/1001.0; 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
41 t2 = t1*t1; in En()
44 f = t2-10.0*r*t1*s+20.0*t5*t1-10.0*t5*r*s+t9; in En()
48 t2 = t1*t1; 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()
55 t2 = t1*t1; in En()
58 f = t2*t1-21.0*r*t2*s+105.0*t6*t2-175.0*t6*r*t1*s+105.0*t11*t1-21.0*t11 in En()
63 t2 = t1*s; 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()
73 t2 = t1*t1; 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
40 t2 = s+r-1.0; in Bn()
41 t4 = t2*t2; in Bn()
42 f = t1+6.0/5.0*t*t2+3.0/10.0*t4; in Bn()
49 t2 = r-1.0; in Bn()
50 t4 = t2*t2; in Bn()
51 f = t1+6.0/7.0*t2*s+t4/7; in Bn()
70 t2 = t*t; 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()
77 t2 = r-1.0; 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
36 t2 = r*s; 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()
43 t2 = t1*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()
53 t2 = t1*t1; 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()
[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
78 t2 = s+r-1.0; in BnDrv()
79 t3 = t*t2; in BnDrv()
80 t4 = t2*t2; in BnDrv()
88 t2 = s+r-1.0; in BnDrv()
89 t3 = t*t2; in BnDrv()
90 t4 = t2*t2; in BnDrv()
98 t2 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0; 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 | modeShapeDrv.c | 30 double rs,rst2,P1P2,t,P1,P2,dP1,dP2,t1,t2; in F_modeShapeTriDrv() local
|
| /phasta/phSolver/incompressible/ |
| H A D | e3lhs.f | 45 dimension t1(npro,3), t2(npro,3), t3(npro,3),
110 t2(:,1) = rmu * shg(:,b,1)
111 t2(:,2) = rmu * shg(:,b,2)
112 t2(:,3) = rmu * shg(:,b,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)
144 & + t2(:,1) * shg(:,aa,1)
147 & + t2(:,2) * shg(:,aa,2)
150 & + t2(:,3) * shg(:,aa,3)
[all …]
|
| H A D | e3source.f | 178 t2 = x**2
179 t3 = t2**2
181 t7 = t2*x
187 t30 = t3*t2
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)
|
| /phasta/phSolver/common/ |
| H A D | blowerControl.f90 | 101 real*8 :: t0, t1, t2, t3, t4, vmax, vmin local
128 t2 = t1 + blower(iBlower)%t_fullOn
129 t3 = t2 + blower(iBlower)%t_fallTime
138 elseif(t1 <= t .and. t < t2) then !high
140 elseif(t2 <= t .and. t < t3) then !falling edge
141 vmagref = (vmax - vmin)*(t - t3)/(t2 - t3) + vmin
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 59 c$$$ t2 = t1**2
61 c$$$ t6 = t2*x*t5
65 c$$$ t12 = t2*t5
75 c$$$ t43 = nu*t2
79 c$$$ t50 = t2*t22*t5
82 c$$$ t65 = t2**2
85 c$$$ t77 = t2*t1*t5
99 c$$$ t2 = t1**2
100 c$$$ t3 = nu*t2
110 c$$$ t33 = t2*t32
[all …]
|
| H A D | e3source.fRHSonly | 59 c$$$ t2 = t1**2
61 c$$$ t6 = t2*x*t5
65 c$$$ t12 = t2*t5
75 c$$$ t43 = nu*t2
79 c$$$ t50 = t2*t22*t5
82 c$$$ t65 = t2**2
85 c$$$ t77 = t2*t1*t5
99 c$$$ t2 = t1**2
100 c$$$ t3 = nu*t2
110 c$$$ t33 = t2*t32
[all …]
|
| /phasta/phSolver/AMG/ |
| H A D | ramg_paratools.f | 504 real(kind=8),dimension(nshg,3) :: t2,t2a local
516 call fLesSparseApG(colm,rowp,lhsP,t1,t2,nshg,nnz_tot)
517 call commIn(t2,ilwork,3,iper,iBC,BC)
519 call fMtxVdimVecMult(t2,diag,t2a,3,4,3,3,nshg)
520 call fMtxVdimVecMult(t2a,diag,t2,3,4,3,3,nshg)
524 t2b(:,1:3) = -t2(:,1:3)
561 real(kind=8),dimension(nshg,3) :: t2,t2a local
570 call fLesSparseApG(colm,rowp,lhsP,t1,t2,nshg,nnz_tot)
571 call commIn(t2,ilwork,3,iper,iBC,BC)
573 call fMtxVdimVecMult(t2,diag,t2a,3,4,3,3,nshg)
[all …]
|