#include <geomObjects.h>

Inheritance diagram for geom::Drum2D:

Collaboration diagram for geom::Drum2D:

Public Member Functions
	Drum2D ()
	Constructor.

	Drum2D (double r, double w, util::Point x=util::Point(0., 0., 0.), util::Point a=util::Point(1., 0., 0.), std::string description="")
	Constructor.

Drum2D &	operator= (const Drum2D &other)
	Assignment operator.

void	transform (const util::Point &translation, const double &scale, const double &angle, const util::Point &axis, const util::Point *rotationPoint) override
	Similarity about pivot \(\mathbf p\) (default: old center `d_x`), then rigid displacement \(\mathbf t\) = `translation:` \(\mathbf y = \mathbf p + s\,\mathbf R(\mathbf x-\mathbf p) + \mathbf t\). When \(\theta=0\) and \(s=1\), \(\mathbf y = \mathbf x + \mathbf t\) for every point and the center.

double	volume () const override
	Computes the volume (area in 2d, length in 1d) of object.

util::Point	center () const override
	Computes the center of object.

std::pair< util::Point, util::Point >	box () const override
	Computes the bounding box of object.

std::pair< util::Point, util::Point >	box (const double &tol) const override
	Computes the bounding box of object.

double	inscribedRadius () const override
	Computes the radius of biggest circle/sphere completely within the object.

double	boundingRadius () const override
	Computes the radius of smallest circle/sphere such that object can be fit into it.

std::string	printStr (int nt, int lvl) const override
	Returns the string containing printable information about the object.

void	print (int nt, int lvl) const override
	Prints the information about the object.

void	print () const override
	Prints the information about the object.

Interaction with point
bool	isInside (const util::Point &x) const override
	Checks if point is inside this object.

bool	isOutside (const util::Point &x) const override
	Checks if point is outside of this object.

bool	isNear (const util::Point &x, const double &tol) const override
	Checks if point is within given distance of this object.

bool	isNearBoundary (const util::Point &x, const double &tol, const bool &within) const override
	cons

bool	doesIntersect (const util::Point &x) const override
	Checks if point lies exactly on the boundary.

Interaction with box
bool	isInside (const std::pair< util::Point, util::Point > &box) const override
	Checks if box is completely inside.

bool	isOutside (const std::pair< util::Point, util::Point > &box) const override
	Checks if box is outside of the object.

bool	isNear (const std::pair< util::Point, util::Point > &box, const double &tol) const override
	Checks if box is within given distance of this object.

bool	doesIntersect (const std::pair< util::Point, util::Point > &box) const override
	Checks if box intersects this object.

Public Member Functions inherited from geom::GeomObject
	GeomObject (std::string name="", std::string description="")
	Constructor.

virtual bool	wallContactQuery (const util::Point &x, WallContactHit &hit) const
	Closest-point / signed-gap query for analytical wall contact. Default: unsupported (returns false, clears `hit`).

Data Fields
std::vector< util::Point >	d_vertices
	Vertices.

util::Point	d_x
	Center.

double	d_w
	Half width of neck.

double	d_r
	Distance between center and the farthest vertex.

util::Point	d_a
	Axis: defined as the vector pointing from center to the first vertex v3 v2.

Data Fields inherited from geom::GeomObject
const std::string	d_name
	name of object

const std::string	d_description
	Further description of object.

std::vector< std::string >	d_tags
	Tags/attributes about the object.

Detailed Description

Defines Drum2D.

        v3   o-------------------o   v2
             \                   /
              \                 /
               \               /
          v4    o      +      o v1
               /       c       \
              /                 \
             /                   \
       v5   o---------------------o  v6

w = distance between c and v1 = half-width of neck r = distance between c and v2 theta = pi/3 = angle between c-v2 and c-v1 a = axis = unit vector from c to v1

Definition at line 1757 of file geomObjects.h.

Constructor & Destructor Documentation

◆ Drum2D() [1/2]

geom::Drum2D::Drum2D ( )

inline

Constructor.

Definition at line 1791 of file geomObjects.h.

      : GeomObject("drum2d", ""),
        d_r(0.),
        d_w(0.),
        d_a(util::Point(1., 0., 0.)),
        d_x(util::Point()),
        d_vertices(std::vector<util::Point>(6, util::Point())) {
    };

