Home
last modified time | relevance | path

Searched refs:r (Results 1 – 25 of 25) sorted by relevance

/phasta/phSolver/compressible/
H A Di3ldu.f1 subroutine i3LDU (Diag, r, code) argument
44 dimension Diag(nshg,nsymdf), r(nshg,nflow)
105 r(:,1) = Diag(:, 1) * r(:,1)
106 r(:,2) = Diag(:, 3) * ( r(:,2)
107 & - r(:,1) * Diag(:, 2) )
108 r(:,3) = Diag(:, 6) * ( r(:,3)
109 & - r(:,1) * Diag(:, 4)
110 & - r(:,2) * Diag(:, 5) )
111 r(:,4) = Diag(:,10) * ( r(:,4)
112 & - r(:,1) * Diag(:, 7)
[all …]
H A Di3lu.f1 subroutine i3LU (Diag, r, code) argument
34 dimension Diag(nshg,nflow,nflow), r(nshg,nflow)
105 r(:,2) = r(:,2) - Diag(:,2,1) * r(:,1)
106 r(:,3) = r(:,3) - Diag(:,3,1) * r(:,1)
107 & - Diag(:,3,2) * r(:,2)
108 r(:,4) = r(:,4) - Diag(:,4,1) * r(:,1)
109 & - Diag(:,4,2) * r(:,2)
110 & - Diag(:,4,3) * r(:,3)
111 r(:,5) = r(:,5) - Diag(:,5,1) * r(:,1)
112 & - Diag(:,5,2) * r(:,2)
[all …]
H A De3source.f278 & r (npro), g (npro),
327 r(:) = sclrm(:)/Stilde(:)/(saKappa*dist2w(:))**2
329 r(:) = 1.0d32
333 g = r + saCw2*(r**6-r)
399 gp = one+saCw2*(6*r**5 - one)
/phasta/shapeFunction/src/
H A DFn.c18 double Fn(int i, int j, double r, double s) { in Fn() argument
32 f = r-1.0/3.0; in Fn()
39 f = r*s-r/4-s/4+1.0/14.0; in Fn()
42 t1 = r*r; in Fn()
43 f = t1-3.0/4.0*r+3.0/28.0; in Fn()
51 f = r*t1-3.0/5.0*r*s+r/15-t1/5+2.0/15.0*s-1.0/60.0; in Fn()
54 t1 = r*r; in Fn()
55 f = t1*s-t1/5-3.0/5.0*r*s+2.0/15.0*r+s/15-1.0/60.0; in Fn()
58 t1 = r*r; in Fn()
59 f = r*t1-6.0/5.0*t1+2.0/5.0*r-1.0/30.0; in Fn()
[all …]
H A DBn.c19 double Bn(int i, int j, int k, double r, double s, double t) { in Bn() argument
30 f = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0; in Bn()
33 f = s+2.0/5.0*r-2.0/5.0; in Bn()
36 f = r-2.0/7.0; in Bn()
40 t2 = s+r-1.0; in Bn()
45 f = (s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); in Bn()
49 t2 = r-1.0; in Bn()
54 f = (r-2.0/9.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); in Bn()
57 f = (r-2.0/9.0)*(s+2.0/5.0*r-2.0/5.0); in Bn()
60 t1 = r*r; in Bn()
[all …]
H A DFnDrv.c18 int FnDrv(int i,int j, double r, double s, double drv[2]) { in FnDrv() argument
43 drv[1] = r-1.0/4.0; in FnDrv()
46 drv[0] = 2.0*r-3.0/4.0; in FnDrv()
57 drv[1] = 2.0*r*s-3.0/5.0*r-2.0/5.0*s+2.0/15.0; in FnDrv()
60 t3 = r*r; in FnDrv()
61 drv[0] = 2.0*r*s-2.0/5.0*r-3.0/5.0*s+2.0/15.0; in FnDrv()
62 drv[1] = t3-3.0/5.0*r+1.0/15.0; in FnDrv()
65 t1 = r*r; in FnDrv()
66 drv[0] = 3.0*t1-12.0/5.0*r+2.0/5.0; in FnDrv()
77 drv[1] = 3.0*r*t1-2.0*r*s+3.0/11.0*r-t1/2+4.0/11.0*s-3.0/55.0; in FnDrv()
[all …]
H A DEn.c17 double En(int ip, double r, double s) { in En() argument
27 f = s-r; in En()
31 t3 = r*r; in En()
32 f = t1-3.0*r*s+t3; in En()
36 t4 = r*r; in En()
37 f = t1*s-6.0*r*t1+6.0*t4*s-t4*r; in En()
42 t5 = r*r; 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()
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()
[all …]
H A DBnDrv.c19 int BnDrv(int i,int j, int k, double r, double s, double t, double drv[3]) { in BnDrv() argument
46 t1 = 6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0; in BnDrv()
49 drv[2] = 2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0; in BnDrv()
52 drv[0] = 2.0/7.0*t+6.0/7.0*s+8.0/21.0*r-8.0/21.0; in BnDrv()
53 drv[1] = t+4.0/3.0*s+6.0/7.0*r-6.0/7.0; in BnDrv()
54 drv[2] = s+2.0/7.0*r-2.0/7.0; in BnDrv()
57 drv[0] = 6.0/7.0*s+2.0/7.0*r-2.0/7.0; in BnDrv()
58 drv[1] = 2.0*s+6.0/7.0*r-6.0/7.0; in BnDrv()
62 drv[0] = t+2.0/3.0*s+4.0/3.0*r-22.0/27.0; in BnDrv()
63 drv[1] = 2.0/3.0*r-4.0/27.0; in BnDrv()
[all …]
H A DEnDrv.c17 int EnDrv(int ip, double r, double s, double drv[2]) { in EnDrv() argument
31 drv[0] = -3.0*s+2.0*r; in EnDrv()
32 drv[1] = 2.0*s-3.0*r; in EnDrv()
36 t2 = r*s; in EnDrv()
37 t3 = r*r; in EnDrv()
44 t3 = r*t1; in EnDrv()
45 t4 = r*r; in EnDrv()
47 t6 = t4*r; in EnDrv()
54 t4 = r*t1*s; in EnDrv()
55 t5 = r*r; in EnDrv()
[all …]
H A DLagrange.c24 double r=par[0],s=par[1],t=par[2],t8,t4; in TetQuadLag() local
25 t8 = 1.0-r-s-t; in TetQuadLag()
26 N[0] = (2.0*r-1.0)*r; in TetQuadLag()
29 N[3] = (1.0-2.0*r-2.0*s-2.0*t)*t8; in TetQuadLag()
30 N[4] = 4.0*r*s; in TetQuadLag()
32 N[6] = 4.0*t*r; in TetQuadLag()
33 N[7] = 4.0*r*t8; in TetQuadLag()
37 t4 = -3.0+4.0*r+4.0*s+4.0*t; in TetQuadLag()
38 dN[0][0] = 4.0*r-1.0; in TetQuadLag()
51 dN[4][1] = 4.0*r; in TetQuadLag()
[all …]
H A DblendDrv.c55 double r,s; in F_edgeBlendTriDrv() local
58 r = L[0] ; s = L[1] ; in F_edgeBlendTriDrv()
64 drv[1] = -2.0*r; in F_edgeBlendTriDrv()
68 drv[1] = -2.0+2.0*r+4.0*s; in F_edgeBlendTriDrv()
71 drv[0] = 4.0*r-2.0+2.0*s; in F_edgeBlendTriDrv()
72 drv[1] = 2.0*r; in F_edgeBlendTriDrv()
96 double r=L[0],s=L[1],t=L[2]; in R_edgeBlendTetDrv() local
102 drv[1] = -2.0*r; in R_edgeBlendTetDrv()
106 drv[0] = -2.0+4.0*r+2.0*s+2.0*t; in R_edgeBlendTetDrv()
107 drv[1] = 2.0*r; in R_edgeBlendTetDrv()
[all …]
H A DshapeFuncInternals.h55 double En(int ip, double r, double s);
56 int EnDrv(int ip, double r, double s, double *drv);
57 double Fn(int i, int j, double r, double s);
58 int FnDrv(int i,int j, double r, double s, double drv[2]);
59 double Bn(int i, int j, int k, double r, double s, double t);
60 int BnDrv(int i,int j, int k, double r, double s, double t, double drv[3]);
/phasta/phSolver/common/
H A DGaussLegendreSimplex.c9 double *r, double *s, double *t, double *J);
12 double *r, double *s,double *J);
57 double *r, double *s, double *t, double *J) { in brickToTet() argument
59 *r = 0.5e0*(1.0e0+xi); in brickToTet()
60 r1 = 1.0e0-(*r); in brickToTet()
62 rs1 = 1.0e0-(*r)-(*s); in brickToTet()
67 void quadToTri(double xi,double eta,double *r, double *s, double *J) { in quadToTri() argument
69 *r = 0.5e0*(1.0e0+xi); in quadToTri()
70 r1 = 1.0e0-(*r); in quadToTri()
H A DphComm.c5 int r; in phcomm_rank() local
6 MPI_Comm_rank(MPI_COMM_WORLD, &r); in phcomm_rank()
7 return r; in phcomm_rank()
H A Dfillsparse.f87 integer f, g, r, s, t local
109 r = (aa-1)*nflow !starting index of the ath node in EGmass
121 lhsK(t+f,k) = lhsK(t+f,k) + EGmass(e,r+f,s+g)
149 integer f, g, r, s, t local
H A Dinput.config132 Radial Cooling Sponge Begins at r: 7.5
/phasta/phSolver/AMG/
H A Dramg_cheby.f9 subroutine ramg_cheby_apply(u,v,r,level,colm,rowp,lhsK,lhsP, argument
26 real(kind=8),intent(in),dimension(amg_nshg(level)) :: v,r local
51 cy1 = r
55 cy1 = v-1.1*3.0/beta*(pAux-r)
60 cz1 = r-pAux
79 epdk = r*invtheta!/theta
84 epdk = r-pAux
96 epdk = dtemp1*epdk+dtemp2*(r-pAux)
H A Dramg_smoother.f59 & r,u,v)
66 real(kind=8),intent(in),dimension(Anshg) :: v,r local
84 u(i) = damp_jacobi*u(i) + r(i)
97 & r,u,v,clevel,pflag,
108 real(kind=8),intent(in),dimension(anshg) :: v,r local
133 tmp = tmp-r
168 & r,u,v,fwdbck,clevel,
179 real(kind=8),intent(in),dimension(anshg) :: v,r local
214 & r,u,v,clevel,pflag,ilwork,BC,iBC,iper)
233 u(ki) = r(ki) - tmp
[all …]
H A Dramg_mls.f563 subroutine ramg_mls_apply(u,v,r,level,colm,rowp,lhsK,lhsP, argument
580 real(kind=8),intent(in),dimension(amg_nshg(level)) :: v,r local
590 aux = -r
594 aux = aux-r
610 subroutine ramg_mls_apply_fwd(u,v,r,level,colm,rowp,lhsK,lhsP, argument
627 real(kind=8),intent(in),dimension(amg_nshg(level)) :: v,r local
635 pAux = r
651 subroutine ramg_mls_apply_post(u,v,r,level,colm,rowp,lhsK,lhsP, argument
668 real(kind=8),intent(in),dimension(amg_nshg(level)) :: v,r local
676 call ramg_mls_apply_fwd(y1,v,r,level,colm,rowp,lhsK,lhsP,
H A Dramg_ITAI.f122 integer :: i,j,k,m,p,q,r,n,s,p0 local
/phasta/phSolver/incompressible/
H A De3source.f17 real*8 r, Stheta, dpdz, rP5 local
26 r = sqrt( xx(iel,1)**2 + xx(iel,2)**2)
27 rP5 = (r/bigR)**5
31 src(iel,1) = -xx(iel,2)/r * Stheta
32 src(iel,2) = xx(iel,1)/r * Stheta
253 real*8 chi, chiP3, fv1, fv2, st, r, local
403 r = zero
406 r = tmp / st
407 r = min( max(r, -8.0d0), 8.0d0)
410 rP5 = r * (r * r) ** 2
[all …]
/phasta/phSolver/testIC-DSIM/
H A Dgeom.sms161 1 24+e+h+r-p
163 1 24-r+i-k-s
174 0 34-r-v-n+x
H A Dgeom.sdm160 1 24+e+h+r-p
162 1 24-r+i-k-s
/phasta/phSolver/incompressible/e3source/
H A De3source.fRHSonly17 real*8 r, Stheta, dpdz, rP5
26 r = sqrt( xx(iel,1)**2 + xx(iel,2)**2)
27 rP5 = (r/bigR)**5
31 src(iel,1) = -xx(iel,2)/r * Stheta
32 src(iel,2) = xx(iel,1)/r * Stheta
253 real*8 chi, chiP3, fv1, fv2, st, r,
365 r = zero
368 r = tmp / st
369 r = min( max(r, -8.0d0), 8.0d0)
372 rP5 = r * (r * r) ** 2
[all …]
H A De3source.fimplicit17 real*8 r, Stheta, dpdz, rP5
26 r = sqrt( xx(iel,1)**2 + xx(iel,2)**2)
27 rP5 = (r/bigR)**5
31 src(iel,1) = -xx(iel,2)/r * Stheta
32 src(iel,2) = xx(iel,1)/r * Stheta
253 real*8 chi, chiP3, fv1, fv2, st, r,
379 r = zero
382 r = tmp / st
383 r = min( max(r, -8.0d0), 8.0d0)
386 rP5 = r * (r * r) ** 2
[all …]