| 57 |  |  | 
| 58 |  | #ifdef IS_MPI | 
| 59 |  |  | 
| 60 | < | enum direction { | 
| 61 | < | Row = 0, | 
| 62 | < | Column = 1 | 
| 60 | > | enum communicatorType { | 
| 61 | > | Global = 0, | 
| 62 | > | Row = 1, | 
| 63 | > | Column = 2 | 
| 64 |  | }; | 
| 65 |  |  | 
| 66 |  | template<typename T> | 
| 79 |  | template<> const MPI::Datatype MPITraits<Mat3x3d>::datatype = MPI_REALTYPE; | 
| 80 |  | template<> const int MPITraits<Mat3x3d>::dim = 9; | 
| 81 |  |  | 
| 82 | < | template<direction D, typename T> | 
| 82 | > | template<communicatorType D, typename T> | 
| 83 |  | class Communicator { | 
| 84 |  | public: | 
| 85 |  |  | 
| 99 |  | rowIndex_ = myRank / nColumns; | 
| 100 |  | columnIndex_ = myRank % nColumns; | 
| 101 |  |  | 
| 102 | < | if (D == Row) { | 
| 102 | > | switch(D) { | 
| 103 | > | case Row : | 
| 104 |  | myComm = MPI::COMM_WORLD.Split(rowIndex_, 0); | 
| 105 | < | } else { | 
| 105 | > | break; | 
| 106 | > | case Column: | 
| 107 |  | myComm = MPI::COMM_WORLD.Split(columnIndex_, 0); | 
| 108 | + | break; | 
| 109 | + | case Global: | 
| 110 | + | myComm = MPI::COMM_WORLD.Split(myRank, 0); | 
| 111 |  | } | 
| 112 |  |  | 
| 113 |  | int nCommProcs = myComm.Get_size(); |