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

don't compute in the halo region unnecessarily

Browse Source
This commit is contained in:
James E McClure 2015-06-26 14:15:17 -04:00
parent febb0beaaf
commit ead3fd59c2
1 changed files with 9 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
tests/lbpm_color_simulator.cpp
View File

@ -943,6 +943,15 @@ int main(int argc, char **argv)
if (rank==0) printf ("Copying phase ID to device \n");
char *ID;
AllocateDeviceMemory((void **) &ID, N); // Allocate device memory
// Don't compute in the halo
for (k=0;k<Nz;k++){
for (j=0;j<Ny;j++){
for (i=0;i<Nx;i++){
n = k*Nx*Ny+j*Nx+i;
if (i==0 || i==Nx-1 || j==0 || j=Ny-1 || k==0 || k == Nz-1) id[n] = 0;
}
}
}
// Copy to the device
CopyToDevice(ID, id, N);
//...........................................................................
@ -1766,16 +1775,6 @@ int main(int argc, char **argv)
fwrite(Averages.Press.get(),8,N,PRESS);
fclose(PRESS);
char *dvcid;
dvcid = new char [N];
CopyToHost(dvcid,ID,N*sizeof(char));
sprintf(LocalRankFilename,"%s%s","DVCID.",LocalRankString);
FILE *DVCID;
DVCID = fopen(LocalRankFilename,"wb");
fwrite(dvcid,1,N,DVCID);
fclose(DVCID);
// ****************************************************
MPI_Barrier(MPI_COMM_WORLD);
MPI_Finalize();