◆ Drum2D() [2/2]

geom::Drum2D::Drum2D	(	double	r,
		double	w,
		util::Point	x = `util::Point(0., 0., 0.)`,
		util::Point	a = `util::Point(1., 0., 0.)`,
		std::string	description = `""`
	)

inline

Constructor.

Parameters

r	Distance between center and farthest vertex of drum2d structure
w	Half of the distance between two vertices in the neck (half width of neck)
x	Center point
a	Axis vector that will be rotated and scaled by r to get the farthest four vertices of structure
description	Description of object (e.g., further classification or any tag)

Definition at line 1810 of file geomObjects.h.

      : GeomObject("drum2d", description),
        d_r(r),
        d_w(w),
        d_a(a),
        d_x(x),
        d_vertices(std::vector<util::Point>(6, util::Point())) {
      // generate vertices in the plane perpendicular to rotate_axis (z for axis a in xy)
      const util::Point rotate_axis(0., 0., 1.);
 
      d_vertices[0] = d_x + d_w * d_a;
      d_vertices[3] = d_x - d_w * d_a;
 
      d_vertices[1] = d_x + d_r * util::rotate(d_a, M_PI / 3., rotate_axis);
      /* Option 2 used (v1 - k*z_hat), which is not in the drum plane for a in xy and
       * breaks coplanarity (e.g. z != 0 on some corners) -> degenerate nodal volumes. */
      d_vertices[2] = d_x + d_r * util::rotate(d_a, 2. * M_PI / 3., rotate_axis);
 
      d_vertices[4] = d_x + d_r * util::rotate(-1. * d_a, M_PI / 3., rotate_axis);
      d_vertices[5] = d_x + d_r * util::rotate(-1. * d_a, 2. * M_PI / 3., rotate_axis);
    };

References d_a, d_r, d_vertices, d_w, d_x, and util::rotate().

Here is the call graph for this function:

Member Function Documentation

◆ boundingRadius()

double geom::Drum2D::boundingRadius ( ) const

overridevirtual

Computes the radius of smallest circle/sphere such that object can be fit into it.

Returns: Radius Radius of bounding circle/sphere

Reimplemented from geom::GeomObject.

Definition at line 935 of file geomObjects.cpp.

                                        {
 
      return d_r;
    }

References d_r.

◆ box() [1/2]

std::pair< util::Point, util::Point > geom::Drum2D::box ( ) const

overridevirtual

Computes the bounding box of object.

Returns: Pair Left-bottom-back and right-top-front corner points of box

Reimplemented from geom::GeomObject.

Definition at line 917 of file geomObjects.cpp.

                                                    {
 
      return box(0.);
    }

References box().

Referenced by box(), doesIntersect(), isInside(), isNear(), isNear(), and isOutside().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ box() [2/2]

std::pair< util::Point, util::Point > geom::Drum2D::box ( const double & tol ) const

overridevirtual

Computes the bounding box of object.

Parameters

tol	Tolerance/padding used in creating bounding box

Returns: Pair Left-bottom-back and right-top-front corner points of box

Reimplemented from geom::GeomObject.

Definition at line 922 of file geomObjects.cpp.

                                                                                       {
 
      auto p1 = d_x - util::Point(d_r + tol, d_r + tol, d_x[2] + tol);
      auto p2 = d_x + util::Point(d_r + tol, d_r + tol, d_x[2] + tol);
      return {p1, p2};
    }

References d_r, and d_x.

◆ center()

util::Point geom::Drum2D::center ( ) const

overridevirtual

Computes the center of object.

Returns: Point Coordinates of center

Reimplemented from geom::GeomObject.

Definition at line 913 of file geomObjects.cpp.

                                   {
      return d_x;
    }

References d_x.

◆ doesIntersect() [1/2]

bool geom::Drum2D::doesIntersect ( const std::pair< util::Point, util::Point > & box ) const

overridevirtual

Checks if box intersects this object.

Parameters

box Box

Returns: True True if intersects

Reimplemented from geom::GeomObject.

Definition at line 1025 of file geomObjects.cpp.

                                                            {
 
      // need to check all four corner points
      for (auto p: geom::getCornerPoints(2, box))
        if (this->isInside(p))
          return true;
 
      return false;
    }

References box(), geom::getCornerPoints(), and isInside().

