IntersectionOfWalls.h

Go to the documentation of this file.

// This file is part of the MercuryDPM project (https://www.mercurydpm.org).
// Copyright (c), The MercuryDPM Developers Team. All rights reserved.
// License: BSD 3-Clause License; see the LICENSE file in the root directory.
 
#ifndef INTERSECTIONOFWALLS_H
#define INTERSECTIONOFWALLS_H
 
#include <vector>
#include "BaseWall.h"
#include "InfiniteWall.h"
#include "Math/Vector.h"
 
class IntersectionOfWalls : public BaseWall
{
public:
    struct normalAndPosition
    {
        Vec3D normal;
        Vec3D position;
    };
    
    IntersectionOfWalls();
    
    IntersectionOfWalls(const IntersectionOfWalls& other);
    
    IntersectionOfWalls(const std::vector<normalAndPosition>& walls, const ParticleSpecies* species);
    
    ~IntersectionOfWalls() override;
    
    IntersectionOfWalls& operator=(const IntersectionOfWalls& other);
    
    IntersectionOfWalls* copy() const override;
    
    void clear();
    
    void setSpecies(const ParticleSpecies* species);
    
    void setHandler(WallHandler* wallHandler) override;
    
    unsigned int getNumberOfObjects();
    
    void addObject(Vec3D normal, Vec3D point);
    
    void addObject(Quaternion orientation, Vec3D position);
    
    void add3PointObject(Vec3D PointA, Vec3D PointB, Vec3D PointC);
    
    void setPointsAndLines(unsigned int n);
    
    void addTetraSTL(Vec3D PointA, Vec3D PointB, Vec3D PointC, Vec3D WallNormal, Mdouble Thickness, int wallidentifier);
    
    void addTetra(const Vec3D& PointA, const Vec3D& PointB, const Vec3D& PointC, Mdouble& Thickness);
    
    
    void addPlate(const Vec3D& PointA, const Vec3D& PointB, const Vec3D& PointC, const Vec3D& WallNormal,
                  const Mdouble& Thickness, int wallidentifier);
    
    MERCURYDPM_DEPRECATED
    void addObject(Vec3D normal, Mdouble position);
    
    void createOpenPrism(std::vector<Vec3D> points, Vec3D prismAxis);
    
    void createPrism(std::vector<Vec3D> points, Vec3D prismAxis);
    
    void createOpenPrism(std::vector<Vec3D> points);
    
    void createPrism(std::vector<Vec3D> points);
    
    bool getDistanceAndNormal(const BaseParticle& p, Mdouble& distance, Vec3D& normal_return) const override;
    
    bool getDistanceAndNormal(const Vec3D& position, Mdouble wallInteractionRadius, Mdouble& distance,
                              Vec3D& normal_return) const;
    
//    void move(const Vec3D& move) override;
    
    void read(std::istream& is) override;
    
    void write(std::ostream& os) const override;
    
    std::string getName() const override;
    
    void writeVTK(VTKContainer& vtk) const override;
 
//    void rotate(const Vec3D& rotate) override;
 
protected:
    std::vector<InfiniteWall> wallObjects_;
 
private:
    std::vector<Vec3D> A_;
    
    std::vector<Vec3D> AB_;
 
protected:
    std::vector<Vec3D> C_;
};
 
#endif