| /phasta/phSolver/incompressible/ |
| H A D | e3ql.f | 25 dimension yl(npro,nshl,ndof), dwl(npro,nshl), 26 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 27 & xl(npro,nenl,nsd), sgn(npro,nshl), 28 & ql(npro,nshl,idflx), xmudmi(npro,ngauss) 33 & g3yi(npro,ndof), shg(npro,nshl,nsd), 36 & rminv(npro,nshl,nshl), 37 & qrl(npro,nshl,nsd*nsd) 39 dimension qdi(npro,nsd*nsd), shape(npro,nshl), 40 & shdrv(npro,nsd,nshl), indx(nshl), 41 & rmass(npro,nshl,nshl) [all …]
|
| H A D | asigmr.f | 24 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 25 & ien(npro,nshl), 30 dimension yl(npro,nshl,ndofl), acl(npro,nshl,ndofl), 32 & rl(npro,nshl,nflow), 33 & ql(npro,nshl,idflx) 35 dimension xKebe(npro,9,nshl,nshl), 36 & xGoC(npro,4,nshl,nshl) 38 dimension rlsl(npro,nshl,6) 41 real*8 lStsVec(npro,nshl,nResDims) 44 dimension sgn(npro,nshl) [all …]
|
| H A D | e3.f | 40 dimension yl(npro,nshl,ndof), 41 & acl(npro,nshl,ndof), 42 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 44 & rl(npro,nshl,nflow), ql(npro,nshl,idflx) 47 dimension xKebe(npro,9,nshl,nshl), xGoC(npro,4,nshl,nshl) 52 & g3yi(npro,ndof), shg(npro,nshl,nsd), 58 & xmudmi(npro,ngauss), sgn(npro,nshl), 59 & shpfun(npro,nshl), shdrv(npro,nsd,nshl), 64 dimension rlsl(npro,nshl,6), rlsli(npro,6) 66 real*8 rerrl(npro,nshl,6) [all …]
|
| H A D | asbmfg.f | 21 & shpb(nshl,ngaussb), 22 & shglb(nsd,nshl,ngaussb), 23 & ienb(npro,nshl), materb(npro), 27 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd), 28 & rl(npro,nshl,nflow), sgn(npro,nshl), 29 & ul(npro,nshl,nsd), acl(npro,nshl,ndofl) 31 dimension xKebe(npro,9,nshl,nshl) 94 & shpb(nshl,*), 95 & shglb(nsd,nshl,*), 96 & ienb(npro,nshl), materb(npro), [all …]
|
| H A D | e3q.f | 24 dimension yl(npro,nshl,ndof), dwl(npro,nenl), 25 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 27 & ql(npro,nshl,idflx), rmassl(npro,nshl), 33 & g3yi(npro,nflow), shg(npro,nshl,nsd), 39 dimension sgn(npro,nshl), shape(npro,nshl), 40 & shdrv(npro,nsd,nshl), shpsum(npro) 111 do i=1,nshl 133 do i=1,nshl 142 do i=1,nshl 167 do n = 1, nshl [all …]
|
| H A D | asiq.f | 28 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 29 & ien(npro,nshl), dwl(npro,nenl), 32 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 33 & ql(npro,nshl,idflx), rmassl(npro,nshl), 36 dimension sgn(npro,nshl) 41 do i=1,nshl 96 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 97 & ien(npro,nshl), dwl(npro,nenl), 100 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 101 & ql(npro,nshl,nsd), rmassl(npro,nshl) [all …]
|
| H A D | errsmooth.f | 40 nshl = lcblk(10,iblk) 50 allocate (tmpshp(nshl,MAXQPT)) 51 allocate (tmpshgl(nsd,nshl,MAXQPT)) 53 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 54 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 135 & shp(nshl,maxsh), 136 & shgl(nsd,nshl,maxsh), 137 & ien(npro,nshl), 142 dimension rerrl(npro,nshl,10), xl(npro,nenl,nsd), 143 & rerrsml(npro,nshl,10), rmassl(npro,nshl) [all …]
|
| H A D | asbflx.f | 20 & shpb(nshl,ngaussb), 21 & shglb(nsd,nshl,ngaussb), 22 & ienb(npro,nshl), 27 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd), 28 & rl(npro,nshl,nflow), sgn(npro,nshl), 29 & flhsl(npro,nshl,1), fnrml(npro,nshl,nsd), 31 & ul(npro,nshl,nsd), acl(npro,nshl,ndofl) 32 real*8 dwl(npro,nshl) 34 dimension xKebe(npro,9,nshl,nshl) 85 do n = 1, nshl
|
| H A D | hessian.f | 33 nshl = lcblk(10,iblk) 42 & shp(lcsyst,1:nshl,:), 43 & shgl(lcsyst,:,1:nshl,:), 62 nshl = lcblk(10,iblk) 72 & shp(lcsyst,1:nshl,:), 73 & shgl(lcsyst,:,1:nshl,:), 93 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 94 & ien(npro,nshl), gradu(nshg,9), 95 & shdrv(npro,nsd,nshl), shape( npro, nshl ), 98 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), [all …]
|
| H A D | genlmass.f | 31 nshl = lcblk(10,iblk) 40 allocate (tmpshp(nshl,MAXQPT)) 41 allocate (tmpshgl(nsd,nshl,MAXQPT)) 42 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 43 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 75 & shp(nshl,maxsh), shgl(nsd,nshl,ngauss), 78 integer ien(npro,nshl) 82 & sgn(npro,nshl), shape(npro,nshl), 83 & locmass(npro,nshl), shg(npro,nshl,nsd), 85 & dxidx(npro,nsd,nsd), shdrv(npro,nsd,nshl) [all …]
|
| /phasta/phSolver/common/ |
| H A D | cmass.f | 34 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT), 35 & em(npro,nshl,nshl), 38 dimension shape(npro,nshl), shdrv(npro,nsd,nshl), 39 & sgn(npro,nshl), dxidx(npro,nsd,nsd), 40 & shg(npro,nshl,nsd), 59 do na = 1, nshl 60 do nb = 1, nshl 87 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT), 88 & shpf(nshl,MAXQPT), shglf(nsd,nshl,MAXQPT), 89 & em(npro,nshl,nshl), eml(npro,nshl), [all …]
|
| H A D | hierarchic.f | 17 dimension ien(npro,nshl), sgn(npro,nshl) 19 do i=nenl+1,nshl 37 dimension shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 38 & sgn(npro,nshl), shape(npro,nshl), 39 & shdrv(npro,nsd,nshl) 49 do i=nenl+1,nshl 67 dimension shp(nshl,ngaussb), shgl(nsd,nshl,ngaussb), 68 & sgn(npro,nshl), shape(npro,nshl), 69 & shdrv(npro,nsd,nshl) 79 do i=nenl+1,nshl [all …]
|
| H A D | fillsparse.c | 24 int nshl = shpdat.nshl; in fillsparsecpetscs() local 25 double* mb = (double*) malloc(sizeof(double)*nshl*nshl); //block to insert in fillsparsecpetscs() 27 PetscInt* locat = (PetscInt*) malloc(sizeof(PetscInt)*nshl); in fillsparsecpetscs() 30 for(aa=0;aa<nshl;aa++) locat[aa]=ieng[e+npro*aa]-1; in fillsparsecpetscs() 32 for (i=0; i<nshl; i++) { // fill up Ke with respective egmass in fillsparsecpetscs() 33 for (j=0; j<nshl; j++) { in fillsparsecpetscs() 34 mb[nshl*i + j] = EGmass[e + npro*(i + nshl*j)]; in fillsparsecpetscs() 39 petsc_nshl = (PetscInt) nshl; in fillsparsecpetscs() 48 int nshl = shpdat.nshl; in fillsparsecpetscc() local 53 PetscInt* locat = (PetscInt*) malloc(sizeof(PetscInt)*nshl); in fillsparsecpetscc() [all …]
|
| H A D | errorsmooth.f | 39 nshl = lcblk(10,iblk) 49 allocate (tmpshp(nshl,MAXQPT)) 50 allocate (tmpshgl(nsd,nshl,MAXQPT)) 52 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 53 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 134 & shp(nshl,maxsh), 135 & shgl(nsd,nshl,maxsh), 136 & ien(npro,nshl), 141 dimension rerrl(npro,nshl,10), xl(npro,nenl,nsd), 142 & rerrsml(npro,nshl,10), rmassl(npro,nshl) [all …]
|
| H A D | fillsparse.f | 8 real*8 xKebe(npro,9,nshl,nshl), xGoC(npro,4,nshl,nshl) 9 integer ien(npro,nshl), col(nshg+1), row(nshg*nnz) 16 integer iens(npro,nshl) 27 do aa = 1, nshl ! loop over the local equation numbers 32 do b = 1, nshl ! local variable number tangent respect 82 integer ien(npro,nshl), col(nshg+1), row(nnz*nshg) 91 integer iens(npro,nshl) 102 do aa = 1, nshl !loop over matrix block column index 110 do b = 1, nshl 144 real*8 EGmass(npro, nshl, nshl) !only contains term (5,5) [all …]
|
| H A D | local.f | 23 dimension global(nshg,n), rlocal(npro,nshl,n), 24 & ien(npro,nshl), ientmp(npro,nshl) 45 do i = 1, nshl 54 gbytes = gbytes + n*nshl*npro 68 do i = 1, nshl 79 sbytes = sbytes + n*nshl*npro 80 flops = flops + n*nshl*npro 94 do i = 1, nshl 137 & ien(npro,nshl) 356 dimension global(nshg,n), rlocal(npro,nshl,n), [all …]
|
| /phasta/phSolver/compressible/ |
| H A D | asigmr.f | 21 & shp(nshl,MAXQPT), 22 & shgl(nsd,nshl,MAXQPT), 23 & ien(npro,nshl), 29 dimension ycl(npro,nshl,ndofl), acl(npro,nshl,ndof), 30 & xl(npro,nenl,nsd), ytargetl(npro,nshl,nflow), 31 & rl(npro,nshl,nflow), rml(npro,nshl,nflow), 32 & BDiagl(npro,nshl,nflow,nflow), 33 & ql(npro,nshl,idflx) 36 dimension sgn(npro,nshl), EGmass(npro,nedof,nedof) 38 dimension rlsl(npro,nshl,6) [all …]
|
| H A D | elmgmr.f | 77 nshl = lcblk(10,iblk) 87 allocate (tmpshp(nshl,MAXQPT)) 88 allocate (tmpshgl(nsd,nshl,MAXQPT)) 90 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 91 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 137 nshl = lcblk(10,iblk) 148 allocate (tmpshp(nshl,MAXQPT)) 149 allocate (tmpshgl(nsd,nshl,MAXQPT)) 151 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 152 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) [all …]
|
| H A D | elmgmrpetsc.f | 91 nshl = lcblk(10,iblk) 101 allocate (tmpshp(nshl,MAXQPT)) 102 allocate (tmpshgl(nsd,nshl,MAXQPT)) 104 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 105 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 152 nshl = lcblk(10,iblk) 170 allocate (tmpshp(nshl,MAXQPT)) 171 allocate (tmpshgl(nsd,nshl,MAXQPT)) 172 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 173 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) [all …]
|
| H A D | asbmfg.f | 16 & shpb(nshl,ngaussb), 17 & shglb(nsd,nshl,ngaussb), 18 & ienb(npro,nshl), materb(npro), 22 dimension ycl(npro,nshl,ndofl), xlb(npro,nenl,nsd), 23 & rl(npro,nshl,nflow), 24 & rml(npro,nshl,nflow), 25 & EGmass(npro, nshl, nshl) 27 dimension sgn(npro,nshl) 83 & shpb(nshl,maxsh), 84 & shglb(nsd,nshl,maxsh), [all …]
|
| H A D | e3ql.f | 24 dimension ycl(npro,nshl,ndof), 25 & shp(nshl,ngauss), 26 & shgl(nsd,nshl,ngauss), 27 & xl(npro,nenl,nsd), sgn(npro,nshl), 28 & ql(npro,nshl,idflx), xmudmi(npro,ngauss) 33 & g3yi(npro,nflow), shg(npro,nshl,nsd), 39 & con(npro), rminv(npro,nshl,nshl), 40 & qrl(npro,nshl,(nflow-1)*nsd) 42 dimension qdi(npro,nsd*(nflow-1)), shape(npro,nshl), 43 & shdrv(npro,nsd,nshl), indx(nshl), [all …]
|
| H A D | elmmfg.f | 76 nshl = lcblk(10,iblk) 86 allocate (tmpshp(nshl,MAXQPT)) 87 allocate (tmpshgl(nsd,nshl,MAXQPT)) 89 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 90 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 135 nshl = lcblk(10,iblk) 145 allocate (tmpshp(nshl,MAXQPT)) 146 allocate (tmpshgl(nsd,nshl,MAXQPT)) 148 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 149 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) [all …]
|
| H A D | asiflx.f | 19 & ien(npro,nshl), 22 dimension ycl(npro,nshl,ndof), acl(npro,nshl,ndof), 24 & rl(npro,nshl,nflow), BDiagl(npro,nshl,nflow,nflow) 26 dimension ql(npro,nshl,(nflow-1)*nsd) 29 dimension sgn(npro,nshl) 30 dimension EGmass(npro,nedof,nedof), rml(npro,nshl,nflow) 32 dimension rlsl(npro,nshl,6) 63 call e3 (ycl, ycl, acl, shp(lcsyst,1:nshl,:), 64 & shgl(lcsyst,:,1:nshl,:), xl, rl, rml,
|
| H A D | asimfg.f | 19 & shp(nshl,MAXQPT), 20 & shgl(nsd,nshl,MAXQPT), 21 & ien(npro,nshl), 27 dimension ycl(npro,nshl,ndofl), acl(npro,nshl,ndofl), 29 & rl(npro,nshl,nflow), rml(npro,nshl,nflow), 30 & BDiagl(npro,nshl,nflow,nflow), 31 & ql(npro,nshl,idflx) 33 dimension rlsl(npro,nshl,6) 35 dimension sgn(npro,nshl) 37 real*8 rerrl(npro,nshl,6), rerr(nshg,10)
|
| H A D | asires.f | 18 & shp(nshl,ngauss), 19 & shgl(nsd,nshl,ngauss), 20 & ien(npro,nshl), mater(npro), 23 dimension yl(npro,nshl,nflow), ycl(npro,nshl,ndofl), 24 & xl(npro,nenl,nsd), acl(npro,nshl,nflow), 25 & rml(npro,nshl,nflow), ql(npro,nshl,(nflow-1)*nsd) 28 dimension sgn(npro,nshl) 30 dimension rlsl(npro,nshl,6) 32 real*8 rerrl(npro,nshl,6)
|