| /phasta/phSolver/incompressible/ |
| H A D | lesSparse.f | 205 & p, q, nNodes,
213 real*8 pLhs(4,nnz_tot), p(nNodes), q(nNodes,3)
221 q(i,1) = 0
222 q(i,2) = 0
223 q(i,3) = 0
235 q(j,1) = q(j,1) - pLhs(1,k) * pisave
236 q(j,2) = q(j,2) - pLhs(2,k) * pisave
237 q(j,3) = q(j,3) - pLhs(3,k) * pisave
253 1 p, q, nNodes,
264 real*8 p(nNodes,4), q(nNodes,3)
[all …]
|
| H A D | filters.f | 1446 & q(nshg), ilwork(nlwork),
1513 & pp, q, nshg,
1516 call commIn(q, ilwork, 1,iper,iBC,BC) ! masters=masters+slaves
1573 pap = pap + p(n) * q(n)
1591 rhs(:) = rhs(:) - alpha * q(:)
|
| /phasta/phSolver/compressible/ |
| H A D | sparseap.f | 21 real*8 p(nshg,nflow), q(nshg,nflow) local
44 q=zero
87 q(i,1) = q(i,1) + tmp1
88 q(i,2) = q(i,2) + tmp2
89 q(i,3) = q(i,3) + tmp3
90 q(i,4) = q(i,4) + tmp4
91 q(i,5) = q(i,5) + tmp5
97 p = q
|
| H A D | bflux.f | 56 & q(nshg,ndof), qq(nshg,8),
70 q(:,1:ndof)=y(:,1:ndof)
72 if (irs .ge. 1) call restar ('out ', q,ac)
80 q(:,1) = ro * q(:,1)
81 q(:,2) = vel * q(:,2)
82 q(:,3) = vel * q(:,3)
84 & q(:,4) = vel * q(:,4)
85 q(:,nflow) = temper * q(:,nflow)
136 write (iout, 2000) n, (q(n,i), i=1,ndof)
141 write (iout) q
[all …]
|
| /phasta/phSolver/common/ |
| H A D | restar.f | 1 subroutine restar (code, q, ac) argument
28 dimension q(nshg,ndof),ac(nshg,ndof)
39 q(:,1:3)=qold(:,2:4)
40 q(:,4)=qold(:,1)
41 if(ndof.gt.4) q(:,5:ndof)=qold(:,5:ndof)
74 qold(:,2:4) = q(:,1:3)
75 qold(:,1) = q(:,4)
76 if(ndof.gt.4) qold(:,5:ndof) = q(:,5:ndof)
|
| H A D | rwvelb.f | 1 subroutine rwvelb (code, q, ifail) argument
15 dimension q(nfath,nflow)
56 read (irstin,err=877,end=877) q
87 write (irstou) q
|
| /phasta/phSolver/AMG/ |
| H A D | ramg_ITAI.f | 14 integer :: i,j,k,p,q,mnnz local
24 q=b_rowp(p)
25 if (itemp(q).ne.i) then
26 itemp(q)=i
46 integer :: i,j,k,p,q,mnnz local
60 q=b_rowp(p)
61 if (itemp(q).ne.i) then
62 itemp(q)=i
64 c_rowp(mnnz) = q
87 integer :: i,j,k,p,q local
[all …]
|
| H A D | ramg_extract.f | 48 integer :: i,j,k,m,n,p,q,ki,kj,ni,nj,qq
288 q = lhsGProwp%p(kj)
289 if (p.eq.q) then
291 k = q
|
| H A D | ramg_smoother.f | 116 integer i,j,k,p,q local
145 q = amg_paramap(clevel)%p(i)
148 if ((p.gt.0).or.(p.ne.q)) then ! master boundary or extended
|
| H A D | ramg_paratools.f | 486 subroutine ramg_PPEAp(q,p, argument
494 real(kind=8),intent(inout),dimension(nshg) :: q local
532 call fMtxVdimVecMult(t1a,amg_ppeDiag(1)%p,q,1,1,1,1,nshg)
542 subroutine ramg_PPEAps(q,p,diag, argument
550 real(kind=8),intent(inout),dimension(nshg) :: q local
585 call fMtxVdimVecMult(t1a,amg_ppeDiag(1)%p,q,1,1,1,1,nshg)
|
| H A D | ramg_mls.f | 59 integer :: p,q,nsize,step,asize,gsize local
239 integer :: p,q,nsize,step,asize,gsize local
|
| H A D | ramg_ggb.f | 42 integer :: p,q,nsize,step,asize,gsize local
|
| H A D | ramg_coarse.f | 38 integer :: i,j,k,m,n,p,q local
|
| /phasta/phSolver/testIC-DSIM/ |
| H A D | geom.sms | 160 1 24+d+p+q-f
162 1 24-q+s-j-g
173 0 34-q-x-l-u
|
| H A D | geom.sdm | 159 1 24+d+p+q-f
161 1 24-q+s-j-g
|