Fix some tests by explicitly using make_array_view().

diff --git a/tests/integrators/advection_01.cc b/tests/integrators/advection_01.cc
index 55f68d1b4c02fb2f3a856c594e1f502c3f2b2982..add2bb55c165682b99831963454e7e93d0cf7acf 100644 (file)
--- a/tests/integrators/advection_01.cc
+++ b/tests/integrators/advection_01.cc
@@ -63,7 +63,7 @@ test_cell(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_values(u, indices, uval, true);
+        fev.get_function_values(u, indices, make_array_view(uval), true);
         cell_residual(w, fev, uval, vel);
         M.vmult(v, u);
         w.add(-1., v);
@@ -118,7 +118,7 @@ test_boundary(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_values(u, indices, uval, true);
+        fev.get_function_values(u, indices, make_array_view(uval), true);
         upwind_value_residual(w, fev, uval, null_val, vel);
         M.vmult(v, u);
         w.add(-1., v);
@@ -191,7 +191,7 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
         u1(i1) = 1.;
         w1     = 0.;
         w2     = 0.;
-        fev1.get_function_values(u1, indices1, u1val, true);
+        fev1.get_function_values(u1, indices1, make_array_view(u1val), true);
         upwind_face_residual(w1, w2, fev1, fev2, u1val, nullval, vel);
         M11.vmult(v1, u1);
         w1.add(-1., v1);
@@ -213,7 +213,7 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
           }
         w1 = 0.;
         w2 = 0.;
-        fev2.get_function_values(u1, indices2, u1val, true);
+        fev2.get_function_values(u1, indices2, make_array_view(u1val), true);
         upwind_face_residual(w1, w2, fev1, fev2, nullval, u1val, vel);
         M12.vmult(v1, u1);
         w1.add(-1., v1);

diff --git a/tests/integrators/divergence_01.cc b/tests/integrators/divergence_01.cc
index 2d49e7b1accac321a7f53abdf9e84d5080c69a49..e497de263df919b3af6bb961ccce6dbb844085cf 100644 (file)
--- a/tests/integrators/divergence_01.cc
+++ b/tests/integrators/divergence_01.cc
@@ -66,7 +66,7 @@ test_cell(const FEValuesBase<dim> &fev, const FEValuesBase<dim> &fes)
       u    = 0.;
       u(i) = 1.;
       w    = 0.;
-      fev.get_function_gradients(u, indices, ugrad, true);
+      fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
       cell_residual<dim>(w, fes, ugrad);
       Md.vmult(v, u);
       w.add(-1., v);
@@ -123,7 +123,7 @@ test_boundary(const FEValuesBase<dim> &fev, const FEValuesBase<dim> &fes)
       u    = 0.;
       u(i) = 1.;
       w    = 0.;
-      fev.get_function_values(u, indices, uval, true);
+      fev.get_function_values(u, indices, make_array_view(uval), true);
       u_dot_n_residual(w, fev, fes, uval);
       M.vmult(v, u);
       w.add(-1., v);

diff --git a/tests/integrators/elasticity_01.cc b/tests/integrators/elasticity_01.cc
index 65a993619e176c43ee88f43a050d636068c1a11f..b29c6b7e9b567bf7ac481afc0b1e44a567af64e2 100644 (file)
--- a/tests/integrators/elasticity_01.cc
+++ b/tests/integrators/elasticity_01.cc
@@ -58,7 +58,7 @@ test_cell(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_gradients(u, indices, ugrad, true);
+        fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
         cell_residual<dim>(w, fev, ugrad);
         M.vmult(v, u);
         w.add(-1., v);
@@ -100,8 +100,8 @@ test_boundary(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_values(u, indices, uval, true);
-        fev.get_function_gradients(u, indices, ugrad, true);
+        fev.get_function_values(u, indices, make_array_view(uval), true);
+        fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
         nitsche_residual<dim>(w, fev, uval, ugrad, null_val, 17);
         M.vmult(v, u);
         w.add(-1., v);
@@ -168,8 +168,11 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
         u1(i1) = 1.;
         w1     = 0.;
         w2     = 0.;
-        fev1.get_function_values(u1, indices1, u1val, true);
-        fev1.get_function_gradients(u1, indices1, u1grad, true);
+        fev1.get_function_values(u1, indices1, make_array_view(u1val), true);
+        fev1.get_function_gradients(u1,
+                                    indices1,
+                                    make_array_view(u1grad),
+                                    true);
         ip_residual<dim>(
           w1, w2, fev1, fev2, u1val, u1grad, nullval, nullgrad, 17);
         M11.vmult(v1, u1);
@@ -180,8 +183,11 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
 
         w1 = 0.;
         w2 = 0.;
-        fev2.get_function_values(u1, indices2, u1val, true);
-        fev2.get_function_gradients(u1, indices2, u1grad, true);
+        fev2.get_function_values(u1, indices2, make_array_view(u1val), true);
+        fev2.get_function_gradients(u1,
+                                    indices2,
+                                    make_array_view(u1grad),
+                                    true);
         ip_residual<dim>(
           w1, w2, fev1, fev2, nullval, nullgrad, u1val, u1grad, 17);
         M12.vmult(v1, u1);

diff --git a/tests/integrators/elasticity_02.cc b/tests/integrators/elasticity_02.cc
index 851d363b0e294a268073d80de2ddd6e88ded29a8..647fead845dcee1206d42fd52e22a653476a2843 100644 (file)
--- a/tests/integrators/elasticity_02.cc
+++ b/tests/integrators/elasticity_02.cc
@@ -62,8 +62,8 @@ test_boundary(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_values(u, indices, uval, true);
-        fev.get_function_gradients(u, indices, ugrad, true);
+        fev.get_function_values(u, indices, make_array_view(uval), true);
+        fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
         nitsche_tangential_residual<dim>(w, fev, uval, ugrad, null_val, 17);
         M.vmult(v, u);
         w.add(-1., v);

diff --git a/tests/integrators/grad_div_01.cc b/tests/integrators/grad_div_01.cc
index 378f499a9a483072b165efe8fd08a556188db6fc..f9f04de5c5fb6664bdebf95d83d40ec147d29c2b 100644 (file)
--- a/tests/integrators/grad_div_01.cc
+++ b/tests/integrators/grad_div_01.cc
@@ -58,7 +58,7 @@ test_cell(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_gradients(u, indices, ugrad, true);
+        fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
         cell_residual<dim>(w, fev, ugrad);
         M.vmult(v, u);
         w.add(-1., v);
@@ -100,8 +100,8 @@ test_boundary(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_values(u, indices, uval, true);
-        fev.get_function_gradients(u, indices, ugrad, true);
+        fev.get_function_values(u, indices, make_array_view(uval), true);
+        fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
         nitsche_residual<dim>(w, fev, uval, ugrad, null_val, 17);
         M.vmult(v, u);
         w.add(-1., v);
@@ -168,8 +168,11 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
         u1(i1) = 1.;
         w1     = 0.;
         w2     = 0.;
-        fev1.get_function_values(u1, indices1, u1val, true);
-        fev1.get_function_gradients(u1, indices1, u1grad, true);
+        fev1.get_function_values(u1, indices1, make_array_view(u1val), true);
+        fev1.get_function_gradients(u1,
+                                    indices1,
+                                    make_array_view(u1grad),
+                                    true);
         ip_residual<dim>(
           w1, w2, fev1, fev2, u1val, u1grad, nullval, nullgrad, 17);
         M11.vmult(v1, u1);
@@ -180,8 +183,11 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
 
         w1 = 0.;
         w2 = 0.;
-        fev2.get_function_values(u1, indices2, u1val, true);
-        fev2.get_function_gradients(u1, indices2, u1grad, true);
+        fev2.get_function_values(u1, indices2, make_array_view(u1val), true);
+        fev2.get_function_gradients(u1,
+                                    indices2,
+                                    make_array_view(u1grad),
+                                    true);
         ip_residual<dim>(
           w1, w2, fev1, fev2, nullval, nullgrad, u1val, u1grad, 17);
         M12.vmult(v1, u1);

diff --git a/tests/integrators/laplacian_01.cc b/tests/integrators/laplacian_01.cc
index 82a8b9dc676427bd90a66874970794e0a3a48da6..38c60acfe59f124052b25fbe0125387a6bc19e0a 100644 (file)
--- a/tests/integrators/laplacian_01.cc
+++ b/tests/integrators/laplacian_01.cc
@@ -57,7 +57,7 @@ test_cell(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_gradients(u, indices, ugrad, true);
+        fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
         cell_residual<dim>(w, fev, ugrad);
         M.vmult(v, u);
         w.add(-1., v);
@@ -106,8 +106,8 @@ test_boundary(const FEValuesBase<dim> &fev)
         u    = 0.;
         u(i) = 1.;
         w    = 0.;
-        fev.get_function_values(u, indices, uval, true);
-        fev.get_function_gradients(u, indices, ugrad, true);
+        fev.get_function_values(u, indices, make_array_view(uval), true);
+        fev.get_function_gradients(u, indices, make_array_view(ugrad), true);
         nitsche_residual<dim>(w, fev, uval, ugrad, null_val, 17);
         M.vmult(v, u);
         w.add(-1., v);
@@ -181,8 +181,11 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
         u1(i1) = 1.;
         w1     = 0.;
         w2     = 0.;
-        fev1.get_function_values(u1, indices1, u1val, true);
-        fev1.get_function_gradients(u1, indices1, u1grad, true);
+        fev1.get_function_values(u1, indices1, make_array_view(u1val), true);
+        fev1.get_function_gradients(u1,
+                                    indices1,
+                                    make_array_view(u1grad),
+                                    true);
         ip_residual<dim>(
           w1, w2, fev1, fev2, u1val, u1grad, nullval, nullgrad, 17);
         M11.vmult(v1, u1);
@@ -213,8 +216,11 @@ test_face(const FEValuesBase<dim> &fev1, const FEValuesBase<dim> &fev2)
 
         w1 = 0.;
         w2 = 0.;
-        fev2.get_function_values(u1, indices2, u1val, true);
-        fev2.get_function_gradients(u1, indices2, u1grad, true);
+        fev2.get_function_values(u1, indices2, make_array_view(u1val), true);
+        fev2.get_function_gradients(u1,
+                                    indices2,
+                                    make_array_view(u1grad),
+                                    true);
         ip_residual<dim>(
           w1, w2, fev1, fev2, nullval, nullgrad, u1val, u1grad, 17);
         M12.vmult(v1, u1);