ImmersedBoundaries.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _IMMERSED_BOUNDARIES_H
#define _IMMERSED_BOUNDARIES_H
 
#ifndef Real
#define Real double 
#endif
 
#include <vector>
 
class Vec3;
 
namespace utl {
  class Point;
}
 
class ElementBlock;
 
typedef std::vector<Real>      RealArray;
typedef std::vector<RealArray> Real2DMat;
typedef std::vector<Real2DMat> Real3DMat;
typedef std::vector<utl::Point> PointVec;
 
 
namespace Immersed 
{
  class Geometry
  {
  protected:
    Geometry() {}
  public:
    virtual ~Geometry() {}
 
    virtual double Alpha(double X, double Y, double Z) const { return 0.0; }
    virtual double Alpha(const Vec3& X) const;
 
    virtual ElementBlock* tesselate() const { return nullptr; }
  };
 
  bool getQuadraturePoints(const Geometry& geo,
                           const std::vector<PointVec>& elmCorner,
                           int max_depth, int p,
                           Real3DMat& quadPoints, ElementBlock* grid = nullptr);
 
  bool getQuadraturePoints(const Geometry& geo,
                           const utl::Point& X1, const utl::Point& X2,
                           const utl::Point& X3, const utl::Point& X4,
                           int max_depth, int nGauss,
                           RealArray& GP1, RealArray& GP2,
                           RealArray& GPw, ElementBlock* grid = nullptr);
 
  bool getQuadraturePoints(const Geometry& geo,
                           const utl::Point& X1, const utl::Point& X2,
                           const utl::Point& X3, const utl::Point& X4,
                           const utl::Point& X5, const utl::Point& X6,
                           const utl::Point& X7, const utl::Point& X8,
                           int max_depth, int nGauss,
                           RealArray& GP1, RealArray& GP2, RealArray& GP3,
                           RealArray& GPw);
 
  enum Stab
  {
    NO_STAB = 0,
    ALL_INTERFACES = 1,
    SUBDIV_INTERFACES = 2
  };
 
  extern int stabilization; 
 
  extern bool plotCells; 
}
 
#endif