#include <geomObjects.h>

Inheritance diagram for geom::Line:

Collaboration diagram for geom::Line:

Public Member Functions
	Line ()
	Constructor.

	Line (util::Point x1, util::Point x2, std::string description="")
	Constructor.

	Line (double L, util::Point x=util::Point(), std::string description="")
	Constructor.

	Line (const Line &other)
	Copy constructor.

Line &	operator= (const Line &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.

bool	wallContactQuery (const util::Point &x, WallContactHit &hit) const override

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
TODO Implement methods
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.

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

util::Point	d_x
	Center.

double	d_L
	Length of line.

double	d_r
	Radius of bounding circle.

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 Line.

Definition at line 522 of file geomObjects.h.

Constructor & Destructor Documentation

◆ Line() [1/4]

geom::Line::Line ( )

inline

Constructor.

Definition at line 540 of file geomObjects.h.

      : GeomObject("line", ""),
        d_r(0.),
        d_L(0.),
        d_x(util::Point()),
        d_vertices({util::Point(), util::Point()}) {
    };

◆ Line() [2/4]

geom::Line::Line	(	util::Point	x1,
		util::Point	x2,
		std::string	description = `""`
	)

inline

Constructor.

Parameters

x1	Left-bottom corner point
x2	Right-top corner point
description	Description of object (e.g., further classification or any tag)

Definition at line 555 of file geomObjects.h.

      : GeomObject("line", description),
        d_L((x1 - x2).length()),
        d_r(0.5 * (x1 - x2).length()),
        d_x(0.5 * (x1 + x2)),
        d_vertices({x1, x2}) {
    };

◆ Line() [3/4]

geom::Line::Line	(	double	L,
		util::Point	x = `util::Point()`,
		std::string	description = `""`
	)

inline

Constructor.

Parameters

L	Length of a line
x	Center point
description	Description of object (e.g., further classification or any tag)

Definition at line 570 of file geomObjects.h.

      : GeomObject("line", description),
        d_L(L),
        d_r(0.5 * L),
        d_x(x),
        d_vertices({
          x + util::Point(-0.5 * L, 0., 0.),
          x + util::Point(0.5 * L, 0., 0.)
        }) {
    };

◆ Line() [4/4]

geom::Line::Line ( const Line & other )

inline

Copy constructor.

Parameters

other Object to copy from

Definition at line 586 of file geomObjects.h.

      : GeomObject(other.d_name, other.d_description),
        d_r(other.d_r),
        d_L(other.d_L),
        d_x(other.d_x),
        d_vertices(other.d_vertices) {
      d_tags = other.d_tags;
    }

References geom::GeomObject::d_tags.

Member Function Documentation

◆ boundingRadius()

double geom::Line::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 136 of file geomObjects.cpp.

                                      {
 
      return d_r;
    }

References d_r.

◆ box() [1/2]

std::pair< util::Point, util::Point > geom::Line::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 120 of file geomObjects.cpp.

                                                  {
 
      return box(0.);
    }

References box().

Referenced by box(), and printStr().

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::Line::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 125 of file geomObjects.cpp.

                                                                                     {
 
      return {d_vertices[0] - tol, d_vertices[1] + tol};
    }

References d_vertices.

◆ center()

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

overridevirtual

Computes the center of object.

Returns: Point Coordinates of center

Reimplemented from geom::GeomObject.

Definition at line 116 of file geomObjects.cpp.

                                 {
      return d_x;
    }

References d_x.

◆ doesIntersect() [1/2]

bool geom::Line::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 221 of file geomObjects.cpp.

                                                            {
 
      return false;
    }

◆ doesIntersect() [2/2]

bool geom::Line::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 197 of file geomObjects.cpp.

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

References isNearBoundary().

Here is the call graph for this function:

◆ inscribedRadius()

double geom::Line::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 131 of file geomObjects.cpp.

                                       {
 
      return d_r;
    }

References d_r.

◆ isInside() [1/2]

bool geom::Line::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 202 of file geomObjects.cpp.

                                                            {
 
      return false;
    }

◆ isInside() [2/2]

bool geom::Line::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 141 of file geomObjects.cpp.

                                                {
 
      auto da = (d_vertices[1] - d_vertices[0]) / d_L;
      auto db = x - d_vertices[0];
      double dot = db * da;
 
      if (util::isLess(dot, 0.) or util::isGreater(dot, d_L))
        return false;
 
 
      auto dx = db - dot * da;
 
      return util::isLess(dx.length(), 1.0E-10);
    }

References d_L, d_vertices, util::isGreater(), and util::isLess().

Referenced by isOutside().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ isNear() [1/2]

bool geom::Line::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 214 of file geomObjects.cpp.

                                     {
 
      return true;
    }

◆ isNear() [2/2]

bool geom::Line::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 161 of file geomObjects.cpp.

                                                               {
 
      auto da = (d_vertices[1] - d_vertices[0]) / d_L;
      auto db = x - d_vertices[0];
      double dot = db * da;
 
      if (util::isLess(dot, 0.) or util::isGreater(dot, d_L))
        return false;
 
      auto dx = db - dot * da;
 
      return util::isLess(dx.length(), tol);
    }

References d_L, d_vertices, util::isGreater(), and util::isLess().

Referenced by isNearBoundary().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ isNearBoundary()

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

overridevirtual

cons

Reimplemented from geom::GeomObject.

Definition at line 176 of file geomObjects.cpp.

                                                             {
 
      // check if particle is inside the object
      if (!isNear(x, within ? 0. : tol))
        return false;
 
      auto da = (d_vertices[1] - d_vertices[0]) / d_L;
      auto db = x - d_vertices[0];
      double dot = db * da;
 
      if (util::isLess(dot, 0.) or util::isGreater(dot, tol)
          or util::isGreater(dot, d_L) or util::isLess(dot, d_L - tol))
        return false;
 
      auto dx = db - dot * da;
 
      return util::isLess(dx.length(), tol);
    }

References d_L, d_vertices, util::isGreater(), util::isLess(), and isNear().

Referenced by doesIntersect().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ isOutside() [1/2]

bool geom::Line::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 208 of file geomObjects.cpp.

                                                            {
 
      return true;
    }

◆ isOutside() [2/2]

bool geom::Line::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 156 of file geomObjects.cpp.

                                                 {
 
      return !isInside(x);
    }

References isInside().

Here is the call graph for this function:

◆ operator=()

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

inline

Assignment operator.

Parameters

other Object to copy from

Returns: Reference to this object

Definition at line 600 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_L = other.d_L;
        d_x = other.d_x;
        d_vertices = other.d_vertices;
      }
      return *this;
    }

