TriangleMeshWall.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 MERCURYDPM_TRIANGLEMESHWALL_H
#define MERCURYDPM_TRIANGLEMESHWALL_H
 
#include "BaseWall.h"
#include "TriangleWall.h"
 
class TriangleMeshWall : public BaseWall
{
public:
    
    TriangleMeshWall() = default;
    
    TriangleMeshWall(const TriangleMeshWall& other);
    
    TriangleMeshWall(const std::vector<Vec3D>& points, const std::vector<std::array<unsigned, 3>>& cells, const ParticleSpecies* species = nullptr);
 
    TriangleMeshWall(const Vec3D& P0, const Vec3D& P1, const Vec3D& P, Mdouble resolutionU, Mdouble resolutionV,
                     const ParticleSpecies* species = nullptr, bool periodicInU = false, bool periodicInV = false);
    
    ~TriangleMeshWall() override = default;
    
    TriangleMeshWall& operator=(const TriangleMeshWall& other);
    
    TriangleMeshWall* copy() const override;
    
    void read(std::istream& is) override;
    
    void write(std::ostream& os) const override;
    
    std::string getName() const override;
    
    void set(const std::vector<Vec3D>& points, const std::vector<std::array<unsigned, 3>>& cells);
    
    bool getDistanceAndNormal(const BaseParticle& p, Mdouble& distance, Vec3D& normal_return) const override;
    
    BaseInteraction* getInteractionWith(BaseParticle* p, unsigned timeStamp, InteractionHandler* interactionHandler) override;
    
    void resetForceTorque(int numberOfOMPthreads) override;
    
    void writeVTK(VTKContainer& vtk) const override;
    
    void setSpecies(const ParticleSpecies* species);
    
    void moveVertex(unsigned index, const Vec3D& dP);
    
    void removeTriangle(unsigned index, bool removeFreeVertex = true);
    
    void refineTriangle(unsigned index, unsigned numberOfTimes = 1);
    
    void setPeriodicCompanions(const std::vector<std::pair<unsigned, unsigned>>& periodicCompanions);
    
    void moveVerticesToMatchVolume(std::vector<Vec3D> displacements, Mdouble targetVolume, int maxNumRecursiveCalls = 15);
    
    Mdouble getVolumeTetrahedron(const Vec3D& a, const Vec3D& b, const Vec3D& c, const Vec3D& d);
    
    void addToMesh(TriangleMeshWall mesh);
 
    void createParallelogramMesh(const Vec3D& P0, const Vec3D& P1, const Vec3D& P2, Mdouble resolutionU, Mdouble resolutionV, bool periodicInU = false, bool periodicInV = false);
 
    void createFourPointMesh(const Vec3D& P0, const Vec3D& P1, const Vec3D& P2, const Vec3D& P3, int numSegmentsU, int numSegmentsV);
    
    bool isLocal(Vec3D& min, Vec3D& max) const override;
    
    void setPosition(const Vec3D& position) override;
    void setOrientation(const Quaternion& orientation) override;
    void move(const Vec3D& move) override;
    void rotate(const Vec3D& angularVelocity) override;
   
protected:
    struct Triangle
    {
        Triangle() = default;
        
        Triangle(const TriangleWall& w, std::array<unsigned, 3> vi) : wall(w), vertexIndices(vi) {};
        
        TriangleWall wall;
        
        std::array<unsigned, 3> vertexIndices{};
    };
    
    struct Vertex
    {
        Vertex() = default;
        
        explicit Vertex(Vec3D P) : position(P) {};
        
        Vec3D position;
        
        std::vector<std::pair<unsigned, unsigned>> triangleIndices;
    };
 
protected:
    std::vector<Triangle> triangles_;
    std::vector<Vertex> vertices_;
    
private:
    std::vector<std::pair<unsigned, unsigned>> periodicCompanions_;
    
    std::vector<unsigned> getPeriodicCompanions(unsigned index);
    
    Vec3D boundingBoxMin_, boundingBoxMax_;
    Vec3D boundingBoxMinGlobal_, boundingBoxMaxGlobal_;
    
    bool isWithinBoundingBox(const Vec3D& position, Mdouble radius) const;
    
    void updateBoundingBox();
    
    void updateBoundingBoxGlobal();
 
    int getNumberOfSegmentsAndResolution(Mdouble length, Mdouble& resolution);
 
    void addZigZagDiagonalCells(std::vector<std::array<unsigned, 3>>& cells, int numU, int numV, int i, int j);
};
 
 
#endif //MERCURYDPM_TRIANGLEMESHWALL_H