From f0f761f8303bd23f4fcd440c16224c10708a9d0a Mon Sep 17 00:00:00 2001 From: James McClure Date: Sun, 23 Mar 2014 07:55:37 -0400 Subject: [PATCH] trying to fix nan --- tests/lb2_Color_wia_mpi.cpp | 15 +-------------- 1 file changed, 1 insertion(+), 14 deletions(-) diff --git a/tests/lb2_Color_wia_mpi.cpp b/tests/lb2_Color_wia_mpi.cpp index 34e03955..cb915512 100644 --- a/tests/lb2_Color_wia_mpi.cpp +++ b/tests/lb2_Color_wia_mpi.cpp @@ -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;