Skip to content
Snippets Groups Projects
Commit 680fdbac authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

use the correct type for global coordinates

[[Imported from SVN: r6497]]
parent d88986de
No related branches found
No related tags found
No related merge requests found
...@@ -247,13 +247,13 @@ int main (int argc, char *argv[]) try ...@@ -247,13 +247,13 @@ int main (int argc, char *argv[]) try
// Interpret TargetSpace as isometrically embedded into an R^m, because this is // Interpret TargetSpace as isometrically embedded into an R^m, because this is
// how the corresponding Sobolev spaces are defined. // how the corresponding Sobolev spaces are defined.
BlockVector<TargetSpace::TangentVector> difference(referenceSolution.size()); BlockVector<TargetSpace::CoordinateType> difference(referenceSolution.size());
for (int j=0; j<referenceSolution.size(); j++) for (int j=0; j<referenceSolution.size(); j++)
difference[j] = solution[j].globalCoordinates() - referenceSolution[j].globalCoordinates(); difference[j] = solution[j].globalCoordinates() - referenceSolution[j].globalCoordinates();
H1SemiNorm< BlockVector<TargetSpace::TangentVector> > h1Norm(laplace); H1SemiNorm< BlockVector<TargetSpace::CoordinateType> > h1Norm(laplace);
H1SemiNorm< BlockVector<TargetSpace::TangentVector> > l2Norm(massMatrix); H1SemiNorm< BlockVector<TargetSpace::CoordinateType> > l2Norm(massMatrix);
// Compute max-norm difference // Compute max-norm difference
std::cout << "Level: " << i-1 std::cout << "Level: " << i-1
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment