From: Timo Heister Date: Wed, 1 Jun 2016 11:06:36 +0000 (+0100) Subject: use compute_global_error X-Git-Tag: v8.5.0-rc1~990^2~2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7eddafaa7a2c3a9b6265c96ba2eec8128fd4c8e5;p=dealii.git use compute_global_error --- diff --git a/examples/step-56/step-56.cc b/examples/step-56/step-56.cc index 4618406b77..912fd5a018 100644 --- a/examples/step-56/step-56.cc +++ b/examples/step-56/step-56.cc @@ -1019,7 +1019,9 @@ namespace Step56 VectorTools::L2_norm, &velocity_mask); - const double Velocity_L2_error = difference_per_cell.l2_norm(); + const double Velocity_L2_error = VectorTools::compute_global_error(triangulation, + difference_per_cell, + VectorTools::L2_norm); VectorTools::integrate_difference (dof_handler, solution, @@ -1029,7 +1031,9 @@ namespace Step56 VectorTools::L2_norm, &pressure_mask); - const double Pressure_L2_error = difference_per_cell.l2_norm(); + const double Pressure_L2_error = VectorTools::compute_global_error(triangulation, + difference_per_cell, + VectorTools::L2_norm); VectorTools::integrate_difference (dof_handler, solution, @@ -1039,7 +1043,9 @@ namespace Step56 VectorTools::H1_norm, &velocity_mask); - const double Velocity_H1_error = difference_per_cell.l2_norm(); + const double Velocity_H1_error = VectorTools::compute_global_error(triangulation, + difference_per_cell, + VectorTools::H1_norm); std::cout << std::endl << " Velocity L2 Error: " << Velocity_L2_error