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

Debugging tests

Browse Source
This commit is contained in:
James E McClure
2018-01-26 15:19:50 -05:00
parent eda88946f9
commit 056ba3a1ae
3 changed files with 9 additions and 40 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
common/ScaLBL.h
View File

@@ -257,7 +257,7 @@ public:
void SendHalo(double *data); void SendHalo(double *data);
void RecvHalo(double *data); void RecvHalo(double *data);
void RegularLayout(IntArray map, double *data, DoubleArray regdata); void RegularLayout(IntArray map, double *data, DoubleArray regdata);
// Routines to set boundary conditions // Routines to set boundary conditions
void Color_BC_z(int *Map, double *Phi, double *Den, double vA, double vB); void Color_BC_z(int *Map, double *Phi, double *Den, double vA, double vB);
void Color_BC_Z(int *Map, double *Phi, double *Den, double vA, double vB); void Color_BC_Z(int *Map, double *Phi, double *Den, double vA, double vB);

42
tests/TestColorBubble.cpp
View File

@@ -521,50 +521,18 @@ int main(int argc, char **argv)
double PoreVel = sum*iVol_global; double PoreVel = sum*iVol_global;
if (rank==0) printf("Velocity = %f \n",PoreVel); if (rank==0) printf("Velocity = %f \n",PoreVel);
double *PHASE; /*
double *PHASE;
PHASE= new double [Nx*Ny*Nz]; PHASE= new double [Nx*Ny*Nz];
SIZE=Nx*Ny*Nz*sizeof(double); SIZE=Nx*Ny*Nz*sizeof(double);
ScaLBL_CopyToHost(&PHASE[0],&Phi[0],SIZE); ScaLBL_CopyToHost(&PHASE[0],&Phi[0],SIZE);
FILE *OUTFILE; FILE *OUTFILE;
OUTFILE = fopen("Phase.raw","wb"); sprintf(LocalRankFilename,"Phase.%05i.raw",rank);
OUTFILE = fopen(LocalRankFilename,"wb");
fwrite(PHASE,8,N,OUTFILE); fwrite(PHASE,8,N,OUTFILE);
fclose(OUTFILE); fclose(OUTFILE);
*/
double *DENA, *DENB, *TMPDAT;
SIZE=Np*sizeof(double);
TMPDAT = new double [Np];
DENA= new double [Nx*Ny*Nz];
DENB= new double [Nx*Ny*Nz];
ScaLBL_CopyToHost(&TMPDAT[0],&Den[0],SIZE);
ScaLBL_Comm.RegularLayout(Map,TMPDAT,DENA);
ScaLBL_CopyToHost(&TMPDAT[0],&Den[Np],SIZE);
ScaLBL_Comm.RegularLayout(Map,TMPDAT,DENB);
FILE *AFILE;
sprintf(LocalRankFilename,"na.%05i.raw",rank);
AFILE = fopen(LocalRankFilename,"wb");
fwrite(DENA,8,N,AFILE);
fclose(AFILE);
FILE *BFILE;
BFILE = fopen("nB.raw","wb");
fwrite(DENB,8,N,BFILE);
fclose(BFILE);
double *CG;
CG= new double [3*Np];
ScaLBL_CopyToHost(&CG[0],&ColorGrad[0],3*SIZE);
for (int idx=0; idx<Np; idx++){
double C=CG[idx]*CG[idx]+CG[Np+idx]*CG[Np+idx]+CG[2*Np+idx]*CG[2*Np+idx];
TMPDAT[idx]=C;
}
ScaLBL_Comm.RegularLayout(Map,TMPDAT,DENB);
FILE *CGFILE;
CGFILE = fopen("cgrad.raw","wb");
fwrite(DENB,8,N,CGFILE);
fclose(CGFILE);
} }
// **************************************************** // ****************************************************

5
tests/TestColorSquareTube.cpp
View File

@@ -483,7 +483,8 @@ int main(int argc, char **argv)
} }
} }
double *PHASE; /*
double *PHASE;
PHASE= new double [Nx*Ny*Nz]; PHASE= new double [Nx*Ny*Nz];
SIZE=Nx*Ny*Nz*sizeof(double); SIZE=Nx*Ny*Nz*sizeof(double);
ScaLBL_CopyToHost(&PHASE[0],&Phi[0],SIZE); ScaLBL_CopyToHost(&PHASE[0],&Phi[0],SIZE);
@@ -531,7 +532,7 @@ int main(int argc, char **argv)
CGFILE = fopen(LocalRankFilename,"wb"); CGFILE = fopen(LocalRankFilename,"wb");
fwrite(DENB,8,N,CGFILE); fwrite(DENB,8,N,CGFILE);
fclose(CGFILE); fclose(CGFILE);
*/
} }
// **************************************************** // ****************************************************