Commit a0dcb435 authored by frey_m's avatar frey_m
Browse files

EllipticDomain: change for loop order

parent 34744fe7
...@@ -84,8 +84,8 @@ void EllipticDomain::compute(Vector_t hr, NDIndex<3> localId) { ...@@ -84,8 +84,8 @@ void EllipticDomain::compute(Vector_t hr, NDIndex<3> localId) {
* grid points per plane --> otherwise we might * grid points per plane --> otherwise we might
* get not unique global indices in the Tpetra::CrsMatrix * get not unique global indices in the Tpetra::CrsMatrix
*/ */
for (y = 0; y < nr_m[1]; ++y) { for (x = 0; x < nr_m[0]; ++x) {
for (x = 0; x < nr_m[0]; ++x) { for (y = 0; y < nr_m[1]; ++y) {
if (isInside(x, y, 1)) { if (isInside(x, y, 1)) {
idxMap_m[toCoordIdx(x, y)] = idx; idxMap_m[toCoordIdx(x, y)] = idx;
coordMap_m[idx++] = toCoordIdx(x, y); coordMap_m[idx++] = toCoordIdx(x, y);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment