Home
last modified time | relevance | path

Searched refs:i (Results 1 – 25 of 224) sorted by relevance

123456789

/phasta/phSolver/common/
H A Dget_a_not_hex.f10 do i = 1, nsd
12 anot(:,1,i) = pt125*(xc(:,1,i)+xc(:,2,i)+xc(:,3,i)+xc(:,4,i)
13 & +xc(:,5,i)+xc(:,6,i)+xc(:,7,i)+xc(:,8,i))
15 anot(:,2,i) = pt125*(-xc(:,1,i)+xc(:,2,i)+xc(:,3,i)-xc(:,4,i)
16 & -xc(:,5,i)+xc(:,6,i)+xc(:,7,i)-xc(:,8,i))
18 anot(:,3,i) = pt125*(-xc(:,1,i)-xc(:,2,i)+xc(:,3,i)+xc(:,4,i)
19 & -xc(:,5,i)-xc(:,6,i)+xc(:,7,i)+xc(:,8,i))
21 anot(:,4,i) = pt125*(-xc(:,1,i)-xc(:,2,i)-xc(:,3,i)-xc(:,4,i)
22 & +xc(:,5,i)+xc(:,6,i)+xc(:,7,i)+xc(:,8,i))
24 anot(:,5,i) = pt125*(xc(:,1,i)-xc(:,2,i)+xc(:,3,i)-xc(:,4,i)
[all …]
H A Dgenscale.f55 do i = 1, nelint
57 yintl(i,:,j+1) = y(ien2D(i,:),j)
59 yintl(i,:,1) = y(ien2D(i,:),4)
62 yintl(i,:,j) = y(ien2D(i,:),j)
72 do i=1,nfint
73 velbarR(i,:)=0
74 do j=1,imax(i)+1
75 call shptet(ipord,xsinfin(i,j,:),shapef(:),shgradl(:,:))
77 velbarR(i,:)=velbarR(i,:)
78 & + yintl(elcnfin(i,j),k,:)*shapef(k)
[all …]
H A Drenum.f49 do i=1,numnp
50 if ((x(i,1).eq.0).and.(x(i,2).eq.0)) then
52 nula(j) = i
54 xcyl(i,2) = atan2(x(i,2),x(i,1))
66 do i=2,npin
67 if((x(nen1(i),1).eq.0).and.(x(nen1(i),2).eq.0)) then
70 if(thmin.gt.xcyl(nen1(i),2)) thmin = xcyl(nen1(i),2)
71 if(thmax.lt.xcyl(nen1(i),2)) thmax = xcyl(nen1(i),2)
76 do i=1,j
77 xcyl(nula(i),2) = thmin
[all …]
H A Ddtn.f25 do i=1,nshg
26 itype=ifeature(i)
27 if(btest(iBC(i),13)) then
29 tmp(j)=y(i,5+j)
39 BC(i,6+j)=-tmp(j)
66 integer i local
78 do i=1,nscalar
79 tmp(i) = 0.0d0
90 do i=1,nscalar
91 tmp(i) = D(i) * ( tmp(i) - 0.0 ) / distance * units
[all …]
H A Dgenshpb.f30 do i=1,nintb(lcsyst)
31 call shpTet(ipord,Qptb(1,1:3,i),shpb(1,:,i),
32 & shglb(1,:,:,i))
39 do i=1,nintb(lcsyst)
40 call shpHex (ipord, Qptb(2,1:3,i),shpb(2,:,i),
41 & shglb(2,:,:,i))
46 do i=1,nintb(lcsyst)
47 call shp6w (ipord,Qptb(3,1:3,i),
48 & shpb(3,:,i),shglb(3,:,:,i))
53 do i=1,nintb(lcsyst)
[all …]
/phasta/phSolver/incompressible/
H A Dftools.f27 do i = 1, nshl
28 i0 = (nsd) * (i - 1)
29 Diagl(:,i,1) = xKebe(1:npro,i0+1,i0+1)
30 Diagl(:,i,2) = xKebe(1:npro,i0+2,i0+2)
31 Diagl(:,i,3) = xKebe(1:npro,i0+3,i0+3)
37 do i = 1, nshl
38 Diagl(:,i,4) = xGoC(1:npro,nn2+i,i)
60 integer i, j local
68 do i = 1, n
69 c(i,1) = a(i,1) * b(i,1)
[all …]
H A Dbc3lhs.f83 do i = 1, nshl
84 xKebe(iel,iadj1,i,inod) = xKebe(iel,iadj1,i,inod)
85 & - BC(in,4) * xKebe(iel,irem1,i,inod)
86 xKebe(iel,iadj2,i,inod) = xKebe(iel,iadj2,i,inod)
87 & - BC(in,4) * xKebe(iel,irem2,i,inod)
88 xKebe(iel,iadj3,i,inod) = xKebe(iel,iadj3,i,inod)
89 & - BC(in,4) * xKebe(iel,irem3,i,inod)
103 do i = 1, nshl
104 xKebe(iel,iadj1,i,inod) = xKebe(iel,iadj1,i,inod)
105 & - BC(in,5) * xKebe(iel,irem1,i,inod)
[all …]
H A De3q.f111 do i=1,nshl
112 ql(:,i,1 ) = ql(:,i,1 )+ shape(:,i)*WdetJ*qdi(:,1 )
113 ql(:,i,2 ) = ql(:,i,2 )+ shape(:,i)*WdetJ*qdi(:,2 )
114 ql(:,i,3 ) = ql(:,i,3 )+ shape(:,i)*WdetJ*qdi(:,3 )
116 ql(:,i,4 ) = ql(:,i,4 )+ shape(:,i)*WdetJ*qdi(:,4 )
117 ql(:,i,5 ) = ql(:,i,5 )+ shape(:,i)*WdetJ*qdi(:,5 )
118 ql(:,i,6 ) = ql(:,i,6 )+ shape(:,i)*WdetJ*qdi(:,6 )
120 ql(:,i,7 ) = ql(:,i,7 )+ shape(:,i)*WdetJ*qdi(:,7 )
121 ql(:,i,8 ) = ql(:,i,8 )+ shape(:,i)*WdetJ*qdi(:,8 )
122 ql(:,i,9 ) = ql(:,i,9 )+ shape(:,i)*WdetJ*qdi(:,9 )
[all …]
H A Dfilters.f115 do i=1,nshl
116 u1 = u1 + shp(i,intp)*yl(:,i,2)
117 u2 = u2 + shp(i,intp)*yl(:,i,3)
118 u3 = u3 + shp(i,intp)*yl(:,i,4)
120 S11 = S11 + shg(:,i,1)*yl(:,i,2)
121 S22 = S22 + shg(:,i,2)*yl(:,i,3)
122 S33 = S33 + shg(:,i,3)*yl(:,i,4)
124 S12 = S12 + shg(:,i,2)*yl(:,i,2)
125 & +shg(:,i,1)*yl(:,i,3)
126 S13 = S13 + shg(:,i,3)*yl(:,i,2)
[all …]
H A Dstats.f96 integer i local
101 do i = 1, nshg
102 res = stsVec(i,:)
109 stsReg(i) = 0
111 stsReg(i) = (res(1) + res(2) + res(3)) / 1000.
118 do i = 1, nshg
119 res = stsVec(i,:)
120 reg = stsReg(i)
133 stsMInv(i,1) = det * (r1 * r2 - r4 * r4)
134 stsMInv(i,2) = det * (r0 * r2 - r5 * r5)
[all …]
H A DadvLES.f235 do i = 1, npro
237 shpfun(i,n) = shp(n,intp)
611 i = iper(j)
612 if (i .ne. j) then
613 fres(i,:) = fres(i,:) + fres(j,:)
618 i = iper(j)
619 if (i .ne. j) then
659 do i = 1, 21
660 call sparseCG (fres(:,i), pfres(:,i), lhsG,
668 do i = 1, 21
[all …]
H A De3qGradV.f62 do i=1,nshl
63 ql(:,i,1 ) = ql(:,i,1 )+ shape(:,i)*WdetJ*g1yi(:,2 ) ! du/dx
64 ql(:,i,2 ) = ql(:,i,2 )+ shape(:,i)*WdetJ*g2yi(:,2 ) ! du/dy
65 ql(:,i,3 ) = ql(:,i,3 )+ shape(:,i)*WdetJ*g3yi(:,2 ) ! du/dz
67 ql(:,i,4 ) = ql(:,i,4 )+ shape(:,i)*WdetJ*g1yi(:,3 ) ! dv/dx
68 ql(:,i,5 ) = ql(:,i,5 )+ shape(:,i)*WdetJ*g2yi(:,3 ) ! dv/dy
69 ql(:,i,6 ) = ql(:,i,6 )+ shape(:,i)*WdetJ*g3yi(:,3 ) ! dv/dz
71 ql(:,i,7 ) = ql(:,i,7 )+ shape(:,i)*WdetJ*g1yi(:,4 ) ! dw/dx
72 ql(:,i,8 ) = ql(:,i,8 )+ shape(:,i)*WdetJ*g2yi(:,4 ) ! dw/dy
73 ql(:,i,9 ) = ql(:,i,9 )+ shape(:,i)*WdetJ*g3yi(:,4 ) ! dw/dz
/phasta/phSolver/compressible/
H A De3wmlt.f74 do i = 1, nshl
75 rl(:,i,1) = rl(:,i,1) + WdetJ * (
76 & shg(:,i,1) * ri(:, 1) + shg(:,i,2) * ri(:, 6)
77 & + shg(:,i,3) * ri(:,11) )
78 rl(:,i,2) = rl(:,i,2) + WdetJ * (
79 & shg(:,i,1) * ri(:, 2) + shg(:,i,2) * ri(:, 7)
80 & + shg(:,i,3) * ri(:,12) )
81 rl(:,i,3) = rl(:,i,3) + WdetJ * (
82 & shg(:,i,1) * ri(:, 3) + shg(:,i,2) * ri(:, 8)
83 & + shg(:,i,3) * ri(:,13) )
[all …]
H A Dbc3lhs.f57 do i = 1, nedof
58 EGmass(iel,i,i1) = zero
59 EGmass(iel,i1,i) = zero
74 do i = 1, nedof
75 EGmass(iel,i3,i) = EGmass(iel,i3,i)
76 & - BC(in,4) * EGmass(iel,i2,i)
77 EGmass(iel,i4,i) = EGmass(iel,i4,i)
78 & - BC(in,5) * EGmass(iel,i2,i)
80 do i = 1, nedof
81 EGmass(iel,i,i3) = EGmass(iel,i,i3)
[all …]
H A De3juel.f84 do i = 1, nshl
85 rml(:,i,1) = rml(:,i,1) + fddt * (yl(:,i,1)+ub(:,1))
86 rml(:,i,2) = rml(:,i,2) + fddt * (yl(:,i,2)+ub(:,2))
87 rml(:,i,3) = rml(:,i,3) + fddt * (yl(:,i,3)+ub(:,3))
88 rml(:,i,4) = rml(:,i,4) + fddt * (yl(:,i,4)+ub(:,4))
89 rml(:,i,5) = rml(:,i,5) + fddt * (yl(:,i,5)+ub(:,5))
103 do i = 1, nshl
104 yl(:,i,1) = fact*(acl(:,i,1)+ub(:,1))
105 yl(:,i,2) = fact*(acl(:,i,2)+ub(:,2))
106 yl(:,i,3) = fact*(acl(:,i,3)+ub(:,3))
[all …]
H A De3q.f153 do i=1,nshl
154 ql(:,i,1 ) = ql(:,i,1 )+ shape(:,i)*WdetJ*qdi(:,1 )
155 ql(:,i,2 ) = ql(:,i,2 )+ shape(:,i)*WdetJ*qdi(:,2 )
156 ql(:,i,3 ) = ql(:,i,3 )+ shape(:,i)*WdetJ*qdi(:,3 )
157 ql(:,i,4 ) = ql(:,i,4 )+ shape(:,i)*WdetJ*qdi(:,4 )
158 ql(:,i,5 ) = ql(:,i,5 )+ shape(:,i)*WdetJ*qdi(:,5 )
159 ql(:,i,6 ) = ql(:,i,6 )+ shape(:,i)*WdetJ*qdi(:,6 )
160 ql(:,i,7 ) = ql(:,i,7 )+ shape(:,i)*WdetJ*qdi(:,7 )
161 ql(:,i,8 ) = ql(:,i,8 )+ shape(:,i)*WdetJ*qdi(:,8 )
162 ql(:,i,9 ) = ql(:,i,9 )+ shape(:,i)*WdetJ*qdi(:,9 )
[all …]
H A Dasaugmr.f32 do i = 1, nflow*nshl, nflow
34 ubBgl(:,i ) = ubBgl(:,i )
35 & + EGmass(:,i ,j ) * uBrgl(:,j )
36 & + EGmass(:,i ,j+1) * uBrgl(:,j+1)
37 & + EGmass(:,i ,j+2) * uBrgl(:,j+2)
38 & + EGmass(:,i ,j+3) * uBrgl(:,j+3)
39 & + EGmass(:,i ,j+4) * uBrgl(:,j+4)
41 ubBgl(:,i+1) = ubBgl(:,i+1)
42 & + EGmass(:,i+1,j ) * uBrgl(:,j )
43 & + EGmass(:,i+1,j+1) * uBrgl(:,j+1)
[all …]
/phasta/svLS/
H A DINCOMMU.f50 INTEGER i, j, k, s, e, ierr, nTasks, tF, stat(MPI_STATUS_SIZE) local
61 i = tF - 1
62 i = cS(i)%ptr + cS(i)%n - 1
63 ALLOCATE(rTmp(dof,i))
66 DO i=1, nTasks
67 IF (cS(i)%tag .NE. 0) THEN
68 s = cS(i)%ptr
69 e = s + cS(i)%n - 1
70 IF (i .LT. tF) THEN
71 CALL MPI_IRECV(rTmp(:,s:e), cS(i)%n*dof, mpreal, i-1,
[all …]
H A DSPARMUL.f53 INTEGER i, j, l, col, s, e local
58 DO i=1, nNo
59 DO j=rowPtr(1,i), rowPtr(2,i)
60 KU(1,i) = KU(1,i) + K(1,j)*U(1,colPtr(j))
64 DO i=1, nNo
65 DO j=rowPtr(1,i), rowPtr(2,i)
67 KU(1,i) = KU(1,i) + K(1,j)*U(1,col) + K(2,j)*U(2,col)
68 KU(2,i) = KU(2,i) + K(3,j)*U(1,col) + K(4,j)*U(2,col)
72 DO i=1, nNo
73 DO j=rowPtr(1,i), rowPtr(2,i)
[all …]
H A DGMRES.f55 INTEGER i, j, k, l local
101 DO i=1, ls%sD
104 2 Val, u(:,:,i), u(:,:,i+1))
107 2 face, u(:,:,i), u(:,:,i+1))
110 unCondU = u(:,:,i+1)
112 2 commu, face, unCondU, u(:,:,i+1))
114 DO j=1, i
115 h(j,i) = DOTV(dof, mynNo, commu, u(:,:,i+1), u(:,:,j))
116 u(:,:,i+1) = u(:,:,i+1) - h(j,i)*u(:,:,j)
118 h(i+1,i) = NORMV(dof, mynNo, commu, u(:,:,i+1))
[all …]
H A DNORM.f49 INTEGER i, ierr local
55 DO i=1, nNo
56 NORMV = NORMV + U(1,i)*U(1,i)
59 DO i=1, nNo
60 NORMV = NORMV + U(1,i)*U(1,i) + U(2,i)*U(2,i)
63 DO i=1, nNo
64 NORMV = NORMV+ U(1,i)*U(1,i) + U(2,i)*U(2,i) + U(3,i)*U(3,i)
67 DO i=1, nNo
68 NORMV = NORMV + U(1,i)*U(1,i) + U(2,i)*U(2,i)
69 2 + U(3,i)*U(3,i) + U(4,i)*U(4,i)
[all …]
H A DDOT.f49 INTEGER i, ierr local
55 DO i=1, nNo
56 DOTV = DOTV + U(1,i)*V(1,i)
59 DO i=1, nNo
60 DOTV = DOTV + U(1,i)*V(1,i) + U(2,i)*V(2,i)
63 DO i=1, nNo
64 DOTV = DOTV + U(1,i)*V(1,i) + U(2,i)*V(2,i) +U(3,i)*V(3,i)
67 DO i=1, nNo
68 DOTV = DOTV + U(1,i)*V(1,i) + U(2,i)*V(2,i)
69 2 + U(3,i)*V(3,i) + U(4,i)*V(4,i)
[all …]
/phasta/converterIO/
H A DconverterN2O.cc45 int i, j, k, N_restart_integer, N_restart_double; in main() local
144 for ( i = 0; i < N_restart_double; i++ ) in main()
146 WriteLockD[i]=0; in main()
147 Dfield[i] = new double*[nppp]; in main()
149 paraD[i] = new int*[nppp]; in main()
151 fieldNameD[i] = new char[128]; in main()
152 fileTypeD[i] = new char[128]; in main()
153 dataTypeD[i] = new char[128]; in main()
154 headerTypeD[i] = new char[128]; in main()
158 for ( i = 0; i < N_restart_integer; i++ ) in main()
[all …]
H A DconverterO2N.cc71 int i, j, k, N_geombc_double, N_geombc_integer, N_restart_double; in main() local
251 for ( i = 0; i < nppp; i++ ) in main()
253 Dfield[i] = new double*[N_geombc_double]; // Space for the datablock of double format in main()
254 … paraD[i] = new int*[N_geombc_double]; // Integers in the header of each field of double format in main()
255 Ifield[i] = new int*[N_geombc_integer]; // Space for the datablock of integer format in main()
256 …paraI[i] = new int*[N_geombc_integer]; // Integers in the header of each field of integer format in main()
275 for ( i = 0; i < N_geombc_double; i++ ) in main()
277 WriteLockD[i]=0; in main()
279 fieldNameD[i] = new char[128]; in main()
280 fileTypeD[i] = new char[128]; in main()
[all …]
/phasta/phSolver/AMG/
H A Dramg_coarse.f38 integer :: i,j,k,m,n,p,q local
85 do i=2,numtask
86 subnei(i)=ilwork(itkbeg+3)+1 ! iother+1
87 subcfrev(subnei(i)) = i
99 do i=1,numtask
100 subcfrev(i) = i
127 do i = 1,amg_nshg(level1)
128 p = iabs(amg_paramap(level1)%p(i))
131 if (amg_F(i).eq.1) then
135 do i = 1,numtask
[all …]

123456789