Home
last modified time | relevance | path

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

12

/phasta/phSolver/compressible/
H A De3dc.f276 DC(:,intp)=max(zero,(fact*sqrt(raLS*gnorm))-(rtLS*gnorm))
298 DC(:,intp) = two * rtLS * gnorm
324 DC(:,intp) = min( max(zero,fact * sqrt(raLS * gnorm)-
343 ri ( :,1) = ri ( :,1) + DC(:,intp) * gAgyi( :,1)
344 rmi( :,1) = rmi( :,1) + DC(:,intp) * gAgyi( :,1)
345 ri ( :,2) = ri ( :,2) + DC(:,intp) * gAgyi( :,2)
346 rmi( :,2) = rmi( :,2) + DC(:,intp) * gAgyi( :,2)
347 ri ( :,3) = ri ( :,3) + DC(:,intp) * gAgyi( :,3)
348 rmi( :,3) = rmi( :,3) + DC(:,intp) * gAgyi( :,3)
349 ri ( :,4) = ri ( :,4) + DC(:,intp) * gAgyi( :,4)
[all …]
H A De3ql.f54 do intp = 1, ngauss
55 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
144 tmp = Qwt(lcsyst,intp)
H A Dgetdiff.f178 xmut = xmudmi(:,intp)
376 xmut(i)=xmudmi(i,intp)
379 xmut(i)=xi*xmut(i)+(one-xi)*xmudmi(1,intp)
H A De3.f106 do intp = 1, ngauss
110 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
226 if((intp.eq.1).and.(ierrcalc.eq.1).and.(nitr.eq.iter)) then
401 do intp = 1, ngauss
405 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
H A De3bvar.f182 WdetJb = (1 - Qwtb(lcsyst,intp)) / (four*temp)
190 WdetJb = Qwtb(lcsyst,intp) / temp
192 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
539 WdetJb = (1 - Qwtb(lcsyst,intp)) / (four*temp)
541 WdetJb = Qwtb(lcsyst,intp)/ temp
543 WdetJb =Qwtb(lcsyst,intp) / (four*temp)
H A De3b.f111 do intp = 1, ngaussb
115 if (Qwtb(lcsyst,intp) .eq. zero) cycle ! precaution
468 do intp = 1, ngaussb
472 if (Qwtb(lcsyst,intp) .eq. zero) cycle ! precaution
/phasta/phSolver/common/
H A Dasithf.f30 do intp = 1, ngaussf
43 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
44 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
45 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
46 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
47 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
48 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
49 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
50 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
51 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
[all …]
H A Dgetstrl.f25 do intp = 1, ngauss
37 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
38 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
39 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
40 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
41 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
42 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
43 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
44 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
45 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
[all …]
H A Dbardmc.f701 do intp = 1, ngaussf
713 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
714 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
715 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
716 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
717 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
718 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
719 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
720 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
721 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
[all …]
H A Dhierarchic.f43 shape(:,i) = shp(i,intp)
45 shdrv(:,j,i) = shgl(j,i,intp)
50 shape(:,i) = sgn(:,i) * shp(i,intp)
52 shdrv(:,j,i) = shgl(j,i,intp)*sgn(:,i)
73 shape(:,i) = shp(i,intp)
75 shdrv(:,j,i) = shgl(j,i,intp)
80 shape(:,i) = sgn(:,i) * shp(i,intp)
82 shdrv(:,j,i) = shgl(j,i,intp)*sgn(:,i)
336 do intp = 1, npts
345 yvl(:,intp,j) = yvl(:,intp,j) + shape(:,i) * ycl(:,i,j)
[all …]
H A Daveprep.f83 do intp = 1, ngauss
86 xmudmi(iel,intp) = cdelsq( lfathe(iel,intp) )*
87 & strl(iel,intp)
127 do intp = 1, ngauss
132 fresli(:,j) = fresli(:,j) + shp(i,intp)*fresl(:,i,j)
216 integer intp local
222 do intp = 1, ngauss
230 yint = yint + xl(nel,n,2)*shp(n,intp)
231 xint = xint + xl(nel,n,1)*shp(n,intp)
232 zint = zint + xl(nel,n,3)*shp(n,intp)
[all …]
H A Df3lhs.f75 do intp = 1, ngaussb
85 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
87 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
89 WdetJb = Qwtb(lcsyst,intp) / (two*temp)
91 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
93 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
95 WdetJb = Qwtb(lcsyst,intp) / (two*temp)
H A Dcmass.f47 do intp = 1, ngauss ! (these are in common.h)
52 & shape, shdrv, intp)
111 do intp = 1, ngaussf ! (these are in common.h)
116 & shape, shdrv, intp)
141 do intp = 1, ngauss ! (these are in common.h)
150 & shape, shdrv, intp)
228 do intp = 1, ngaussf ! (these are in common.h)
233 & shape, shdrv, intp)
258 do intp = 1, ngauss ! (these are in common.h)
267 & shape, shdrv, intp)
[all …]
H A Dasbnabi.f115 do intp = 1, ngaussb
123 WdetJb(:) = Qwtb(lcsyst,intp) * Wfactor(:)
271 do intp = 1, ngaussb
273 WdetJb(:) = Qwtb(lcsyst,intp) * Wfactor(:)
H A Derrorsmooth.f173 do intp = 1, ngauss
174 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
/phasta/phSolver/incompressible/
H A Dfilters.f32 do intp = 1, ngaussf ! Loop over quadrature points
44 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
45 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
46 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
47 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
48 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
49 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
50 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
51 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
52 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
[all …]
H A De3qGradV.f42 do intp = 1, ngauss
43 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
H A De3.f83 do intp = 1, ngauss
85 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
216 do intp = 1, ngauss
218 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
H A DadvLES.f84 do intp = 1, ngaussf
95 pres(:) = pres(:) + shp(n,intp) * yl(:,n,1)
96 u1(:) = u1(:) + shp(n,intp) * yl(:,n,2)
97 u2(:) = u2(:) + shp(n,intp) * yl(:,n,3)
98 u3(:) = u3(:) + shp(n,intp) * yl(:,n,4)
106 aci(:,1) = aci(:,1) + shp(n,intp) * acl(:,n,2)
107 aci(:,2) = aci(:,2) + shp(n,intp) * acl(:,n,3)
108 aci(:,3) = aci(:,3) + shp(n,intp) * acl(:,n,4)
120 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
121 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
[all …]
H A De3q.f57 do intp = 1, ngauss
58 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
243 do intp = 1, ngauss
244 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
H A De3ql.f52 do intp = 1, ngauss
90 tmp = Qwt(lcsyst,intp)
215 do intp = 1, ngauss
241 tmp = Qwt(lcsyst,intp)
H A Dhessian.f128 do intp = 1, ngauss
130 if ( Qwt( lcsyst, intp ) .eq. zero ) cycle
203 do intp = 1, ngauss
205 if ( Qwt( lcsyst, intp ) .eq. zero ) cycle
351 WdetJ = Qwt(lcsyst,intp)/ tmp
H A De3bvar.f234 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
236 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
238 WdetJb = Qwtb(lcsyst,intp) / (two*temp)
240 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
242 WdetJb = Qwtb(lcsyst,intp) / (four*temp)
244 WdetJb = Qwtb(lcsyst,intp) / (two*temp)
382 if (intp.eq.ngaussb) then ! do this only for the last gauss point
682 if (intp.eq.ngaussb) then ! do this only for the last gauss point
871 WdetJb = (1 - Qwtb(lcsyst,intp)) / (four*temp)
873 WdetJb = Qwtb(lcsyst,intp) / temp
[all …]
H A Derrsmooth.f174 do intp = 1, ngauss
175 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
H A Dgetdiff.f86 rmu = rmu + xmudmi(:,intp)
116 xmut(i)=xmudmi(i,intp)
119 xmut(i)=xi*xmut(i)+(one-xi)*(xmudmi(1,intp)+datmat(1,2,2))

12