Correcting warning message

# Conflicts:
#	fem/bilinearform.cpp
#	fem/bilinearform.hpp

fem/bilinearform.cpp

@@ -273,6 +273,32 @@ void BilinearForm::AddBdrFaceIntegrator(BilinearFormIntegrator *bfi,
    boundary_face_integs_marker.Append(&bdr_marker);
 }
 
+void BilinearForm::ElementWiseInnerProduct(const GridFunction &x,
+                                           const GridFunction &y,
+                                           Vector &v)
+{
+   MFEM_ASSERT(v.Size() == fes->GetNE(),
+               "Incorrect size for result vector");
+
+   v = 0.0;
+
+   Vector loc_x, loc_y;
+   DenseMatrix elmat;
+
+   Array<int> vdofs;
+
+   for (int i = 0; i < fes->GetNE(); i++)
+   {
+      fes->GetElementVDofs(i, vdofs);
+      x.GetSubVector(vdofs, loc_x);
+      y.GetSubVector(vdofs, loc_y);
+
+      ComputeElementMatrix(i, elmat);
+
+      v[i] = elmat.InnerProduct(loc_x, loc_y);
+   }
+}
+
 void BilinearForm::ComputeElementMatrix(int i, DenseMatrix &elmat) const
 {
    if (element_matrices)

fem/bilinearform.hpp

@@ -334,6 +334,14 @@ public:
    real_t InnerProduct(const Vector &x, const Vector &y) const
    { return mat->InnerProduct (x, y); }
 
+   /// @brief Compute \f$ v_e = y_e^T M_e x_e\f$
+   ///
+   /// @note Recomputes element matrices with each invocation unless
+   ///       the AssemblyLevel is ELEMENT. This can be computationally
+   ///       expensive.
+   void ElementWiseInnerProduct(const GridFunction &x, const GridFunction &y,
+                                Vector &v);
+
    /** @brief Returns a pointer to (approximation) of the matrix inverse:
        $ M^{-1} $ (currently returns NULL) */
    MatrixInverse *Inverse() const override;