Remove deprecated functions in fe_interface_values.h

diff --git a/doc/news/changes/incompatibilities/20230629DanielArndt-2 b/doc/news/changes/incompatibilities/20230629DanielArndt-2
new file mode 100644 (file)
index 0000000..a5f3eda
--- /dev/null
+++ b/doc/news/changes/incompatibilities/20230629DanielArndt-2
@@ -0,0 +1,3 @@
+Removed: Deprecated functions in the FEInterfaceValues nampespace have been removed.
+<br>
+(Daniel Arndt, 2023/06/29)

diff --git a/include/deal.II/fe/fe_interface_values.h b/include/deal.II/fe/fe_interface_values.h
index e032f72d35d7b50850137d7e95e04c3236cea03f..dc6684e2298b91f9143cb69b918591d595b079c3 100644 (file)
--- a/include/deal.II/fe/fe_interface_values.h
+++ b/include/deal.II/fe/fe_interface_values.h
@@ -198,16 +198,6 @@ namespace FEInterfaceViews
     jump_in_values(const unsigned int interface_dof_index,
                    const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the jump_in_values() function instead.
-     */
-    DEAL_II_DEPRECATED
-    value_type
-    jump(const unsigned int interface_dof_index,
-         const unsigned int q_point) const;
-
     /**
      * Return the jump of the gradient $\jump{nabla u}$ on the interface for
      * the shape
@@ -222,16 +212,6 @@ namespace FEInterfaceViews
     jump_in_gradients(const unsigned int interface_dof_index,
                       const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the jump_in_gradients() function instead.
-     */
-    DEAL_II_DEPRECATED
-    gradient_type
-    jump_gradient(const unsigned int interface_dof_index,
-                  const unsigned int q_point) const;
-
     /**
      * Return the jump in the gradient $\jump{\nabla u}=\nabla u_{\text{cell0}}
      * - \nabla u_{\text{cell1}}$ on the interface for the shape function @p
@@ -246,16 +226,6 @@ namespace FEInterfaceViews
     jump_in_hessians(const unsigned int interface_dof_index,
                      const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the jump_in_hessians() function instead.
-     */
-    DEAL_II_DEPRECATED
-    hessian_type
-    jump_hessian(const unsigned int interface_dof_index,
-                 const unsigned int q_point) const;
-
     /**
      * Return the jump in the third derivative $\jump{\nabla^3 u} = \nabla^3
      * u_{\text{cell0}} - \nabla^3 u_{\text{cell1}}$ on the interface for the
@@ -270,16 +240,6 @@ namespace FEInterfaceViews
     jump_in_third_derivatives(const unsigned int interface_dof_index,
                               const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the jump_in_third_derivatives() function instead.
-     */
-    DEAL_II_DEPRECATED
-    third_derivative_type
-    jump_3rd_derivative(const unsigned int interface_dof_index,
-                        const unsigned int q_point) const;
-
     /** @} */
 
     /**
@@ -301,26 +261,6 @@ namespace FEInterfaceViews
     average_of_values(const unsigned int interface_dof_index,
                       const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the average_of_values() function instead.
-     */
-    DEAL_II_DEPRECATED
-    value_type
-    average_value(const unsigned int interface_dof_index,
-                  const unsigned int q_point) const;
-
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the average_of_values() function instead.
-     */
-    DEAL_II_DEPRECATED
-    value_type
-    average(const unsigned int interface_dof_index,
-            const unsigned int q_point) const;
-
     /**
      * Return the average of the gradient $\average{\nabla u}$ on the interface
      * for the shape
@@ -335,16 +275,6 @@ namespace FEInterfaceViews
     average_of_gradients(const unsigned int interface_dof_index,
                          const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the average_of_gradients() function instead.
-     */
-    DEAL_II_DEPRECATED
-    gradient_type
-    average_gradient(const unsigned int interface_dof_index,
-                     const unsigned int q_point) const;
-
     /**
      * Return the average of the Hessian $\average{\nabla^2 u} =
      * \frac{1}{2}\nabla^2 u_{\text{cell0}} + \frac{1}{2} \nabla^2
@@ -360,16 +290,6 @@ namespace FEInterfaceViews
     average_of_hessians(const unsigned int interface_dof_index,
                         const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the average_of_hessians() function instead.
-     */
-    DEAL_II_DEPRECATED
-    hessian_type
-    average_hessian(const unsigned int interface_dof_index,
-                    const unsigned int q_point) const;
-
     /** @} */
 
     /**
@@ -830,16 +750,6 @@ namespace FEInterfaceViews
     jump_in_values(const unsigned int interface_dof_index,
                    const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the jump_in_values() function instead.
-     */
-    DEAL_II_DEPRECATED
-    value_type
-    jump(const unsigned int interface_dof_index,
-         const unsigned int q_point) const;
-
     /**
      * Return the jump of the gradient (a tensor of rank 2) $\jump{\nabla
      * \mathbf{u}}$ on the interface for the shape
@@ -899,16 +809,6 @@ namespace FEInterfaceViews
     jump_in_third_derivatives(const unsigned int interface_dof_index,
                               const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the jump_in_third_derivatives() function instead.
-     */
-    DEAL_II_DEPRECATED
-    third_derivative_type
-    jump_3rd_derivative(const unsigned int interface_dof_index,
-                        const unsigned int q_point) const;
-
     /** @} */
 
     /**
@@ -929,16 +829,6 @@ namespace FEInterfaceViews
     average_of_values(const unsigned int interface_dof_index,
                       const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the average_of_values() function instead.
-     */
-    DEAL_II_DEPRECATED
-    value_type
-    average(const unsigned int interface_dof_index,
-            const unsigned int q_point) const;
-
     /**
      * Return the average of the gradient (a tensor of rank 2) $\average{\nabla
      * \mathbf{u}}$ on the interface for the shape
@@ -952,16 +842,6 @@ namespace FEInterfaceViews
     average_of_gradients(const unsigned int interface_dof_index,
                          const unsigned int q_point) const;
 
-    /**
-     * The same as above.
-     *
-     * @deprecated Use the average_of_gradients() function instead.
-     */
-    DEAL_II_DEPRECATED
-    gradient_type
-    average_gradient(const unsigned int interface_dof_index,
-                     const unsigned int q_point) const;
-
     /**
      * Return the average of the Hessian $\average{\nabla^2 u} =
      * \frac{1}{2}\nabla^2 u_{\text{cell0}} + \frac{1}{2} \nabla^2
@@ -1903,17 +1783,6 @@ public:
                        const unsigned int q_point,
                        const unsigned int component = 0) const;
 
-  /**
-   * The same as above.
-   *
-   * @deprecated Use the jump_in_shape_values() function instead.
-   */
-  DEAL_II_DEPRECATED
-  double
-  jump(const unsigned int interface_dof_index,
-       const unsigned int q_point,
-       const unsigned int component = 0) const;
-
   /**
    * Return the jump in the gradient $\jump{\nabla u}=\nabla u_{\text{cell0}} -
    * \nabla u_{\text{cell1}}$ on the interface for the shape function @p
@@ -1932,17 +1801,6 @@ public:
                           const unsigned int q_point,
                           const unsigned int component = 0) const;
 
-  /**
-   * The same as above.
-   *
-   * @deprecated Use the jump_in_shape_gradients() function instead.
-   */
-  DEAL_II_DEPRECATED
-  Tensor<1, spacedim>
-  jump_gradient(const unsigned int interface_dof_index,
-                const unsigned int q_point,
-                const unsigned int component = 0) const;
-
   /**
    * Return the jump in the Hessian $\jump{\nabla^2 u} = \nabla^2
    * u_{\text{cell0}} - \nabla^2 u_{\text{cell1}}$ on the interface for the
@@ -1962,17 +1820,6 @@ public:
                          const unsigned int q_point,
                          const unsigned int component = 0) const;
 
-  /**
-   * The same as above.
-   *
-   * @deprecated Use the jump_in_shape_hessians() function instead.
-   */
-  DEAL_II_DEPRECATED
-  Tensor<2, spacedim>
-  jump_hessian(const unsigned int interface_dof_index,
-               const unsigned int q_point,
-               const unsigned int component = 0) const;
-
   /**
    * Return the jump in the third derivative $\jump{\nabla^3 u} = \nabla^3
    * u_{\text{cell0}} - \nabla^3 u_{\text{cell1}}$ on the interface for the
@@ -1991,17 +1838,6 @@ public:
                                 const unsigned int q_point,
                                 const unsigned int component = 0) const;
 
-  /**
-   * The same as above.
-   *
-   * @deprecated Use the jump_in_shape_3rd_derivatives() function instead.
-   */
-  DEAL_II_DEPRECATED
-  Tensor<3, spacedim>
-  jump_3rd_derivative(const unsigned int interface_dof_index,
-                      const unsigned int q_point,
-                      const unsigned int component = 0) const;
-
   /**
    * @}
    */
@@ -2029,17 +1865,6 @@ public:
                           const unsigned int q_point,
                           const unsigned int component = 0) const;
 
