physics_post_processor.h

#ifndef __PHYSICS_POST__
#define __PHYSICS_POST__

#include <deal.II/numerics/data_postprocessor.h>

#include "physics/physics.h"
#include "physics/model.h"

namespace PHiLiP {
namespace Postprocess {

template <int dim>
class PostprocessorFactory
{
public:
    static std::unique_ptr< dealii::DataPostprocessor<dim> > create_Postprocessor(const Parameters::AllParameters *const parameters_input);
};


template <int dim, int nstate>
class PhysicsPostprocessor : public dealii::DataPostprocessor<dim>
{
public:
    explicit PhysicsPostprocessor (const Parameters::AllParameters *const parameters_input);

    std::shared_ptr < Physics::ModelBase<dim, nstate, double > > model;

    std::shared_ptr < Physics::PhysicsBase<dim, nstate, double > > physics;

    virtual void evaluate_vector_field (const dealii::DataPostprocessorInputs::Vector<dim> &inputs, std::vector<dealii::Vector<double>> &computed_quantities) const override;
    virtual void evaluate_scalar_field (const dealii::DataPostprocessorInputs::Scalar<dim> &inputs, std::vector<dealii::Vector<double>> &computed_quantities) const override;
    virtual std::vector<std::string> get_names () const override;
    virtual std::vector<dealii::DataComponentInterpretation::DataComponentInterpretation> get_data_component_interpretation () const override;
    virtual dealii::UpdateFlags get_needed_update_flags () const override;
};

} // Postprocess namespace
} // PHiLiP namespace

#endif