]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Update tests
authorJean-Paul Pelteret <jppelteret@gmail.com>
Sat, 31 Jul 2021 10:46:39 +0000 (12:46 +0200)
committerJean-Paul Pelteret <jppelteret@gmail.com>
Wed, 4 Aug 2021 17:15:20 +0000 (19:15 +0200)
24 files changed:
tests/feinterface/face_orientation_01.cc
tests/feinterface/fe_interface_values_02.cc
tests/feinterface/fe_interface_values_02.output
tests/feinterface/fe_interface_values_03.cc
tests/feinterface/fe_interface_values_03.output
tests/feinterface/fe_interface_values_04.cc
tests/feinterface/fe_interface_values_04.output
tests/feinterface/fe_interface_values_06.cc
tests/feinterface/fe_interface_values_06.output
tests/feinterface/fe_interface_values_07.cc
tests/feinterface/fe_interface_values_07.output
tests/feinterface/fe_interface_values_08.cc
tests/feinterface/fe_interface_values_08.output
tests/feinterface/fe_interface_values_09.cc
tests/feinterface/fe_interface_values_09.output
tests/feinterface/fe_interface_view_01.cc
tests/feinterface/fe_interface_view_02.cc
tests/feinterface/step-12.cc
tests/meshworker/scratch_data_08.cc
tests/simplex/matrix_free_03.cc
tests/simplex/poisson_02.cc
tests/simplex/step-12.cc
tests/simplex/step-12a.cc
tests/simplex/step-74.cc

index beecec988d317647d10b7adcd92928865b010c34..3abc7ed73f303c710bdf5145544da0d18fe01b9a 100644 (file)
@@ -88,7 +88,8 @@ test()
               for (unsigned int q = 0; q < fiv.get_quadrature().size(); ++q)
                 for (unsigned int i = 0; i < fiv.n_current_interface_dofs();
                      ++i)
-                  Assert(std::abs(fiv.jump(i, q)) <= 1e-12, ExcInternalError());
+                  Assert(std::abs(fiv.jump_in_shape_values(i, q)) <= 1e-12,
+                         ExcInternalError());
             }
 
       deallog << "Twist " << twist << "*pi/2 seems to work just fine..."
index 689f2be36dbd2632aa73c7de030984c08f68b5f7..251359905d743022c62768cbbee766d4387f2bcc 100644 (file)
@@ -14,7 +14,8 @@
 // ---------------------------------------------------------------------
 
 
-// evaluate jump(), average(), shape_value() of FEInterfaceValues
+// evaluate jump_in_shape_values(), average_of_shape_values(), shape_value() of
+// FEInterfaceValues
 
 #include <deal.II/base/quadrature_lib.h>
 
@@ -114,16 +115,16 @@ test(unsigned int fe_degree)
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.jump(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "jump(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.jump_in_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "jump_in_shape_values(): " << cell_vector << std::endl;
 
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.average(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "average(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.average_of_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "average_of_shape_values(): " << cell_vector << std::endl;
       }
 }
 
index b0023c0ce52e1668ef6299b1a34821843befbdab..fe3840c85501f13500b3fbebf4ac5195c4c505ba 100644 (file)
@@ -2,20 +2,20 @@
 DEAL::FE_DGQ<2>(0)
 DEAL::shape_value(true): 1.00000 0.00000
 DEAL::shape_value(false): 0.00000 1.00000
-DEAL::jump(): 1.00000 -1.00000
-DEAL::average(): 0.500000 0.500000
+DEAL::jump_in_shape_values(): 1.00000 -1.00000
+DEAL::average_of_shape_values(): 0.500000 0.500000
 DEAL::FE_DGQ<2>(1)
 DEAL::shape_value(true): 0.00000 0.500000 0.00000 0.500000 0.00000 0.00000 0.00000 0.00000
 DEAL::shape_value(false): 0.00000 0.00000 0.00000 0.00000 0.500000 0.00000 0.500000 0.00000
-DEAL::jump(): 0.00000 0.500000 0.00000 0.500000 -0.500000 0.00000 -0.500000 0.00000
-DEAL::average(): 0.00000 0.250000 0.00000 0.250000 0.250000 0.00000 0.250000 0.00000
+DEAL::jump_in_shape_values(): 0.00000 0.500000 0.00000 0.500000 -0.500000 0.00000 -0.500000 0.00000
+DEAL::average_of_shape_values(): 0.00000 0.250000 0.00000 0.250000 0.250000 0.00000 0.250000 0.00000
 DEAL::FE_DGQ<3>(0)
 DEAL::shape_value(true): 1.00000 0.00000
 DEAL::shape_value(false): 0.00000 1.00000
-DEAL::jump(): 1.00000 -1.00000
-DEAL::average(): 0.500000 0.500000
+DEAL::jump_in_shape_values(): 1.00000 -1.00000
+DEAL::average_of_shape_values(): 0.500000 0.500000
 DEAL::FE_DGQ<3>(1)
 DEAL::shape_value(true): 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
 DEAL::shape_value(false): 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000
-DEAL::jump(): 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 -0.250000 0.00000 -0.250000 0.00000 -0.250000 0.00000 -0.250000 0.00000
-DEAL::average(): 0.00000 0.125000 0.00000 0.125000 0.00000 0.125000 0.00000 0.125000 0.125000 0.00000 0.125000 0.00000 0.125000 0.00000 0.125000 0.00000
+DEAL::jump_in_shape_values(): 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 -0.250000 0.00000 -0.250000 0.00000 -0.250000 0.00000 -0.250000 0.00000
+DEAL::average_of_shape_values(): 0.00000 0.125000 0.00000 0.125000 0.00000 0.125000 0.00000 0.125000 0.125000 0.00000 0.125000 0.00000 0.125000 0.00000 0.125000 0.00000
index b241abdf1b0c25e5a09cb0be006a155239596ec1..b0823148d7788fc016d6dbeed337fa5f8875aa01 100644 (file)
@@ -14,8 +14,8 @@
 // ---------------------------------------------------------------------
 
 
-// evaluate jump(), average(), shape_value() of FEInterfaceValues on an adaptive
-// mesh
+// evaluate jump_in_shape_values(), average_of_shape_values(), shape_value() of
+// FEInterfaceValues on an adaptive mesh
 
 #include <deal.II/base/quadrature_lib.h>
 
@@ -133,16 +133,16 @@ test(unsigned int fe_degree)
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.jump(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "jump(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.jump_in_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "jump_in_shape_values(): " << cell_vector << std::endl;
 
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.average(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "average(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.average_of_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "average_of_shape_values(): " << cell_vector << std::endl;
       }
 }
 
index 477b2dea180aa31139be6d2dc8a091cd467045b1..cfc64df654dd5ec0595172de27d8d7ba78109819 100644 (file)
@@ -5,8 +5,8 @@ DEAL::  idx: 0 global: 0 dof indices: -1 | 0
 DEAL::  idx: 1 global: 2 dof indices: 0 | -1
 DEAL::shape_value(true): 0.00000 0.500000
 DEAL::shape_value(false): 0.500000 0.00000