References d_L, d_r, geom::GeomObject::d_tags, d_vertices, and d_x.

◆ print() [1/2]

void geom::Line::print ( ) const

inlineoverridevirtual

Prints the information about the object.

Reimplemented from geom::GeomObject.

Definition at line 753 of file geomObjects.h.

753{ print(0, 0); };

geom::Line::print

void print() const override

Prints the information about the object.

Definition geomObjects.h:753

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::Line::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 746 of file geomObjects.h.

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

References printStr().

Here is the call graph for this function:

◆ printStr()

std::string geom::Line::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 227 of file geomObjects.cpp.

                                                  {
 
      auto tabS = util::io::getTabS(nt);
 
      std::ostringstream oss;
 
      oss << tabS << "------- Line --------" << std::endl << std::endl;
      oss << tabS << "Name = " << d_name << std::endl;
      oss << tabS << "Length = " << d_L << std::endl;
      oss << tabS << "Point 1 = " << d_vertices[0].printStr(0, lvl) << std::endl;
      oss << tabS << "Point 2 = " << d_vertices[1].printStr(0, lvl) << std::endl;
      oss << std::endl;
 
      if (lvl > 0)
        oss << tabS << "Bounding box: "
            << util::io::printBoxStr(box(0.), nt + 1);
 
      if (lvl == 0)
        oss << std::endl;
 
      return oss.str();
    }

