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

refactoring lbpm_color_simulator

Browse Source
This commit is contained in:
James E McClure 2015-06-24 17:24:39 -04:00
parent 6e0913d5f9
commit d301d8b14e
1 changed files with 3 additions and 41 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

44
tests/lbpm_color_simulator.cpp
View File

@ -1127,6 +1127,9 @@ int main(int argc, char **argv)
ColorBC_outlet(Phi,Den,A_even,A_odd,B_even,B_odd,Nx,Ny,Nz);
}
ComputePressureD3Q19(ID,f_even,f_odd,Pressure,Nx,Ny,Nz);
ComputeVelocityD3Q19(ID,f_even,f_odd,Velocity,Nx,Ny,Nz);
//...........................................................................
// Copy the phase indicator field for the earlier timestep
DeviceBarrier();
@ -1762,47 +1765,6 @@ int main(int argc, char **argv)
fwrite(Averages.Press.get(),8,N,PRESS);
fclose(PRESS);
/* sprintf(LocalRankFilename,"%s%s","dPdt.",LocalRankString);
FILE *SPEED;
SPEED = fopen(LocalRankFilename,"wb");
fwrite(dPdt.get(),8,N,SPEED);
fclose(SPEED);
sprintf(LocalRankFilename,"%s%s","DenA.",LocalRankString);
FILE *DENA;
DENA = fopen(LocalRankFilename,"wb");
fwrite(&Den[0],8,N,DENA);
fclose(DENA);
sprintf(LocalRankFilename,"%s%s","DenB.",LocalRankString);
FILE *DENB;
DENB = fopen(LocalRankFilename,"wb");
fwrite(&Den[N],8,N,DENB);
fclose(DENB);
sprintf(LocalRankFilename,"%s%s","GradMag.",LocalRankString);
FILE *GRAD;
GRAD = fopen(LocalRankFilename,"wb");
for (k=0; k<Nz; k++){
for (j=0; j<Ny; j++){
for (i=0; i<Nx; i++){
Phase(i,j,k) = sqrt(Phase_x(i,j,k)*Phase_x(i,j,k) + Phase_y(i,j,k)*Phase_y(i,j,k) + Phase_z(i,j,k)*Phase_z(i,j,k));
}
}
}
fwrite(Phase.get(),8,N,GRAD);
fclose(GRAD);
*/
/* double *DensityValues;
DensityValues = new double [2*N];
CopyToHost(DensityValues,Copy,2*N*sizeof(double));
FILE *PHASE;
PHASE = fopen(LocalRankFilename,"wb");
fwrite(DensityValues,8,2*N,PHASE);
fclose(PHASE);
*/ //************************************************************************/
// ****************************************************
MPI_Barrier(MPI_COMM_WORLD);
MPI_Finalize();