Here is the call graph for this function:

◆ doesIntersect() [2/2]

bool geom::Drum2D::doesIntersect ( const util::Point & x ) const

overridevirtual

Checks if point lies exactly on the boundary.

Parameters

x Point

Returns: True True if it lies on the boundary

Reimplemented from geom::GeomObject.

Definition at line 992 of file geomObjects.cpp.

                                                       {
 
      return isNearBoundary(x, 1.0E-8, false);
    }

References isNearBoundary().

Here is the call graph for this function:

◆ inscribedRadius()

double geom::Drum2D::inscribedRadius ( ) const

overridevirtual

Computes the radius of biggest circle/sphere completely within the object.

Returns: Radius Radius of inscribed circle/sphere

Reimplemented from geom::GeomObject.

Definition at line 930 of file geomObjects.cpp.

                                         {
 
      return d_w;
    }

References d_w.

Referenced by isInside(), and isNearBoundary().

Here is the caller graph for this function:

◆ isInside() [1/2]

bool geom::Drum2D::isInside ( const std::pair< util::Point, util::Point > & box ) const

overridevirtual

Checks if box is completely inside.

Parameters

box Box

Returns: True True if box lies inside

Reimplemented from geom::GeomObject.

Definition at line 997 of file geomObjects.cpp.

                                                            {
 
      for (auto p: geom::getCornerPoints(2, box))
        if (!this->isInside(p))
          return false;
 
      return true;
    }

References box(), geom::getCornerPoints(), and isInside().

Here is the call graph for this function:

◆ isInside() [2/2]

bool geom::Drum2D::isInside ( const util::Point & x ) const

overridevirtual

Checks if point is inside this object.

Parameters

x Point

Returns: True If point lies inside

Reimplemented from geom::GeomObject.

Definition at line 940 of file geomObjects.cpp.

                                                  {
 
      if ((x - d_x).length() > d_r)
        return false;
 
      if ((x - d_x).length() < inscribedRadius())
        return true;
 
      // rotate axis to get orthogonal axis
      auto ortho_axis = util::rotate(d_a, M_PI * 0.5, util::Point(0., 0., 1.));
 
      //
      //                                   + v2
      //                                  /
      //                                /    x
      //                              /
      //                            /
      //                     o----+v1
      //
      auto ox = x - d_x;
      double angle_ox_ov1 = std::acos(std::abs(d_a.dot(ox)) / ox.length());
      double max_length = d_w + angle_ox_ov1 * (d_r - d_w) / (M_PI / 3.);
 
      return ox.length() <= max_length;
    }

References d_a, d_r, d_w, d_x, util::Point::dot(), inscribedRadius(), and util::rotate().

Referenced by doesIntersect(), isInside(), isOutside(), and isOutside().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ isNear() [1/2]

bool geom::Drum2D::isNear	(	const std::pair< util::Point, util::Point > &	box,
		const double &	tol
	)		const

overridevirtual

Checks if box is within given distance of this object.

Parameters

box	Box
tol	Tolerance used in checking the nearness

Returns: True True if box is inside within the tol distance

Reimplemented from geom::GeomObject.

Definition at line 1018 of file geomObjects.cpp.

                                     {
 
      return geom::areBoxesNear(this->box(), box, tol, 2);
    }

References geom::areBoxesNear(), and box().

Here is the call graph for this function:

◆ isNear() [2/2]

bool geom::Drum2D::isNear	(	const util::Point &	x,
		const double &	tol
	)		const

overridevirtual

Checks if point is within given distance of this object.

Parameters

x	Point
tol	Tolerance used in checking the nearness

Returns: True True if within the tol distance

Reimplemented from geom::GeomObject.

Definition at line 970 of file geomObjects.cpp.

                                                                 {
 
      // get a bigger box containing this object
      auto bbox = box(tol);
 
      return geom::isPointInsideBox(x, 2, bbox);
    }

References box(), and geom::isPointInsideBox().

Here is the call graph for this function:

◆ isNearBoundary()

bool geom::Drum2D::isNearBoundary	(	const util::Point &	x,
		const double &	tol,
		const bool &	within
	)		const

overridevirtual

cons

Reimplemented from geom::GeomObject.