-DEAL::jump(): -0.500000 0.500000
-DEAL::average(): 0.250000 0.250000
+DEAL::jump_in_shape_values(): -0.500000 0.500000
+DEAL::average_of_shape_values(): 0.250000 0.250000
 DEAL::FE_DGQ<2>(1)
 DEAL::qpoint 0: 1.00000 0.105662
 DEAL::qpoint 1: 1.00000 0.394338
@@ -20,16 +20,16 @@ DEAL::  idx: 6 global: 10 dof indices: 2 | -1
 DEAL::  idx: 7 global: 11 dof indices: 3 | -1
 DEAL::shape_value(true): 0.00000 0.00000 0.00000 0.00000 0.00000 0.250000 0.00000 0.250000
 DEAL::shape_value(false): 0.375000 0.00000 0.125000 0.00000 0.00000 0.00000 0.00000 0.00000
-DEAL::jump(): -0.375000 0.00000 -0.125000 0.00000 0.00000 0.250000 0.00000 0.250000
-DEAL::average(): 0.187500 0.00000 0.0625000 0.00000 0.00000 0.125000 0.00000 0.125000
+DEAL::jump_in_shape_values(): -0.375000 0.00000 -0.125000 0.00000 0.00000 0.250000 0.00000 0.250000
+DEAL::average_of_shape_values(): 0.187500 0.00000 0.0625000 0.00000 0.00000 0.125000 0.00000 0.125000
 DEAL::FE_DGQ<3>(0)
 DEAL::qpoint 0: 1.00000 0.250000 0.250000
 DEAL::  idx: 0 global: 0 dof indices: -1 | 0
 DEAL::  idx: 1 global: 2 dof indices: 0 | -1
 DEAL::shape_value(true): 0.00000 0.250000
 DEAL::shape_value(false): 0.250000 0.00000
-DEAL::jump(): -0.250000 0.250000
-DEAL::average(): 0.125000 0.125000
+DEAL::jump_in_shape_values(): -0.250000 0.250000
+DEAL::average_of_shape_values(): 0.125000 0.125000
 DEAL::FE_DGQ<3>(1)
 DEAL::qpoint 0: 1.00000 0.105662 0.105662
 DEAL::qpoint 1: 1.00000 0.394338 0.105662
@@ -53,5 +53,5 @@ DEAL::  idx: 14 global: 22 dof indices: 6 | -1
 DEAL::  idx: 15 global: 23 dof indices: 7 | -1
 DEAL::shape_value(true): 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.0625000 0.00000 0.0625000 0.00000 0.0625000 0.00000 0.0625000
 DEAL::shape_value(false): 0.140625 0.00000 0.0468750 0.00000 0.0468750 0.00000 0.0156250 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
-DEAL::jump(): -0.140625 0.00000 -0.0468750 0.00000 -0.0468750 0.00000 -0.0156250 0.00000 0.00000 0.0625000 0.00000 0.0625000 0.00000 0.0625000 0.00000 0.0625000
-DEAL::average(): 0.0703125 0.00000 0.0234375 0.00000 0.0234375 0.00000 0.00781250 0.00000 0.00000 0.0312500 0.00000 0.0312500 0.00000 0.0312500 0.00000 0.0312500
+DEAL::jump_in_shape_values(): -0.140625 0.00000 -0.0468750 0.00000 -0.0468750 0.00000 -0.0156250 0.00000 0.00000 0.0625000 0.00000 0.0625000 0.00000 0.0625000 0.00000 0.0625000
+DEAL::average_of_shape_values(): 0.0703125 0.00000 0.0234375 0.00000 0.0234375 0.00000 0.00781250 0.00000 0.00000 0.0312500 0.00000 0.0312500 0.00000 0.0312500 0.00000 0.0312500
index bfb923052f5834ac2d6326195fff6e619455186c..c92f4e481423d561cef5040aaffefe1878d921c8 100644 (file)
@@ -14,8 +14,8 @@
 // ---------------------------------------------------------------------
 
 
-// evaluate jump(), average(), shape_value() of FEInterfaceValues for
-// at_boundary()
+// evaluate jump_in_shape_values(), average_of_shape_values(), shape_value() of
+// FEInterfaceValues for at_boundary()
 
 #include <deal.II/base/quadrature_lib.h>
 
@@ -81,14 +81,16 @@ test(unsigned int fe_degree)
   cell_vector = 0.0;
   for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
     for (unsigned int i = 0; i < n_dofs; ++i)
-      cell_vector(i) += fiv.jump(i, qpoint) * fiv.get_JxW_values()[qpoint];
-  deallog << "jump(): " << cell_vector << std::endl;
+      cell_vector(i) +=
+        fiv.jump_in_shape_values(i, qpoint) * fiv.get_JxW_values()[qpoint];
+  deallog << "jump_in_shape_values(): " << cell_vector << std::endl;
 
   cell_vector = 0.0;
   for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
     for (unsigned int i = 0; i < n_dofs; ++i)
-      cell_vector(i) += fiv.average(i, qpoint) * fiv.get_JxW_values()[qpoint];
-  deallog << "average(): " << cell_vector << std::endl;
+      cell_vector(i) +=
+        fiv.average_of_shape_values(i, qpoint) * fiv.get_JxW_values()[qpoint];
+  deallog << "average_of_shape_values(): " << cell_vector << std::endl;
 }
 
 
index 3da4157a11198015fbf259b4bd6bb848ed49ed87..1e25c01888bac573886fac61acfbbec69871da2b 100644 (file)
@@ -2,20 +2,20 @@
 DEAL::FE_DGQ<2>(0)
 DEAL::shape_value(true): 1.00000
 DEAL::shape_value(false): 0.00000
-DEAL::jump(): 1.00000
-DEAL::average(): 1.00000
+DEAL::jump_in_shape_values(): 1.00000
+DEAL::average_of_shape_values(): 1.00000
 DEAL::FE_DGQ<2>(1)
 DEAL::shape_value(true): 0.500000 0.500000 0.00000 0.00000
 DEAL::shape_value(false): 0.00000 0.00000 0.00000 0.00000
-DEAL::jump(): 0.500000 0.500000 0.00000 0.00000
-DEAL::average(): 0.500000 0.500000 0.00000 0.00000
+DEAL::jump_in_shape_values(): 0.500000 0.500000 0.00000 0.00000
+DEAL::average_of_shape_values(): 0.500000 0.500000 0.00000 0.00000
 DEAL::FE_DGQ<3>(0)
 DEAL::shape_value(true): 1.00000
 DEAL::shape_value(false): 0.00000
-DEAL::jump(): 1.00000
-DEAL::average(): 1.00000
+DEAL::jump_in_shape_values(): 1.00000
+DEAL::average_of_shape_values(): 1.00000
 DEAL::FE_DGQ<3>(1)
 DEAL::shape_value(true): 0.250000 0.250000 0.00000 0.00000 0.250000 0.250000 0.00000 0.00000
 DEAL::shape_value(false): 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
-DEAL::jump(): 0.250000 0.250000 0.00000 0.00000 0.250000 0.250000 0.00000 0.00000
-DEAL::average(): 0.250000 0.250000 0.00000 0.00000 0.250000 0.250000 0.00000 0.00000
+DEAL::jump_in_shape_values(): 0.250000 0.250000 0.00000 0.00000 0.250000 0.250000 0.00000 0.00000
+DEAL::average_of_shape_values(): 0.250000 0.250000 0.00000 0.00000 0.250000 0.250000 0.00000 0.00000
index edcdda30ada01421ff6a939d5e46d1795afe3fb2..6edda1755eb68160205256bb0d96b230423103b5 100644 (file)
@@ -14,8 +14,8 @@
 // ---------------------------------------------------------------------
 
 
