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

Updated BasicSimulator

Browse Source
This commit is contained in:
James McClure
2014-07-26 13:58:34 -04:00
parent b3f53f35e4
commit 2421289b5d
1 changed files with 7 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
tests/BasicSimulator.cpp
View File

@@ -265,7 +265,6 @@ int main(int argc, char **argv)
printf("das = %f \n", das);
printf("dbs = %f \n", dbs);
printf("Value of phi at solid surface = %f \n", phi_s);
printf("Distance to phi = 0.0: %f \n", xIntPos);
printf("gamma_{wn} = %f \n", 5.796*alpha);
// printf("cos theta_c = %f \n", 1.05332*Ps);
printf("Initial Force(x) = %f \n", Fx);
@@ -274,6 +273,12 @@ int main(int argc, char **argv)
printf("Sub-domain size = %i x %i x %i\n",Nz,Nz,Nz);
printf("Parallel domain size = %i x %i x %i\n",nprocx,nprocy,nprocz);
printf("********************************************************\n");
if (rank==0) printf("Number of paired forces = %i \n", SimCount);
if (rank==0) printf("Percent difference for force pair = %i \n", SimDelta);
if (rank==0) printf("********************************************************\n");
}
InitializeRanks( rank, nprocx, nprocy, nprocz, iproc, jproc, kproc,
@@ -289,11 +294,6 @@ int main(int argc, char **argv)
int N = Nx*Ny*Nz;
int dist_mem_size = N*sizeof(double);
if (rank==0) printf("Number of paired forces = %i \n", SimCount);
if (rank==0) printf("Percent difference for force pair = %i \n", SimDelta);
if (rank==0) printf("********************************************************\n");
//.......................................................................
if (rank == 0) printf("Read input media... \n");
//.......................................................................
@@ -1462,6 +1462,7 @@ int main(int argc, char **argv)
Fy += SimDelta*Fy;
Fz += SimDelta*Fz;
}
printf("Simulating {Fx,Fy,Fz} = %f, %f, %f", Fx,Fy,Fz);
//************ MAIN ITERATION LOOP ***************************************/
MPI_Barrier(MPI_COMM_WORLD);