GeometryHelpers.h

Go to the documentation of this file.

#ifndef included_AMP_GeometryHelpers
#define included_AMP_GeometryHelpers
 
#include "AMP/utils/MeshPoint.h"
 
#include <array>
#include <vector>
 
 
namespace AMP::Mesh {
using Point = AMP::Mesh::MeshPoint<double>;
} // namespace AMP::Mesh
 
namespace AMP::Geometry::GeometryHelpers {
 
 
using Point1D = std::array<double, 1>;
using Point2D = std::array<double, 2>;
using Point3D = std::array<double, 3>;
 
 
Point2D map_logical_circle( double R, int method, double x, double y );
 
 
Point2D map_circle_logical( double R, int method, double x, double y );
 
 
Point2D map_logical_poly( int N, double R, double x, double y );
 
 
Point2D map_poly_logical( int N, double R, double x, double y );
 
 
std::vector<Point2D> get_poly_vertices( int N, double R );
 
 
Point3D map_logical_sphere( double R, double x, double y, double z );
 
 
Point3D map_sphere_logical( double R, double x, double y, double z );
 
 
Point3D map_logical_sphere_surface( int method, double R, double x, double y );
 
Point2D map_sphere_surface_logical( int method, double R, double x, double y, double z );
 
 
Point3D map_logical_shell( double r1, double r2, double x, double y, double z );
 
 
Point3D map_shell_logical( double r1, double r2, double x, double y, double z );
 
 
double distance( const Point3D &x1, const Point3D &x2 );
 
 
double
distanceToLine( const Point2D &pos, const Point2D &ang, const Point2D &v1, const Point2D &v2 );
 
 
double
distanceToLine( const Point3D &pos, const Point3D &ang, const Point3D &v1, const Point3D &v2 );
 
 
template<std::size_t NDIM>
double distanceToBox( const std::array<double, NDIM> &pos,
                      const std::array<double, NDIM> &ang,
                      const std::array<double, NDIM> &lb,
                      const std::array<double, NDIM> &ub );
 
 
double
distanceToPlane( const Point3D &n, const Point3D &p0, const Point3D &pos, const Point3D &ang );
 
 
double distanceToCircle( double r, const Point2D &pos, const Point2D &ang );
 
 
double distanceToCylinder( double r, double h, const Point3D &pos, const Point3D &ang );
 
 
double distanceToTube( double r1, double r2, double h, const Point3D &pos, const Point3D &ang );
 
 
double distanceToSphere( double r, const Point3D &pos, const Point3D &ang );
 
 
double distanceToCone( const Point3D &V, double theta, const Point3D &pos, const Point3D &ang );
 
 
double
distanceToTriangle( const std::array<Point2D, 3> &tri, const Point2D &pos, const Point2D &ang );
 
 
double
distanceToTriangle( const std::array<Point3D, 3> &tri, const Point3D &pos, const Point3D &ang );
 
 
double
distanceToTetrahedron( const std::array<Point3D, 4> &tet, const Point3D &pos, const Point3D &ang );
 
 
Point3D normalToQuadrilateral( const std::array<Point3D, 4> &quad );
 
double distanceToQuadrilateral( const std::array<Point2D, 4> &quad,
                                const Point2D &pos,
                                const Point2D &ang );
 
 
double distanceToQuadrilateral( const std::array<Point3D, 4> &quad,
                                const Point3D &pos,
                                const Point3D &ang );
 
 
double
distanceToCircularFrustum( double rb, double rt, double h, const Point3D &pos, const Point3D &ang );
 
 
template<int NP, int NDIM>
std::array<double, NP> barycentric( const std::array<std::array<double, NDIM>, NP> &x,
                                    const std::array<double, NDIM> &p );
 
 
Point2D normal( const Point2D &a, const Point2D &b );
 
 
Point3D normal( const Point3D &a, const Point3D &b, const Point3D &c );
 
 
Point2D nearest( const Point2D &v1, const Point2D &v2, const Point2D &p );
 
 
Point3D nearest( const Point3D &v1, const Point3D &v2, const Point3D &p );
 
 
Point3D nearest( const std::array<Point3D, 3> &v, const Point3D &p );
 
 
Point3D nearest( const std::array<Point3D, 4> &v, const Point3D &p );
 
 
std::vector<AMP::Mesh::Point> subdivide( const std::array<AMP::Mesh::Point, 3> &v, double res );
 
 
std::vector<Point3D> sampleLine( const std::array<Point3D, 2> &v, double d0, bool interior = true );
 
 
std::vector<Point3D> sampleTri( const std::array<Point3D, 3> &v, double d0, bool interior = true );
 
 
std::vector<Point3D> sampleQuad( const std::array<Point3D, 4> &v, double d0, bool interior = true );
 
 
std::vector<Point3D> sampleTet( const std::array<Point3D, 4> &v, double d0, bool interior = true );
 
 
std::vector<int> assignRanks( const std::vector<double> &x, int N_ranks );
 
std::vector<int> assignRanks( const std::vector<Point1D> &x, int N_ranks );
 
std::vector<int> assignRanks( const std::vector<Point2D> &x, int N_ranks );
 
std::vector<int> assignRanks( const std::vector<Point3D> &x, int N_ranks );
 
 
} // namespace AMP::Geometry::GeometryHelpers
 
#endif