Definition at line 979 of file geomObjects.cpp.

                                                               {
 
      if ((x - d_x).length() > d_r + tol)
        return false;
 
      if ((x - d_x).length() < this->inscribedRadius() - tol)
        return false;
 
      return true;
    }

References d_r, d_x, and inscribedRadius().

Referenced by doesIntersect().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ isOutside() [1/2]

bool geom::Drum2D::isOutside ( const std::pair< util::Point, util::Point > & box ) const

overridevirtual

Checks if box is outside of the object.

Parameters

box Box

Returns: True True if box lies outside

Reimplemented from geom::GeomObject.

Definition at line 1007 of file geomObjects.cpp.

                                                            {
 
      bool intersect = false;
      for (auto p: geom::getCornerPoints(2, box))
        if (!intersect)
          intersect = this->isInside(p);
 
      return !intersect;
    }

References box(), geom::getCornerPoints(), and isInside().

Here is the call graph for this function:

◆ isOutside() [2/2]

bool geom::Drum2D::isOutside ( const util::Point & x ) const

overridevirtual

Checks if point is outside of this object.

Parameters

x Point

Returns: True If point lies outside

Reimplemented from geom::GeomObject.

Definition at line 966 of file geomObjects.cpp.

                                                   {
      return !isInside(x);
    }

References isInside().

Here is the call graph for this function:

◆ operator=()

Drum2D & geom::Drum2D::operator= ( const Drum2D & other )

inline

Assignment operator.

Parameters

other Object to copy from

Returns: Reference to this object

Definition at line 1839 of file geomObjects.h.

                                           {
      if (this != &other) {
        // Copy base class members
        d_tags = other.d_tags;
        
        // Copy member variables
        d_r = other.d_r;
        d_w = other.d_w;
        d_a = other.d_a;
        d_x = other.d_x;
        d_vertices = other.d_vertices;
      }
      return *this;
    }

References d_a, d_r, geom::GeomObject::d_tags, d_vertices, d_w, and d_x.

◆ print() [1/2]

void geom::Drum2D::print ( ) const

inlineoverridevirtual

Prints the information about the object.

Reimplemented from geom::GeomObject.

Definition at line 1985 of file geomObjects.h.

1985{ print(0, 0); };

geom::Drum2D::print

void print() const override

Prints the information about the object.

Definition geomObjects.h:1985

References print().

Referenced by print().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ print() [2/2]

void geom::Drum2D::print	(	int	nt,
		int	lvl
	)		const

inlineoverridevirtual

Prints the information about the object.

Parameters

nt	Number of tabs to append before printing
lvl	Information level (higher means more information)

Reimplemented from geom::GeomObject.

Definition at line 1978 of file geomObjects.h.

                                               {
      std::cout << printStr(nt, lvl);
    };

References printStr().

Here is the call graph for this function:

◆ printStr()

std::string geom::Drum2D::printStr	(	int	nt,
		int	lvl
	)		const

overridevirtual

Returns the string containing printable information about the object.

Parameters

nt	Number of tabs to append before printing
lvl	Information level (higher means more information)

Returns: string String containing printable information about the object

Reimplemented from geom::GeomObject.

Definition at line 1036 of file geomObjects.cpp.

                                                    {
 
      auto tabS = util::io::getTabS(nt);
 
      std::ostringstream oss;
 
      oss << tabS << "------- Drum2D --------" << std::endl << std::endl;
      oss << tabS << "Name = " << d_name << std::endl;
      oss << tabS << "Radius = " << d_r << std::endl;
      oss << tabS << "Neck half-width = " << d_w << std::endl;
      oss << tabS << "Center = " << d_x.printStr(0, lvl) << std::endl;
      oss << tabS << "Axis = " << d_a.printStr(0, lvl) << std::endl;
      oss << tabS << "Vertices = " << util::io::printStr(d_vertices, lvl) <<
          std::endl;
      oss << std::endl;
 
      if (lvl == 0)
        oss << std::endl;
 
      return oss.str();
    }

References d_a, geom::GeomObject::d_name, d_r, d_vertices, d_w, d_x, util::io::getTabS(), util::io::printStr(), and util::Point::printStr().

Referenced by print().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ transform()

void geom::Drum2D::transform	(	const util::Point &	translation,
		const double &	scale,
		const double &	angle,
		const util::Point &	axis,
		const util::Point *	rotationPoint
	)

inlineoverridevirtual

