Skip to content
Snippets Groups Projects
Commit c2531957 authored by Robert Klöfkorn's avatar Robert Klöfkorn
Browse files

fixed Verschlimmbesserung.

[[Imported from SVN: r4021]]
parent adc042fd
No related branches found
No related tags found
No related merge requests found
......@@ -59,14 +59,10 @@ void checkIntersectionIterator(const GridType& grid) {
FieldVector<ctype, Geometry::coorddimension> center(0);
for (int j=0; j<intersectionGlobal.corners(); j++)
center += intersectionGlobal[j];
center /= (double) intersectionGlobal.corners();
std::cout << center << " center \n";
// The geometry center in local coordinates
FieldVector<ctype, Geometry::mydimension> localCenter = intersectionGlobal.local(center);
std::cout << localCenter << " localcenter\n";
// Back to global coordinates to check for correctness
FieldVector<ctype, Geometry::coorddimension> worldCenter = intersectionGlobal.global(localCenter);
if ((center-worldCenter).infinity_norm() > 1e-6)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment