amr.hpp

Go to the documentation of this file.

#ifndef AMR_HPP
#define AMR_HPP 1

#include "computational_cell_2d.hpp"
#include "extensive.hpp"
#include "../common/equation_of_state.hpp"
#include "OuterBoundary.hpp"
#include "../../tessellation/tessellation.hpp"
#include "../../tessellation/ConvexHull.hpp"
#include "clipper.hpp"
#include "../test_2d/main_loop_2d.hpp"
#include "../../tessellation/polygon_overlap_area.hpp"
#include <boost/scoped_ptr.hpp>
#include "interpolations/LinearGaussImproved.hpp"

class AMRCellUpdater
{
public:

        virtual ComputationalCell ConvertExtensiveToPrimitve(const Extensive& extensive,const EquationOfState& eos,
                double volume,ComputationalCell const& old_cell,TracerStickerNames const& tracerstickernames) const = 0;

        virtual ~AMRCellUpdater(void);
};

class AMRExtensiveUpdater
{
public:

        virtual Extensive ConvertPrimitveToExtensive(const ComputationalCell& cell, const EquationOfState& eos,
                double volume,TracerStickerNames const& tracerstickernames) const = 0;

        virtual ~AMRExtensiveUpdater(void);
};

class SimpleAMRExtensiveUpdater : public AMRExtensiveUpdater
{
public:
        Extensive ConvertPrimitveToExtensive(const ComputationalCell& cell, const EquationOfState& eos,
                double volume, TracerStickerNames const& tracerstickernames) const;
};

class SimpleAMRCellUpdater : public AMRCellUpdater
{
private:
        const vector<string>  toskip_;
public:
  explicit SimpleAMRCellUpdater(vector<string> toskip);

        ComputationalCell ConvertExtensiveToPrimitve(const Extensive& extensive, const EquationOfState& eos,
                double volume, ComputationalCell const& old_cell,TracerStickerNames const& tracerstickernames) const;
};

class SimpleAMRExtensiveUpdaterSR : public AMRExtensiveUpdater
{
public:
        Extensive ConvertPrimitveToExtensive(const ComputationalCell& cell, const EquationOfState& eos,
                double volume, TracerStickerNames const& tracerstickernames) const;
};

class SimpleAMRCellUpdaterSR : public AMRCellUpdater
{
private:
        const double G_;
        const vector<string>  toskip_;
public:
        explicit SimpleAMRCellUpdaterSR(double G,vector<string> toskip);

        ComputationalCell ConvertExtensiveToPrimitve(const Extensive& extensive, const EquationOfState& eos,
                double volume, ComputationalCell const& old_cell, TracerStickerNames const& tracerstickernames) const;
};

class CellsToRemove
{
public:
        virtual std::pair<vector<size_t>,vector<double> > ToRemove(Tessellation const& tess,
                vector<ComputationalCell> const& cells,double time,
                TracerStickerNames const& tracerstickernames)const=0;

        virtual ~CellsToRemove(void);
};

class CellsToRefine
{
public:
        virtual vector<size_t> ToRefine(Tessellation const& tess, vector<ComputationalCell> const& cells, double time,
                TracerStickerNames const& tracerstickernames)const = 0;
        
        virtual ~CellsToRefine(void);
};

class AMR : public Manipulate
{
protected:

#ifdef RICH_MPI
#endif // RICH_MPI
  void GetNewPoints
  (vector<size_t> const& ToRefine,
   Tessellation const& tess,
   vector<std::pair<size_t, Vector2D> > &NewPoints, 
   vector<Vector2D> &Moved,
   OuterBoundary const& obc
#ifdef RICH_MPI
          , vector<Vector2D> const& proc_chull
#endif
          )const;
#ifdef RICH_MPI

  vector<size_t> RemoveNearBoundaryPoints(vector<size_t> const&ToRemove,
          Tessellation const& tess,vector<double> &merits)const;
#endif

