reduction for scalar minkowski

analysis/Minkowski.cpp

@@ -234,6 +234,14 @@ void Minkowski::ComputeScalar(const DoubleArray Field, const double isovalue)
 			}
 		}
 	}
+	MPI_Barrier(Dm->Comm);
+	// Phase averages
+	MPI_Allreduce(&Vi,&Vi_global,1,MPI_DOUBLE,MPI_SUM,Dm->Comm);
+	MPI_Allreduce(&Xi,&Xi_global,1,MPI_DOUBLE,MPI_SUM,Dm->Comm);
+	MPI_Allreduce(&Ai,&Ai_global,1,MPI_DOUBLE,MPI_SUM,Dm->Comm);
+	MPI_Allreduce(&Ji,&Ji_global,1,MPI_DOUBLE,MPI_SUM,Dm->Comm);
+	MPI_Barrier(Dm->Comm);
+
 }
 
 /*

analysis/Minkowski.h

@@ -60,6 +60,8 @@ public:
 	double euler,Kn,Jn,An;
 	double euler_global,Kn_global,Jn_global,An_global;
 	double Ai,Ji,Xi,Vi;
+	double Ai_global,Ji_global,Xi_global,Vi_global;
+
 	//...........................................................................
 	int Nx,Ny,Nz;
 	IntArray PhaseID;	// Phase ID array (solid=0, non-wetting=1, wetting=2)