Home
last modified time | relevance | path

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

1234

/phasta/phSolver/AMG/
H A Dramg_extract.f48 integer :: i,j,k,m,n,p,q,ki,kj,ni,nj,qq
69 deallocate(ramg_flowDiag%p)
75 allocate(amg_paramap(1)%p(nshg))
76 allocate(amg_paraext(1)%p(nshg))
97 allocate(ramg_flowDiag%p(nshg,4))
109 ramg_flowDiag%p(:,i) = mflowDiag(:,i)
118 call ramg_initBCflag(amg_paramap(1)%p,ilwork,BC,iBC,iper)
127 amg_paraext(1)%p = amg_paramap(1)%p
129 if (amg_paramap(1)%p(i).ne.(myrank+1)) then
132 do m=lhsGPcolm%p(i),lhsGPcolm%p(i+1)-1
[all …]
H A Dramg_ITAI.f14 integer :: i,j,k,p,q,mnnz local
23 do p=b_colm(j),b_colm(j+1)-1
24 q=b_rowp(p)
46 integer :: i,j,k,p,q,mnnz local
59 do p=b_colm(j),b_colm(j+1)-1
60 q=b_rowp(p)
87 integer :: i,j,k,p,q local
99 do p = b_colm(j),b_colm(j+1)-1
100 q = b_rowp(p)
101 c_mtx(itemp(q)) = c_mtx(itemp(q)) + rtemp*b_mtx(p)
[all …]
H A Dramg_paratools.f31 integer :: numtask,i,j,k,m,p,ki,kj,krowp,ierr local
57 allocate(lhsGP%p(redun,annz_tot),stat=ierr)
58 allocate(lhsGProwp%p(annz_tot))
59 allocate(lhsGPcolm%p(nshg+1))
60 lhsGP%p(:,:) = alhsP(:,:)
61 lhsGProwp%p(:) = arowp(:)
62 lhsGPcolm%p(:) = acolm(:)
81 allocate(sub_nnz%p(numtask))
82 allocate(sub_nnz2%p(numtask))
83 allocate(sub_nshg%p(numtask))
[all …]
H A Dramg_control.f50 call ramg_direct_LU(amg_A_colm(ramg_levelx)%p,
51 & amg_A_rowp(ramg_levelx)%p,
52 & amg_A_lhs(ramg_levelx)%p,amg_nshg(ramg_levelx),
78 integer :: i,p,p2 variable
100 p = 1
101 if (maxstopsign) p = 0
102 call MPI_AllReduce(p,p2,1,MPI_INTEGER,MPI_SUM,
119 deallocate(amg_A_colm(1)%p)
120 deallocate(amg_A_rowp(1)%p)
121 deallocate(amg_A_lhs(1)%p)
[all …]
H A Dramg_coarse.f38 integer :: i,j,k,m,n,p,q local
105 & amg_A_colm(level1)%p,amg_A_rowp(level1)%p,
106 & amg_A_lhs(level1)%p,amg_paramap(level1)%p,
128 p = iabs(amg_paramap(level1)%p(i))
129 p = subcfrev(p)
130 subnei(p) = subnei(p) + 1
132 subcf(p) = subcf(p) + 1
146 allocate(CF_map(level1)%p(amg_nshg(level1)))
147 allocate(CF_revmap(level1)%p(amg_nshg(level1)))
152 CF_map(level1)%p(k) = i
[all …]
H A Dramg_tools.f36 integer :: i,j,k,p local
44 do k=I_cf_colm(level1)%p(i),I_cf_colm(level1)%p(i+1)-1
45 j = I_cf_rowp(level1)%p(k)
46 vf(i) = vf(i) + VC(j)*I_cf(level1)%p(k)
72 integer :: i,j,k,p local
80 do k=I_fc_colm(level1)%p(i),I_fc_colm(level1)%p(i+1)-1
81 j = I_fc_rowp(level1)%p(k)
82 VC(i) = VC(i) + vf(j)*I_fc(level1)%p(k)
117 integer :: i,j,k,p local
157 integer :: i,j,k,p,mi,mj local
[all …]
H A Dramg_data.f6 real*8, pointer :: p(:) component
10 real*8, pointer :: p(:,:) component
14 real*8, pointer :: p(:,:,:) component
18 integer, pointer :: p(:) component
22 integer, pointer :: p(:,:) component
26 integer, pointer :: p(:,:,:) component
/phasta/phSolver/compressible/
H A Dsparseap.f1 subroutine SparseAp(iper, ilwork, iBC, col, row, lhsK, p) argument
21 real*8 p(nshg,nflow), q(nshg,nflow) local
28 call commu (p, ilwork, nflow , 'out')
34 p(:,j)=p(iper(:),j)
40 & call rotabc(p(1,2), iBC, 'out')
59 tmp1 = tmp1 + lhsK(1 ,k)*p(j,1)
60 1 + lhsK(6 ,k)*p(j,2)
61 2 + lhsK(11,k)*p(j,3)
62 3 + lhsK(16,k)*p(j,4)
63 4 + lhsK(21,k)*p(j,5)
[all …]
H A Delmgmr.f96 & mien(iblk)%p, mxmudmi(iblk)%p,
155 & x, mxmudmi(iblk)%p,
157 & tmpshgl, mien(iblk)%p,
158 & mmat(iblk)%p, res,
165 call bc3LHS (iBC, BC, mien(iblk)%p,
210 & mienb(iblk)%p, mmatb(iblk)%p,
211 & miBCB(iblk)%p, mBCB(iblk)%p,
378 & mien(iblk)%p, mxmudmi(iblk)%p,
442 & x, mxmudmi(iblk)%p,
444 & tmpshgl, mien(iblk)%p,
[all …]
H A Delmgmrpetsc.f110 & mien(iblk)%p, mxmudmi(iblk)%p,
176 & x, mxmudmi(iblk)%p,
178 & tmpshgl, mien(iblk)%p,
179 & mmat(iblk)%p, res,
187 call bc3LHS (iBC, BC, mien(iblk)%p,
195 call fillsparsecpetscc( mieng(iblk)%p, EGmass, lhsP)
246 & mienb(iblk)%p, mmatb(iblk)%p,
247 & miBCB(iblk)%p, mBCB(iblk)%p,
405 & mien(iblk)%p,
406 & mmat(iblk)%p, rest,
[all …]
/phasta/shapeFunction/src/
H A DmodeShape.c21 double E_modeShape(int p, double *L) { in E_modeShape() argument
24 if( p < 2 ) in E_modeShape()
27 return En(p-2,L[0],L[1]); in E_modeShape()
30 double F_modeShapeTri(int p, int i, double *L) { in F_modeShapeTri() argument
35 if( p < 3 || i < 0 || i > p-3 ) in F_modeShapeTri()
39 for(alpha=0; alpha <= p-3; alpha++) { in F_modeShapeTri()
40 for(beta=0; beta <= p-3; beta++) { in F_modeShapeTri()
41 if( alpha+beta == p-3 ) { in F_modeShapeTri()
59 double F_modeShapeQuad(int p, int i, double *L) { in F_modeShapeQuad() argument
64 double R_modeShapeTet(int p, int i, double *L) { in R_modeShapeTet() argument
[all …]
H A DmodeShapeDrv.c21 int E_modeShapeDrv(int p, double *L, double drv[2]) { in E_modeShapeDrv() argument
25 return EnDrv(p-2,L[0],L[1],drv); in E_modeShapeDrv()
28 int F_modeShapeTriDrv(int p, int i, double *L, double mdrv[2]) { in F_modeShapeTriDrv() argument
34 if( p < 3 || i < 0 || i > p-3 ) in F_modeShapeTriDrv()
38 for(alpha=0; alpha <= p-3; alpha++) { in F_modeShapeTriDrv()
39 for(beta=0; beta <= p-3; beta++) { in F_modeShapeTriDrv()
40 if( alpha+beta == p-3 ) { in F_modeShapeTriDrv()
58 int F_modeShapeQuadDrv(int p, int i, double *L, double mdrv[2]) { in F_modeShapeQuadDrv() argument
62 int R_modeShapeTetDrv(int p, int i, double *L, double mdrv[3]) { in R_modeShapeTetDrv() argument
68 if( p < 4 || i < 0 || i > (((p-2)*(p-3)/2)-1) ) in R_modeShapeTetDrv()
[all …]
H A DshapeFuncInternals.h42 int F_modeShapeTriDrv(int p, int i, double *L, double mdrv[2]);
43 int R_modeShapeTetDrv(int p, int i, double *L, double mdrv[3]);
44 int E_modeShapeDrv(int p, double *L, double drv[2]);
46 double E_modeShape(int p, double *L);
47 double F_modeShapeTri(int p, int i, double *L);
48 double F_modeShapeQuad(int p, int i, double *L);
49 double F_modeShapeQuad(int p, int i, double *L);
50 double R_modeShapeTet(int p, int i, double *L);
51 double R_modeShapeHex(int p, int i, double *L);
H A DuniformP.c8 int TetShapeAndDrv(int p,double par[3],double N[],double dN[][3]) { in TetShapeAndDrv() argument
14 if(p<1) in TetShapeAndDrv()
35 if( p > 1 ) { in TetShapeAndDrv()
43 for(ip=2; ip <= p; ip++) { in TetShapeAndDrv()
57 if( p > 2 ) { in TetShapeAndDrv()
66 for(ip=3; ip <= p; ip++) { in TetShapeAndDrv()
83 if( p > 3 ) { in TetShapeAndDrv()
85 for(ip=4; ip <= p; ip++) { in TetShapeAndDrv()
112 int TriShapeAndDrv(int p,double par[2],double N[],double dN[][2]){ in TriShapeAndDrv() argument
116 if(p > 2) /* not supported */ in TriShapeAndDrv()
[all …]
/phasta/phSolver/incompressible/
H A DlesSparse.f205 & p, q, nNodes,
213 real*8 pLhs(4,nnz_tot), p(nNodes), q(nNodes,3)
230 pisave = p(i)
253 1 p, q, nNodes,
264 real*8 p(nNodes,4), q(nNodes,3)
284 pisave = p(i,4)
289 1 + kLhs(1,k) * p(j,1)
290 2 + kLhs(4,k) * p(j,2)
291 3 + kLhs(7,k) * p(j,3)
293 1 + kLhs(2,k) * p(j,1)
[all …]
H A Delmgmr.f105 & mien(iblk)%p,
111 & mien(iblk)%p, mxmudmi(iblk)%p,
170 & x, mxmudmi(iblk)%p,
173 & mien(iblk)%p,
182 & call bc3lhs (iBC, BC,mien(iblk)%p, xKebe)
183 call fillsparseI (mien(iblk)%p,
294 & mienb(iblk)%p, mmatb(iblk)%p,
295 & miBCB(iblk)%p, mBCB(iblk)%p,
308 & call bc3lhs (iBC, BC,mienb(iblk)%p, xKebe)
309 call fillsparseI (mienb(iblk)%p,
[all …]
/phasta/phSolver/common/
H A Dpointer.f12 real*8, pointer :: p(:) component
16 real*8, pointer :: p(:,:) component
20 real*8, pointer :: p(:,:,:) component
24 integer, pointer :: p(:) component
28 integer, pointer :: p(:,:) component
32 integer*8, pointer :: p(:,:) component
36 integer, pointer :: p(:,:,:) component
H A Dtopo_shapedefs.h5 extern "C" double phi(int p, double x);
6 extern "C" double phiDrv(int p,double x);
7 extern "C" int HexShapeAndDrv(int p, double par[3], double N[], double
9 extern "C" int WedgeShapeAndDrv(int p, double Inputpar[3], double N[], double
11 extern "C" int PyrShapeAndDrv (int p, double Inputpar[3], double N[], double
43 extern "C" int mesh_edge(double xi1,int gOrd[3], int p, double* entfn,
45 extern "C" int quad_face(double xi[3], int gOrd[3], int p, double*
47 extern "C" int hex_regn(double xi[3], int p, double*
H A Dnewshape.cc43 double phi(int p, double x) in phi() argument
47 PHI = LP(p,x)-LP(p-2,x); in phi()
48 PHI = PHI / (2*p-1); in phi()
53 double phiDrv(int p, double x) in phiDrv() argument
56 Phidrv = LP(p-1,x); in phiDrv()
292 int mesh_edge(double xi1, int gOrd[3], int p,double* entfn,double** edrv) in mesh_edge() argument
294 int nem = p-1; in mesh_edge()
323 int quad_face(double xi[3], int gOrd[3], int p, double* entfn, double** edrv) in quad_face() argument
337 if(p > 3){ in quad_face()
339 nfm = (p-2)*(p-3)/2; in quad_face()
[all …]
H A DgenblkPosix.f74 allocate (mmat(nelblk)%p(npro))
76 allocate (mien(nelblk)%p(npro,nshl))
77 allocate (mxmudmi(nelblk)%p(npro,maxsh))
79 allocate (mienG(nelblk)%p(1,1))
81 allocate (mienG(nelblk)%p(npro,nshl))
92 & mater, mien(nelblk)%p,
93 & mienG(nelblk)%p,
94 & mmat(nelblk)%p)
H A DgenblkSyncIO.f117 allocate (mmat(nelblk)%p(npro))
119 allocate (mien(nelblk)%p(npro,nshl))
120 allocate (mxmudmi(nelblk)%p(npro,maxsh))
122 allocate (mienG(nelblk)%p(1,1))
124 allocate (mienG(nelblk)%p(npro,nshl))
135 & mater, mien(nelblk)%p,
136 & mienG(nelblk)%p,
137 & mmat(nelblk)%p)
H A Dgenblk.f126 allocate (mmat(nelblk)%p(npro))
128 allocate (mien(nelblk)%p(npro,nshl))
129 allocate (mxmudmi(nelblk)%p(npro,maxsh))
131 allocate (mienG(nelblk)%p(1,1))
133 allocate (mienG(nelblk)%p(npro,nshl))
144 & mater, mien(nelblk)%p,
145 & mienG(nelblk)%p,
146 & mmat(nelblk)%p)
H A DgenbkbPosix.f133 allocate (mienb(nelblb)%p(npro,nshl))
134 allocate (miBCB(nelblb)%p(npro,ndiBCB))
135 allocate (mBCB(nelblb)%p(npro,nshlb,ndBCB))
136 allocate (mmatb(nelblb)%p(npro))
142 & materb, mienb(nelblb)%p,
143 & miBCB(nelblb)%p, mBCB(nelblb)%p,
144 & mmatb(nelblb)%p)
173 2100 format(2x,i5,1p,1x,6e12.4)
H A Dgenbkb.f188 allocate (mienb(nelblb)%p(npro,nshl))
189 allocate (miBCB(nelblb)%p(npro,ndiBCB))
190 allocate (mBCB(nelblb)%p(npro,nshlb,ndBCB))
191 allocate (mmatb(nelblb)%p(npro))
197 & materb, mienb(nelblb)%p,
198 & miBCB(nelblb)%p, mBCB(nelblb)%p,
199 & mmatb(nelblb)%p)
219 2100 format(2x,i5,1p,1x,6e12.4)
H A DgenbkbSyncIO.f174 allocate (mienb(nelblb)%p(npro,nshl))
175 allocate (miBCB(nelblb)%p(npro,ndiBCB))
176 allocate (mBCB(nelblb)%p(npro,nshlb,ndBCB))
177 allocate (mmatb(nelblb)%p(npro))
183 & materb, mienb(nelblb)%p,
184 & miBCB(nelblb)%p, mBCB(nelblb)%p,
185 & mmatb(nelblb)%p)
215 2100 format(2x,i5,1p,1x,6e12.4)

1234