LinearGaussImproved.hpp

Go to the documentation of this file.

#ifndef LINEAR_GAUSS_IMPROVED
#define LINEAR_GAUSS_IMPROVED 1

#include "../../common/equation_of_state.hpp"
#include "../spatial_reconstruction.hpp"
#include <cmath>
#include "../../../misc/universal_error.hpp"
#include "../GhostPointGenerator.hpp"
#include "../../../misc/serializable.hpp"

class LinearGaussImproved : public SpatialReconstruction
{
public:

        LinearGaussImproved
                (EquationOfState const& eos,
                GhostPointGenerator const& ghost,
                bool slf = true,
                double delta_v = 0.2,
                double theta = 0.5,
                double delta_P = 0.7,
                const vector<string>& calc_tracers =
                vector<string>(),
                string skip_key=string(),
                bool SR=false);

        void operator() (const Tessellation& tess,const vector<ComputationalCell>& cells,double time,
                vector<pair<ComputationalCell, ComputationalCell> > &res,TracerStickerNames const& tracerstickersnames,CacheData const& cd)const;

        ComputationalCell Interp(ComputationalCell const& cell, size_t cell_index, Vector2D const& cm, Vector2D const& target)const;

        vector<Slope>& GetSlopes(void)const;

        vector<Slope>& GetSlopesUnlimited(void)const;

private:
        EquationOfState const& eos_;
        GhostPointGenerator const& ghost_;
        mutable vector<Slope> rslopes_;
        mutable vector<Slope> naive_rslopes_;
        const bool slf_;
        const double shockratio_;
        const double diffusecoeff_;
        const double pressure_ratio_;
        const vector<string> calc_tracers_;
        const string skip_key_;
        mutable vector<size_t> to_skip_;
        bool SR_;

        LinearGaussImproved
                (const LinearGaussImproved& origin);

        LinearGaussImproved& operator=
                (const LinearGaussImproved& origin);

};
#endif //LINEAR_GAUSS_IMPROVED