1 // Copyright (C) 2006 Mathias Froehlich - Mathias.Froehlich@web.de
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Library General Public
5 // License as published by the Free Software Foundation; either
6 // version 2 of the License, or (at your option) any later version.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Library General Public License for more details.
13 // You should have received a copy of the GNU General Public License
14 // along with this program; if not, write to the Free Software
15 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
27 // Avoid "_center not initialized" warnings.
28 # pragma GCC diagnostic ignored "-Wuninitialized"
33 * Do not initialize _center to save unneeded initialization time.
34 * Fix 'may be used uninitialized' warnings locally instead
36 // _center(0.0, 0.0, 0.0),
39 SGSphere(const SGVec3<T>& center, const T& radius) :
44 explicit SGSphere(const SGSphere<S>& sphere) :
45 _center(sphere.getCenter()),
46 _radius(sphere.getRadius())
50 // Restore warning settings.
51 # pragma GCC diagnostic warning "-Wuninitialized"
54 const SGVec3<T>& getCenter() const
56 void setCenter(const SGVec3<T>& center)
59 const T& getRadius() const
61 void setRadius(const T& radius)
64 { return _radius*_radius; }
70 { return 0 <= _radius; }
75 /// Return true if this is inside sphere
76 bool inside(const SGSphere<T>& sphere) const
83 T dist = sphere.getRadius() - getRadius();
84 return distSqr(getCenter(), sphere.getCenter()) <= dist*dist;
87 void expandBy(const SGVec3<T>& v)
95 T dist2 = distSqr(_center, v);
96 if (dist2 <= getRadius2())
100 T newRadius = T(0.5)*(_radius + dist);
101 _center += ((newRadius - _radius)/dist)*(v - _center);
105 void expandBy(const SGSphere<T>& s)
111 _center = s.getCenter();
112 _radius = s.getRadius();
116 T dist = length(_center - s.getCenter());
117 if (dist <= SGLimits<T>::min()) {
118 _radius = SGMisc<T>::max(_radius, s._radius);
123 if (dist + s.getRadius() <= _radius)
126 // new one includes all
127 if (dist + _radius <= s.getRadius()) {
128 _center = s.getCenter();
129 _radius = s.getRadius();
133 T newRadius = T(0.5)*(_radius + dist + s.getRadius());
134 T ratio = (newRadius - _radius) / dist;
137 _center[0] += ratio*(s._center[0] - _center[0]);
138 _center[1] += ratio*(s._center[1] - _center[1]);
139 _center[2] += ratio*(s._center[2] - _center[2]);
142 void expandBy(const SGBox<T>& box)
148 _center = box.getCenter();
149 _radius = T(0.5)*length(box.getSize());
153 SGVec3<T> boxCenter = box.getCenter();
155 for (unsigned i = 0; i < 3; ++i) {
156 if (_center[i] < boxCenter[i])
157 corner[i] = box.getMax()[i];
159 corner[i] = box.getMin()[i];
169 /// Output to an ostream
170 template<typename char_type, typename traits_type, typename T>
172 std::basic_ostream<char_type, traits_type>&
173 operator<<(std::basic_ostream<char_type, traits_type>& s,
174 const SGSphere<T>& sphere)
176 return s << "center = " << sphere.getCenter()
177 << ", radius = " << sphere.getRadius();