Changeset c8c9f0 for test/interfaces/interface_sinus.cpp

Timestamp:

Apr 19, 2013, 8:41:55 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

8b2e07

Parents:

6f7f90

Message:

Added error output to sine testing interface.

File:

• test/interfaces/interface_sinus.cpp (modified) (1 diff)

test/interfaces/interface_sinus.cpp

@@ -59 +59 @@
 void InterfaceSinus::ExportSolution(Grid& grid)
 {
+  Index i;
+  vmg_float err_1 = 0.0;
+  vmg_float err_2 = 0.0;
+  vmg_float err_inf = 0.0;
+
+  const Index begin_local = grid.Global().LocalBegin() - grid.Local().HaloSize1();
+
+  for (i.X()=grid.Local().Begin().X(); i.X()<grid.Local().End().X(); ++i.X())
+    for (i.Y()=grid.Local().Begin().Y(); i.Y()<grid.Local().End().Y(); ++i.Y())
+      for (i.Z()=grid.Local().Begin().Z(); i.Z()<grid.Local().End().Z(); ++i.Z()) {
+        const Vector pos = grid.Extent().MeshWidth() * static_cast<Vector>(begin_local + i);
+        const vmg_float err = std::abs(grid.GetVal(i) - std::sin(sine_factor * pos.X()) * std::sin(sine_factor * pos.Y()) * std::sin(sine_factor * pos.Z()));
+        err_1 += err;
+        err_2 += err * err;
+        err_inf = std::max(err_inf, err);
+      }
+
+  err_1 = grid.Extent().MeshWidth().Product() * err_1;
+  err_2 = std::sqrt(grid.Extent().MeshWidth().Product() * err_2);
+  err_inf = grid.Extent().MeshWidth().Product() * err_inf;
+
+  std::cout << std::scientific << "Error L1-Norm:   " << err_1 << std::endl
+            << "Error L2-Norm:   " << err_2 << std::endl
+            << "Error Inf-Norm: " << err_inf << std::endl;
 }