TransformIF.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
*      _______              __
*     / ___/ /  ___  __ _  / /  ___
*    / /__/ _ \/ _ \/  V \/ _ \/ _ \
*    \___/_//_/\___/_/_/_/_.__/\___/
*    Please refer to Copyright.txt, in Chombo's root directory.
*/
#endif

#ifndef _TRANSFORMIF_H_
#define _TRANSFORMIF_H_

#include "MayDay.H"
#include "RealVect.H"
#include "Vector.H"

#include "BaseIF.H"

#include "NamespaceHeader.H"


class TransformIF: public BaseIF
{
public:

  TransformIF(const BaseIF& a_impFunc);

  TransformIF(const TransformIF& a_inputIF);

  virtual ~TransformIF();


  void translate(const RealVect& a_trans);


  void scale(const Real& a_scale);


  void scale(const RealVect& a_scale);


  void rotate(const Real&     a_angle,
              const RealVect& a_point = RealVect::Zero,
              const RealVect& a_axis  = BASISREALV(0));


  void rotate(const RealVect& a_axis1,
              const RealVect& a_axis2,
              const RealVect& a_point = RealVect::Zero);


  virtual Real value(const RealVect& a_point) const;

  Real value(const IndexTM<Real,GLOBALDIM>& a_point) const;

  virtual IndexTM<Real,GLOBALDIM> normal(const IndexTM<Real,GLOBALDIM>& a_point) const;

  virtual Vector< IndexTM<Real,GLOBALDIM> > gradNormal(const IndexTM<Real,GLOBALDIM>& a_point) const;

  virtual BaseIF* newImplicitFunction() const;

  virtual bool fastIntersection(const RealVect& a_lo,
                                const RealVect& a_hi) const ;

  virtual GeometryService::InOut InsideOutside(const RealVect& a_lo,
                                               const RealVect& a_hi) const ;

protected:

  TransformIF(const BaseIF& a_impFunc,
              const Real    a_transform[SpaceDim+1][SpaceDim+1],
              const Real    a_invTransform[SpaceDim+1][SpaceDim+1]);

  // transformation (in homogeneous coordinates)
  Real m_transform[SpaceDim+1][SpaceDim+1];

  // inverse transformation (in homogeneous coordinates)
  Real m_invTransform[SpaceDim+1][SpaceDim+1];

  BaseIF* m_impFunc; // implicit function to transform

  void vectorMultiply(RealVect&       m_outPoint,
                      const Real      m_intrans[SpaceDim+1][SpaceDim+1],
                      const RealVect& m_inPoint) const;

  void vectorMultiply(IndexTM<Real,GLOBALDIM>&       m_outPoint,
                      const Real                     m_intrans[SpaceDim+1][SpaceDim+1],
                      const IndexTM<Real,GLOBALDIM>& m_inPoint) const;

  void matrixIdentity(Real m_trans[SpaceDim+1][SpaceDim+1]) const;

  void matrixMultiply(Real       m_outTrans[SpaceDim+1][SpaceDim+1],
                      const Real m_intrans1[SpaceDim+1][SpaceDim+1],
                      const Real m_intrans2[SpaceDim+1][SpaceDim+1]) const;

  void matrixTranslate(Real            m_trans[SpaceDim+1][SpaceDim+1],
                       const RealVect& m_translate) const;

  void matrixScale(Real            m_trans[SpaceDim+1][SpaceDim+1],
                   const RealVect& m_scale) const;

private:
  TransformIF()
  {
    MayDay::Error("TransformIF uses strong construction");
  }

  void operator=(const TransformIF& a_inputIF)
  {
    MayDay::Error("TransformIF doesn't allow assignment");
  }
};

#include "NamespaceFooter.H"
#endif

---