ADTypes.hpp

#ifndef __AD_TYPES__
#define __AD_TYPES__

#include <Sacado.hpp>
#include <CoDiPack/include/codi.hpp>
#include <deal.II/differentiation/ad/sacado_math.h>
#include <deal.II/differentiation/ad/sacado_number_types.h>
#include <deal.II/differentiation/ad/sacado_product_types.h>

namespace PHiLiP {
using FadType = Sacado::Fad::DFad<double>; 
using FadFadType = Sacado::Fad::DFad<FadType>; 

static constexpr int dimForwardAD = 1; 
static constexpr int dimReverseAD = 1; 

using codi_FadType = codi::RealForwardGen<double, codi::Direction<double,dimForwardAD>>; 
//using codi_FadType = codi::RealForwardGen<double, codi::DirectionVar<double>>;

using codi_JacobianComputationType = codi::RealReverseIndexVec<dimReverseAD>; 
using codi_HessianComputationType  = codi::RealReversePrimalIndexGen< codi::RealForwardVec<dimForwardAD>,
                                                  codi::Direction< codi::RealForwardVec<dimForwardAD>, dimReverseAD>
                                                >; 

//using RadFadType = Sacado::Rad::ADvar<FadType>; ///< Sacado AD type that allows 2nd derivatives.
//using RadFadType = codi_JacobianComputationType; ///< Reverse only mode that only allows Jacobian computation.
using RadType = codi_JacobianComputationType; 
using RadFadType = codi_HessianComputationType ; 
} // PHiLiP namespace

#endif