ASMs2DIB.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ASM_S2D_IB_H
#define _ASM_S2D_IB_H
 
#include "ASMs2D.h"
 
namespace Immersed { class Geometry; }
 
 
class ASMs2DIB : public ASMs2D
{
  class Intersected : public ASMs2D::InterfaceChecker
  {
    bool            myAll;   
    bool            alsoSW;  
  public:
    explicit Intersected(const ASMs2DIB& pch, bool all = false, bool sw = false)
      : InterfaceChecker(pch), myAll(all), alsoSW(sw) {}
    virtual ~Intersected() {}
    virtual short int hasContribution(int, int I, int J, int = -1) const;
  };
 
public:
  ASMs2DIB(unsigned char n_s = 2, unsigned char n_f = 1, int max_depth = 5);
  ASMs2DIB(const ASMs2DIB& patch, unsigned char n_f = 0);
  virtual ~ASMs2DIB();
 
  virtual void addHole(double R, double Xc, double Yc);
  virtual void addHole(double R, double X1, double Y1, double X2, double Y2);
 
  virtual bool setGeometry(RealFunc* f, double power, double threshold);
 
  virtual ElementBlock* immersedGeometry(char* name) const;
 
  virtual void getNoIntPoints(size_t& nPt, size_t& nIPt);
 
  virtual bool generateFEMTopology();
 
  bool isIntersected(int iel, bool checkIfInDomainOnly = false) const;
 
  using ASMs2D::integrate;
  virtual bool integrate(Integrand& integrand,
                         GlobalIntegral& glbInt, const TimeDomain& time);
 
  virtual void filterResults(Matrix& field, const ElementBlock* grid) const;
 
private:
  Immersed::Geometry* myGeometry; 
  ElementBlock*       myLines;    
 
  Real3DMat quadPoints; 
  int       maxDepth;   
};
 
#endif