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

Updates to track interface speed using signed distance function instead of phase density field

Browse Source
This commit is contained in:
James E McClure
2015-01-28 11:18:48 -05:00
parent 7274795794
commit 434bdf374e
1 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
tests/lb2_Color_wia_mpi.cpp
View File

@@ -1110,8 +1110,6 @@ int main(int argc, char **argv)
Averages.UpdateSolid(); Averages.UpdateSolid();
//....................................................................... //.......................................................................
//......................................................................
//************************************************************************* //*************************************************************************
// Compute the phase indicator field and reset Copy, Den // Compute the phase indicator field and reset Copy, Den
//************************************************************************* //*************************************************************************
@@ -1665,7 +1663,8 @@ int main(int argc, char **argv)
//........................................................................... //...........................................................................
// Copy the phase indicator field for the earlier timestep // Copy the phase indicator field for the earlier timestep
DeviceBarrier(); DeviceBarrier();
CopyToHost(Averages.Phase_tplus.data,Phi,N*sizeof(double)); CopyToHost(Averages.Phase.data,Phi,N*sizeof(double));
Averages.ColorToSignedDistance(beta,Averages.Phase.data,Averages.Phase_tplus.data);
//........................................................................... //...........................................................................
} }
if (timestep%1000 == 0){ if (timestep%1000 == 0){
@@ -1688,9 +1687,11 @@ int main(int argc, char **argv)
// Copy the phase indicator field for the later timestep // Copy the phase indicator field for the later timestep
DeviceBarrier(); DeviceBarrier();
CopyToHost(Averages.Phase_tminus.data,Phi,N*sizeof(double)); CopyToHost(Averages.Phase_tminus.data,Phi,N*sizeof(double));
Averages.ColorToSignedDistance(beta,Averages.Phase_tminus.data,Averages.Phase_tminus.data);
//.................................................................... //....................................................................
// The following need to be called each time new averages are computed // The following need to be called each time new averages are computed
Averages.Initialize(); Averages.Initialize();
Averages.ColorToSignedDistance(beta,Averages.Phase.data,Averages.SDn.data);
Averages.UpdateMeshValues(); Averages.UpdateMeshValues();
Averages.ComputeLocal(); Averages.ComputeLocal();
Averages.Reduce(); Averages.Reduce();