| 2 |  | * @file Communicator.hpp | 
| 3 |  | * @author Charles Vardeman <cvardema.at.nd.edu> | 
| 4 |  | * @date 08/18/2010 | 
| 5 | – | * @time 11:56am | 
| 5 |  | * @version 1.0 | 
| 6 |  | * | 
| 7 |  | * @section LICENSE | 
| 41 |  | * | 
| 42 |  | * [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). | 
| 43 |  | * [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). | 
| 44 | < | * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). | 
| 45 | < | * [4]  Vardeman & Gezelter, in progress (2009). | 
| 44 | > | * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). | 
| 45 | > | * [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010). | 
| 46 | > | * [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). | 
| 47 |  | */ | 
| 48 |  |  | 
| 49 |  | #ifndef PARALLEL_COMMUNICATOR_HPP | 
| 71 |  | static int Length() { return 1; }; | 
| 72 |  | }; | 
| 73 |  |  | 
| 74 | < | template<> inline MPI::Datatype MPITraits<int>::Type() { return MPI_INT; } | 
| 75 | < | template<> inline MPI::Datatype MPITraits<RealType>::Type() { return MPI_REALTYPE; } | 
| 74 | > | template<> inline MPI::Datatype MPITraits<int>::Type() { return MPI::INT; } | 
| 75 | > | template<> inline MPI::Datatype MPITraits<RealType>::Type() { return MPI::REALTYPE; } | 
| 76 |  |  | 
| 77 |  | template<class T, unsigned int Dim> | 
| 78 |  | class MPITraits< Vector<T, Dim> > { | 
| 88 |  | static int Length() {return 3;} | 
| 89 |  | }; | 
| 90 |  |  | 
| 91 | < | template<class T, unsigned int Row, unsigned int Col> | 
| 92 | < | class MPITraits< RectMatrix<T, Row, Col> > { | 
| 91 | > | template<class T, unsigned int R, unsigned int C> | 
| 92 | > | class MPITraits< RectMatrix<T, R, C> > { | 
| 93 |  | public: | 
| 94 |  | static MPI::Datatype Type() { return MPITraits<T>::Type(); } | 
| 95 | < | static int Length() {return Row * Col;} | 
| 95 | > | static int Length() {return R * C;} | 
| 96 |  | }; | 
| 97 |  |  | 
| 98 |  | template<class T> | 
| 119 |  | if (nProc % i == 0) nColumns = i; | 
| 120 |  | } | 
| 121 |  |  | 
| 122 | < | int nRows = nProc / nColumns; | 
| 122 | > | // int nRows = nProc / nColumns; | 
| 123 |  | rowIndex_ = myRank / nColumns; | 
| 124 |  | columnIndex_ = myRank % nColumns; | 
| 125 |  |  | 
| 149 |  | class Plan { | 
| 150 |  | public: | 
| 151 |  |  | 
| 152 | < | Plan<T>(MPI::Intracomm comm, int nObjects) { | 
| 153 | < | myComm = comm; | 
| 152 | > | Plan<T>(MPI::Intracomm comm, int nObjects) : myComm(comm) { | 
| 153 |  | int nCommProcs = myComm.Get_size(); | 
| 154 |  |  | 
| 155 |  | counts.resize(nCommProcs, 0); |