-// evaluate jump(), average(), shape_value() of FEInterfaceValues on an adaptive
-// mesh for a continuous element
+// evaluate jump_in_shape_values(), average_of_shape_values(), shape_value() of
+// FEInterfaceValues on an adaptive mesh for a continuous element
 
 #include <deal.II/base/quadrature_lib.h>
 
@@ -134,16 +134,16 @@ test(unsigned int fe_degree)
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.jump(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "jump(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.jump_in_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "jump_in_shape_values(): " << cell_vector << std::endl;
 
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.average(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "average(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.average_of_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "average_of_shape_values(): " << cell_vector << std::endl;
       }
 }
 
index 6488c7fc32d79f1dc47562a28f6a12b6b83a93b5..1bd6762b95e1bdc40dfc00dcb12e66cec171b598 100644 (file)
@@ -11,8 +11,8 @@ DEAL::  idx: 5 global: 7 dof indices: 2 | -1
 DEAL::  idx: 6 global: 8 dof indices: 3 | -1
 DEAL::shape_value(true): 0.250000 0.00000 0.00000 0.00000 0.00000 0.00000 0.250000
 DEAL::shape_value(false): 0.375000 0.00000 0.125000 0.00000 0.00000 0.00000 0.00000
-DEAL::jump(): -0.125000 0.00000 -0.125000 0.00000 0.00000 0.00000 0.250000
-DEAL::average(): 0.312500 0.00000 0.0625000 0.00000 0.00000 0.00000 0.125000
+DEAL::jump_in_shape_values(): -0.125000 0.00000 -0.125000 0.00000 0.00000 0.00000 0.250000
+DEAL::average_of_shape_values(): 0.312500 0.00000 0.0625000 0.00000 0.00000 0.00000 0.125000
 DEAL::FE_Q<3>(1)
 DEAL::qpoint 0: 1.00000 0.105662 0.105662
 DEAL::qpoint 1: 1.00000 0.394338 0.105662
@@ -35,5 +35,5 @@ DEAL::  idx: 13 global: 17 dof indices: 5 | -1
 DEAL::  idx: 14 global: 18 dof indices: 7 | -1
 DEAL::shape_value(true): 0.0625000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.0625000 0.0625000 0.0625000
 DEAL::shape_value(false): 0.140625 0.00000 0.0468750 0.00000 0.0468750 0.00000 0.0156250 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
-DEAL::jump(): -0.0781250 0.00000 -0.0468750 0.00000 -0.0468750 0.00000 -0.0156250 0.00000 0.00000 0.00000 0.00000 0.00000 0.0625000 0.0625000 0.0625000
-DEAL::average(): 0.101562 0.00000 0.0234375 0.00000 0.0234375 0.00000 0.00781250 0.00000 0.00000 0.00000 0.00000 0.00000 0.0312500 0.0312500 0.0312500
+DEAL::jump_in_shape_values(): -0.0781250 0.00000 -0.0468750 0.00000 -0.0468750 0.00000 -0.0156250 0.00000 0.00000 0.00000 0.00000 0.00000 0.0625000 0.0625000 0.0625000
+DEAL::average_of_shape_values(): 0.101562 0.00000 0.0234375 0.00000 0.0234375 0.00000 0.00781250 0.00000 0.00000 0.00000 0.00000 0.00000 0.0312500 0.0312500 0.0312500
index f38bb06ecf5245d3dd15d3dc5ad78d689f1ab961..10ede02e2ce6ae9d7042f83a2e19e820e1058560 100644 (file)
@@ -14,8 +14,8 @@
 // ---------------------------------------------------------------------
 
 
-// evaluate jump(), average(), shape_value() of FEInterfaceValues
-// like fe_interface_values_02, but for a continuous element
+// evaluate jump_in_shape_values(), average_of_shape_values(), shape_value() of
+// FEInterfaceValues like fe_interface_values_02, but for a continuous element
 
 #include <deal.II/base/quadrature_lib.h>
 
@@ -116,16 +116,16 @@ test(unsigned int fe_degree)
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.jump(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "jump(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.jump_in_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "jump_in_shape_values(): " << cell_vector << std::endl;
 
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) +=
-              fiv.average(i, qpoint) * fiv.get_JxW_values()[qpoint];
-        deallog << "average(): " << cell_vector << std::endl;
+            cell_vector(i) += fiv.average_of_shape_values(i, qpoint) *
+                              fiv.get_JxW_values()[qpoint];
+        deallog << "average_of_shape_values(): " << cell_vector << std::endl;
       }
 }
 
index 2cae4be57f66e6e46c1dc1be407e41de84f10c87..b30862e2b338f4254ad75d73af1c9da2bd823cbe 100644 (file)
@@ -2,10 +2,10 @@
 DEAL::FE_Q<2>(1)
 DEAL::shape_value(true): 0.00000 0.500000 0.00000 0.500000 0.00000 0.00000
 DEAL::shape_value(false): 0.00000 0.500000 0.00000 0.500000 0.00000 0.00000
-DEAL::jump(): 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
-DEAL::average(): 0.00000 0.500000 0.00000 0.500000 0.00000 0.00000
+DEAL::jump_in_shape_values(): 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
+DEAL::average_of_shape_values(): 0.00000 0.500000 0.00000 0.500000 0.00000 0.00000
 DEAL::FE_Q<3>(1)
 DEAL::shape_value(true): 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.00000 0.00000 0.00000
 DEAL::shape_value(false): 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.00000 0.00000 0.00000
-DEAL::jump(): 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
-DEAL::average(): 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.00000 0.00000 0.00000
+DEAL::jump_in_shape_values(): 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
+DEAL::average_of_shape_values(): 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.250000 0.00000 0.00000 0.00000 0.00000
index cd448e3379c4510b7b688c1d1f811c2746fe5d6d..d751b91c58fbfb8a609c8f3b533739b42f9a1efc 100644 (file)
@@ -14,7 +14,8 @@
 // ---------------------------------------------------------------------
 
 
-// evaluate average_gradient average_hessian jump_gradient of FEInterfaceValues
+// evaluate average_gradient average_of_shape_hessians( jump_gradient of
+// FEInterfaceValues
 
 #include <deal.II/base/quadrature_lib.h>
 
@@ -97,23 +98,24 @@ test(const FiniteElement<dim> &fe)
     cell_vector = 0.0;
     for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
       for (unsigned int i = 0; i < n_dofs; ++i)
-        cell_vector(i) +=
-          fiv.average_gradient(i, qpoint).norm() * fiv.get_JxW_values()[qpoint];
-    deallog << "average_gradient.norm(): " << cell_vector << std::endl;
+        cell_vector(i) += fiv.average_of_shape_gradients(i, qpoint).norm() *
+                          fiv.get_JxW_values()[qpoint];
+    deallog << "average_of_shape_gradients.norm(): " << cell_vector
+            << std::endl;
 
     cell_vector = 0.0;
     for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
       for (unsigned int i = 0; i < n_dofs; ++i)
