sim_fdm_snapshot.hpp

/* sim_fdm_snapshot.hpp */
/* $Id$ */
#ifndef __sim_fdm_snapshot_hpp
#define __sim_fdm_snapshot_hpp

#include <vector>
#include <s2proc/sim_service.hpp>
#include <s2proc/sim_point_2d.hpp>

namespace sim
{
  namespace fdm
  {
    class Snapshot
    {
    public:
      
      class Area
      {
      public:
    
    Area() {/* nothing to do*/}
    
    Area(const Point2D &top_left_coords,
         const std::vector<CoordT> &dx_vec,
         const std::vector<CoordT> &dz_vec);
    
    Area(const Point2D &top_left_coords,
         const Point2D &bottom_right_coords,
         const Point2D &coord_steps);
    
    const Point2D& topLeftCoords() const {return m_top_left_coords;}
    
    const std::vector<CoordT>& stepVec(const CoordsMask &coord) const;
    
      protected:
    
    Point2D m_top_left_coords;    
    std::vector<CoordT> m_dx_vec; 
    std::vector<CoordT> m_dz_vec; 
      };
      
      Snapshot();
      
      Snapshot(const Area &area);
      
      
      std::vector<AmplT>& operator()(const CoordsMask &coord);
      
      const std::vector<AmplT>& operator()(const CoordsMask &coord) const;
      
      
      int id(const Point2D &coords) const;
      
      
      AmplT* operator()(const size_t x_id,
            const size_t z_id,
            const CoordsMask &coord);
      
      const AmplT* operator()(const size_t x_id,
                  const size_t z_id,
                  const CoordsMask &coord) const;
      
      const Area& area() const {return m_area;}
      
      
      static bool test_to_file(const Snapshot &snapshot, const std::string &fname);
      
    private:
      
      Area m_area;                     
      std::vector<AmplT> m_u_snapshot; 
      std::vector<AmplT> m_w_snapshot; 
    };
    
  }; // finite differences modeling namespace
  
}; // simulation namespace

#endif /* sim_fdm_snapshot.hpp */