--- trunk/OOPSE/libmdtools/mpiSimulation.cpp	2003/04/01 16:50:14	441
+++ trunk/OOPSE/libmdtools/mpiSimulation.cpp	2003/08/20 22:23:34	708
@@ -4,7 +4,6 @@
 #include <cstring>
 #include <cmath>
 #include <mpi.h>
-#include <mpi++.h>
 
 #include "mpiSimulation.hpp"
 #include "simError.h"
@@ -20,7 +19,7 @@ mpiSimulation::mpiSimulation(SimInfo* the_entryPlug)
   entryPlug = the_entryPlug;
   mpiPlug = new mpiSimData;
   
-  mpiPlug->numberProcessors = MPI::COMM_WORLD.Get_size();
+  MPI_Comm_size(MPI_COMM_WORLD, &(mpiPlug->numberProcessors) );
   mpiPlug->myNode = worldRank;
 
   MolToProcMap = new int[entryPlug->n_mol];
@@ -67,7 +66,7 @@ int* mpiSimulation::divideLabor( void ){
   int nmol_global, natoms_global;
   int local_index, index;
   int smallDiff, bigDiff;
-  int baseSeed = BASE_SEED;
+  int baseSeed = entryPlug->getSeed();
 
   int testSum;
 
@@ -83,6 +82,7 @@ int* mpiSimulation::divideLabor( void ){
   mpiPlug->nSRIGlobal = entryPlug->n_SRI;
   mpiPlug->nMolGlobal = entryPlug->n_mol;
 
+  
   myRandom = new randomSPRNG( baseSeed );
 
   a = 3.0 * (double)mpiPlug->nMolGlobal / (double)mpiPlug->nAtomsGlobal;
@@ -203,32 +203,32 @@ int* mpiSimulation::divideLabor( void ){
 
     // Spray out this nonsense to all other processors:
 
-    MPI::COMM_WORLD.Bcast(MolToProcMap, mpiPlug->nMolGlobal, 
-                          MPI_INT, 0);
+    MPI_Bcast(MolToProcMap, mpiPlug->nMolGlobal, 
+	      MPI_INT, 0, MPI_COMM_WORLD);
 
-    MPI::COMM_WORLD.Bcast(AtomToProcMap, mpiPlug->nAtomsGlobal, 
-                          MPI_INT, 0);
+    MPI_Bcast(AtomToProcMap, mpiPlug->nAtomsGlobal, 
+	      MPI_INT, 0, MPI_COMM_WORLD);
 
-    MPI::COMM_WORLD.Bcast(MolComponentType, mpiPlug->nMolGlobal, 
-                          MPI_INT, 0);
+    MPI_Bcast(MolComponentType, mpiPlug->nMolGlobal, 
+	      MPI_INT, 0, MPI_COMM_WORLD);
 
-    MPI::COMM_WORLD.Bcast(AtomsPerProc, mpiPlug->numberProcessors,
-                          MPI_INT, 0);    
+    MPI_Bcast(AtomsPerProc, mpiPlug->numberProcessors,
+	      MPI_INT, 0, MPI_COMM_WORLD);    
   } else {
 
     // Listen to your marching orders from processor 0:
     
-    MPI::COMM_WORLD.Bcast(MolToProcMap, mpiPlug->nMolGlobal, 
-                          MPI_INT, 0);
+    MPI_Bcast(MolToProcMap, mpiPlug->nMolGlobal, 
+	      MPI_INT, 0, MPI_COMM_WORLD);
     
-    MPI::COMM_WORLD.Bcast(AtomToProcMap, mpiPlug->nAtomsGlobal, 
-                          MPI_INT, 0);
+    MPI_Bcast(AtomToProcMap, mpiPlug->nAtomsGlobal, 
+	      MPI_INT, 0, MPI_COMM_WORLD);
 
-    MPI::COMM_WORLD.Bcast(MolComponentType, mpiPlug->nMolGlobal, 
-                          MPI_INT, 0);
+    MPI_Bcast(MolComponentType, mpiPlug->nMolGlobal, 
+	      MPI_INT, 0, MPI_COMM_WORLD);
     
-    MPI::COMM_WORLD.Bcast(AtomsPerProc, mpiPlug->numberProcessors,
-                          MPI_INT, 0);
+    MPI_Bcast(AtomsPerProc, mpiPlug->numberProcessors,
+	      MPI_INT, 0, MPI_COMM_WORLD);
 
 
   }
@@ -250,8 +250,10 @@ int* mpiSimulation::divideLabor( void ){
     }
   }
 
-  MPI::COMM_WORLD.Allreduce(&nmol_local,&nmol_global,1,MPI_INT,MPI_SUM);
-  MPI::COMM_WORLD.Allreduce(&natoms_local,&natoms_global,1,MPI_INT,MPI_SUM);
+  MPI_Allreduce(&nmol_local,&nmol_global,1,MPI_INT,MPI_SUM, 
+		MPI_COMM_WORLD);
+  MPI_Allreduce(&natoms_local,&natoms_global,1,MPI_INT,
+		MPI_SUM, MPI_COMM_WORLD);
   
   if( nmol_global != entryPlug->n_mol ){
     sprintf( painCave.errMsg,