-        cell_vector(i) +=
-          fiv.average_hessian(i, qpoint).norm() * fiv.get_JxW_values()[qpoint];
-    deallog << "average_hessian.norm(): " << cell_vector << std::endl;
+        cell_vector(i) += fiv.average_of_shape_hessians(i, qpoint).norm() *
+                          fiv.get_JxW_values()[qpoint];
+    deallog << "average_of_shape_hessians.norm(): " << cell_vector << std::endl;
 
     cell_vector = 0.0;
     for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
       for (unsigned int i = 0; i < n_dofs; ++i)
-        cell_vector(i) +=
-          fiv.jump_gradient(i, qpoint).norm() * fiv.get_JxW_values()[qpoint];
-    deallog << "jump_gradient.norm(): " << cell_vector << std::endl;
+        cell_vector(i) += fiv.jump_in_shape_gradients(i, qpoint).norm() *
+                          fiv.get_JxW_values()[qpoint];
+    deallog << "jump_in_shape_gradients.norm(): " << cell_vector << std::endl;
   }
 
 
@@ -135,23 +137,26 @@ test(const FiniteElement<dim> &fe)
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) += fiv.average_gradient(i, qpoint).norm() *
+            cell_vector(i) += fiv.average_of_shape_gradients(i, qpoint).norm() *
                               fiv.get_JxW_values()[qpoint];
-        deallog << "average_gradient.norm(): " << cell_vector << std::endl;
+        deallog << "average_of_shape_gradients.norm(): " << cell_vector
+                << std::endl;
 
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) += fiv.average_hessian(i, qpoint).norm() *
+            cell_vector(i) += fiv.average_of_shape_hessians(i, qpoint).norm() *
                               fiv.get_JxW_values()[qpoint];
-        deallog << "average_hessian.norm(): " << cell_vector << std::endl;
+        deallog << "average_of_shape_hessians.norm(): " << cell_vector
+                << std::endl;
 
         cell_vector = 0.0;
         for (unsigned int qpoint = 0; qpoint < q_points.size(); ++qpoint)
           for (unsigned int i = 0; i < n_dofs; ++i)
-            cell_vector(i) += fiv.jump_gradient(i, qpoint).norm() *
+            cell_vector(i) += fiv.jump_in_shape_gradients(i, qpoint).norm() *
                               fiv.get_JxW_values()[qpoint];
-        deallog << "jump_gradient.norm(): " << cell_vector << std::endl;
+        deallog << "jump_in_shape_gradients.norm(): " << cell_vector
+                << std::endl;
       }
 }
 
index f3ef235513cf5ed29f61a9807ebd2300988c52b8..7c4e78a3c6ee472e55430fe77db8c64369485c94 100644 (file)
@@ -1,50 +1,50 @@
 
 DEAL::FE_Q<2>(1)
-DEAL::average_gradient.norm(): 1.14783 0.500000 1.14783 0.500000
-DEAL::average_hessian.norm(): 1.41421 1.41421 1.41421 1.41421
-DEAL::jump_gradient.norm(): 1.14783 0.500000 1.14783 0.500000
-DEAL::average_gradient.norm(): 0.250000 1.00000 0.250000 1.00000 0.250000 0.250000
-DEAL::average_hessian.norm(): 0.707107 0.00000 0.707107 0.00000 0.707107 0.707107
-DEAL::jump_gradient.norm(): 0.500000 1.00000 0.500000 1.00000 0.500000 0.500000
+DEAL::average_of_shape_gradients.norm(): 1.14783 0.500000 1.14783 0.500000
+DEAL::average_of_shape_hessians.norm(): 1.41421 1.41421 1.41421 1.41421
+DEAL::jump_in_shape_gradients.norm(): 1.14783 0.500000 1.14783 0.500000
+DEAL::average_of_shape_gradients.norm(): 0.250000 1.00000 0.250000 1.00000 0.250000 0.250000
+DEAL::average_of_shape_hessians.norm(): 0.707107 0.00000 0.707107 0.00000 0.707107 0.707107
+DEAL::jump_in_shape_gradients.norm(): 0.500000 1.00000 0.500000 1.00000 0.500000 0.500000
 DEAL::FE_DGQ<2>(0)
-DEAL::average_gradient.norm(): 0.00000
-DEAL::average_hessian.norm(): 0.00000
-DEAL::jump_gradient.norm(): 0.00000
-DEAL::average_gradient.norm(): 0.00000 0.00000
-DEAL::average_hessian.norm(): 0.00000 0.00000
-DEAL::jump_gradient.norm(): 0.00000 0.00000
+DEAL::average_of_shape_gradients.norm(): 0.00000
+DEAL::average_of_shape_hessians.norm(): 0.00000
+DEAL::jump_in_shape_gradients.norm(): 0.00000
+DEAL::average_of_shape_gradients.norm(): 0.00000 0.00000
+DEAL::average_of_shape_hessians.norm(): 0.00000 0.00000
+DEAL::jump_in_shape_gradients.norm(): 0.00000 0.00000
 DEAL::FE_DGQ<2>(1)
-DEAL::average_gradient.norm(): 1.14783 0.500000 1.14783 0.500000
-DEAL::average_hessian.norm(): 1.41421 1.41421 1.41421 1.41421
-DEAL::jump_gradient.norm(): 1.14783 0.500000 1.14783 0.500000
-DEAL::average_gradient.norm(): 0.250000 0.573917 0.250000 0.573917 0.573917 0.250000 0.573917 0.250000
-DEAL::average_hessian.norm(): 0.707107 0.707107 0.707107 0.707107 0.707107 0.707107 0.707107 0.707107
-DEAL::jump_gradient.norm(): 0.500000 1.14783 0.500000 1.14783 1.14783 0.500000 1.14783 0.500000
+DEAL::average_of_shape_gradients.norm(): 1.14783 0.500000 1.14783 0.500000
+DEAL::average_of_shape_hessians.norm(): 1.41421 1.41421 1.41421 1.41421
+DEAL::jump_in_shape_gradients.norm(): 1.14783 0.500000 1.14783 0.500000
+DEAL::average_of_shape_gradients.norm(): 0.250000 0.573917 0.250000 0.573917 0.573917 0.250000 0.573917 0.250000
+DEAL::average_of_shape_hessians.norm(): 0.707107 0.707107 0.707107 0.707107 0.707107 0.707107 0.707107 0.707107
+DEAL::jump_in_shape_gradients.norm(): 0.500000 1.14783 0.500000 1.14783 1.14783 0.500000 1.14783 0.500000
 DEAL::FE_DGQ<2>(2)
