localenergy.hh

#ifndef DUNE_GFE_LOCALENERGY_HH
#define DUNE_GFE_LOCALENERGY_HH

#include <vector>

namespace Dune {

namespace GFE {

/** \brief Base class for energies defined by integrating over one grid element */
template<class Basis, class... TargetSpaces>
class LocalEnergy
{
public:
  using RT = typename std::common_type<typename TargetSpaces::ctype...>::type;

  /** \brief Compute the energy
   *
   * \param localView Local view specifying the current element and the FE space there
   * \param coefficients The coefficients of a FE function on the current element
   */
  virtual RT
  energy (const typename Basis::LocalView& localView,
          const std::vector<TargetSpaces>&... localSolution) const = 0;

  /** Empty virtual default destructor
   *
   * To allow proper destruction of derived classes through a base class pointer
   */
  virtual ~LocalEnergy() = default;

};

}  // namespace GFE

}  // namespace Dune

#endif  // DUNE_GFE_LOCALENERGY_HH