ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/brains/SimCreator.cpp
(Generate patch)

Comparing branches/development/src/brains/SimCreator.cpp (file contents):
Revision 1794 by gezelter, Thu Sep 6 19:44:06 2012 UTC vs.
Revision 1798 by gezelter, Thu Sep 13 14:10:11 2012 UTC

# Line 507 | Line 507 | namespace OpenMD {
507      int nGlobalMols = info->getNGlobalMolecules();
508      std::vector<int> molToProcMap(nGlobalMols, -1); // default to an error condition:
509      
510 <    MPI_Comm_size(MPI_COMM_WORLD, &nProcessors);
510 >    nProcessors = MPI::COMM_WORLD.Get_size();
511      
512      if (nProcessors > nGlobalMols) {
513        sprintf(painCave.errMsg,
# Line 622 | Line 622 | namespace OpenMD {
622        delete myRandom;
623        
624        // Spray out this nonsense to all other processors:
625 <      
626 <      MPI_Bcast(&molToProcMap[0], nGlobalMols, MPI_INT, 0, MPI_COMM_WORLD);
625 >      MPI::COMM_WORLD.Bcast(&molToProcMap[0], nGlobalMols, MPI::INT, 0);
626      } else {
627        
628        // Listen to your marching orders from processor 0:
629 <      
631 <      MPI_Bcast(&molToProcMap[0], nGlobalMols, MPI_INT, 0, MPI_COMM_WORLD);
629 >      MPI::COMM_WORLD.Bcast(&molToProcMap[0], nGlobalMols, MPI::INT, 0);
630      }
631      
632      info->setMolToProcMap(molToProcMap);
# Line 860 | Line 858 | namespace OpenMD {
858      // This would be prettier if we could use MPI_IN_PLACE like the MPI-2
859      // docs said we could.
860      std::vector<int> tmpGroupMembership(info->getNGlobalAtoms(), 0);
861 <    MPI_Allreduce(&globalGroupMembership[0], &tmpGroupMembership[0], nGlobalAtoms,
862 <                  MPI_INT, MPI_SUM, MPI_COMM_WORLD);
861 >    MPI::COMM_WORLD.Allreduce(&globalGroupMembership[0],
862 >                              &tmpGroupMembership[0], nGlobalAtoms,
863 >                              MPI::INT, MPI::SUM);
864      info->setGlobalGroupMembership(tmpGroupMembership);
865   #else
866      info->setGlobalGroupMembership(globalGroupMembership);
# Line 878 | Line 877 | namespace OpenMD {
877      
878   #ifdef IS_MPI
879      std::vector<int> tmpMolMembership(info->getNGlobalAtoms(), 0);
880 +    MPI::COMM_WORLD.Allreduce(&globalMolMembership[0], &tmpMolMembership[0],
881 +                              nGlobalAtoms,
882 +                              MPI::INT, MPI::SUM);
883      
882 –    MPI_Allreduce(&globalMolMembership[0], &tmpMolMembership[0], nGlobalAtoms,
883 –                  MPI_INT, MPI_SUM, MPI_COMM_WORLD);
884 –    
884      info->setGlobalMolMembership(tmpMolMembership);
885   #else
886      info->setGlobalMolMembership(globalMolMembership);
# Line 897 | Line 896 | namespace OpenMD {
896      
897   #ifdef IS_MPI
898      std::vector<int> numIntegrableObjectsPerMol(info->getNGlobalMolecules(), 0);
899 <    MPI_Allreduce(&nIOPerMol[0], &numIntegrableObjectsPerMol[0],
900 <                  info->getNGlobalMolecules(), MPI_INT, MPI_SUM, MPI_COMM_WORLD);
899 >    MPI::COMM_WORLD.Allreduce(&nIOPerMol[0], &numIntegrableObjectsPerMol[0],
900 >                              info->getNGlobalMolecules(), MPI::INT, MPI::SUM);
901   #else
902      std::vector<int> numIntegrableObjectsPerMol = nIOPerMol;
903   #endif    

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines