| /phasta/phSolver/common/ |
| H A D | genshp.f | 23 lcsyst = lcblk(3,iblk) 28 select case ( lcsyst ) 30 maxnint=max(maxnint,nint(lcsyst)) 31 do i=1,nint(lcsyst) 34 shgl(1,:,1:nshl,1:nint(lcsyst)) = 35 & shgl(1,:,1:nshl,1:nint(lcsyst))/two 39 maxnint=max(maxnint,nint(lcsyst)) 40 do i=1,nint(lcsyst) 47 maxnint=max(maxnint,nint(lcsyst)) 48 do i=1,nint(lcsyst) [all …]
|
| H A D | asbnabi.f | 56 if(lcsyst.eq.1) then ! set to curl into element all others out 59 elseif(lcsyst.eq.2) then 62 elseif(lcsyst.eq.3) then 65 elseif(lcsyst.eq.4) then 68 elseif(lcsyst.eq.5) then 71 elseif(lcsyst.eq.6) then 91 if (lcsyst .eq. 1) then 93 elseif (lcsyst .eq. 2) then 95 elseif (lcsyst .eq. 3) then 97 elseif (lcsyst .eq. 4) then [all …]
|
| H A D | genshpb.f | 23 lcsyst = lcblkb(3,iblk) 27 select case ( lcsyst ) 30 do i=1,nintb(lcsyst) 34 shglb(1,:,1:nshl,1:nintb(lcsyst)) = 35 & shglb(1,:,1:nshl,1:nintb(lcsyst))/two 39 do i=1,nintb(lcsyst) 46 do i=1,nintb(lcsyst) 53 do i=1,nintb(lcsyst) 59 do i=1,nintb(lcsyst) 66 do i=1,nintb(lcsyst)
|
| H A D | f3lhs.f | 60 if (lcsyst .eq. 1) then 84 if (lcsyst .eq. 1) then 85 WdetJb = Qwtb(lcsyst,intp) / (four*temp) 86 elseif (lcsyst .eq. 2) then 87 WdetJb = Qwtb(lcsyst,intp) / (four*temp) 88 elseif (lcsyst .eq. 3) then 89 WdetJb = Qwtb(lcsyst,intp) / (two*temp) 90 elseif (lcsyst .eq. 4) then 91 WdetJb = Qwtb(lcsyst,intp) / (four*temp) 92 elseif (lcsyst .eq. 5) then [all …]
|
| H A D | pvsqbi.f | 78 lcsyst = lcblkb(3,iblk) 89 if(lcsyst.eq.3) lcsyst=nenbl 91 if(lcsyst.eq.3 .or. lcsyst.eq.4) then 92 ngaussb = nintb(lcsyst) 94 ngaussb = nintb(lcsyst) 103 tmpshpb(1:nshl,:) = shpb(lcsyst,1:nshl,:)
|
| H A D | getdmc.f | 72 lcsyst = lcblk(3,iblk) 80 ngauss = nint(lcsyst) 81 ngaussf = nintf(lcsyst) 84 & shglf(lcsyst,:,1:nshl,:), 85 & shpf(lcsyst,1:nshl,:),Qwtf(lcsyst,1:ngaussf)) 91 lcsyst = lcblk(3,iblk) 99 ngauss = nint(lcsyst) 100 ngaussf = nintf(lcsyst) 104 & strl(iel:inum,:), shgl(lcsyst,:,1:nshl,:), 105 & shp(lcsyst,1:nshl,:)) [all …]
|
| H A D | aveprep.f | 42 lcsyst = lcblk(3,iblk) 50 ngauss = nint(lcsyst) 52 call getylist( ylist, ifathe(iel:inum,:), shp(lcsyst,1:nshl,:), 318 lcsyst = lcblk(3,iblk) 326 ngauss = nint(lcsyst) 327 ngaussf = nintf(lcsyst) 330 & shglf(lcsyst,:,1:nshl,:), 331 & shpf(lcsyst,1:nshl,:),Qwtf(lcsyst,1:ngaussf)) 369 lcsyst = lcblk(3,iblk) 375 ngauss = nint(lcsyst) [all …]
|
| H A D | filtprep.f | 141 lcsyst = lcblk(3,iblk) 146 if (lcsyst .eq. 1) then ! tets 159 else if (lcsyst .eq. 2) then ! hexes 169 else if (lcsyst .eq. 3) then ! wedges 176 else if (lcsyst .eq. 5) then ! pyramids
|
| H A D | hierarchic.f | 102 if (lcsyst .eq. 1) then 121 else if(lcsyst .eq. 2) then 144 else if(lcsyst .eq. 3) then 160 else if(lcsyst .eq. 4) then 181 else if(lcsyst .eq. 5) then 202 else if(lcsyst .eq. 6) then 224 & ,lcsyst 266 lcsyst = lcblk(3,iblk)
|
| H A D | errorsmooth.f | 36 lcsyst = lcblk(3,iblk) 44 ngauss = nint(lcsyst) 52 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 53 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 174 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
|
| H A D | bardmc.f | 114 lcsyst = lcblk(3,iblk) 122 ngauss = nint(lcsyst) 123 ngaussf = nintf(lcsyst) 126 & shglf(lcsyst,:,1:nshl,:), 127 & shpf(lcsyst,1:nshl,:),Qwtf(lcsyst,1:ngaussf)) 238 lcsyst = lcblk(3,iblk) 246 ngauss = nint(lcsyst) 247 ngaussf = nintf(lcsyst) 250 & fres, hfres, shgl(lcsyst,:,1:nshl,:), 251 & shp(lcsyst,1:nshl,:),Qwtf(lcsyst,1:ngaussf)) [all …]
|
| /phasta/phSolver/incompressible/ |
| H A D | elmgmr.f | 87 lcsyst = lcblk(3,iblk) 95 ngauss = nint(lcsyst) 103 & shp(lcsyst,1:nshl,:), 104 & shgl(lcsyst,:,1:nshl,:), 109 & shp(lcsyst,1:nshl,:), 110 & shgl(lcsyst,:,1:nshl,:), 145 lcsyst = lcblk(3,iblk) 154 ngauss = nint(lcsyst) 166 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 167 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) [all …]
|
| H A D | bflux.f | 70 lcsyst = lcblkb(3,iblk) 112 lcsyst = lcblk(3,iblk) 114 ngauss = nint(lcsyst) 142 & shp(lcsyst,1:nshl,:), 143 & shgl(lcsyst,:,1:nshl,:), 166 lcsyst = lcblkb(3,iblk) 173 if(lcsyst.eq.3) lcsyst=nenbl 175 if(lcsyst.eq.3 .or. lcsyst.eq.4) then 176 ngaussb = nintb(lcsyst) 178 ngaussb = nintb(lcsyst) [all …]
|
| H A D | advLES.f | 567 lcsyst = lcblk(3,iblk) 575 ngauss = nint(lcsyst) 576 ngaussf = nintf(lcsyst) 579 & shglf(lcsyst,:,1:nshl,:), 580 & shpf(lcsyst,1:nshl,:),Qwtf(lcsyst,1:ngaussf)) 588 & shgl(lcsyst,:,1:nshl,:), shp(lcsyst,1:nshl,:), 589 & shglf(lcsyst,:,1:nshl,:), shpf(lcsyst,1:nshl,:), em, 590 & Qwtf(lcsyst,1:ngaussf)) 675 lcsyst = lcblk(3,iblk) 683 ngauss = nint(lcsyst) [all …]
|
| H A D | hessian.f | 30 lcsyst = lcblk(3,iblk) 38 ngauss = nint(lcsyst) 42 & shp(lcsyst,1:nshl,:), 43 & shgl(lcsyst,:,1:nshl,:), 59 lcsyst = lcblk(3,iblk) 67 ngauss = nint(lcsyst) 72 & shp(lcsyst,1:nshl,:), 73 & shgl(lcsyst,:,1:nshl,:), 130 if ( Qwt( lcsyst, intp ) .eq. zero ) cycle 205 if ( Qwt( lcsyst, intp ) .eq. zero ) cycle [all …]
|
| H A D | e3bvar.f | 197 if(lcsyst.eq.1) then ! set to curl into element all others out 200 elseif(lcsyst.eq.2) then 203 elseif(lcsyst.eq.3) then 206 elseif(lcsyst.eq.4) then 209 elseif(lcsyst.eq.5) then 212 elseif(lcsyst.eq.6) then 233 if (lcsyst .eq. 1) then 234 WdetJb = Qwtb(lcsyst,intp) / (four*temp) 235 elseif (lcsyst .eq. 2) then 236 WdetJb = Qwtb(lcsyst,intp) / (four*temp) [all …]
|
| H A D | errsmooth.f | 37 lcsyst = lcblk(3,iblk) 45 ngauss = nint(lcsyst) 53 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 54 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 175 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
|
| H A D | genlmass.f | 28 lcsyst = lcblk(3,iblk) 35 ngauss = nint(lcsyst) 42 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 43 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 110 if (Qwt(lcsyst,intp) .eq. zero) cycle ! precaution
|
| H A D | e3b.f | 92 if(lcsyst.eq.3.or.lcsyst.eq.4) then 93 ngaussb = nintb(lcsyst) 95 ngaussb = nintb(lcsyst) 398 if(lcsyst.eq.3.or.lcsyst.eq.4) then 399 ngaussb = nintb(lcsyst) 401 ngaussb = nintb(lcsyst)
|
| /phasta/phSolver/compressible/ |
| H A D | elmgmr.f | 74 lcsyst = lcblk(3,iblk) 82 ngauss = nint(lcsyst) 90 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 91 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 134 lcsyst = lcblk(3,iblk) 143 ngauss = nint(lcsyst) 151 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 152 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 185 lcsyst = lcblkb(3,iblk) 194 if(lcsyst.eq.3) lcsyst=nenbl [all …]
|
| H A D | elmgmrpetsc.f | 88 lcsyst = lcblk(3,iblk) 96 ngauss = nint(lcsyst) 104 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 105 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 149 lcsyst = lcblk(3,iblk) 158 ngauss = nint(lcsyst) 172 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 173 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 221 lcsyst = lcblkb(3,iblk) 230 if(lcsyst.eq.3) lcsyst=nenbl [all …]
|
| H A D | elmmfg.f | 73 lcsyst = lcblk(3,iblk) 81 ngauss = nint(lcsyst) 89 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 90 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 132 lcsyst = lcblk(3,iblk) 141 ngauss = nint(lcsyst) 148 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 149 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 178 lcsyst = lcblkb(3,iblk) 187 if(lcsyst.eq.3) lcsyst=nenbl [all …]
|
| H A D | itrres.f | 59 lcsyst = lcblk(3,iblk) 67 ngauss = nint(lcsyst) 75 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 76 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:) 112 lcsyst = lcblkb(3,iblk) 121 if(lcsyst.eq.3) lcsyst=nenbl 122 ngaussb = nintb(lcsyst) 127 tmpshpb(1:nshl,:) = shpb(lcsyst,1:nshl,:) 128 tmpshglb(:,1:nshl,:) = shglb(lcsyst,:,1:nshl,:)
|
| H A D | bflux.f | 202 lcsyst = lcblk(3,iblk) 210 ngauss = nint(lcsyst) 243 lcsyst = lcblkb(3,iblk) 252 if(lcsyst.eq.3) lcsyst=nenbl 253 ngaussb = nintb(lcsyst) 260 & shpb(lcsyst,1:nshl,:), 261 & shglb(lcsyst,:,1:nshl,:),
|
| H A D | e3bvar.f | 157 if (lcsyst .eq. 4) then 165 elseif (lcsyst .eq. 1) then 181 if (lcsyst .eq. 3) then 182 WdetJb = (1 - Qwtb(lcsyst,intp)) / (four*temp) 183 elseif (lcsyst .eq. 4) then 190 WdetJb = Qwtb(lcsyst,intp) / temp 192 WdetJb = Qwtb(lcsyst,intp) / (four*temp) 514 if (lcsyst .eq. 4) then 522 elseif (lcsyst .eq. 1) then 538 if (lcsyst .eq. 3) then [all …]
|