Commit f5a27b63 authored by Christof Metzger-Kraus's avatar Christof Metzger-Kraus

ArbitraryDomain only works when partitionend in z-direction; throw exception otherwise

parent f72d1842
......@@ -111,8 +111,18 @@ MGPoissonSolver::MGPoissonSolver ( PartBunch &beam,
throw OpalException("MGPoissonSolver::MGPoissonSolver",
"Geometry not known");
}
} else
} else {
NDIndex<3> localId = layout_m->getLocalNDIndex();
if (localId[0].length() != domain_m[0].length() ||
localId[1].length() != domain_m[1].length()) {
throw OpalException("ArbitraryDomain::compute",
"The class ArbitraryDomain only works with parallelization\n"
"in z-direction.\n"
"Please set PARFFTX=FALSE, PARFFTY=FALSE, PARFFTT=TRUE in \n"
"the definition of the field solver in the input file.\n");
}
bp = new ArbitraryDomain(currentGeometry, orig_nr_m, hr_m, interpl);
}
Map = 0;
A = Teuchos::null;
......
......@@ -1807,6 +1807,7 @@ Change orientation if diff is:
*gmsg << "* Triangle barycent built done." << endl;
*gmsg << *this << endl;
Ippl::Comm->barrier();
IpplTimings::stopTimer (Tinitialize_m);
}
......
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