| /phasta/phSolver/incompressible/ |
| H A D | bc3global.f | 36 globMas(i0+1,i0+1) = one
47 globMas(i0+2,i0+2) = one
55 globMas(i0+3,i0+3) = one
64 globMas(i0+2,i0+2) = one
67 globMas(i0+3,i0+3) = one
75 globMas(i0+4,i0+4) = one
84 globMas(i0+2,i0+2) = one
87 globMas(i0+4,i0+4) = one
96 globMas(i0+3,i0+3) = one
99 globMas(i0+4,i0+4) = one
[all …]
|
| H A D | itrPC.f | 18 almi = one
19 alfi = one
20 gami = one
23 almi = (three-rhoinf(itseq))/(one+rhoinf(itseq))/two
24 alfi = one/(one+rhoinf(itseq))
44 y=y-(one-alfi)*Delt(1)*acold
52 Dtgl = one / Delt(itseq) ! caution: inverse of time step
76 fct = (gami-one)/gami
91 y = yold + (one - gami)/Dtgl * acold
101 fct1 = alfi/(one-alfi)
[all …]
|
| H A D | e3stab.f | 62 rhoinv=one/rho
98 dtsi=one/dts
101 tauM = one/fact
130 tauC=tauM*(one+tauM*rmu*rmu)
131 tauC=one/tauC
158 dt=one/Dtgl
181 tauM=one/fact ! turn it right side up.
190 tauC=tauM*(one+tauM*rmu*rmu)
191 tauC=one/tauC
217 dt=one/Dtgl
[all …]
|
| H A D | itrdrv.f | 435 Dtgl = one / Delt(1)
1258 POnly(:)=one ! one to get area
1876 tfact=one/istepsinybar
1887 ybar(:,1) = tfact*yold(:,1) + (one-tfact)*ybar(:,1)
1888 ybar(:,2) = tfact*yold(:,2) + (one-tfact)*ybar(:,2)
1889 ybar(:,3) = tfact*yold(:,3) + (one-tfact)*ybar(:,3)
1890 ybar(:,4) = tfact*yold(:,4) + (one-tfact)*ybar(:,4)
1892 & yold(:,3)**2) + (one-tfact)*ybar(:,5)
1894 & (one-tfact)*ybar(:,6)
1896 & (one-tfact)*ybar(:,7)
[all …]
|
| H A D | timeseries.f | 145 if ((abs(zi0(e,1)).lt.(one+tolpt)).and.
146 & (abs(zi0(e,2)).lt.(one+tolpt)).and.
147 & (abs(zi0(e,3)).lt.(one+tolpt))) then ! got the element
194 if ((zi0(e,1)+zi0(e,2)+zi0(e,3)).lt.(one+tolpt).and.
195 …& zi0(e,1).lt.(one +tolpt).and. !should not be necessary; the limit in the tet …
197 & zi0(e,2).lt.(one +tolpt).and. !should not be necessary
199 & zi0(e,3).lt.(one +tolpt).and. !should not be necessary
|
| H A D | asiqGradV.f | 40 sgn(:,i) = -one
42 sgn(:,i) = one
|
| /phasta/phSolver/compressible/ |
| H A D | itrPC.f | 17 almi = one
18 alfi = one
19 gami = one
22 almi = (three-rhoinf(itseq))/(one+rhoinf(itseq))/two
23 alfi = one/(one+rhoinf(itseq))
32 & y=y-(one-alfi)*Delt(1)*CFLfl(1)*acold
42 Dtgl = one / Delt(itseq) ! caution: inverse of time step
72 ac = acold*(one-almi/gami)
79 y = yold+alfi/Dtgl*acold*(one-gami)
85 ac = acold*(one-almi)
[all …]
|
| H A D | getthm.f | 125 prop_blend(i) = 0.5*(one + Sclr(i)/epsilon_ls +
128 prop_blend(i) = one
168 alfap = one / T
169 betaT = one / pres
180 prop_blend(i) = 0.5*(one + Sclr(i)/epsilon_ls +
183 prop_blend(i) = one
201 tmpg = one / T ! alfap for gas phase
205 tmpg = one / pres ! betaT for gas phase
279 prop_blend(i) = 0.5*(one + Sclr(i)/epsilon_ls +
282 prop_blend(i) = one
[all …]
|
| H A D | bc3lhs.f | 61 EGmass(iel,i1,i1) = one
90 EGmass(iel,i2,i2) = one
113 EGmass(iel,i3,i3) = one
136 EGmass(iel,i2,i2) = one
137 EGmass(iel,i3,i3) = one
160 EGmass(iel,i4,i4) = one
183 EGmass(iel,i2,i2) = one
184 EGmass(iel,i4,i4) = one
207 EGmass(iel,i3,i3) = one
208 EGmass(iel,i4,i4) = one
[all …]
|
| H A D | bc3bdg.f | 60 BDiag(:,1,1) = one
74 BDiag(:,1,1) = one
115 BDiag(:,2,2) = one
153 BDiag(:,3,3) = one
186 BDiag(:,3,3) = one
187 BDiag(:,2,2) = one
225 BDiag(:,4,4) = one
258 BDiag(:,4,4) = one
259 BDiag(:,2,2) = one
292 BDiag(:,4,4) = one
[all …]
|
| H A D | e3eig1.f | 57 fact = one / u
62 Rcs(:,1,1) = one
86 & + sign(one, Rcs(:,1,2)*Rcs(:,1,3)) * Rcs(:,1,3)
88 & - sign(one, Rcs(:,1,2)*Rcs(:,1,3)) * Rcs(:,1,3)
89 Rcs(:,2,3) = sign(one, Rcs(:,1,2)*Rcs(:,1,3)) *
96 fact = one / sqrt(fact)
173 fact = (one / sqrt( two * rho * T )) / c
199 fact = (one / sqrt( cp * rho )) / T
292 t = sign(one, t) / ( abs(t) + sqrt(one + t**2) )
296 Rcs = one / sqrt( one + t**2 )
|
| H A D | e3source.f | 309 fd=one-tanh((8.0000000000000000d0*rd)**3)
322 fv2 = one - chi/(one+chi*fv1)
343 fwog = ((one+saCw3**6)/(g**6+saCw3**6))**sixth
356 if(one.eq.zero) then
362 saCb1Scale = one
367 srcrat = saCb1Scale*saCb1*(one-ft2)*Stilde
370 srcrat=saCb1*(one-ft2)*Stilde
389 fv2p = (chi**2)*fv1p-(one/rnu)
390 fv2p = fv2p/(one+chi*fv1)**2
399 gp = one+saCw2*(6*r**5 - one)
[all …]
|
| H A D | i3ldu.f | 52 Diag(:, 1) = one / sqrt (Diag(:, 1))
56 Diag(:, 3) = one / sqrt (Diag(:, 3))
63 Diag(:, 6) = one / sqrt (Diag(:, 6))
75 Diag(:,10) = one / sqrt (Diag(:,10))
92 Diag(:,15) = one / sqrt (Diag(:,15))
|
| H A D | i3lu.f | 43 Diag(:,1,1) = one / Diag(:,1,1)
54 Diag(:,2,2) = one / Diag(:,2,2)
70 Diag(:,3,3) = one / Diag(:,3,3)
82 Diag(:,4,4) = one / Diag(:,4,4)
93 Diag(:,5,5) = one / Diag(:,5,5)
|
| H A D | timeseries.f | 145 if ((abs(zi0(e,1)).lt.(one+tolpt)).and.
146 & (abs(zi0(e,2)).lt.(one+tolpt)).and.
147 & (abs(zi0(e,3)).lt.(one+tolpt))) then ! got the element
194 if ((zi0(e,1)+zi0(e,2)+zi0(e,3)).lt.(one+tolpt).and.
195 …& zi0(e,1).lt.(one +tolpt).and. !should not be necessary; the limit in the tet …
197 & zi0(e,2).lt.(one +tolpt).and. !should not be necessary
199 & zi0(e,3).lt.(one +tolpt).and. !should not be necessary
|
| /phasta/phSolver/common/test/ |
| H A D | phIOwrite.f | 8 integer, target :: rank, ierror, one, ppf, peers, fish, nfiles variable
22 one = 1
33 call syncio_setup_write(nfiles, one, ppf, handle(1))
42 call phio_writeheader(handle(i), phrase, c_loc(fish), one, one,
45 & one, dataDbl, iotype)
54 & one, dataDbl, iotype)
56 & one, dataDbl, iotype)
|
| H A D | phIOwrite.cc | 23 int one = 1; in main() local
31 syncio_setup_write(nfiles, one, ppf, &(file[0])); in main()
35 phio_writeheader(file[i], phrase, &fish, &one, &one, type, iotype); in main()
36 phio_writedatablock(file[i], phrase, &(fishWeight[i]), &one, type, iotype); in main()
43 phio_readheader(file[i], phrase, &(numFish[i]), &one, type, iotype); in main()
44 phio_readdatablock(file[i], phrase, &(fishWeight[i]), &one, type, iotype); in main()
|
| H A D | phIOwriteReadZeroSz.cc | 28 int one = 1; in main() local
37 syncio_setup_write(nfiles, one, ppf, &(file[0])); in main()
44 phio_writeheader(file[i], phrase, &zero, &one, &zero, type, iotype); in main()
56 phio_readheader(file[i], phrase, &numFish, &one, type, iotype); in main()
|
| H A D | phIOwriteFields.cc | 30 int one = 1; in main() local
52 phio_writeheader(file[i], str, &blockEntries, &one, &blockEntries, type, iotype); in main()
65 phio_readheader(file[i], phrase, &numFish, &one, type, iotype); in main()
|
| H A D | phIOreadheader.cc | 26 int one = 1; in main() local
32 &one, "integer", iotype); in main()
|
| /phasta/phSolver/common/ |
| H A D | perprep.f | 33 one=1.00000000000
34 rcount=one
38 rcount(i) = rcount(i) + one
42 if(rcount(k).ne.one) then
43 rcount(k)=one/rcount(k)
|
| H A D | asbnabi.f | 86 temp = one / sqrt ( temp1**2 + temp2**2 + temp3**2 )
92 Wfactor(:) = one / (four*temp(:))
94 Wfactor(:) = one / (four*temp(:))
96 Wfactor(:) = one / (two*temp(:))
98 Wfactor(:) = one / (four*temp(:))
100 Wfactor(:) = one / (four*temp(:))
102 Wfactor(:) = one / (two*temp(:))
244 temp = one / sqrt ( temp1**2 + temp2**2 + temp3**2 )
248 Wfactor(:) = one / (four*temp(:))
250 Wfactor(:) = one / (four*temp(:))
[all …]
|
| H A D | getvel.f | 99 rinvsons = one/nsons ! division is expensive
106 denom=max(one*(lstep),one)
108 tavef=one/denom
120 velbar(:,:)=tavef*velft(:,:)+(one-tavef)*velbar(:,:)
|
| H A D | elem-search.f | 79 if (zi0(e,1).lt.(one+tolpt).and.
84 if (zi0(e,2).lt.(one+tolpt).and.
89 if (zi0(e,3).lt.(one+tolpt).and.
94 if ((1-zi0(e,1)-zi0(e,2)-zi0(e,3)).lt.(one+tolpt).and.
|
| H A D | proces.f | 137 & +(one-xi)*bcinterp(iupper-1,4:6))
139 & +(one-xi)*bcinterp(iupper-1,2))
141 & +(one-xi)*bcinterp(iupper-1,3))
147 & +(one-xi)*bcinterp(iupper-1,j+1))
|