Brick.hxx

Go to the documentation of this file.

#ifndef _BRICK_HXX_
#define _BRICK_HXX_

#include "MureHeaders.hxx"
#include "Cylinder.hxx"
#include <string>
using namespace std;


class Brick : public Shape
{
 public :

    Brick(double HalfX=0.01,double HalfY=0.01,double HalfZ=0.01, double Signe=-1);
    

    Brick(double *Origin, double *Ux, double *Uy, double *Uz, double HalfX=0.01,double HalfY=0.01,double HalfZ=0.01, double Signe=-1);
    

    Brick(const char* Infinite, double HalfX=0.01,double HalfY=0.01, double Signe=-1);
    unsigned IsBrick() {return 1;} 

    
   double GetHalfX() {return fHalfSize[0];} 
   double GetHalfY() {return fHalfSize[1];} 
   double GetHalfZ() {return fHalfSize[2];} 
   const double* GetOrigin() {return fOrigin;} 
   double* GetUx() {return fUx;} 
   double* GetUy() {return fUy;} 
   double* GetUz() {return fUz;} 
  

    void SetPeriodicBoundary(bool state=true,string TopBottPlaneBC="mirror");
   

    virtual void Translate(double dx,double dy,double dz,int count=0);

    virtual void Rotate(double phi,double theta=0,double psi=0,double *center=0,int count=0);

    unsigned Inclus(Shape_ptr OtherShape); 
    unsigned Disjoint(Shape_ptr OtherShape); 
     

    bool PointInShape(double *P,int border);
    
    string  Print(); 
    void    DefineSurface(); 
    
    double MinDistanceFrom2OppositeFace(double *N,int whichplane,double *S);
    
    Shape_ptr Clone(int count=0); 
    void Copy(Shape_ptr B); 
    
    virtual int NumberOfShapes(){return 6;}
 private:
    Brick(double HalfX,double HalfY,double HalfZ, double Signe,bool Infinite); 
  void Init(double HalfX,double HalfY,double HalfZ); 

    unsigned CornerInBrick(Brick *brick,int all,int border);

    unsigned CornerInCylinder(Cylinder *cyl,int all,int border);
    
    double   Dist2Segment(Cylinder *cyl,double Segment[6][2],double C0[6][2],int &segmentNum);
    void FindBrickProjection(Cylinder *cyl,double Segment[6][2],double C0[6][2],double Origin[2]);
    
    unsigned GeneralIntersectionOf2Brick(Brick *brick);

    unsigned LineCrossPlane(double *X,double *N,double *P, double *U,double *D);
    unsigned CylinderCrossBrick(Cylinder *Cyl);
    unsigned PointInSegment(double *X,double *A,double *B,int border);
    
    unsigned TwoFaces3Segments(Brick *brick,double *N,int whichSize,double *C0,double C1[3][3]);
    
    unsigned TwoFaces2Segments(Brick *brick,double *N,int whichSize,double *C0,double C1[2][3]);

    double GetVirtualSphereRadius(); 
    double* GetVirtualSphereCenter(){return fOrigin;} 
    
    void CalculNot(); 
    
    double MinDist1Point2Face(double *U0,double *U1,double *U2,double *P0,double *S,int whichplane);
    double MinDist1Point2FaceInfinite(double *U0,double *U1,double *P0,double *S,int whichplane);

 private:
    Surface* DefinePlane(double *N,double D);

    double fHalfSize[3];
    double fOrigin[3]; 
    double fUx[3]; 
    double fUy[3]; 
    double fUz[3]; 

};

#endif