MeshUtilities.h

Go to the documentation of this file.

#ifndef included_AMP_MeshUtilities
#define included_AMP_MeshUtilities
 
 
#include "AMP/geometry/Geometry.h"
#include "AMP/mesh/Mesh.h"
#include "AMP/mesh/MeshElement.h"
#include "AMP/utils/Array.h"
#include "AMP/utils/MeshPoint.h"
#include "AMP/utils/kdtree2.h"
 
 
namespace AMP::Mesh {
 
 
std::tuple<std::vector<Point>, std::vector<MeshElementID>> sample( const Mesh &mesh, double dx );
 
 
std::vector<Point> sample( const MeshElement &elem, double dx );
 
 
Array<double> volumeOverlap( const AMP::Geometry::Geometry &geom, const std::vector<int> &N );
 
 
class ElementFinder
{
public:
    ElementFinder() : d_pos_hash( 0 ), d_dist( 0 ) {}
 
    ElementFinder( std::shared_ptr<AMP::Mesh::Mesh> mesh );
 
    ElementFinder( std::shared_ptr<AMP::Mesh::Mesh> mesh, AMP::Mesh::MeshIterator it );
 
    ElementFinder( const ElementFinder & ) = delete;
 
    ElementFinder( ElementFinder && ) = default;
 
    ElementFinder &operator=( ElementFinder && ) = default;
 
    std::pair<std::unique_ptr<MeshElement>, Point> nearest( const Point &x ) const;
 
    double distance( const Point &pos, const Point &dir ) const;
 
    inline const auto &getTree() const { return d_tree; }
 
    void update() const;
 
private:
    void initialize() const;
 
private:
    std::shared_ptr<AMP::Mesh::Mesh> d_mesh;
    mutable uint64_t d_pos_hash;
    AMP::Mesh::MeshIterator d_elements;
    mutable double d_dist;
    mutable kdtree2<3, AMP::Mesh::MeshElementID> d_tree;
};
 
 
} // namespace AMP::Mesh
 
 
#endif