| 58 | 
  | 
#ifdef IS_MPI | 
| 59 | 
  | 
 | 
| 60 | 
  | 
  enum direction { | 
| 61 | 
< | 
    I = 0, | 
| 62 | 
< | 
    J = 1 | 
| 61 | 
> | 
    Row = 0, | 
| 62 | 
> | 
    Column = 1 | 
| 63 | 
  | 
  }; | 
| 64 | 
  | 
     | 
| 65 | 
  | 
  template<typename T>  | 
| 79 | 
  | 
  template<> const int MPITraits<Mat3x3d>::dim = 9; | 
| 80 | 
  | 
   | 
| 81 | 
  | 
  template<direction D, typename T> | 
| 82 | 
< | 
  class Comm {  | 
| 82 | 
> | 
  class Communicator {  | 
| 83 | 
  | 
  public:  | 
| 84 | 
  | 
     | 
| 85 | 
< | 
    Comm<D, T>(int nObjects) { | 
| 85 | 
> | 
    Communicator<D, T>(int nObjects) { | 
| 86 | 
  | 
       | 
| 87 | 
  | 
      int nProc = MPI::COMM_WORLD.Get_size(); | 
| 88 | 
  | 
      int myRank = MPI::COMM_WORLD.Get_rank(); | 
| 98 | 
  | 
      rowIndex_ = myRank / nColumns;       | 
| 99 | 
  | 
      columnIndex_ = myRank % nColumns; | 
| 100 | 
  | 
 | 
| 101 | 
< | 
      if (D == I) { | 
| 101 | 
> | 
      if (D == Row) { | 
| 102 | 
  | 
        myComm = MPI::COMM_WORLD.Split(rowIndex_, 0); | 
| 103 | 
  | 
      } else { | 
| 104 | 
  | 
        myComm = MPI::COMM_WORLD.Split(columnIndex_, 0); |