-  /**
-   * The same as above.
-   *
-   * @deprecated Use the average_of_shape_values() function instead.
-   */
-  DEAL_II_DEPRECATED
-  double
-  average(const unsigned int interface_dof_index,
-          const unsigned int q_point,
-          const unsigned int component = 0) const;
-
   /**
    * Return the average of the gradient $\average{\nabla u} = \frac{1}{2}\nabla
    * u_{\text{cell0}} + \frac{1}{2} \nabla u_{\text{cell1}}$ on the interface
@@ -2058,17 +1883,6 @@ public:
                              const unsigned int q_point,
                              const unsigned int component = 0) const;
 
-  /**
-   * The same as above.
-   *
-   * @deprecated Use the average_of_shape_gradients() function instead.
-   */
-  DEAL_II_DEPRECATED
-  Tensor<1, spacedim>
-  average_gradient(const unsigned int interface_dof_index,
-                   const unsigned int q_point,
-                   const unsigned int component = 0) const;
-
   /**
    * Return the average of the Hessian $\average{\nabla^2 u} =
    * \frac{1}{2}\nabla^2 u_{\text{cell0}} + \frac{1}{2} \nabla^2
@@ -2088,17 +1902,6 @@ public:
                             const unsigned int q_point,
                             const unsigned int component = 0) const;
 
-  /**
-   * The same as above.
-   *
-   * @deprecated Use the average_of_shape_hessians() function instead.
-   */
-  DEAL_II_DEPRECATED
-  Tensor<2, spacedim>
-  average_hessian(const unsigned int interface_dof_index,
-                  const unsigned int q_point,
-                  const unsigned int component = 0) const;
-
   /**
    * @}
    */
