MathBrick.hxx

Go to the documentation of this file.

/*
  This file is part of MURE,
  Copyright (C) 2007-2021 MURE developers.
 
  MURE is free software: you can redistribute it and/or modify
  it under the terms of the GNU Lesser General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  MURE is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU Lesser General Public License for more details.
 
  You should have received a copy of the GNU Lesser General Public License
  along with MURE.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef _MathBRICK_HXX_
#define _MathBRICK_HXX_
 
#include <string>
#include <algorithm>
#include <iosfwd>
#include <vector>
 
#include "Shape.hxx"
#include "MathCylinder.hxx"
#include "MCSource.hxx"
 
class MathCylinder;
 
using namespace std;
 
 
class MathBrick : public Shape
{
  public :
 
    MathBrick(double HalfX = 0.01, double HalfY = 0.01, double HalfZ = 0.01, int Sign = - 1);
 
 
    MathBrick(double *Origin, double *Ux, double *Uy, double *Uz, double HalfX = 0.01, double HalfY = 0.01, double HalfZ = 0.01, int Sign = - 1);
 
 
    MathBrick(const char *Infinite, double HalfX = 0.01, double HalfY = 0.01, int Sign = - 1);
    Shape_ptr Clone(int count = 0) override; 
    void Copy(Shape_ptr B) override; 
 
    unsigned IsBrick() override
    {
        return 1;   
    }
    double GetHalfX()
    {
        return fHalfSize[0];   
    }
    double GetHalfY()
    {
        return fHalfSize[1];   
    }
    double GetHalfZ()
    {
        return fHalfSize[2];   
    }
    vector < double > GetOrigin()
    {
        return fOrigin;   
    }
    vector < double > GetUx()
    {
        return fUx;   
    }
    vector < double > GetUy()
    {
        return fUy;   
    }
    vector < double > GetUz()
    {
        return fUz;   
    }
 
 
    void SetPeriodicBoundary(bool state = true, string TopBottPlaneBC = "mirror") override;
 
 
    void Translate(double dx, double dy, double dz, int count = 0) override;
 
    void Rotate(double phi, double theta = 0, double psi = 0, double *center = nullptr, int count = 0) override;
 
    unsigned IsIncluded(Shape_ptr OtherShape) override; 
    unsigned IsDisjoint(Shape_ptr OtherShape) override; 
 
 
    bool PointInShape(double *P, int border) override;
 
 
    double MinDistanceFrom2OppositeFace(vector < double > N, int whichplane, double *S);
 
    void CalculNot() override; 
 
    int NumberOfShapes() override
    {
        if(fInfinite) return 4 ;
        return 6;
    }
 
  protected:
 
    MathBrick(double HalfX, double HalfY, double HalfZ, int Sign, bool Infinite); 
 
    MathBrick *CreateNewInstance(int sign = 1) override; 
 
    void Init(double HalfX, double HalfY, double HalfZ); 
 
    unsigned CornerInBrick(MathBrick *brick, int all, int border);
 
 
    unsigned CornerInCylinder(MathCylinder *cyl, int all, int border);
 
    double   Dist2Segment(MathCylinder *cyl, double Segment[6][2], double C0[6][2], int &segmentNum);
    void FindBrickProjection(MathCylinder *cyl, double Segment[6][2], double C0[6][2], double Origin[2]);
 
    unsigned GeneralIntersectionOf2Brick(MathBrick *brick);
    unsigned LineCrossPlane(double *X, vector < double > N, double *P, vector < double > U, double *D);
    unsigned CylinderCrossBrick(MathCylinder *Cyl);
    unsigned PointInSegment(double *X, double *A, double *B, int border);
 
    unsigned TwoFaces3Segments(MathBrick *brick, vector < double > N, int whichSize, double *C0, double C1[3][3]);
 
    unsigned TwoFaces2Segments(MathBrick *brick, vector < double > N, int whichSize, double *C0, double C1[2][3]);
 
    double GetVirtualSphereRadius() override; 
    double *GetVirtualSphereCenter() override; 
 
 
    double MinDist1Point2Face(vector < double > U0, vector < double > U1, vector < double > U2, double *P0, double *S, int whichplane);
    double MinDist1Point2FaceInfinite(vector < double > U0, vector < double > U1, double *P0, double *S, int whichplane);
 
    double fHalfSize[3];
    vector < double > fOrigin; 
    vector < double > fUx; 
    vector < double > fUy; 
    vector < double > fUz; 
 
};
 
#endif