-DEAL::average_gradient.norm(): 1.52420 0.860663 0.215166 3.17925 2.66667 0.666667 1.52420 0.860663 0.215166
-DEAL::average_hessian.norm(): 7.17398 7.68564 2.12446 12.5704 13.4538 4.36931 7.17398 7.68564 2.12446
-DEAL::jump_gradient.norm(): 1.52420 0.860663 0.215166 3.17925 2.66667 0.666667 1.52420 0.860663 0.215166
-DEAL::average_gradient.norm(): 0.107583 0.430331 0.762102 0.333333 1.33333 1.58962 0.107583 0.430331 0.762102 0.762102 0.430331 0.107583 1.58962 1.33333 0.333333 0.762102 0.430331 0.107583
-DEAL::average_hessian.norm(): 1.06223 3.84282 3.58699 2.18466 6.72690 6.28519 1.06223 3.84282 3.58699 3.58699 3.84282 1.06223 6.28519 6.72690 2.18466 3.58699 3.84282 1.06223
-DEAL::jump_gradient.norm(): 0.215166 0.860663 1.52420 0.666667 2.66667 3.17925 0.215166 0.860663 1.52420 1.52420 0.860663 0.215166 3.17925 2.66667 0.666667 1.52420 0.860663 0.215166
+DEAL::average_of_shape_gradients.norm(): 1.52420 0.860663 0.215166 3.17925 2.66667 0.666667 1.52420 0.860663 0.215166
+DEAL::average_of_shape_hessians.norm(): 7.17398 7.68564 2.12446 12.5704 13.4538 4.36931 7.17398 7.68564 2.12446
+DEAL::jump_in_shape_gradients.norm(): 1.52420 0.860663 0.215166 3.17925 2.66667 0.666667 1.52420 0.860663 0.215166
+DEAL::average_of_shape_gradients.norm(): 0.107583 0.430331 0.762102 0.333333 1.33333 1.58962 0.107583 0.430331 0.762102 0.762102 0.430331 0.107583 1.58962 1.33333 0.333333 0.762102 0.430331 0.107583
+DEAL::average_of_shape_hessians.norm(): 1.06223 3.84282 3.58699 2.18466 6.72690 6.28519 1.06223 3.84282 3.58699 3.58699 3.84282 1.06223 6.28519 6.72690 2.18466 3.58699 3.84282 1.06223
+DEAL::jump_in_shape_gradients.norm(): 0.215166 0.860663 1.52420 0.666667 2.66667 3.17925 0.215166 0.860663 1.52420 1.52420 0.860663 0.215166 3.17925 2.66667 0.666667 1.52420 0.860663 0.215166
 DEAL::FE_Q<3>(1)
-DEAL::average_gradient.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
-DEAL::average_hessian.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
-DEAL::jump_gradient.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
-DEAL::average_gradient.norm(): 0.125000 0.761802 0.125000 0.761802 0.125000 0.761802 0.125000 0.761802 0.125000 0.125000 0.125000 0.125000
-DEAL::average_hessian.norm(): 0.538675 1.41421 0.538675 1.41421 0.538675 1.41421 0.538675 1.41421 0.538675 0.538675 0.538675 0.538675
-DEAL::jump_gradient.norm(): 0.250000 0.500000 0.250000 0.500000 0.250000 0.500000 0.250000 0.500000 0.250000 0.250000 0.250000 0.250000
+DEAL::average_of_shape_gradients.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
+DEAL::average_of_shape_hessians.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
+DEAL::jump_in_shape_gradients.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
+DEAL::average_of_shape_gradients.norm(): 0.125000 0.761802 0.125000 0.761802 0.125000 0.761802 0.125000 0.761802 0.125000 0.125000 0.125000 0.125000
+DEAL::average_of_shape_hessians.norm(): 0.538675 1.41421 0.538675 1.41421 0.538675 1.41421 0.538675 1.41421 0.538675 0.538675 0.538675 0.538675
+DEAL::jump_in_shape_gradients.norm(): 0.250000 0.500000 0.250000 0.500000 0.250000 0.500000 0.250000 0.500000 0.250000 0.250000 0.250000 0.250000
 DEAL::FE_DGQ<3>(0)
-DEAL::average_gradient.norm(): 0.00000
-DEAL::average_hessian.norm(): 0.00000
-DEAL::jump_gradient.norm(): 0.00000
-DEAL::average_gradient.norm(): 0.00000 0.00000
-DEAL::average_hessian.norm(): 0.00000 0.00000
-DEAL::jump_gradient.norm(): 0.00000 0.00000
+DEAL::average_of_shape_gradients.norm(): 0.00000
+DEAL::average_of_shape_hessians.norm(): 0.00000
+DEAL::jump_in_shape_gradients.norm(): 0.00000
+DEAL::average_of_shape_gradients.norm(): 0.00000 0.00000
+DEAL::average_of_shape_hessians.norm(): 0.00000 0.00000
+DEAL::jump_in_shape_gradients.norm(): 0.00000 0.00000
 DEAL::FE_DGQ<3>(1)
-DEAL::average_gradient.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
-DEAL::average_hessian.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
-DEAL::jump_gradient.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
-DEAL::average_gradient.norm(): 0.125000 0.405739 0.125000 0.405739 0.125000 0.405739 0.125000 0.405739 0.405739 0.125000 0.405739 0.125000 0.405739 0.125000 0.405739 0.125000
-DEAL::average_hessian.norm(): 0.538675 0.905750 0.538675 0.905750 0.538675 0.905750 0.538675 0.905750 0.905750 0.538675 0.905750 0.538675 0.905750 0.538675 0.905750 0.538675
-DEAL::jump_gradient.norm(): 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
+DEAL::average_of_shape_gradients.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
+DEAL::average_of_shape_hessians.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
+DEAL::jump_in_shape_gradients.norm(): 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
+DEAL::average_of_shape_gradients.norm(): 0.125000 0.405739 0.125000 0.405739 0.125000 0.405739 0.125000 0.405739 0.405739 0.125000 0.405739 0.125000 0.405739 0.125000 0.405739 0.125000
+DEAL::average_of_shape_hessians.norm(): 0.538675 0.905750 0.538675 0.905750 0.538675 0.905750 0.538675 0.905750 0.905750 0.538675 0.905750 0.538675 0.905750 0.538675 0.905750 0.538675
+DEAL::jump_in_shape_gradients.norm(): 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000 0.811479 0.250000
index dd114c7ce74746da0caff02ba8b1d3054eb45050..6478c5a5b3de63137282776705bf078b4b3d738a 100644 (file)
@@ -74,16 +74,17 @@ print_norm_of_average_over_quadrature_points(const FEInterfaceValues<dim> &fiv)
   cell_vector = 0.0;
   for (unsigned int qpoint = 0; qpoint < fiv.n_quadrature_points; ++qpoint)
     for (unsigned int i = 0; i < n_dofs; ++i)
-      cell_vector(i) +=
-        fiv.jump_hessian(i, qpoint).norm() * fiv.get_JxW_values()[qpoint];
-  deallog << "jump_hessian.norm(): " << cell_vector << std::endl;
+      cell_vector(i) += fiv.jump_in_shape_hessians(i, qpoint).norm() *
+                        fiv.get_JxW_values()[qpoint];
+  deallog << "jump_in_shape_hessians.norm(): " << cell_vector << std::endl;
 
   cell_vector = 0.0;
   for (unsigned int qpoint = 0; qpoint < fiv.n_quadrature_points; ++qpoint)
     for (unsigned int i = 0; i < n_dofs; ++i)
-      cell_vector(i) += fiv.jump_3rd_derivative(i, qpoint).norm() *
+      cell_vector(i) += fiv.jump_in_shape_3rd_derivatives(i, qpoint).norm() *
                         fiv.get_JxW_values()[qpoint];
