burgers_rewienski.h

#ifndef __BURGERS_REWIENSKI__
#define __BURGERS_REWIENSKI__

#include "burgers.h"

namespace PHiLiP {
namespace Physics {
template <int dim, int nspecies, int nstate, typename real>
class BurgersRewienski : public Burgers <dim, nspecies, nstate, real>
{
protected:
    // For overloading the virtual functions defined in PhysicsBase
    using PhysicsBase<dim,nspecies,nstate,real>::dissipative_flux;
    using PhysicsBase<dim,nspecies,nstate,real>::source_term;
    using PhysicsBase<dim,nspecies,nstate,real>::boundary_face_values;
public:
    BurgersRewienski(
            const Parameters::AllParameters *const                    parameters_input,
            const double                                              rewienski_a,
            const double                                              rewienski_b,
            const bool                                                rewienski_manufactured_solution,
            const bool                                                convection,
            const bool                                                diffusion,
            const dealii::Tensor<2,3,double>                          input_diffusion_tensor = Parameters::ManufacturedSolutionParam::get_default_diffusion_tensor(),
            std::shared_ptr< ManufacturedSolutionFunction<dim,nspecies,real> > manufactured_solution_function = nullptr);

    const double rewienski_a;

    const double rewienski_b;

    const bool rewienski_manufactured_solution;

    std::array<real,nstate> source_term (
            const dealii::Point<dim,real> &pos,
            const std::array<real,nstate> &solution,
            const real current_time) const override;

    void boundary_face_values (
            const int /*boundary_type*/,
            const dealii::Point<dim, real> &/*pos*/,
            const dealii::Tensor<1,dim,real> &/*normal*/,
            const std::array<real,nstate> &/*soln_int*/,
            const std::array<dealii::Tensor<1,dim,real>,nstate> &/*soln_grad_int*/,
            std::array<real,nstate> &/*soln_bc*/,
            std::array<dealii::Tensor<1,dim,real>,nstate> &/*soln_grad_bc*/) const override;

};
} // Physics namespace
} // PHiLiP namespace

#endif