Similarity about pivot \(\mathbf p\) (default: old center d_x), then rigid displacement \(\mathbf t\) = translation: \(\mathbf y = \mathbf p + s\,\mathbf R(\mathbf x-\mathbf p) + \mathbf t\). When \(\theta=0\) and \(s=1\), \(\mathbf y = \mathbf x + \mathbf t\) for every point and the center.

Parameters

translation	Rigid translation vector \(\mathbf t\) (added after rotate+scale about \(\mathbf p\))
scale	Uniform scale factor
angle	Rotation angle (radians)
axis	Axis of rotation
rotationPoint	If non-null, use as \(\mathbf p\); if null, use old `d_x` as pivot.

Reimplemented from geom::GeomObject.

Definition at line 1854 of file geomObjects.h.

                                                                                                       {
      // Scale the dimensions
      d_r *= scale;
      d_w *= scale;
 
      const util::Point c0 = d_x;
      for (auto &v : d_vertices) {
        v = mapSimilarity(v, c0, translation, scale, angle, axis, rotationPoint);
      }
      d_x = mapSimilarity(c0, c0, translation, scale, angle, axis, rotationPoint);
      d_a = util::rotate(d_a, angle, axis);
    }  

References d_a, d_r, d_vertices, d_w, d_x, geom::mapSimilarity(), and util::rotate().

Here is the call graph for this function:

◆ volume()

double geom::Drum2D::volume ( ) const

overridevirtual

Computes the volume (area in 2d, length in 1d) of object.

Returns: Volume Volume of object

Reimplemented from geom::GeomObject.

Definition at line 908 of file geomObjects.cpp.

                                {
 
      return (2. * d_r * d_r - d_r * (d_r - 2. * d_w)) * std::sin(M_PI / 3.);
    }

References d_r, and d_w.

Field Documentation

◆ d_a

util::Point geom::Drum2D::d_a

Axis: defined as the vector pointing from center to the first vertex v3 v2.

+

                 +         o           +
                v4         x            v1

          +                                +
          v5                               v6

Axis is a unit vector from x to v1

Definition at line 1785 of file geomObjects.h.

Referenced by Drum2D(), isInside(), operator=(), printStr(), and transform().

◆ d_r

double geom::Drum2D::d_r

Distance between center and the farthest vertex.

Definition at line 1769 of file geomObjects.h.

Referenced by boundingRadius(), box(), Drum2D(), isInside(), isNearBoundary(), operator=(), printStr(), transform(), and volume().

◆ d_vertices

std::vector<util::Point> geom::Drum2D::d_vertices

Vertices.

Definition at line 1760 of file geomObjects.h.

Referenced by Drum2D(), mesh_gen::anonymous_namespace{builtinGmshGeometry.cpp}::drum2dPolygonBoundary(), operator=(), printStr(), and transform().

◆ d_w

double geom::Drum2D::d_w

Half width of neck.

Definition at line 1766 of file geomObjects.h.

Referenced by Drum2D(), inscribedRadius(), isInside(), operator=(), printStr(), transform(), and volume().

◆ d_x

util::Point geom::Drum2D::d_x

Center.

Definition at line 1763 of file geomObjects.h.

Referenced by box(), center(), Drum2D(), isInside(), isNearBoundary(), operator=(), printStr(), and transform().

The documentation for this class was generated from the following files:

/home/user/project/src/geom/geomObjects.h
/home/user/project/src/geom/geomObjects.cpp

Public Member Functions

Data Fields

Detailed Description

Constructor & Destructor Documentation

◆ Drum2D() [1/2]

◆ Drum2D() [2/2]

Member Function Documentation

◆ boundingRadius()

◆ box() [1/2]

◆ box() [2/2]

◆ center()

◆ doesIntersect() [1/2]

◆ doesIntersect() [2/2]

◆ inscribedRadius()

◆ isInside() [1/2]

◆ isInside() [2/2]

◆ isNear() [1/2]

◆ isNear() [2/2]

◆ isNearBoundary()

◆ isOutside() [1/2]

◆ isOutside() [2/2]

◆ operator=()

◆ print() [1/2]

◆ print() [2/2]

◆ printStr()

◆ transform()

◆ volume()

Field Documentation

◆ d_a

◆ d_r

◆ d_vertices

◆ d_w

◆ d_x