@@ -3043,17 +2846,6 @@ FEInterfaceValues<dim, spacedim>::jump_in_shape_values(
 
 
 
-template <int dim, int spacedim>
-double
-FEInterfaceValues<dim, spacedim>::jump(const unsigned int interface_dof_index,
-                                       const unsigned int q_point,
-                                       const unsigned int component) const
-{
-  return jump_in_shape_values(interface_dof_index, q_point, component);
-}
-
-
-
 template <int dim, int spacedim>
 double
 FEInterfaceValues<dim, spacedim>::average_of_shape_values(
@@ -3084,18 +2876,6 @@ FEInterfaceValues<dim, spacedim>::average_of_shape_values(
 
 
 
-template <int dim, int spacedim>
-double
-FEInterfaceValues<dim, spacedim>::average(
-  const unsigned int interface_dof_index,
-  const unsigned int q_point,
-  const unsigned int component) const
-{
-  return average_of_shape_values(interface_dof_index, q_point, component);
-}
-
-
-
 template <int dim, int spacedim>
 Tensor<1, spacedim>
 FEInterfaceValues<dim, spacedim>::average_of_shape_gradients(
@@ -3126,18 +2906,6 @@ FEInterfaceValues<dim, spacedim>::average_of_shape_gradients(
 
 
 
-template <int dim, int spacedim>
-Tensor<1, spacedim>
-FEInterfaceValues<dim, spacedim>::average_gradient(
-  const unsigned int interface_dof_index,
-  const unsigned int q_point,
-  const unsigned int component) const
-{
-  return average_of_shape_gradients(interface_dof_index, q_point, component);
-}
-
-
-
 template <int dim, int spacedim>
 Tensor<2, spacedim>
 FEInterfaceValues<dim, spacedim>::average_of_shape_hessians(
@@ -3168,18 +2936,6 @@ FEInterfaceValues<dim, spacedim>::average_of_shape_hessians(
 
 
 
-template <int dim, int spacedim>
-Tensor<2, spacedim>
-FEInterfaceValues<dim, spacedim>::average_hessian(
-  const unsigned int interface_dof_index,
-  const unsigned int q_point,
-  const unsigned int component) const
-{
-  return average_of_shape_hessians(interface_dof_index, q_point, component);
-}
-
-
-
 template <int dim, int spacedim>
 Tensor<1, spacedim>
 FEInterfaceValues<dim, spacedim>::jump_in_shape_gradients(
@@ -3210,18 +2966,6 @@ FEInterfaceValues<dim, spacedim>::jump_in_shape_gradients(
 
 
 
-template <int dim, int spacedim>
-Tensor<1, spacedim>
-FEInterfaceValues<dim, spacedim>::jump_gradient(
-  const unsigned int interface_dof_index,
-  const unsigned int q_point,
-  const unsigned int component) const
-{
-  return jump_in_shape_gradients(interface_dof_index, q_point, component);
-}
-
-
-
 template <int dim, int spacedim>
 Tensor<2, spacedim>
 FEInterfaceValues<dim, spacedim>::jump_in_shape_hessians(
@@ -3252,18 +2996,6 @@ FEInterfaceValues<dim, spacedim>::jump_in_shape_hessians(
 
 
 
-template <int dim, int spacedim>
-Tensor<2, spacedim>
-FEInterfaceValues<dim, spacedim>::jump_hessian(
-  const unsigned int interface_dof_index,
-  const unsigned int q_point,
-  const unsigned int component) const
-{
-  return jump_in_shape_hessians(interface_dof_index, q_point, component);
-}
-
-
-
 template <int dim, int spacedim>
 Tensor<3, spacedim>
 FEInterfaceValues<dim, spacedim>::jump_in_shape_3rd_derivatives(
@@ -3294,18 +3026,6 @@ FEInterfaceValues<dim, spacedim>::jump_in_shape_3rd_derivatives(
 
 
 
-template <int dim, int spacedim>
-Tensor<3, spacedim>
-FEInterfaceValues<dim, spacedim>::jump_3rd_derivative(
-  const unsigned int interface_dof_index,
-  const unsigned int q_point,
-  const unsigned int component) const
-{
-  return jump_in_shape_3rd_derivatives(interface_dof_index, q_point, component);
-}
-
-
-
 template <int dim, int spacedim>
 template <class InputVector>
 void
@@ -3540,16 +3260,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Scalar<dim, spacedim>::value_type
-  Scalar<dim, spacedim>::jump(const unsigned int interface_dof_index,
-                              const unsigned int q_point) const
-  {
-    return jump_in_values(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Scalar<dim, spacedim>::value_type
   Scalar<dim, spacedim>::average_of_values(
@@ -3580,16 +3290,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Scalar<dim, spacedim>::value_type
-  Scalar<dim, spacedim>::average(const unsigned int interface_dof_index,
-                                 const unsigned int q_point) const
-  {
-    return average_of_values(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Scalar<dim, spacedim>::gradient_type
   Scalar<dim, spacedim>::average_of_gradients(
@@ -3618,16 +3318,6 @@ namespace FEInterfaceViews
   }
 
 
-  template <int dim, int spacedim>
-  typename Scalar<dim, spacedim>::gradient_type
-  Scalar<dim, spacedim>::average_gradient(
-    const unsigned int interface_dof_index,
-    const unsigned int q_point) const
-  {
-    return average_of_gradients(interface_dof_index, q_point);
-  }
-
-
 
   template <int dim, int spacedim>
   typename Scalar<dim, spacedim>::gradient_type
@@ -3658,16 +3348,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Scalar<dim, spacedim>::gradient_type
-  Scalar<dim, spacedim>::jump_gradient(const unsigned int interface_dof_index,
-                                       const unsigned int q_point) const
-  {
-    return jump_in_gradients(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Scalar<dim, spacedim>::hessian_type
   Scalar<dim, spacedim>::average_of_hessians(
@@ -3697,16 +3377,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Scalar<dim, spacedim>::hessian_type
-  Scalar<dim, spacedim>::average_hessian(const unsigned int interface_dof_index,
-                                         const unsigned int q_point) const
-  {
-    return average_of_hessians(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Scalar<dim, spacedim>::third_derivative_type
   Scalar<dim, spacedim>::jump_in_third_derivatives(
@@ -3735,17 +3405,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Scalar<dim, spacedim>::third_derivative_type
-  Scalar<dim, spacedim>::jump_3rd_derivative(
-    const unsigned int interface_dof_index,
-    const unsigned int q_point) const
-  {
-    return jump_in_third_derivatives(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Scalar<dim, spacedim>::hessian_type
   Scalar<dim, spacedim>::jump_in_hessians(
@@ -3775,16 +3434,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Scalar<dim, spacedim>::hessian_type
-  Scalar<dim, spacedim>::jump_hessian(const unsigned int interface_dof_index,
-                                      const unsigned int q_point) const
-  {
-    return jump_in_hessians(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   template <class InputVector>
   void
@@ -4176,16 +3825,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Vector<dim, spacedim>::value_type
-  Vector<dim, spacedim>::jump(const unsigned int interface_dof_index,
-                              const unsigned int q_point) const
-  {
-    return jump_in_values(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Vector<dim, spacedim>::value_type
   Vector<dim, spacedim>::average_of_values(
@@ -4216,16 +3855,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Vector<dim, spacedim>::value_type
-  Vector<dim, spacedim>::average(const unsigned int interface_dof_index,
-                                 const unsigned int q_point) const
-  {
-    return average_of_values(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Vector<dim, spacedim>::gradient_type
   Vector<dim, spacedim>::average_of_gradients(
@@ -4255,17 +3884,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Vector<dim, spacedim>::gradient_type
-  Vector<dim, spacedim>::average_gradient(
-    const unsigned int interface_dof_index,
-    const unsigned int q_point) const
-  {
-    return average_of_gradients(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   typename Vector<dim, spacedim>::gradient_type
   Vector<dim, spacedim>::jump_in_gradients(
@@ -4411,17 +4029,6 @@ namespace FEInterfaceViews
 
 
 
-  template <int dim, int spacedim>
-  typename Vector<dim, spacedim>::third_derivative_type
-  Vector<dim, spacedim>::jump_3rd_derivative(
-    const unsigned int interface_dof_index,
-    const unsigned int q_point) const
-  {
-    return jump_in_third_derivatives(interface_dof_index, q_point);
-  }
-
-
-
   template <int dim, int spacedim>
   template <class InputVector>
   void

diff --git a/tests/feinterface/stokes.cc b/tests/feinterface/stokes.cc
index 883c28e279563c9787386563093a2b8f3d49b827..8ba022a6735d04753b3e6c13e57795f20ac0d8cc 100644 (file)
--- a/tests/feinterface/stokes.cc
+++ b/tests/feinterface/stokes.cc
@@ -665,28 +665,28 @@ namespace StokesTests
                 (
                   // - nu {\nabla u}n . [v] (consistency)
                   -nu *
-                    (fe_fv[velocities].average_gradient(j, point) *
+                    (fe_fv[velocities].average_of_gradients(j, point) *
                      normals[point]) *
-                    fe_fv[velocities].jump(i, point)
+                    fe_fv[velocities].jump_in_values(i, point)
 
                   // - nu [u] . {\nabla v}n  (symmetry) // NIPG: use +
-                  - nu * fe_fv[velocities].jump(j, point) *
-                      (fe_fv[velocities].average_gradient(i, point) *
+                  - nu * fe_fv[velocities].jump_in_values(j, point) *
+                      (fe_fv[velocities].average_of_gradients(i, point) *
                        normals[point])
 
                   // nu sigma [u].[v] (penalty)
                   + nu * penalty *
-                      scalar_product(fe_fv[velocities].jump(j, point),
-                                     fe_fv[velocities].jump(i, point))
+                      scalar_product(fe_fv[velocities].jump_in_values(j, point),
+                                     fe_fv[velocities].jump_in_values(i, point))
 
                   // {p} ([v].n)
-                  + fe_fv[pressure].average(j, point) *
-                      scalar_product(fe_fv[velocities].jump(i, point),
+                  + fe_fv[pressure].average_of_values(j, point) *
+                      scalar_product(fe_fv[velocities].jump_in_values(i, point),
                                      normals[point])
 
                   // {q} ([u].n)
-                  + fe_fv[pressure].average(i, point) *
-                      scalar_product(fe_fv[velocities].jump(j, point),
+                  + fe_fv[pressure].average_of_values(i, point) *
+                      scalar_product(fe_fv[velocities].jump_in_values(j, point),
                                      normals[point])) *
                 JxW[point];
         }

diff --git a/tests/meshworker/scratch_data_08.cc b/tests/meshworker/scratch_data_08.cc
index 6ee4288b10134e508d052c543f619739093d6f1a..09b17d7577957487b23e034f968b0b0755027414 100644 (file)
--- a/tests/meshworker/scratch_data_08.cc
+++ b/tests/meshworker/scratch_data_08.cc
@@ -209,12 +209,13 @@ test()
       for (unsigned int i = 0; i < n_dofs; ++i)
         for (unsigned int j = 0; j < n_dofs; ++j)
           {
-            face_matrix(i, j) +=
-              (-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];
+            face_matrix(i, j) += (-fev.jump_in_shape_gradients(i, q) * n[q] *
+                                    fev.average_of_shape_values(j, q) -
+                                  fev.average_of_shape_values(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];
           }
   };