diff --git a/tests/BlobAnalyzeParallel.cpp b/tests/BlobAnalyzeParallel.cpp index 2f132e87..dd519028 100644 --- a/tests/BlobAnalyzeParallel.cpp +++ b/tests/BlobAnalyzeParallel.cpp @@ -38,7 +38,7 @@ inline void WriteBlobStates(TwoPhase TCAT, double D, double porosity){ pw = TCAT.paw_global; aws = TCAT.aws; // Compute the averages over the entire non-wetting phsae - for (a=0; a<(int)TCAT.BlobAverages.size(1); a++){ + for (a=0; a0; a--){ + for (a=TCAT.nblobs_global-1; a>0; a--){ // Subtract the features one-by-one vol_n -= TCAT.BlobAverages(0,a); pan -= TCAT.BlobAverages(2,a)*TCAT.BlobAverages(0,a);