74 const Teuchos::ParameterList &pl_):
79 actual_procCoords(NULL),
81 actual_machine_extent(NULL),
82 is_transformed(false),
85 actual_machine_extent = machine_extent =
new int[networkDim];
87 actual_machine_extent = machine_extent;
90 actual_procCoords = procCoords =
new pcoord_t *[networkDim];
91 for (
int i = 0; i < networkDim; ++i) {
92 procCoords[i] =
new pcoord_t[this->
numRanks];
93 memset(procCoords[i], 0,
sizeof(pcoord_t) * this->
numRanks);
96 pcoord_t *xyz =
new pcoord_t[networkDim];
98 for (
int i = 0; i < networkDim; i++)
99 procCoords[i][this->
myRank] = xyz[i];
104 gatherMachineCoordinates(comm);
106 const Teuchos::ParameterEntry *pe2 =
107 this->pl->getEntryPtr(
"Machine_Optimization_Level");
111 int optimization_level;
112 optimization_level = pe2->getValue<
int>(&optimization_level);
114 if (optimization_level == 1) {
115 is_transformed =
true;
116 this->networkDim = 3;
117 procCoords =
new pcoord_t * [networkDim];
118 for(
int i = 0; i < networkDim; ++i) {
119 procCoords[i] =
new pcoord_t[this->
numRanks] ;
122 for (
int i = 0; i < this->
numRanks; ++i) {
123 procCoords[0][i] = this->actual_procCoords[0][i] * 8;
124 int yordinal = this->actual_procCoords[1][i];
125 procCoords[1][i] = yordinal/2 * (16 + 8) + (yordinal %2) * 8;
126 int zordinal = this->actual_procCoords[2][i];
127 procCoords[2][i] = zordinal * 5 + (zordinal / 8) * 3;
129 int mx = this->machine_extent[0];
130 int my = this->machine_extent[1];
131 int mz = this->machine_extent[2];
134 this->machine_extent =
new int[networkDim];
135 this->machine_extent[0] = mx * 8;
136 this->machine_extent[1] = my/2 * (16 + 8) + (my %2) * 8;
137 this->machine_extent[2] = mz * 5 + (mz / 8) * 3;
139 std::cout <<
"Transforming the coordinates" << std::endl;
142 else if(optimization_level >= 3) {
143 is_transformed =
true;
144 this->networkDim = 6;
145 procCoords =
new pcoord_t * [networkDim];
146 for(
int i = 0; i < networkDim; ++i) {
147 procCoords[i] =
new pcoord_t[this->
numRanks] ;
152 this->machine_extent =
new int[networkDim];
154 this->machine_extent[0] =
155 ceil (
int (this->actual_machine_extent[0]) / 2.0) * 64 ;
156 this->machine_extent[3] = 2 * 8 ;
157 this->machine_extent[1] =
158 ceil(
int (this->actual_machine_extent[1]) / 2.0) * 8 * 2400;
159 this->machine_extent[4] = 2 * 8;
160 this->machine_extent[2] =
161 ceil((
int (this->actual_machine_extent[2])) / 8.0) * 160;
162 this->machine_extent[5] = 8 * 5;
164 for (
int k = 0; k < this->
numRanks ; k++) {
176 procCoords[0][k] = (int (this->actual_procCoords[0][k]) / 2) * 64;
179 procCoords[3][k] = (int (this->actual_procCoords[0][k]) % 2) * 8 ;
186 (int (this->actual_procCoords[1][k]) / 2) * 8 * 2400;
188 procCoords[4][k] = (int (this->actual_procCoords[1][k]) % 2) * 8;
199 ((int (this->actual_procCoords[2][k])) / 8) * 160;
202 ((int (this->actual_procCoords[2][k])) % 8) * 5;
205 else if(optimization_level == 2) {
207 is_transformed =
true;
208 this->networkDim = 6;
209 procCoords =
new pcoord_t * [networkDim];
210 for(
int i = 0; i < networkDim; ++i) {
211 procCoords[i] =
new pcoord_t[this->
numRanks] ;
216 this->machine_extent =
new int[networkDim];
218 this->machine_extent[0] =
219 ceil(
int (this->actual_machine_extent[0]) / 3.0) * 128 ;
220 this->machine_extent[3] = 3 * 8 ;
221 this->machine_extent[1] =
222 ceil(
int (this->actual_machine_extent[1]) / 2.0) * 8 * 2400;
223 this->machine_extent[4] = 2 * 8;
224 this->machine_extent[2] =
225 ceil((
int (this->actual_machine_extent[2])) / 8.0) * 160;
226 this->machine_extent[5] = 8 * 5;
229 for (
int k = 0; k < this->
numRanks ; k++) {
240 procCoords[0][k] = (int (this->actual_procCoords[0][k]) / 3) * 128;
243 procCoords[3][k] = (int (this->actual_procCoords[0][k]) % 3) * 8 ;
250 (int (this->actual_procCoords[1][k]) / 2) * 8 * 2400;
252 procCoords[4][k] = (int (this->actual_procCoords[1][k]) % 2) * 8;
256 ((int (this->actual_procCoords[2][k])) / 8) * 160;
258 procCoords[5][k] = ((int (this->actual_procCoords[2][k])) % 8) * 5;