l2projection.hh

#ifndef __DUNE_ACFEM_L2PROJECTION_HH__
#define __DUNE_ACFEM_L2PROJECTION_HH__
 
#include "functionals/functionals.hh"
#include "../models/basicmodels.hh"
#include "../algorithms/ellipticfemscheme.hh"
 
namespace Dune {
 
  namespace ACFem {
 
    template<class Functional, class DiscreteFunction,
             std::enable_if_t<(IsLinearFunctional<Functional>::value
                               && IsDiscreteFunction<DiscreteFunction>::value
      ), int> = 0>
    void l2Projection(const Functional& phi, DiscreteFunction& result)
    {
      // skip the initialization, just do a quick loop over "result"
      // in order to determine that all values are finite
      for (auto&& dof : dofs(result)) {
        if (!std::isfinite(dof)) {
          dof = 0.;
 
        }
      }
 
      ellipticFemScheme(result, massModel(result), phi, "acfem.l2projection").solve();
    }
 
    template<class GridFunction, class DiscreteFunction,
             std::enable_if_t<(IsWrappableByConstLocalFunction<GridFunction>::value
                               && IsDiscreteFunction<DiscreteFunction>::value
      ), int> = 0>
    void l2Projection(GridFunction&& fct, DiscreteFunction& result)
    {
      // skip the initialization, just do a quick loop over "result"
      // in order to determine that all values are finite
      for (auto&& dof : dofs(result)) {
        if (!std::isfinite(dof)) {
          dof = 0.;
 
        }
      }
 
      ellipticFemScheme(result, massModel(result) - fct, "acfem.l2projection").solve();
    }
 
 
 
  } // ACFem::
 
} // Dune::
 
#endif // __DUNE_ACFEM_L2PROJECTION_HH__