waveprop/pmlTransformation.hh

Go to the documentation of this file.

#ifndef PMLTransformation_hh
#define PMLTransformation_hh

#include "basics.hh"
#include "toolbox.hh"
#include "formula.hh"
#include "geometry.hh"
#include "space.hh"
namespace concepts {

  // ************************************************* FrmE_PMLTransformation **

  template<uint dim>
  class FrmE_PMLTransformation : 
    public ElementFormula< typename GeneralMapping<Cmplx,dim>::Type >
  {
  public:
    FrmE_PMLTransformation(const Formula<Real>& sigma,
                           const Formula<Real>& sigmaDeriv);

    virtual ~FrmE_PMLTransformation();

    void addSector(const Real2d vtx0, const Real2d vtx1, 
                   const Real2d vtx2, const Real2d vtx3);
    void addSector(const Sequence<Real2d> vtx);
    
    typename GeneralMapping<Cmplx,dim>::Type operator() 
      (const ElementWithCell<Real>& elm, const Real p, 
       const Real t = 0.0) const;
    typename GeneralMapping<Cmplx,dim>::Type operator() 
      (const ElementWithCell<Real>& elm, const Real2d& p,
       const Real t = 0.0) const;
    typename GeneralMapping<Cmplx,dim>::Type operator() 
      (const ElementWithCell<Real>& elm, const Real3d& p,
       const Real t = 0.0) const;
    virtual FrmE_PMLTransformation<dim>* clone() const;

  private:
    Sequence<InverseVertexQuadSector2d> mapFromSectorToRefSector_;
    std::auto_ptr<const Formula<Real> > sigma_,sigmaDeriv_;
    mutable const ElementWithCell<Real>* lastElm_;
    mutable uint lastSector_;
    Real2d localCoord_(const ElementWithCell<Real>& elm,
                       const Real2d& p) const;
    FrmE_PMLTransformation
    (const Formula<Real>* sigma, const Formula<Real>* sigmaDeriv,
     Sequence<InverseVertexQuadSector2d> mapFromSectorToRefSector);

    static const uint OUTSIDEPMLSECTOR_ = -1;
  };
} // namespace concepts

#endif //PMLTransformation_hh