sim_fdm_scheme_iso_it.hpp

/* sim_fdm_scheme_iso_it.hpp */
/* $Id$ */
#ifndef __sim_fdm_scheme_iso_it_hpp
#define __sim_fdm_scheme_iso_it_hpp

#include <s2proc/sim_service.hpp>

namespace sim
{
  namespace fdm
  {
    struct GridIso
    {
    public:
      
      GridIso():
    n(0), m(0)  
      {/* nothing to do */}
      
      GridIso(const int i_n, const int i_m);
      
      int n;
      int m;
      std::vector<CoordT> lambdah;
      std::vector<CoordT> lambdav;
      std::vector<CoordT> muh;
      std::vector<CoordT> muv;
      std::vector<CoordT> rho;
    };
    
    class Snapshot;
    
    class SchemeIsoIt
    {
    public:
      
      SchemeIsoIt(const GridIso &grid,
          const Snapshot &fp,
          const Snapshot &f,
          Snapshot &fn,
          const TimeT &i_dt);
      
      void init(const GridIso &grid,
        const Snapshot &fp,
        const Snapshot &f,
        Snapshot &fn,
        const TimeT &i_dt);
      
      
      SchemeIsoIt& operator++();
      
      void step(const size_t &x_step, const size_t &z_step);
      
      void reset();
      
      
      void do_scheme_full_form();
      
      void do_scheme_short_form();
    
      void do_stacey_top();
      
      void do_stacey_bottom();
      
      void do_stacey_left();
      
      void do_stacey_right();
      
      void do_corners_boundary_conds();
      
      void do_top_left_corner();
      
      void do_top_right_corner();
      
      void do_bottom_left_corner();
      
      void do_bottom_right_corner();
      
      // protected:
      
      size_t cur_ind;
      size_t cur_n;
      size_t cur_m;
      size_t x_size;
      size_t z_size;
      const CoordT *pdx;
      const CoordT *pdz;
      TimeT dt;
      TimeT dt2;
      const AmplT *pup;
      const AmplT *pwp;
      const AmplT *pu;
      const AmplT *pw;
      AmplT *pun;
      AmplT *pwn;
      const CoordT *pden;
      const CoordT *plambdah;
      const CoordT *plambdav;
      const CoordT *pmuh;
      const CoordT *pmuv;
      
    }; // class SchemeIsoIt
    
  }; // finite differences modeling namespace
  
}; // simulation namespace    

#endif /* sim_fdm_scheme_iso_it.hpp */