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

updated TestBlobAnalyze

Browse Source
This commit is contained in:
James E McClure 2015-06-05 15:41:29 -04:00
parent 4ee455ddc7
commit 9b33aa132f
1 changed files with 10 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
tests/TestBlobAnalyze.cpp
View File

@ -234,20 +234,22 @@ int main(int argc, char **argv)
double beta = 0.95;
Averages.SetupCubes(Dm);
Averages.UpdateSolid();
if (rank==0) printf("initializing the system \n");
Averages.Initialize();
if (rank==0) printf("updating mesh \n");
if (rank==0) printf("updating mesh halos \n");
Averages.UpdateMeshValues();
if (rank==0) printf("computing blobs \n");
if (rank==0) printf("computing local averages \n");
Averages.ComputeLocalBlob();
if (rank==0) printf("reducing averages \n");
Averages.Reduce();
if (rank==0) printf("Writing blobs \n");
if (rank==0){
FILE *PHASE;
PHASE = fopen("Phase.00000","wb");
fwrite(Averages.SDn.get(),8,Nx*Ny*Nz,PHASE);
fclose(PHASE);
}
// Write the local blob ids
sprintf(LocalRankFilename,"BlobLabel.%05i",rank);
FILE *BLOBLOCAL = fopen(LocalRankFilename,"wb");
fwrite(GlobalBlobID.get(),4,GlobalBlobID.length(),BLOBLOCAL);
fclose(BLOBLOCAL);
printf("Wrote BlobLabel.%05i \n",rank);
if (rank==0) printf("Sorting averages \n");
// Blobs.Set(Averages.BlobAverages.NBLOBS);