8 #ifndef IMPALGEBRA_SPHERE_D_H
9 #define IMPALGEBRA_SPHERE_D_H
19 IMPALGEBRA_BEGIN_NAMESPACE
28 #if IMP_HAS_CHECKS >= IMP_USAGE
29 radius_ = std::numeric_limits<double>::quiet_NaN();
33 : center_(center), radius_(radius) {
36 double get_radius()
const {
38 "Attempt to use uninitialized sphere.");
41 const VectorD<D> &get_center()
const {
return center_; }
44 double d = (get_center() - o.get_center()).get_magnitude();
45 return (d + o.get_radius() < get_radius());
50 return ((p - center_).get_squared_magnitude() <= get_squared(radius_));
53 { out <<
"(" << spaces_io(center_) <<
": " << get_radius() <<
")"; });
56 VectorD<D> &_access_center() {
return center_; }
57 void _set_radius(
double d) { radius_ = d; }
58 double &operator[](
unsigned int i) {
66 double operator[](
unsigned int i)
const {
76 unsigned int get_dimension()
const {
return center_.get_dimension(); }
84 {
return PI * 4.0 * get_squared(g.get_radius()); },
85 {
return PI * (4.0 / 3.0) * std::pow(g.get_radius(), 3.0); },
86 return BoundingBoxD<D>(g.get_center()) +
89 template <
unsigned int D>
91 return SphereD<D>(get_zero_vector_d<D>(), 1.0);
94 inline SphereD<-1> get_unit_sphere_kd(
unsigned int d) {
104 double d = (a.get_center() - b.get_center()).get_magnitude();
105 return d - a.get_radius() - b.get_radius();
115 double d = (a.get_center() - b.get_center()).get_squared_magnitude();
116 return d - square(a.get_radius()) - square(b.get_radius());
124 double sr = a.get_radius() + b.get_radius();
125 for (
unsigned int i = 0; i < a.get_dimension(); ++i) {
126 double delta = std::abs(a.get_center()[i] - b.get_center()[i]);
127 if (delta >= sr)
return false;
132 #if !defined(SWIG) && !defined(IMP_DOXYGEN)
136 struct SphereSpacesIO {
137 const SphereD<D> &v_;
138 SphereSpacesIO(
const SphereD<D> &v) : v_(v) {}
141 inline std::ostream &operator<<(std::ostream &out, const SphereSpacesIO<D> &s) {
142 for (
unsigned int i = 0; i < s.v_.get_center().get_dimension(); ++i) {
143 out << s.v_.get_center()[i] <<
" ";
145 out << s.v_.get_radius();
156 inline internal::SphereSpacesIO<D> spaces_io(
const SphereD<D> &v) {
157 return internal::SphereSpacesIO<D>(v);
163 template <
class Geometry>
166 template <
class Geometry>
169 template <
class Geometry>
172 template <
class Geometry>
177 VectorD<D> get_vector_geometry(
const SphereD<D> &s) {
178 return s.get_center();
181 IMPALGEBRA_END_NAMESPACE
bool get_contains(const SphereD< D > &o) const
Return true if this sphere contains the other one.
#define IMP_SHOWABLE_INLINE(Name, how_to_show)
Declare the methods needed by an object that can be printed.
VectorD< D > get_zero_vector_kd(int Di)
Return a dynamically sized vector of zeros.
static const double PI
the constant pi
double get_volume(const Cone3D &g)
double get_squared_distance(const VectorD< D > &v1, const VectorD< D > &v2)
Compute the squared distance between two vectors.
#define IMP_INTERNAL_CHECK(expr, message)
An assertion to check for internal errors in IMP. An IMP::ErrorException will be thrown.
Base class for geometric types.
Functions to deal with very common math operations.
A Cartesian vector in D-dimensions.
bool get_interiors_intersect(const SphereD< D > &a, const SphereD< D > &b)
Return true if the two balls bounded by the two spheres intersect.
BoundingBoxD< 3 > get_bounding_box(const Cone3D &g)
#define IMP_VOLUME_GEOMETRY_METHODS_D(Name, name, area, volume, bounding_box)
Implement the needed namespace methods for a geometry type.
A bounding box in D dimensions.
double get_surface_area(const Cone3D &g)
double get_area(const Plane3D &g)
Various useful constants.
bool isnan(const T &a)
Return true if a number is NaN.
double get_distance(const SphereD< D > &a, const SphereD< D > &b)
Return the distance between the two spheres if they are disjoint.
#define IMP_USAGE_CHECK(expr, message)
A runtime test for incorrect usage of a class or method.
bool get_contains(const VectorD< D > &p) const
Return true if the point is in or on the surface of the sphere.
Represent a sphere in D-dimensions.
Various important macros for implementing geometry.
double get_power_distance(const SphereD< D > &a, const SphereD< D > &b)
Return the power distance between the two spheres.