-  deallog << "jump_3rd_derivative.norm(): " << cell_vector << std::endl;
+  deallog << "jump_in_shape_3rd_derivatives.norm(): " << cell_vector
+          << std::endl;
 }
 
 
index bbe33191706eee39c7d71caa6d98a14564b014ff..f68424300c6a27bde09c42eaef20a33209a74f69 100644 (file)
@@ -1,21 +1,21 @@
 
 DEAL::FE_Q<2>(2)
-DEAL::jump_hessian.norm(): 7.17398 2.12446 7.17398 2.12446 12.5704 4.36931 7.68564 7.68564 13.4538
-DEAL::jump_3rd_derivative.norm(): 23.1831 11.8202 23.1831 11.8202 44.4667 20.3528 34.2248 34.2248 63.5828
-DEAL::jump_hessian.norm(): 2.12446 11.0742 2.12446 11.0742 4.36931 14.6059 7.68564 7.68564 13.4538 2.12446 2.12446 4.36931 7.68564 7.68564 13.4538
-DEAL::jump_3rd_derivative.norm(): 11.8202 41.5692 11.8202 41.5692 20.3528 83.1384 34.2248 34.2248 63.5828 11.8202 11.8202 20.3528 34.2248 34.2248 63.5828
+DEAL::jump_in_shape_hessians.norm(): 7.17398 2.12446 7.17398 2.12446 12.5704 4.36931 7.68564 7.68564 13.4538
+DEAL::jump_in_shape_3rd_derivatives.norm(): 23.1831 11.8202 23.1831 11.8202 44.4667 20.3528 34.2248 34.2248 63.5828
+DEAL::jump_in_shape_hessians.norm(): 2.12446 11.0742 2.12446 11.0742 4.36931 14.6059 7.68564 7.68564 13.4538 2.12446 2.12446 4.36931 7.68564 7.68564 13.4538
+DEAL::jump_in_shape_3rd_derivatives.norm(): 11.8202 41.5692 11.8202 41.5692 20.3528 83.1384 34.2248 34.2248 63.5828 11.8202 11.8202 20.3528 34.2248 34.2248 63.5828
 DEAL::FE_DGQ<2>(2)
-DEAL::jump_hessian.norm(): 7.17398 7.68564 2.12446 12.5704 13.4538 4.36931 7.17398 7.68564 2.12446
-DEAL::jump_3rd_derivative.norm(): 23.1831 34.2248 11.8202 44.4667 63.5828 20.3528 23.1831 34.2248 11.8202
-DEAL::jump_hessian.norm(): 2.12446 7.68564 7.17398 4.36931 13.4538 12.5704 2.12446 7.68564 7.17398 7.17398 7.68564 2.12446 12.5704 13.4538 4.36931 7.17398 7.68564 2.12446
-DEAL::jump_3rd_derivative.norm(): 11.8202 34.2248 23.1831 20.3528 63.5828 44.4667 11.8202 34.2248 23.1831 23.1831 34.2248 11.8202 44.4667 63.5828 20.3528 23.1831 34.2248 11.8202
+DEAL::jump_in_shape_hessians.norm(): 7.17398 7.68564 2.12446 12.5704 13.4538 4.36931 7.17398 7.68564 2.12446
+DEAL::jump_in_shape_3rd_derivatives.norm(): 23.1831 34.2248 11.8202 44.4667 63.5828 20.3528 23.1831 34.2248 11.8202
+DEAL::jump_in_shape_hessians.norm(): 2.12446 7.68564 7.17398 4.36931 13.4538 12.5704 2.12446 7.68564 7.17398 7.17398 7.68564 2.12446 12.5704 13.4538 4.36931 7.17398 7.68564 2.12446
+DEAL::jump_in_shape_3rd_derivatives.norm(): 11.8202 34.2248 23.1831 20.3528 63.5828 44.4667 11.8202 34.2248 23.1831 23.1831 34.2248 11.8202 44.4667 63.5828 20.3528 23.1831 34.2248 11.8202
 DEAL::FE_Q<3>(1)
-DEAL::jump_hessian.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
-DEAL::jump_3rd_derivative.norm(): 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949
-DEAL::jump_hessian.norm(): 1.07735 2.15470 1.07735 2.15470 1.07735 2.15470 1.07735 2.15470 1.07735 1.07735 1.07735 1.07735
-DEAL::jump_3rd_derivative.norm(): 2.44949 4.89898 2.44949 4.89898 2.44949 4.89898 2.44949 4.89898 2.44949 2.44949 2.44949 2.44949
+DEAL::jump_in_shape_hessians.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
+DEAL::jump_in_shape_3rd_derivatives.norm(): 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949
+DEAL::jump_in_shape_hessians.norm(): 1.07735 2.15470 1.07735 2.15470 1.07735 2.15470 1.07735 2.15470 1.07735 1.07735 1.07735 1.07735
+DEAL::jump_in_shape_3rd_derivatives.norm(): 2.44949 4.89898 2.44949 4.89898 2.44949 4.89898 2.44949 4.89898 2.44949 2.44949 2.44949 2.44949
 DEAL::FE_DGQ<3>(1)
-DEAL::jump_hessian.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
-DEAL::jump_3rd_derivative.norm(): 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949
-DEAL::jump_hessian.norm(): 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
-DEAL::jump_3rd_derivative.norm(): 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949
+DEAL::jump_in_shape_hessians.norm(): 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
+DEAL::jump_in_shape_3rd_derivatives.norm(): 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949
+DEAL::jump_in_shape_hessians.norm(): 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735 1.81150 1.07735
+DEAL::jump_in_shape_3rd_derivatives.norm(): 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949 2.44949
index 748d7b39cbc32733b39a5f64e1e0a19ab710ac84..ca1541c3423ee2678bbf01965a54fa123dccd4c9 100644 (file)
@@ -79,17 +79,22 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
                   if (fe_iv[single_component].value(false, i, q) != 0)
                     deallog << fe_iv[single_component].value(false, i, q)
                             << "  ";
-                  if (fe_iv[single_component].jump(i, q) != 0)
-                    deallog << fe_iv[single_component].jump(i, q) << "  ";
-                  if (fe_iv[single_component].average(i, q) != 0)
-                    deallog << fe_iv[single_component].average(i, q) << "  ";
-                  if (fe_iv[single_component].jump_gradient(i, q).norm() != 0)
-                    deallog << fe_iv[single_component].jump_gradient(i, q)
+                  if (fe_iv[single_component].jump_in_values(i, q) != 0)
+                    deallog << fe_iv[single_component].jump_in_values(i, q)
                             << "  ";
