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

trying to fix nan

Browse Source
This commit is contained in:
James McClure 2014-03-23 07:55:37 -04:00
parent 199509c56b
commit f0f761f830
1 changed files with 1 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
tests/lb2_Color_wia_mpi.cpp
View File

@ -477,7 +477,7 @@ int main(int argc, char **argv)
// Generate the residual NWP
if (rank==0) printf("Initializing with NWP saturation = %f \n",wp_saturation);
// if (!pBC) GenerateResidual(id,Nx,Ny,Nz,wp_saturation);
if (!pBC) GenerateResidual(id,Nx,Ny,Nz,wp_saturation);
#endif
// Set up MPI communication structurese
@ -2226,19 +2226,6 @@ int main(int argc, char **argv)
van_global(1) = van_global(1) / vol_n_global;
van_global(2) = van_global(2) / vol_n_global;
if (rank==0){
printf("vol_w = %f \n", vol_w);
printf("paw = %f \n", paw);
printf("vol_w_global = %f \n", vol_w_global);
printf("paw_global = %f \n", paw_global);
printf("vaw_global(0) = %f \n", vaw_global(0));
printf("vaw_global(1) = %f \n", vaw_global(1));
printf("vaw_global(2) = %f \n", vaw_global(2));
}
if (rank==25){
printf("rank25: vol_w = %f \n", vol_w);
printf("rank25: paw = %f \n", paw);
}
// Normalize surface averages by the interfacial area
Jwn_global /= awn_global;
efawns_global /= lwns_global;