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

Fixed bug in tests/TestContactAngle

Browse Source
This commit is contained in:
James McClure
2014-12-29 18:33:38 -05:00
parent 9747aff13e
commit eb98f9b793
1 changed files with 13 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
tests/TestContactAngle.cpp
View File

@@ -94,6 +94,19 @@ int main (int argc, char *argv[])
int ncubes = (Nx-2)*(Ny-2)*(Nz-2); // Exclude the "upper" halo int ncubes = (Nx-2)*(Ny-2)*(Nz-2); // Exclude the "upper" halo
IntArray cubeList(3,ncubes); IntArray cubeList(3,ncubes);
pmmc_CubeListFromMesh(cubeList, ncubes, Nx, Ny, Nz); pmmc_CubeListFromMesh(cubeList, ncubes, Nx, Ny, Nz);
int nc=1;
for (k=1; k<Nz-1; k++){
for (j=1; j<Ny-1; j++){
for (i=1; i<Nx-1; i++){
cubeList(0,nc) = i;
cubeList(1,nc) = j;
cubeList(2,nc) = k;
nc++;
}
}
}
ncubes = nc;
//........................................................................... //...........................................................................
double Cx,Cy,Cz; double Cx,Cy,Cz;
double dist1,dist2; double dist1,dist2;