References box(), d_L, geom::GeomObject::d_name, d_vertices, util::io::getTabS(), and util::io::printBoxStr().

Referenced by print().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ transform()

void geom::Line::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 614 of file geomObjects.h.

                                                                                                       {
      // Scale the length
      d_L *= scale;
      d_r *= scale;
 
      const util::Point c0 = d_x;
      // Rotate and translate vertices
      for (auto &v : d_vertices) {
        v = mapSimilarity(v, c0, translation, scale, angle, axis, rotationPoint);
      }
      d_x = mapSimilarity(c0, c0, translation, scale, angle, axis, rotationPoint);
    }

References d_L, d_r, d_vertices, d_x, and geom::mapSimilarity().

Here is the call graph for this function:

◆ volume()

double geom::Line::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 112 of file geomObjects.cpp.

                              {
      return d_L;
    }

References d_L.

◆ wallContactQuery()

bool geom::Line::wallContactQuery	(	const util::Point &	x,
		WallContactHit &	hit
	)		const

overridevirtual

Finite segment wall in 2D: free space is to the left of directed vertex0 → vertex1 (outward normal = rotate tangent 90° CCW).

Reimplemented from geom::GeomObject.

Definition at line 67 of file geomObjects.cpp.

                                                           {
      hit = WallContactHit();
      if (d_vertices.size() < 2)
        return false;
      const util::Point &a = d_vertices[0];
      const util::Point &b = d_vertices[1];
      const util::Point t = b - a;
      const double t2 = t.lengthSq();
      if (!(t2 > 1.e-24))
        return false;
 
      // Outward = left of directed segment in xy (CCW rotate).
      util::Point n(-t.d_y, t.d_x, 0.);
      const double nl = n.length();
      if (!(nl > 1.e-16))
        return false;
      n = n / nl;
 
      double s = ((x - a) * t) / t2;
      if (s < 0.)
        s = 0.;
      else if (s > 1.)
        s = 1.;
      hit.closest = a + s * t;
      const util::Point dx = x - hit.closest;
      if (s > 1.e-12 && s < 1. - 1.e-12) {
        hit.signed_gap = dx * n;
        hit.outward_n = n;
      } else {
        const double dist = dx.length();
        if (!(dist > 1.e-16)) {
          hit.signed_gap = 0.;
          hit.outward_n = n;
        } else {
          hit.outward_n = dx / dist;
          hit.signed_gap = (hit.outward_n * n >= 0.) ? dist : -dist;
          if (hit.outward_n * n < 0.)
            hit.outward_n = n;
        }
      }
      hit.active = true;
      return true;
    }

References geom::WallContactHit::active, geom::WallContactHit::closest, d_vertices, util::Point::d_x, util::Point::d_y, util::Point::length(), util::Point::lengthSq(), geom::WallContactHit::outward_n, and geom::WallContactHit::signed_gap.

Here is the call graph for this function:

Field Documentation

◆ d_L

double geom::Line::d_L

Length of line.

Definition at line 531 of file geomObjects.h.

Referenced by isInside(), isNear(), isNearBoundary(), operator=(), printStr(), transform(), and volume().

◆ d_r

double geom::Line::d_r

Radius of bounding circle.

Definition at line 534 of file geomObjects.h.

Referenced by boundingRadius(), inscribedRadius(), operator=(), and transform().

◆ d_vertices

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

Vertices.

Definition at line 525 of file geomObjects.h.

Referenced by box(), isInside(), isNear(), isNearBoundary(), operator=(), printStr(), transform(), and wallContactQuery().

◆ d_x

util::Point geom::Line::d_x

Center.

Definition at line 528 of file geomObjects.h.

Referenced by center(), operator=(), 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

◆ Line() [1/4]

◆ Line() [2/4]

◆ Line() [3/4]

◆ Line() [4/4]

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()

◆ wallContactQuery()

Field Documentation

◆ d_L

◆ d_r

◆ d_vertices

◆ d_x