fe_values_shape_hessian.h

#ifndef __FE_VALUES_SHAPE_HESSIAN_H__
#define __FE_VALUES_SHAPE_HESSIAN_H__

#include <deal.II/fe/fe_values.h>
#include <deal.II/fe/fe_system.h>

namespace PHiLiP {

template<int dim>
class FEValuesShapeHessian {

public:
    void reinit(const dealii::FEValues<dim,dim> &fe_values_volume, const unsigned int iquad); 


    dealii::Tensor<2,dim,double> shape_hessian_component(
        const unsigned int idof, 
        const unsigned int iquad, 
        const unsigned int istate, 
        const dealii::FESystem<dim,dim> &fe_ref) const;

private:
    dealii::Tensor<2,dim,double> jacobian_inverse;
    
    dealii::Tensor<3,dim,double> derivative_jacobian_inverse_wrt_phys_q;
    
    dealii::Point<dim> ref_point;
};
} // PHiLiP namespace
#endif