MeshGeometry.h

Go to the documentation of this file.

#ifndef included_AMP_MeshGeometry
#define included_AMP_MeshGeometry
 
#include "AMP/geometry/Geometry.h"
#include "AMP/mesh/Mesh.h"
#include "AMP/mesh/MeshUtilities.h"
 
#include <memory>
#include <vector>
 
 
namespace AMP::Geometry {
 
 
class MeshGeometry final : public Geometry
{
public:
    MeshGeometry( std::shared_ptr<AMP::Mesh::Mesh> mesh );
 
    virtual ~MeshGeometry() = default;
 
    MeshGeometry( const MeshGeometry & ) = delete;
 
    MeshGeometry &operator=( const MeshGeometry & ) = delete;
 
    std::string getName() const override { return "MeshGeometry"; }
 
    bool isConvex() const override;
 
    Point nearest( const Point &pos ) const override;
 
    double distance( const Point &pos, const Point &dir ) const override;
 
    bool inside( const Point &pos ) const override;
 
    int NSurface() const override;
 
    int surface( const Point &x ) const override;
 
    Point surfaceNorm( const Point &x ) const override;
 
    Point centroid() const override;
 
    std::pair<Point, Point> box() const override;
 
    double volume() const override;
 
    void displace( const double *x ) override;
 
    std::unique_ptr<AMP::Geometry::Geometry> clone() const override;
 
    const AMP::Mesh::Mesh &getMesh() const { return *d_mesh; }
 
    bool operator==( const Geometry &rhs ) const override;
 
    MeshGeometry( int64_t );
 
protected: // Write/read restart data
    void writeRestart( int64_t ) const override;
 
private:
    void updateCache() const; // Update cached data if underlying mesh has moved
 
private:                                     // Internal data
    std::shared_ptr<AMP::Mesh::Mesh> d_mesh; // Underlying mesh
    std::vector<int> d_surfaceIds;           // Surface ids
    mutable uint64_t d_pos_hash;             // Position hash to update cached data
    mutable bool d_isConvex;                 // Check if the mesh is convex
    mutable double d_volume;                 // Cached value for the volume
    mutable Point d_centroid;                // Cached value for the centroid
    AMP::Mesh::ElementFinder d_find;         // Nearest element finder
    mutable kdtree2<3, bool> d_inside;       // Lookup to find points inside/outside
};
 
 
} // namespace AMP::Geometry
 
 
#endif