local_solution.hpp

#ifndef PHILIP_LOCAL_SOLUTION_HPP
#define PHILIP_LOCAL_SOLUTION_HPP

#include <deal.II/fe/fe_system.h>
#include <deal.II/base/types.h>


namespace PHiLiP {


template <typename real, int dim, int n_components>
class LocalSolution {
   public:
    std::vector<real> coefficients;
    const dealii::FESystem<dim, dim> &finite_element;
    const unsigned int n_dofs;

    explicit LocalSolution(const dealii::FESystem<dim, dim> &finite_element);

    std::vector<std::array<real, n_components>> evaluate_values(const std::vector<dealii::Point<dim>> &unit_points) const;

    std::vector<std::array<dealii::Tensor<1, dim, real>, n_components>> evaluate_reference_gradients(
        const std::vector<dealii::Point<dim>> &unit_points) const;
};

}  // namespace PHiLiP
#endif  // PHILIP_LOCAL_SOLUTION_HPP