  void GetNewPoints2
          (vector<size_t> const& ToRefine,
                  Tessellation const& tess,
                  vector<std::pair<size_t, Vector2D> > &NewPoints,
                  vector<Vector2D> &Moved,
                  OuterBoundary const& obc)const;

public:
        virtual void operator() (hdsim &sim) = 0;
#ifdef RICH_MPI
#endif // RICH_MPI
        virtual void UpdateCellsRefine(Tessellation &tess,
                OuterBoundary const& obc, vector<ComputationalCell> &cells, EquationOfState const& eos,
                vector<Extensive> &extensives, double time,
#ifdef RICH_MPI
                Tessellation const& proctess,
#endif
                TracerStickerNames const& tracerstickernames,
                CacheData const& cd,
                PhysicalGeometry const& pg)const = 0;
#ifdef RICH_MPI
#endif // RICH_MPI
        virtual void UpdateCellsRemove(Tessellation &tess,
                OuterBoundary const& obc, vector<ComputationalCell> &cells, vector<Extensive> &extensives,
                EquationOfState const& eos, double time,
#ifdef RICH_MPI
                Tessellation const& proctess,
#endif
                TracerStickerNames const& tracerstickernames,
                CacheData const& cd,
                PhysicalGeometry const& pg)const = 0;
        virtual ~AMR(void);
};

class ConservativeAMR : public AMR
{
private:
        CellsToRefine const& refine_;
        CellsToRemove const& remove_;
  SimpleAMRCellUpdater scu_;
  SimpleAMRExtensiveUpdater seu_;
  const bool periodic_;
        AMRCellUpdater* cu_;
        AMRExtensiveUpdater* eu_;
        LinearGaussImproved *interp_;

        ConservativeAMR(ConservativeAMR const& other);

        ConservativeAMR& operator=(ConservativeAMR const& other);

public:
        void operator() (hdsim &sim);

        ConservativeAMR
        (CellsToRefine const& refine,
         CellsToRemove const& remove,
                bool periodic = false,
         LinearGaussImproved *slopes = 0,
         AMRCellUpdater* cu=0,
         AMRExtensiveUpdater* eu=0);

        void UpdateCellsRefine(Tessellation &tess,
                OuterBoundary const& obc, vector<ComputationalCell> &cells,EquationOfState const& eos,
                vector<Extensive> &extensives,double time,
#ifdef RICH_MPI
                Tessellation const& proctess,
#endif
                TracerStickerNames const& tracerstickernames,
                CacheData const& cd,
                PhysicalGeometry const& pg)const;

        void UpdateCellsRemove(Tessellation &tess,
                OuterBoundary const& obc, vector<ComputationalCell> &cells, vector<Extensive> &extensives,
                EquationOfState const& eos,double time,
#ifdef RICH_MPI
                Tessellation const& proctess,
#endif
                TracerStickerNames const& tracerstickernames,
                CacheData const& cd,
                PhysicalGeometry const& pg)const;
};

class NonConservativeAMR : public AMR
{
private:
        CellsToRefine const& refine_;
        CellsToRemove const& remove_;
  SimpleAMRCellUpdater scu_;
  SimpleAMRExtensiveUpdater seu_;
  LinearGaussImproved *interp_;
        AMRExtensiveUpdater* eu_;

        NonConservativeAMR(NonConservativeAMR const& other);

        NonConservativeAMR& operator=(NonConservativeAMR const& other);

public:
        void operator() (hdsim &sim);

  NonConservativeAMR
  (CellsToRefine const& refine,
   CellsToRemove const& remove,
          LinearGaussImproved *slopes = 0,
   AMRExtensiveUpdater* eu = 0);

        void UpdateCellsRefine(Tessellation &tess,
                OuterBoundary const& obc, vector<ComputationalCell> &cells, EquationOfState const& eos,
                vector<Extensive> &extensives, double time,
#ifdef RICH_MPI
                Tessellation const& proctess,
#endif
                TracerStickerNames const& tracerstickernames,
                CacheData const& cd,
                PhysicalGeometry const& pg)const;

        void UpdateCellsRemove(Tessellation &tess,
                OuterBoundary const& obc, vector<ComputationalCell> &cells, vector<Extensive> &extensives,
                EquationOfState const& eos, double time,
#ifdef RICH_MPI
                Tessellation const& proctess,
#endif
                TracerStickerNames const& tracerstickernames,
                CacheData const& cd,
                PhysicalGeometry const& pg)const;
};

#endif // AMR_HPP