| /phasta/M2NFixBnd/src/ |
| H A D | cname.f | 8 ii = i 10 i0 = mod (ii,10) 11 ii = (ii - i0) / 10 12 i1 = mod (ii,10) 13 ii = (ii - i1) / 10 14 i2 = mod (ii,10) 15 ii = (ii - i2) / 10 16 i3 = mod (ii,10) 48 ii = i 50 il(0) = mod(ii,10) [all …]
|
| /phasta/M2N/src/ |
| H A D | cname.f | 8 ii = i 10 i0 = mod (ii,10) 11 ii = (ii - i0) / 10 12 i1 = mod (ii,10) 13 ii = (ii - i1) / 10 14 i2 = mod (ii,10) 15 ii = (ii - i2) / 10 16 i3 = mod (ii,10) 48 ii = i 50 il(0) = mod(ii,10) [all …]
|
| /phasta/phSolver/common/ |
| H A D | cname.f | 8 ii = i 10 i0 = mod (ii,10) 11 ii = (ii - i0) / 10 12 i1 = mod (ii,10) 13 ii = (ii - i1) / 10 14 i2 = mod (ii,10) 15 ii = (ii - i2) / 10 16 i3 = mod (ii,10) 48 ii = i 50 il(0) = mod(ii,10) [all …]
|
| H A D | lubksb.f | 11 ii=0 16 if (ii.ne.0)then 17 do j=ii,i-1 21 ii=i
|
| H A D | slpwBC.f | 34 ii=nmax(wn1(1),wn1(2),wn1(3)) 35 if(ii.eq.1)then ! u has the largest constrain, u should be constrained 40 elseif(ii.eq.2)then 45 elseif(ii.eq.3)then 210 ii=nmax(wn1(1),wn1(2),wn1(3)) 211 if(ii.eq.1)then ! u has the largest constrain 213 elseif(ii.eq.2)then ! v has the largest constrain 215 elseif(ii.eq.3)then ! w has the largest constrain
|
| /phasta/AcuStat/src/ |
| H A D | cname.f | 8 ii = i 10 i0 = mod (ii,10) 11 ii = (ii - i0) / 10 12 i1 = mod (ii,10) 13 ii = (ii - i1) / 10 14 i2 = mod (ii,10) 15 ii = (ii - i2) / 10 16 i3 = mod (ii,10) 48 ii = i 50 il(0) = mod(ii,10) [all …]
|
| /phasta/shapeFunction/src/ |
| H A D | parDrv.c | 58 int ind[2],ii; in E_parDrv() local 92 for(ii=0; ii<2; ii++) { in E_parDrv() 93 if(ind[ii] == 0) /* r'= r */ in E_parDrv() 94 drv[0][ii] = 1.0; in E_parDrv() 95 else if(ind[ii] == 1) /* r'= s */ in E_parDrv() 96 drv[1][ii] = 1.0; in E_parDrv() 97 else if(ind[ii] == 2) /* r'= t */ in E_parDrv() 98 drv[2][ii] = 1.0; in E_parDrv() 99 else if(ind[ii] == 3) /* r'= 1-r-s-t */ in E_parDrv() 100 drv[0][ii]=drv[1][ii]=drv[2][ii]=-1.0; in E_parDrv()
|
| /phasta/phSolver/incompressible/ |
| H A D | e3source.f | 482 do ii=1,npro 484 kay(ii) = kay(ii) + shape_funct(ii,jj) * yl(ii,jj,6) 485 epsilon(ii) = epsilon(ii) 486 & + shape_funct(ii,jj) * yl(ii,jj,7) 525 do ii=1,npro 529 sclr_ls(ii) = sclr_ls(ii) 530 & + shape_funct(ii,jj) * yl(ii,jj,6) 534 if (sclr_ls(ii) .gt. epsilon_ls)then 536 sign_levelset(ii) = one 538 elseif (abs(sclr_ls(ii)) .le. epsilon_ls)then [all …]
|
| H A D | soldir.f | 75 do ii=1,4 76 write(543,20) (globMas(i0+ii,j0+kk), kk=1,4)
|
| /phasta/phSolver/compressible/ |
| H A D | timer.f | 41 ii = i 43 if (ii .eq. icd + 2) call error ('timer ', tcode, 0) 50 if (ii .eq. icd) then 72 if (ii .eq. icd+2) then 83 if (ii .ne. icd+1) then 86 icode = ii
|
| H A D | e3b.f | 500 do ii=1,npro 504 sclr_ls(ii) = sclr_ls(ii) 505 & + shape(ii,jj) * ycl(ii,jj,6) 508 if (sclr_ls(ii) .lt. - epsilon_ls)then 510 sign_levelset(ii) = - one 512 elseif (abs(sclr_ls(ii)) .le. epsilon_ls)then 514 sign_levelset(ii) =sclr_ls(ii)/epsilon_ls 515 & + sin(pi*sclr_ls(ii)/epsilon_ls)/pi 517 elseif (sclr_ls(ii) .gt. epsilon_ls) then 519 sign_levelset(ii) = one [all …]
|
| H A D | e3source.f | 437 do ii=1,npro 441 sclr_ls(ii) = sclr_ls(ii) 442 & + shape(ii,jj) * ycl(ii,jj,6) 446 if (sclr_ls(ii) .lt. - epsilon_ls)then 448 sign_levelset(ii) = - one 450 elseif (abs(sclr_ls(ii)) .le. epsilon_ls)then 453 sign_levelset(ii) =sclr_ls(ii)/epsilon_ls 454 & + sin(pi*sclr_ls(ii)/epsilon_ls)/pi 457 elseif (sclr_ls(ii) .gt. epsilon_ls) then 459 sign_levelset(ii) = one [all …]
|
| H A D | e3tau.f | 914 do ii = 1, npro 915 if (Pe(ii).lt.3.0) then 916 rlamtmp(ii,i,1) = rlamtmp(ii,i,1)* 917 & 0.00054*Pe(ii)**2 920 if ((Pe(ii).ge.3).and.(Pe(ii).lt.17.20)) then 921 rlamtmp(ii,i,1) = rlamtmp(ii,i,1)*(0.0114*Pe(ii) 925 if (Pe(ii).ge.17.20) then 926 rlamtmp(ii,i,1) = rlamtmp(ii,i,1)*(1.0/6.0)
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fRHSonly | 444 do ii=1,npro 446 kay(ii) = kay(ii) + shape_funct(ii,jj) * yl(ii,jj,6) 447 epsilon(ii) = epsilon(ii) 448 & + shape_funct(ii,jj) * yl(ii,jj,7) 487 do ii=1,npro 489 do jj = 1, nshl ! first find the value of levelset at point ii 491 sclr_ls(ii) = sclr_ls(ii) 492 & + shape_funct(ii,jj) * yl(ii,jj,6) 496 if (sclr_ls(ii) .gt. epsilon_ls)then 498 sign_levelset(ii) = one [all …]
|
| H A D | e3source.fimplicit | 458 do ii=1,npro 460 kay(ii) = kay(ii) + shape_funct(ii,jj) * yl(ii,jj,6) 461 epsilon(ii) = epsilon(ii) 462 & + shape_funct(ii,jj) * yl(ii,jj,7) 501 do ii=1,npro 503 do jj = 1, nshl ! first find the value of levelset at point ii 505 sclr_ls(ii) = sclr_ls(ii) 506 & + shape_funct(ii,jj) * yl(ii,jj,6) 510 if (sclr_ls(ii) .gt. epsilon_ls)then 512 sign_levelset(ii) = one [all …]
|
| /phasta/phSolver/AMG/ |
| H A D | ramg_tools.f | 841 integer i,j,ii local 850 ii = ncorp_map(i) 852 ii = i 854 write(264,mformat)ii,(rarray(i,j),j=1,rn2) 946 integer ii,jj local 964 ii = ncorp_map(i)!ipmap(i)) 969 write(264,mformat)ii,jj,(alhs(k,p),k=1,redun)
|
| /phasta/converterIO/ |
| H A D | README.md | 32 ii) mod(NP,@)=0
|