OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Minor update

Browse Source
This commit is contained in:
James McClure
2014-11-08 10:31:25 -05:00
parent ca78dc7818
commit ff16995143
2 changed files with 5 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
common/pmmc.h
View File

@@ -4478,17 +4478,17 @@ inline void pmmc_InterfaceSpeed(DoubleArray &dPdt, DoubleArray &P_x, DoubleArray
P.y = 0.33333333333333333*(A.y+B.y+C.y); P.y = 0.33333333333333333*(A.y+B.y+C.y);
P.z = 0.33333333333333333*(A.z+B.z+C.z); P.z = 0.33333333333333333*(A.z+B.z+C.z);
if (temp > 0.0){ if (temp > 0.0){
x = Px.eval(P); x = Px.eval(P);
y = Py.eval(P); y = Py.eval(P);
z = Pz.eval(P); z = Pz.eval(P);
norm = sqrt(x*x+y*y+z*z); norm = sqrt(x*x+y*y+z*z);
if (norm==0.0) norm=1.0; if (norm==0.0) norm=1.0;
zeta = -Pt.eval(P) / norm; zeta = -Pt.eval(P) / norm;
temp = sqrt(temp)/norm;
AvgVel(0) += sqrt(temp)*zeta*x/norm; AvgVel(0) += temp*zeta*x;
AvgVel(1) += sqrt(temp)*zeta*y/norm; AvgVel(1) += temp*zeta*y;
AvgVel(2) += sqrt(temp)*zeta*z/norm; AvgVel(2) += temp*zeta*z;
} }
} }
//............................................................................. //.............................................................................

1
tests/lb2_Color_wia_mpi.cpp
View File

@@ -2193,6 +2193,7 @@ int main(int argc, char **argv)
Gns(0) = Gns(1) = Gns(2) = 0.0; Gns(0) = Gns(1) = Gns(2) = 0.0;
Gns(3) = Gns(4) = Gns(5) = 0.0; Gns(3) = Gns(4) = Gns(5) = 0.0;
vol_w = vol_n =0.0; vol_w = vol_n =0.0;
KGwns = KNwns = 0.0;
Jwn = Kwn = efawns = 0.0; Jwn = Kwn = efawns = 0.0;
trJwn = trawn = trRwn = 0.0; trJwn = trawn = trRwn = 0.0;