This class represents a 3 dimensional sphere shape. More...

#include <Sphere.h>

Inheritance diagram for inet::Sphere:

Public Member Functions
	Sphere (double radius)

const double	getRadius () const

void	setRadius (double radius)

virtual Coord	computeBoundingBoxSize () const override
	Computes the 3 dimensional size of the shapes's bounding box. More...

virtual bool	computeIntersection (const LineSegment &lineSegment, Coord &intersection1, Coord &intersection2, Coord &normal1, Coord &normal2) const override
	Computes the intersection with the given line segment in the shape's coordinate system. More...

Public Member Functions inherited from inet::ShapeBase
	ShapeBase ()

virtual	~ShapeBase ()

Protected Attributes
double	radius

Detailed Description

This class represents a 3 dimensional sphere shape.

Constructor & Destructor Documentation

inet::Sphere::Sphere ( double radius )

                             :
     radius(radius)
 {
 }

Member Function Documentation

virtual Coord inet::Sphere::computeBoundingBoxSize ( ) const

inlineoverridevirtual

Computes the 3 dimensional size of the shapes's bounding box.

Implements inet::ShapeBase.

39 { return Coord(radius * 2, radius * 2, radius * 2); }

inet::Sphere::radius

double radius

Definition: Sphere.h:31

bool inet::Sphere::computeIntersection	(	const LineSegment &	lineSegment,
		Coord &	intersection1,
		Coord &	intersection2,
		Coord &	normal1,
		Coord &	normal2
	)		const

overridevirtual

Computes the intersection with the given line segment in the shape's coordinate system.

Implements inet::ShapeBase.

 {
     // NOTE: based on http://paulbourke.net/geometry/circlesphere/
     Coord point1 = lineSegment.getPoint1();
     Coord point2 = lineSegment.getPoint2();
     double a, b, c;
     double bb4ac;
     Coord dp = point2 - point1;
     a = dp * dp;
     b = 2 * (dp.x * point1.x + dp.y * point1.y + dp.z * point1.z);
     c = point1 * point1 - radius * radius;
     bb4ac = b * b - 4 * a * c;
     if (bb4ac >= 0) {
         double alpha1 = (-b + sqrt(bb4ac)) / (2 * a);
         double alpha2 = (-b - sqrt(bb4ac)) / (2 * a);
         bool inside1 = false;
         if (alpha1 < 0)
             alpha1 = 0;
         else if (alpha1 > 1)
             alpha1 = 1;
         else
             inside1 = true;
         bool inside2 = false;
         if (alpha2 < 0)
             alpha2 = 0;
         else if (alpha2 > 1)
             alpha2 = 1;
         else
             inside2 = true;
         if (alpha1 == alpha2)
             return false;
         else {
             intersection1 = point1 * (1 - alpha1) + point2 * alpha1;
             intersection2 = point1 * (1 - alpha2) + point2 * alpha2;
             if (inside1)
                 normal1 = intersection1 / intersection1.length();
             if (inside2)
                 normal2 = intersection2 / intersection2.length();
             return true;
         }
     }
     else
         return false;
 }