| /phasta/phSolver/incompressible/ |
| H A D | e3lhs.f | 49 integer aa, b local
86 do aa = 1, nshl
87 tmp1 = tsFct * shpfun(:,aa) * shpfun(:,b)
88 tmp2 = tmp1 + t1(:,1) * shpfun(:,aa)
94 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp2
95 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp2
96 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp2
135 aa = b
137 tmp = t3(:,1) * shg(:,aa,1)
138 & + t3(:,2) * shg(:,aa,2)
[all …]
|
| H A D | e3res.f | 59 integer aa local
227 do aa = 1,nshl
231 rl(:,aa,4) = rl(:,aa,4) + WdetJ
232 & * ( shg(:,aa,1) * uBar(:,1)
233 & + shg(:,aa,2) * uBar(:,2)
234 & + shg(:,aa,3) * uBar(:,3) )
238 rl(:,aa,1) = rl(:,aa,1) - WdetJ
239 & * ( shpfun(:,aa) * rNa(:,1)
240 & + shg(:,aa,1) * rGNa(:,1,1)
241 & + shg(:,aa,2) * rGNa(:,1,2)
[all …]
|
| H A D | genlmass.f | 87 integer aa local
127 do aa = 1,nshl
128 locmass(:,aa)= locmass(:,aa)
129 & + shape(:,aa)*shape(:,aa)*WdetJ
150 do aa = 1, nshl
151 temp = temp + locmass(:,aa) !reusing temp to save memory
158 do aa = 1, nshl
159 locmass(:,aa) = locmass(:,aa) * fmstot / temp
|
| H A D | e3bvar.f | 111 integer aa, b local
661 do aa = 1, nshlb
662 tmp1 = tsFctvw * shpb(:,aa) * shpb(:,b)
666 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp1
667 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp1
668 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp1
|
| H A D | e3.f | 67 integer aa local
206 integer aa, b local
|
| /phasta/phSolver/common/ |
| H A D | ludcmp.f | 1 subroutine ludcmp(aa,n,np,indx,d) argument
11 dimension aa(np,np), indx(n), vv(500)
17 if (abs(aa(i,j)).gt.aamax) aamax=abs(aa(i,j))
24 sum=aa(i,j)
26 sum=sum-aa(i,k)*aa(k,j)
28 aa(i,j)=sum
32 sum=aa(i,j)
34 sum=sum-aa(i,k)*aa(k,j)
36 aa(i,j)=sum
45 dum=aa(imax,k)
[all …]
|
| H A D | fillsparse.f | 12 integer aa, b, c, e, i, k, n local
27 do aa = 1, nshl ! loop over the local equation numbers
28 i = ien(e,aa) ! finds the global equation number or
44 lhsK(1,k) = lhsK(1,k) + xKebe(e,1,aa,b)
45 lhsK(2,k) = lhsK(2,k) + xKebe(e,2,aa,b)
46 lhsK(3,k) = lhsK(3,k) + xKebe(e,3,aa,b)
47 lhsK(4,k) = lhsK(4,k) + xKebe(e,4,aa,b)
48 lhsK(5,k) = lhsK(5,k) + xKebe(e,5,aa,b)
49 lhsK(6,k) = lhsK(6,k) + xKebe(e,6,aa,b)
50 lhsK(7,k) = lhsK(7,k) + xKebe(e,7,aa,b)
[all …]
|
| H A D | fillsparse.c | 26 int e,i,j,aa; //following along with fillsparse.f in fillsparsecpetscs() local
30 for(aa=0;aa<nshl;aa++) locat[aa]=ieng[e+npro*aa]-1; in fillsparsecpetscs()
51 int e,i,j,aa; //following along with fillsparse.f in fillsparsecpetscc() local
56 for(aa=0;aa<nshl;aa++) locat[aa]=ieng[e+npro*aa]-1; in fillsparsecpetscc()
|
| H A D | lubksb.f | 1 subroutine lubksb(aa,n,np,indx,bb) argument
9 dimension indx(n),aa(np,np),bb(n)
18 sum=sum-aa(i,j)*bb(j)
28 sum=sum-aa(i,j)*bb(j)
30 bb(i)=sum/aa(i,i)
|
| /phasta/phSolver/compressible/ |
| H A D | e3bvar.f | 83 integer aa local
305 do aa = 1, nshl !TODO: get rid of the intermediary dNadxi
307 dNadxi(:,aa,1) = shglb(:,1,aa) * 1 !would normally be a sum over
308 dNadxi(:,aa,2) = shglb(:,2,aa) * 1 !all nodes, but N = 0 for a /= 1
309 dNadxi(:,aa,3) = shglb(:,3,aa) * 1
312 do aa = 1, nshl
313 dNadx(:,aa,1) = dNadxi(:,aa,1) * dxidxb(:,1,1) +
314 & dNadxi(:,aa,2) * dxidxb(:,2,1) +
315 & dNadxi(:,aa,3) * dxidxb(:,3,1)
316 dNadx(:,aa,2) = dNadxi(:,aa,1) * dxidxb(:,1,2) +
[all …]
|
| H A D | e3b.f | 94 integer aa, b local
233 do aa = 1, nshl
234 dNadn(:,aa) = dNadx(:,aa,1)*bnorm(:,1)
235 & + dNadx(:,aa,2)*bnorm(:,2)
236 & + dNadx(:,aa,3)*bnorm(:,3)
243 do aa = 1,nshl
245 EGmass(:,b, aa) = con * WdetJb * shape(:,b) * dNadn(:,aa)
|
| /phasta/phastaIO/ |
| H A D | phastaIO.cc | 1775 double * aa = ( double * )datablock_header; in writeheader() local
1776 memcpy(PhastaIOActiveFiles[i]->double_chunk, aa, DB_HEADER_SIZE); in writeheader()
1786 int * aa = ( int * )datablock_header; in writeheader() local
1787 memcpy(PhastaIOActiveFiles[i]->int_chunk, aa, DB_HEADER_SIZE); in writeheader()
|