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); |