123 |
|
|
124 |
|
MPI_Status ierr; |
125 |
|
int intObIndex; |
126 |
< |
double zAngle; |
126 |
> |
RealType zAngle; |
127 |
|
|
128 |
|
if (masterNode == 0) { |
129 |
< |
std::map<int, double> zAngData; |
129 |
> |
std::map<int, RealType> zAngData; |
130 |
|
for(int i = 0 ; i < nproc; ++i) { |
131 |
|
if (i == masterNode) { |
132 |
|
for (mol = info_->beginMolecule(mi); mol != NULL; |
138 |
|
|
139 |
|
intObIndex = integrableObject->getGlobalIndex() ; |
140 |
|
zAngle = integrableObject->getZangle(); |
141 |
< |
zAngData.insert(std::pair<int, double>(intObIndex, zAngle)); |
141 |
> |
zAngData.insert(std::pair<int, RealType>(intObIndex, zAngle)); |
142 |
|
} |
143 |
|
} |
144 |
|
|
145 |
|
} else { |
146 |
|
for(int k = 0; k < nIntObjectsInProc[i]; ++k) { |
147 |
|
MPI_Recv(&intObIndex, 1, MPI_INT, i, 0, MPI_COMM_WORLD,&ierr); |
148 |
< |
MPI_Recv(&zAngle, 1, MPI_DOUBLE, i, 0, MPI_COMM_WORLD,&ierr); |
149 |
< |
zAngData.insert(std::pair<int, double>(intObIndex, zAngle)); |
148 |
> |
MPI_Recv(&zAngle, 1, MPI_REALTYPE, i, 0, MPI_COMM_WORLD,&ierr); |
149 |
> |
zAngData.insert(std::pair<int, RealType>(intObIndex, zAngle)); |
150 |
|
} |
151 |
|
} |
152 |
|
|
156 |
|
<< info_->getSnapshotManager()->getCurrentSnapshot()->getTime() |
157 |
|
<< " : omega values at this time\n"; |
158 |
|
|
159 |
< |
std::map<int, double>::iterator l; |
159 |
> |
std::map<int, RealType>::iterator l; |
160 |
|
for (l = zAngData.begin(); l != zAngData.end(); ++l) { |
161 |
|
finalOut << l->second << "\n"; |
162 |
|
} |
172 |
|
intObIndex = integrableObject->getGlobalIndex(); |
173 |
|
zAngle = integrableObject->getZangle(); |
174 |
|
MPI_Send(&intObIndex, 1, MPI_INT, masterNode, 0, MPI_COMM_WORLD); |
175 |
< |
MPI_Send(&zAngle, 1, MPI_DOUBLE, masterNode, 0, MPI_COMM_WORLD); |
175 |
> |
MPI_Send(&zAngle, 1, MPI_REALTYPE, masterNode, 0, MPI_COMM_WORLD); |
176 |
|
} |
177 |
|
} |
178 |
|
} |