-                  if (fe_iv[single_component].average_gradient(i, q).norm() !=
+                  if (fe_iv[single_component].average_of_values(i, q) != 0)
+                    deallog << fe_iv[single_component].average_of_values(i, q)
+                            << "  ";
+                  if (fe_iv[single_component].jump_in_gradients(i, q).norm() !=
                       0)
-                    deallog << fe_iv[single_component].average_gradient(i, q)
-                            << std::endl;
+                    deallog << fe_iv[single_component].jump_in_gradients(i, q)
+                            << "  ";
+                  if (fe_iv[single_component]
+                        .average_of_gradients(i, q)
+                        .norm() != 0)
+                    deallog
+                      << fe_iv[single_component].average_of_gradients(i, q)
+                      << std::endl;
 
 
                   Assert(fe_iv[single_component].value(true, i, q) ==
@@ -98,25 +103,24 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
                   Assert(fe_iv[single_component].value(false, i, q) ==
                            fe_iv.shape_value(false, i, q, c),
                          ExcInternalError());
-                  Assert(fe_iv[single_component].jump(i, q) ==
-                           fe_iv.jump(i, q, c),
+                  Assert(fe_iv[single_component].jump_in_values(i, q) ==
+                           fe_iv.jump_in_shape_values(i, q, c),
                          ExcInternalError());
-                  Assert(fe_iv[single_component].average(i, q) ==
-                           fe_iv.average(i, q, c),
+                  Assert(fe_iv[single_component].average_of_values(i, q) ==
+                           fe_iv.average_of_shape_values(i, q, c),
                          ExcInternalError());
-                  Assert(fe_iv[single_component].average_hessian(i, q) ==
-                           fe_iv.average_hessian(i, q, c),
+                  Assert(fe_iv[single_component].average_of_hessians(i, q) ==
+                           fe_iv.average_of_shape_hessians(i, q, c),
                          ExcInternalError());
-                  Assert(fe_iv[single_component].jump_hessian(i, q) ==
-                           fe_iv.jump_hessian(i, q, c),
+                  Assert(fe_iv[single_component].jump_in_hessians(i, q) ==
+                           fe_iv.jump_in_shape_hessians(i, q, c),
                          ExcInternalError());
 
-                  Assert(fe_iv[single_component].jump_3rd_derivative(i, q) ==
-                           fe_iv.jump_3rd_derivative(i, q, c),
+                  Assert(fe_iv[single_component].jump_in_third_derivatives(i,
+                                                                           q) ==
+                           fe_iv.jump_in_shape_3rd_derivatives(i, q, c),
                          ExcInternalError());
 
-
-
                 } // q loop
 
           } // component loop
index f788af0daab0f322f76e33d1693efd7b0499359e..28d17a8d525cdf5897e9696e26675c4b5ffd51bf 100644 (file)
@@ -89,17 +89,26 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
                       if (fe_iv[vec_components].value(false, i, q).norm() != 0)
                         deallog << fe_iv[vec_components].value(false, i, q)
                                 << "  ";
-                      if (fe_iv[vec_components].jump(i, q).norm() != 0)
-                        deallog << fe_iv[vec_components].jump(i, q) << "  ";
-                      if (fe_iv[vec_components].average(i, q).norm() != 0)
-                        deallog << fe_iv[vec_components].average(i, q) << "  ";
-                      if (fe_iv[vec_components].jump_gradient(i, q).norm() != 0)
-                        deallog << fe_iv[vec_components].jump_gradient(i, q)
-                                << "  ";
-                      if (fe_iv[vec_components].average_gradient(i, q).norm() !=
+                      if (fe_iv[vec_components].jump_in_values(i, q).norm() !=
                           0)
-                        deallog << fe_iv[vec_components].average_gradient(i, q)
-                                << std::endl;
+                        deallog << fe_iv[vec_components].jump_in_values(i, q)
+                                << "  ";
+                      if (fe_iv[vec_components]
+                            .average_of_values(i, q)
+                            .norm() != 0)
+                        deallog << fe_iv[vec_components].average_of_values(i, q)
+                                << "  ";
+                      if (fe_iv[vec_components]
+                            .jump_in_gradients(i, q)
+                            .norm() != 0)
+                        deallog << fe_iv[vec_components].jump_in_gradients(i, q)
+                                << "  ";
+                      if (fe_iv[vec_components]
+                            .average_of_gradients(i, q)
+                            .norm() != 0)
+                        deallog
+                          << fe_iv[vec_components].average_of_gradients(i, q)
+                          << std::endl;
 
                       for (unsigned int d = 0; d < dim; ++d)
                         {
@@ -109,24 +118,28 @@ test(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
                           Assert(fe_iv[vec_components].value(false, i, q)[d] ==
                                    fe_iv.shape_value(false, i, q, c + d),
                                  ExcInternalError());
-                          Assert(fe_iv[vec_components].jump(i, q)[d] ==
-                                   fe_iv.jump(i, q, c + d),
+                          Assert(fe_iv[vec_components].jump_in_values(i,
+                                                                      q)[d] ==
+                                   fe_iv.jump_in_shape_values(i, q, c + d),
                                  ExcInternalError());
-                          Assert(fe_iv[vec_components].average(i, q)[d] ==
-                                   fe_iv.average(i, q, c + d),
-                                 ExcInternalError());
-                          Assert(fe_iv[vec_components].average_hessian(i,
-                                                                       q)[d] ==
-                                   fe_iv.average_hessian(i, q, c + d),
+                          Assert(
+                            fe_iv[vec_components].average_of_values(i, q)[d] ==
+                              fe_iv.average_of_shape_values(i, q, c + d),
+                            ExcInternalError());
+                          Assert(fe_iv[vec_components].average_of_hessians(
+                                   i, q)[d] ==
+                                   fe_iv.average_of_shape_hessians(i, q, c + d),
                                  ExcInternalError());
-                          Assert(fe_iv[vec_components].jump_hessian(i, q)[d] ==
-                                   fe_iv.jump_hessian(i, q, c + d),
+                          Assert(fe_iv[vec_components].jump_in_hessians(i,
+                                                                        q)[d] ==
+                                   fe_iv.jump_in_shape_hessians(i, q, c + d),
                                  ExcInternalError());
 
-                          Assert(fe_iv[vec_components].jump_3rd_derivative(
-                                   i, q)[d] ==
-                                   fe_iv.jump_3rd_derivative(i, q, c + d),
-                                 ExcInternalError());
+                          Assert(
+                            fe_iv[vec_components].jump_in_third_derivatives(
+                              i, q)[d] ==
+                              fe_iv.jump_in_shape_3rd_derivatives(i, q, c + d),
+                            ExcInternalError());
                         }
                     }
               }
index 7554c4374a008c2dc34b690c5c3067a90568d321..badfc44fe9cad6d541587ea28ac26a516058cbee 100644 (file)
@@ -365,7 +365,7 @@ namespace Step12
           for (unsigned int i = 0; i < n_dofs; ++i)
             for (unsigned int j = 0; j < n_dofs; ++j)
               copy_data_face.cell_matrix(i, j) +=
-                fe_facet.jump(i, qpoint)                        // [\phi_i]
+                fe_facet.jump_in_shape_values(i, qpoint)        // [\phi_i]
                 * fe_facet.shape_value((beta_n > 0), j, qpoint) // phi_j^{UP}
                 * beta_n                                        // (\beta . n)
                 * JxW[qpoint];                                  // dx
index 2eaee8f3951ef38a2547babda8efae5a8c306d38..6ee4288b10134e508d052c543f619739093d6f1a 100644 (file)
@@ -210,9 +210,10 @@ test()
         for (unsigned int j = 0; j < n_dofs; ++j)
           {
             face_matrix(i, j) +=
-              (-fev.jump_gradient(i, q) * n[q] * fev.average(j, q) -
-               fev.average(i, q) * fev.jump_gradient(j, q) * n[q] +
-               gh * fev.jump(i, q) * fev.jump(j, q)) *
+              (-fev.jump_in_shape_gradients(i, q) * n[q] * fev.average(j, q) -
+               fev.average(i, q) * fev.jump_in_shape_gradients(j, q) * n[q] +
+               gh * fev.jump_in_shape_values(i, q) *
+                 fev.jump_in_shape_values(j, q)) *
               JxW[q];
           }
   };
index e9062cfee8932e75c3a5a68e71393acb6fcacaa5..45d195b14c71941d8085962ef2d9390b507082b7 100644 (file)
@@ -534,19 +534,21 @@ test(const unsigned int degree)
           for (unsigned int i = 0; i < n_dofs_face; ++i)
             for (unsigned int j = 0; j < n_dofs_face; ++j)
               copy_data_face.cell_matrix(i, j) +=
-                (-diffusion_coefficient *              // - nu
-                   fe_iv.jump(i, point) *              // [v_h]
-                   (fe_iv.average_gradient(j, point) * // ({grad u_h} .
-                    normals[point])                    //  n)
-
-                 - diffusion_coefficient *               // - nu
-                     (fe_iv.average_gradient(i, point) * // (grad v_h .
-                      normals[point]) *                  //  n)
-                     fe_iv.jump(j, point)                // [u_h]
-
-                 + diffusion_coefficient * penalty * // + nu sigma
-                     fe_iv.jump(i, point) *          // [v_h]
-                     fe_iv.jump(j, point)            // [u_h]
+                (-diffusion_coefficient *                 // - nu
+                   fe_iv.jump_in_shape_values(i, point) * // [v_h]
+                   (fe_iv.average_of_shape_gradients(j,
+                                                     point) * // ({grad u_h} .
+                    normals[point])                           //  n)
+
+                 -
+                 diffusion_coefficient *                         // - nu
+                   (fe_iv.average_of_shape_gradients(i, point) * // (grad v_h .
+                    normals[point]) *                            //  n)
+                   fe_iv.jump_in_shape_values(j, point)          // [u_h]
+
+                 + diffusion_coefficient * penalty *        // + nu sigma
+                     fe_iv.jump_in_shape_values(i, point) * // [v_h]
+                     fe_iv.jump_in_shape_values(j, point)   // [u_h]
 
                  ) *
                 JxW[point]; // dx
index efe02a8aa0cf4f0d5f6085eb2c87aa2584826774..257e14330ecce32dcaa713c5e152753c234ab727 100644 (file)
@@ -439,18 +439,19 @@ DGHeat<dim>::assemble_system()
             for (unsigned int j = 0; j < n_dofs; ++j)
               {
                 copy_data_face.cell_matrix(i, j) +=
-                  -normals[qpoint] * fe_iv.average_gradient(i, qpoint) *
-                  fe_iv.jump(j, qpoint) * JxW[qpoint];
+                  -normals[qpoint] *
+                  fe_iv.average_of_shape_gradients(i, qpoint) *
+                  fe_iv.jump_in_shape_values(j, qpoint) * JxW[qpoint];
 
                 copy_data_face.cell_matrix(i, j) +=
-                  -fe_iv.jump(i, qpoint) // \phi_i
-                  * fe_iv.average_gradient(j, qpoint) *
+                  -fe_iv.jump_in_shape_values(i, qpoint) // \phi_i
+                  * fe_iv.average_of_shape_gradients(j, qpoint) *
                   normals[qpoint] // n*\nabla \phi_j
                   * JxW[qpoint];  // dx
 
                 copy_data_face.cell_matrix(i, j) +=
-                  beta * 1. / h * fe_iv.jump(i, qpoint) *
-                  fe_iv.jump(j, qpoint) * JxW[qpoint];
+                  beta * 1. / h * fe_iv.jump_in_shape_values(i, qpoint) *
+                  fe_iv.jump_in_shape_values(j, qpoint) * JxW[qpoint];
               }
           }
       }
index fa7b2d9fbc9550cbb04552db1f040ed6308859f7..9e3ab2fd4b670dd8bc1a4321e9ad0602828c3bb1 100644 (file)
@@ -348,7 +348,7 @@ namespace Step12
           for (unsigned int i = 0; i < n_dofs; ++i)
             for (unsigned int j = 0; j < n_dofs; ++j)
               copy_data_face.cell_matrix(i, j) +=
-                fe_iv.jump(i, qpoint) // [\phi_i]
+                fe_iv.jump_in_shape_values(i, qpoint) // [\phi_i]
                 *
                 fe_iv.shape_value((beta_dot_n > 0), j, qpoint) // phi_j^{upwind}
                 * beta_dot_n                                   // (\beta . n)
index 8b5660d8361b45b43401f15c7fff8de2a2f809c5..3b4dc502385626a59f28f786b499f81b34071115 100644 (file)
@@ -363,7 +363,7 @@ namespace Step12
           for (unsigned int i = 0; i < n_dofs; ++i)
             for (unsigned int j = 0; j < n_dofs; ++j)
               copy_data_face.cell_matrix(i, j) +=
-                fe_iv.jump(i, qpoint) // [\phi_i]
+                fe_iv.jump_in_shape_values(i, qpoint) // [\phi_i]
                 *
                 fe_iv.shape_value((beta_dot_n > 0), j, qpoint) // phi_j^{upwind}
                 * beta_dot_n                                   // (\beta . n)
index ac2f573c8087a357ed20dce50573fe7dd7111620..13d9bde58952e19d7812e098a4f7199e543aed20 100644 (file)
@@ -479,19 +479,21 @@ namespace Step74
           for (unsigned int i = 0; i < n_dofs_face; ++i)
             for (unsigned int j = 0; j < n_dofs_face; ++j)
               copy_data_face.cell_matrix(i, j) +=
-                (-diffusion_coefficient *              // - nu
-                   fe_iv.jump(i, point) *              // [v_h]
-                   (fe_iv.average_gradient(j, point) * // ({grad u_h} .
-                    normals[point])                    //  n)
-
-                 - diffusion_coefficient *               // - nu
-                     (fe_iv.average_gradient(i, point) * // (grad v_h .
-                      normals[point]) *                  //  n)
-                     fe_iv.jump(j, point)                // [u_h]
-
-                 + diffusion_coefficient * penalty * // + nu sigma
-                     fe_iv.jump(i, point) *          // [v_h]
-                     fe_iv.jump(j, point)            // [u_h]
+                (-diffusion_coefficient *                 // - nu
+                   fe_iv.jump_in_shape_values(i, point) * // [v_h]
+                   (fe_iv.average_of_shape_gradients(j,
+                                                     point) * // ({grad u_h} .
+                    normals[point])                           //  n)
+
+                 -
+                 diffusion_coefficient *                         // - nu
+                   (fe_iv.average_of_shape_gradients(i, point) * // (grad v_h .
+                    normals[point]) *                            //  n)
+                   fe_iv.jump_in_shape_values(j, point)          // [u_h]
+
+                 + diffusion_coefficient * penalty *        // + nu sigma
+                     fe_iv.jump_in_shape_values(i, point) * // [v_h]
+                     fe_iv.jump_in_shape_values(j, point)   // [u_h]
 
                  ) *
                 JxW[point]; // dx

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.