TriangleMeshWall.h

Go to the documentation of this file.

//Copyright (c) 2013-2023, The MercuryDPM Developers Team. All rights reserved.
//For the list of developers, see <http://www.MercuryDPM.org/Team>.
//
//Redistribution and use in source and binary forms, with or without
//modification, are permitted provided that the following conditions are met:
//  * Redistributions of source code must retain the above copyright
//    notice, this list of conditions and the following disclaimer.
//  * Redistributions in binary form must reproduce the above copyright
//    notice, this list of conditions and the following disclaimer in the
//    documentation and/or other materials provided with the distribution.
//  * Neither the name MercuryDPM nor the
//    names of its contributors may be used to endorse or promote products
//    derived from this software without specific prior written permission.
//
//THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
//ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
//WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
//DISCLAIMED. IN NO EVENT SHALL THE MERCURYDPM DEVELOPERS TEAM BE LIABLE FOR ANY
//DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
//(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
//LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
//ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
//(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
//SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 
#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