ACFDTD22.hpp

Go to the documentation of this file.

#ifndef ACFDTD22_HPP
#define ACFDTD22_HPP

#include "constant.hpp"

#include "Model.hpp"
#include "HardPlaneWaveSource.hpp"
#include "ObservationRegion.hpp"
#include "PointSensor.hpp"
#include "LineSensor.hpp"
#include "UniqueCoefficientContainer.hpp"

class ACFDTD22 {
private:
  unsigned short **bya, **byb, **bza, **bzb;
  unsigned short **pra, **prb, **prc, **prd, **pre;
  unsigned short **aya, **ayb, **aza, **azb;
  unsigned short **uya, **uyb, **uyc, **uza, **uzb, **uzc;

  void parseModel(const Model &model, const Grid &grid);

protected:
  UniqueCoefficientContainer container; 
  std::string sim_name;
  double f_c; 
  double S; 
  double dt; 
  unsigned nTimeSteps; 
public:
  ACFDTD22();
  
  ACFDTD22(const double frequency, const double cfl_number); 

  virtual void calcCoefficients(const Model &model, const Grid &grid);

  virtual void execute(const Grid &grid, 
                       const CartPosition &sig_loc, 
                       const Function *signal, 
                       const double time_span,
                       const std::vector<Sensor *> &sen_vec,
                       const unsigned negY=BC_PRESS,
                       const unsigned posY=BC_PRESS,
                       const unsigned negZ=BC_PRESS,
                       const unsigned posZ=BC_PRESS);
  
  virtual void execute(const Model &model,
               const Grid  &grid,
               const HardPlaneWaveSource &source,
               const double time_span,
               const std::vector<Sensor *> &sen_vec,
               const unsigned negY=BC_PRESS,
               const unsigned posY=BC_PRESS,
               const unsigned negZ=BC_PRESS,
               const unsigned posZ=BC_PRESS);

  virtual double getDeltaT() const ;
  virtual unsigned getnTimeSteps() const;
  virtual ~ACFDTD22();
};



#endif