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

debugging dfh bubble

Browse Source
This commit is contained in:
James E McClure
2018-05-02 21:12:37 -04:00
parent 1d34a86d8e
commit 8550c57e76
2 changed files with 5 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
cpu/dfh.cpp
View File

@@ -16,9 +16,9 @@ extern "C" void ScaLBL_Gradient_Unpack(double weight, double Cqx, double Cqy, do
value=weight*recvbuf[idx];
if (!(n<0)){
// PARALLEL UPDATE MUST BE DONE ATOMICALLY
grad[n] += Cqx*value;
grad[N+n] += Cqy*value;
grad[2*N+n] += Cqz*value;
grad[n] -= Cqx*value;
grad[N+n] -= Cqy*value;
grad[2*N+n] -= Cqz*value;
}
}
}

3
tests/TestBubbleDFH.cpp
View File

@@ -92,7 +92,7 @@ int main(int argc, char **argv)
int RESTART_INTERVAL=20000;
//int ANALYSIS_)INTERVAL=1000;
int BLOB_ANALYSIS_INTERVAL=1000;
int timestep = 6;
int timestep = 0;
if (rank==0){
//.............................................................
@@ -497,6 +497,7 @@ int main(int argc, char **argv)
err = 1.0;
double sat_w_previous = 1.01; // slightly impossible value!
if (rank==0) printf("Begin timesteps: error tolerance is %f \n", tol);
if (rank==0) printf("ScaLBL_Comm.next=%i, ScaLBL_Comm.first_interior=%i, %i",ScaLBL_Comm.next, ScaLBL_Comm.first_interior, ScaLBL_Comm.last_interior);
if (rank==0){
printf("Analysis intervals: (restart) %i, (TCAT) %i, (blobtracking) %i \n",RESTART_INTERVAL,ANALYSIS_INTERVAL,BLOBID_INTERVAL);
}