Simbody  3.6
Geo_Box.h
Go to the documentation of this file.
1 #ifndef SimTK_SIMMATH_GEO_BOX_H_
2 #define SimTK_SIMMATH_GEO_BOX_H_
3 
4 /* -------------------------------------------------------------------------- *
5  * Simbody(tm): SimTKmath *
6  * -------------------------------------------------------------------------- *
7  * This is part of the SimTK biosimulation toolkit originating from *
8  * Simbios, the NIH National Center for Physics-Based Simulation of *
9  * Biological Structures at Stanford, funded under the NIH Roadmap for *
10  * Medical Research, grant U54 GM072970. See https://simtk.org/home/simbody. *
11  * *
12  * Portions copyright (c) 2011-14 Stanford University and the Authors. *
13  * Authors: Michael Sherman *
14  * Contributors: *
15  * *
16  * Licensed under the Apache License, Version 2.0 (the "License"); you may *
17  * not use this file except in compliance with the License. You may obtain a *
18  * copy of the License at http://www.apache.org/licenses/LICENSE-2.0. *
19  * *
20  * Unless required by applicable law or agreed to in writing, software *
21  * distributed under the License is distributed on an "AS IS" BASIS, *
22  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. *
23  * See the License for the specific language governing permissions and *
24  * limitations under the License. *
25  * -------------------------------------------------------------------------- */
26 
30 #include "SimTKcommon.h"
32 #include "simmath/internal/Geo.h"
33 
34 #include <cassert>
35 #include <cmath>
36 #include <algorithm>
37 
38 namespace SimTK {
39 
40 //==============================================================================
41 // GEO BOX
42 //==============================================================================
47 template <class P>
48 class Geo::Box_ {
49 typedef P RealP;
50 typedef Vec<2,P> Vec2P;
51 typedef Vec<3,P> Vec3P;
52 typedef UnitVec<P,1> UnitVec3P;
53 typedef Mat<3,3,P> Mat33P;
54 typedef Rotation_<P> RotationP;
55 typedef Transform_<P> TransformP;
56 
57 public:
59 Box_() {}
62 Box_(const Vec3P& halfLengths) {setHalfLengths(halfLengths);}
63 
66 Box_& setHalfLengths(const Vec3P& halfLengths) {
67  SimTK_ERRCHK3(halfLengths >= 0, "Geo::Box_::setHalfLengths()",
68  "Half lengths must be nonnegative; got %g,%g,%g.",
69  (double)halfLengths[0],(double)halfLengths[1],(double)halfLengths[2]);
70  h = halfLengths;
71  sortEdges();
72  return *this;
73 }
74 
77 Box_& addToHalfLengths(const Vec3P& incr) {
78  h += incr;
79  SimTK_ERRCHK3(h >= 0, "Geo::Box_::addToHalfLengths()",
80  "Half lengths must be nonnegative but were %g,%g,%g after change.",
81  (double)h[0],(double)h[1],(double)h[2]);
82  sortEdges();
83  return *this;
84 }
85 
88 const Vec3P& getHalfLengths() const {return h;}
89 
91 RealP getOrderedHalfLength(int i) const {
92  SimTK_INDEXCHECK(i, 3, "Geo::Box_::getOrderedHalfLength()");
93  return h[order[i]];
94 }
95 
98  SimTK_INDEXCHECK(i, 3, "Geo::Box_::getOrderedAxis()");
99  return CoordinateAxis(order[i]);
100 }
101 
103 RealP findVolume() const {return 8*h[0]*h[1]*h[2];}
106 RealP findArea() const {return 8*(h[0]*h[1] + h[0]*h[2] + h[1]*h[2]);}
107 
111 bool containsPoint(const Vec3P& pt) const {
112  const Vec3P absPt = pt.abs(); // reflect to first quadrant
113  return absPt <= h;
114 }
115 
121 Vec3P findClosestPointOfSolidBox(const Vec3P& pt, bool& ptWasInside) const {
122  Vec3P c(pt);
123  ptWasInside = true; // tentatively
124  for (int i=0; i<3; ++i) {
125  if (c[i] < -h[i]) {c[i]=-h[i]; ptWasInside=false;}
126  else if (c[i] > h[i]) {c[i]= h[i]; ptWasInside=false;}
127  }
128  return c;
129 }
130 
135 Vec3P findClosestPointOnSurface(const Vec3P& pt, bool& ptWasInside) const {
136  Vec3P c = findClosestPointOfSolidBox(pt, ptWasInside);
137 
138  if (ptWasInside) { // every |c[i]| <= h[i]
139  RealP dToSide = h[0]-std::abs(c[0]); // distance to closest x-face
140  int which=0; RealP minDist=dToSide;
141  dToSide = h[1]-std::abs(c[1]);
142  if (dToSide < minDist) {which=1; minDist=dToSide;}
143  dToSide = h[2]-std::abs(c[2]);
144  if (dToSide < minDist) {which=2; minDist=dToSide;}
145  // Now project the point to the nearest side.
146  c[which] = c[which] < 0 ? -h[which] : h[which];
147  }
148  return c;
149 }
150 
155 RealP findDistanceSqrToPoint(const Vec3P& pt) const {
156  const Vec3P absPt = pt.abs(); // reflect to first quadrant
157  RealP d2 = 0;
158  if (absPt[0] > h[0]) d2 += square(absPt[0]-h[0]);
159  if (absPt[1] > h[1]) d2 += square(absPt[1]-h[1]);
160  if (absPt[2] > h[2]) d2 += square(absPt[2]-h[2]);
161  return d2;
162 }
163 
164 
171 Vec3P findSupportPoint(const Vec3& d) const {
172  // Basically transferring the sign from d to h, but with 0 treated as 1.
173  return Vec3P(d[0]<0 ? -h[0]:h[0], d[1]<0 ? -h[1]:h[1], d[2]<0 ? -h[2]:h[2]);
174 }
175 
180 RealP findDistanceSqrToSphere(const Geo::Sphere_<P>& sphere) const {
181  const Vec3P absCtr = sphere.getCenter().abs(); // reflect to first quadrant
182  const Vec3P grow = h + sphere.getRadius(); // 3 flops
183  RealP d2 = 0;
184  if (absCtr[0] > grow[0]) d2 += square(absCtr[0]-grow[0]);
185  if (absCtr[1] > grow[1]) d2 += square(absCtr[1]-grow[1]);
186  if (absCtr[2] > grow[2]) d2 += square(absCtr[2]-grow[2]);
187  return d2;
188 }
189 
195  const Vec3P absCtr = aab.getCenter().abs(); // reflect to first quadrant
196  const Vec3P grow = h + aab.getHalfLengths();
197  RealP d2 = 0;
198  if (absCtr[0] > grow[0]) d2 += square(absCtr[0]-grow[0]);
199  if (absCtr[1] > grow[1]) d2 += square(absCtr[1]-grow[1]);
200  if (absCtr[2] > grow[2]) d2 += square(absCtr[2]-grow[2]);
201  return d2;
202 }
203 
208 bool intersectsSphere(const Geo::Sphere_<P>& sphere) const {
209  const Vec3P absCtr = sphere.getCenter().abs(); // reflect to first quadrant
210  const RealP r = sphere.getRadius();
211  if (absCtr[0] > h[0]+r) return false;
212  if (absCtr[1] > h[1]+r) return false;
213  if (absCtr[2] > h[2]+r) return false;
214  return true;
215 }
216 
222  const Vec3P absCtr = aab.getCenter().abs(); // reflect to first quadrant
223  const Vec3P& aabh = aab.getHalfLengths();
224  if (absCtr[0] > h[0]+aabh[0]) return false;
225  if (absCtr[1] > h[1]+aabh[1]) return false;
226  if (absCtr[2] > h[2]+aabh[2]) return false;
227  return true;
228 }
229 
240 
250 
251 
255 static int getNumVertices() {return 8;}
256 static int getNumEdges() {return 12;}
257 static int getNumFaces() {return 6;}
258 
268 Vec3P getVertexPos(int vx) const {
269  SimTK_INDEXCHECK(vx,8,"Geo::Box::getVertexPos()");
270  return Vec3P(vx&0x4 ? h[0]:-h[0], vx&0x2 ? h[1]:-h[1], vx&0x1 ? h[2]:-h[2]);
271 }
272 
275 int findSupportVertex(const Vec3P& d) const {
276  int vx = 0;
277  if (d[0] >= 0) vx += 0x4; // see table in getVertexPos().
278  if (d[1] >= 0) vx += 0x2;
279  if (d[2] >= 0) vx += 0x1;
280  return vx;
281 }
282 
285 UnitVec3P getVertexNormal(int vx) const {
286  SimTK_INDEXCHECK(vx,8,"Geo::Box::getVertexNormal()");
287  const RealP c = 1/(RealP)SimTK_SQRT3, nc = -c;
288  return UnitVec3P(Vec3P(vx&0x4 ? c:nc, vx&0x2 ? c:nc, vx&0x1 ? c:nc),true);
289 }
290 
292 Vec3P getEdgeCenter(int ex) const {
293  SimTK_INDEXCHECK(ex,12,"Geo::Box::getEdgeCenter()");
294  int faces[2], which[2];
295  getEdgeFaces(ex, faces, which);
296  Vec3P c(0);
297  for (int i=0; i<2; ++i) {
298  const CoordinateDirection fd = getFaceCoordinateDirection(faces[i]);
299  c[fd.getAxis()] = fd.getDirection() * h[fd.getAxis()];
300  }
301  return c;
302 }
303 
306 UnitVec3P getEdgeNormal(int ex) const {
307  SimTK_INDEXCHECK(ex,12,"Geo::Box::getEdgeNormal()");
308  const RealP oosqrt2 = (RealP)SimTK_OOSQRT2;
309  int faces[2], which[2];
310  getEdgeFaces(ex, faces, which);
311  Vec3P n(0);
312  for (int i=0; i<2; ++i) {
313  const CoordinateDirection fd = getFaceCoordinateDirection(faces[i]);
314  n[fd.getAxis()] = fd.getDirection() * oosqrt2;
315  }
316  return UnitVec3P(n, true);
317 }
318 
319 
320 
324  SimTK_INDEXCHECK(ex,12,"Geo::Box::getEdgeDirection()");
325  static const int axis[12] = {2,1,2,1,0,2,0,0,1,1,2,0};
326  return CoordinateDirection(CoordinateAxis(axis[ex]), 1); // all in + dir
327 }
328 
329 
335 UnitVec3P getEdgeDirection(int ex) const {
336  SimTK_INDEXCHECK(ex,12,"Geo::Box::getEdgeUnitVec()");
337  return UnitVec3P(getEdgeCoordinateDirection(ex));
338 }
339 
342  SimTK_INDEXCHECK(fx,6,"Geo::Box::getFaceDirection()");
343  static const int axis[6] = {0, 1, 2, 0, 1, 2};
344  static const int dir[6] = {-1,-1,-1, 1, 1, 1};
345  return CoordinateDirection(CoordinateAxis(axis[fx]), dir[fx]);
346 }
347 
351 Vec3P getFaceCenter(int fx) const {
352  SimTK_INDEXCHECK(fx,6,"Geo::Box::getFaceCenter()");
354  Vec3P c(0);
355  c[fd.getAxis()] = fd.getDirection() * h[fd.getAxis()];
356  return c;
357 }
358 
361 UnitVec3P getFaceNormal(int fx) const {
362  SimTK_INDEXCHECK(fx,6,"Geo::Box::getFaceNormal()");
363  return UnitVec3P(getFaceCoordinateDirection(fx));
364 }
365 
367 void getFaceVertices(int fx, int v[4]) const {
368  SimTK_INDEXCHECK(fx,6,"Geo::Box::getFaceVertices()");
369  static const int verts[6][4] = {{0,1,3,2},{0,4,5,1},{0,2,6,4},
370  {7,5,4,6},{7,6,2,3},{7,3,1,5}};
371  for (int i=0; i<4; ++i) v[i] = verts[fx][i];
372 }
373 
376 void getVertexFaces(int vx, int f[3], int w[3]) const {
377  SimTK_INDEXCHECK(vx,8,"Geo::Box::getVertexFaces()");
378  static const int faces[8][3] = {{0,1,2},{0,1,5},{0,2,4},{0,4,5},
379  {1,2,3},{1,3,5},{2,3,4},{3,4,5}};
380  static const int which[8][3] = {{0,0,0},{1,3,2},{3,1,2},{2,3,1},
381  {1,3,2},{2,1,3},{2,3,1},{0,0,0}};
382  for (int i=0; i<3; ++i) {f[i]=faces[vx][i]; w[i]=which[vx][i];}
383 }
384 
386 void getEdgeVertices(int ex, int v[2]) const {
387  SimTK_INDEXCHECK(ex,12,"Geo::Box::getEdgeVertices()");
388  static const int verts[12][2] = {{0,1},{1,3},{2,3},{0,2}, // 0-3
389  {0,4},{4,5},{1,5},{2,6}, // 4-7
390  {4,6},{5,7},{6,7},{3,7}}; // 8-11
391  for (int i=0; i<2; ++i) v[i] = verts[ex][i];
392 }
395 void getVertexEdges(int vx, int e[3], int w[3]) const {
396  SimTK_INDEXCHECK(vx,8,"Geo::Box::getVertexEdges()");
397  static const int edges[8][3] = {{0,3,4},{0,1,6},{2,3,7},{1,2,11},
398  {4,5,8},{5,6,9},{7,8,10},{9,10,11}};
399  static const int which[8][3] = {{0,0,0},{1,0,0},{0,1,0},{1,1,0},
400  {1,0,0},{1,1,0},{1,1,0},{1,1,1}};
401  for (int i=0; i<3; ++i) {e[i]=edges[vx][i]; w[i]=which[vx][i];}
402 }
403 
404 
407 void getFaceEdges(int fx, int e[4]) const {
408  SimTK_INDEXCHECK(fx,6,"Geo::Box::getFaceEdges()");
409  static const int edges[6][4] = {{0,1,2, 3},{ 4,5,6, 0},{ 3,7,8,4},
410  {9,5,8,10},{10,7,2,11},{11,1,6,9}};
411  for (int i=0; i<4; ++i) e[i] = edges[fx][i];
412 }
415 void getEdgeFaces(int ex, int f[2], int w[2]) const {
416  SimTK_INDEXCHECK(ex,12,"Geo::Box::getEdgeFaces()");
417  static const int faces[12][2] = {{0,1},{0,5},{0,4},{0,2}, // 0-3
418  {1,2},{1,3},{1,5},{2,4}, // 4-7
419  {2,3},{3,5},{3,4},{4,5}}; // 8-11
420  static const int which[12][2] = {{0,3},{1,1},{2,2},{3,0}, // 0-3
421  {0,3},{1,1},{2,2},{1,1}, // 4-7
422  {2,2},{0,3},{3,0},{3,0}}; // 8-11
423  for (int i=0; i<2; ++i) {f[i] = faces[ex][i]; w[i] = which[ex][i];}
424 }
425 
426 
427 
428 private:
429 // Call this whenever an edge length changes. Each axis will appear once.
430 void sortEdges() {
431  CoordinateAxis shortest = XAxis, longest = ZAxis;
432  if (h[YAxis] < h[shortest]) shortest=YAxis;
433  if (h[ZAxis] < h[shortest]) shortest=ZAxis;
434  if (h[XAxis] > h[longest]) longest=XAxis;
435  if (h[YAxis] > h[longest]) longest=YAxis;
436  order[0] = shortest; order[2] = longest;
437  order[1] = shortest.getThirdAxis(longest); // not shortest or longest
438 }
439 
440 int intersectsOrientedBoxHelper(const OrientedBox_<P>& O,
441  Mat33P& absR_BO,
442  Vec3P& absP_BO) const;
443 
444 Vec3P h; // half-dimensions of the box
445 unsigned char order[3]; // 0,1,2 reordered short to long
446 };
447 
448 
449 
450 //==============================================================================
451 // GEO ALIGNED BOX
452 //==============================================================================
456 template <class P>
458 typedef P RealP;
459 typedef Vec<3,RealP> Vec3P;
460 
461 public:
467 AlignedBox_(const Vec3P& center, const Geo::Box_<P>& box)
468 : center(center), box(box) {}
471 AlignedBox_(const Vec3P& center, const Vec3P& halfLengths)
472 : center(center), box(halfLengths) {}
473 
475 AlignedBox_& setCenter(const Vec3P& center)
476 { this->center=center; return *this; }
477 
479 AlignedBox_& setHalfLengths(const Vec3P& halfLengths)
480 { box.setHalfLengths(halfLengths); return *this; }
481 
483 const Vec3P& getCenter() const {return center;}
485 Vec3P& updCenter() {return center;}
488 const Vec3P& getHalfLengths() const {return box.getHalfLengths();}
489 // no updHalfLengths()
490 const Box_<P>& getBox() const {return box;}
491 Box_<P>& updBox() {return box;}
492 
496 bool containsPoint(const Vec3P& pt_F) const
497 { return box.containsPoint(pt_F - center); } // shift to box frame B
498 
507  const RealP tol = Geo::getDefaultTol<P>();
508  const RealP maxdim = max(getCenter().abs());
509  const RealP maxrad = max(getHalfLengths());
510  const RealP scale = std::max(maxdim, maxrad);
511  const RealP incr = std::max(scale*Geo::getEps<P>(), tol);
512  box.addToHalfLengths(Vec3P(incr));
513  return *this;
514 }
515 
516 private:
517 Vec3P center;
518 Geo::Box_<P> box;
519 };
520 
521 
522 //==============================================================================
523 // GEO ORIENTED BOX
524 //==============================================================================
527 template <class P>
529 typedef P RealP;
530 typedef Vec<3,P> Vec3P;
531 typedef Rotation_<P> RotationP;
532 typedef Transform_<P> TransformP;
533 
534 public:
541 OrientedBox_(const TransformP& X_FB, const Geo::Box_<P>& box)
542 : X_FB(X_FB), box(box) {}
545 OrientedBox_(const TransformP& X_FB, const Vec3P& halfLengths)
546 : X_FB(X_FB), box(halfLengths) {}
547 
548 
550 OrientedBox_& setTransform(const TransformP& newX_FB)
551 { X_FB=newX_FB; return *this; }
552 
554 OrientedBox_& setHalfLengths(const Vec3P& halfLengths)
555 { box.setHalfLengths(halfLengths); return *this; }
556 
557 const Vec3P& getCenter() const {return X_FB.p();}
558 Vec3P& updCenter() {return X_FB.updP();}
559 const RotationP& getOrientation() const {return X_FB.R();}
560 RotationP& updOrientation() {return X_FB.updR();}
561 const TransformP& getTransform() const {return X_FB;}
562 TransformP& updTransform() {return X_FB;}
563 const Vec3P& getHalfLengths() const {return box.getHalfLengths();}
564 // no updHalfLengths()
565 const Box_<P>& getBox() const {return box;}
566 Box_<P>& updBox() {return box;}
567 
568 
572 bool containsPoint(const Vec3P& pt_F) const
573 { return box.containsPoint(~X_FB*pt_F); } // shift to box frame B
574 
583  const RealP tol = Geo::getDefaultTol<P>();
584  const RealP maxdim = max(getCenter().abs());
585  const RealP maxrad = max(getHalfLengths());
586  const RealP scale = std::max(maxdim, maxrad);
587  const RealP incr = std::max(scale*Geo::getEps<P>(), tol);
588  box.addToHalfLengths(Vec3P(incr));
589  return *this;
590 }
591 
592 private:
593 TransformP X_FB;
594 Geo::Box_<P> box;
595 };
596 
597 
598 } // namespace SimTK
599 
600 #endif // SimTK_SIMMATH_GEO_BOX_H_
bool mayIntersectOrientedBox(const Geo::OrientedBox_< P > &ob) const
Given an oriented box whose pose is measured and expressed in the frame of this box, return true if the two boxes may be intersecting.
AlignedBox_(const Vec3P &center, const Geo::Box_< P > &box)
Construct an AlignedBox with the given box shape with the center located as given.
Definition: Geo_Box.h:467
void getVertexFaces(int vx, int f[3], int w[3]) const
Each vertex has three incident faces.
Definition: Geo_Box.h:376
Box_< P > & updBox()
Definition: Geo_Box.h:566
AlignedBox_ & setHalfLengths(const Vec3P &halfLengths)
Change the dimensions of this box.
Definition: Geo_Box.h:479
bool containsPoint(const Vec3P &pt) const
Given a point measured and expressed in the box frame, determine whether it is inside the box (we cou...
Definition: Geo_Box.h:111
bool containsPoint(const Vec3P &pt_F) const
Given a point measured and expressed in the base frame F, determine whether it is strictly contained ...
Definition: Geo_Box.h:496
void getVertexEdges(int vx, int e[3], int w[3]) const
Each vertex has three incident edges.
Definition: Geo_Box.h:395
RealP findDistanceSqrToAlignedBox(const Geo::AlignedBox_< P > &aab) const
Return the square of the distance from this box to an axis-aligned box whose center location is measu...
Definition: Geo_Box.h:194
TransformP & updTransform()
Definition: Geo_Box.h:562
Box_ & setHalfLengths(const Vec3P &halfLengths)
Change the half-dimensions of this box.
Definition: Geo_Box.h:66
TAbs abs() const
Elementwise absolute value; that is, the return value has the same dimension as this Vec but with eac...
Definition: Vec.h:347
CoordinateAxis getOrderedAxis(int i) const
Get axes in order shortest to longest; 0 is shortest, 2 is longest.
Definition: Geo_Box.h:97
static int getNumEdges()
Use bits in the vertex number to pick the signs, with 0=negative, 1=positive:
Definition: Geo_Box.h:256
This is the top-level SimTK namespace into which all SimTK names are placed to avoid collision with o...
Definition: Assembler.h:37
const Vec3P & getHalfLengths() const
Return the half-lengths of this box as a Vec3 from the center to the first quadrant vertex...
Definition: Geo_Box.h:488
void getFaceVertices(int fx, int v[4]) const
A face has four vertices ordered counterclockwise about the face normal.
Definition: Geo_Box.h:367
A CoordinateDirection is a CoordinateAxis plus a direction indicating the positive or negative direct...
Definition: CoordinateAxis.h:244
The Rotation class is a Mat33 that guarantees that the matrix can be interpreted as a legitimate 3x3 ...
Definition: Quaternion.h:40
CoordinateAxis getAxis() const
This is the coordinate axis XAxis, YAxis, or ZAxis contained in this CoordinateDirection. Use getDirection() to determine whether this is the positive or negative direction.
Definition: CoordinateAxis.h:274
bool intersectsOrientedBox(const Geo::OrientedBox_< P > &ob) const
Given an oriented box whose pose is measured and expressed in the frame of this box, return true if the two boxes intersect.
bool containsPoint(const Vec3P &pt_F) const
Given a point measured and expressed in the base frame F, determine whether it is strictly contained ...
Definition: Geo_Box.h:572
A 3d rectangular box aligned with an unspecified frame F and centered at that frame&#39;s origin...
Definition: Geo.h:61
AlignedBox_()
Construct an uninitialized AlignedBox object; the dimensions and location will be garbage...
Definition: Geo_Box.h:464
const CoordinateAxis::ZCoordinateAxis ZAxis
Constant representing the Z coordinate axis; will implicitly convert to the integer 2 when used in a ...
AlignedBox_ & stretchBoundary()
Stretch this box in place by a small amount to ensure that there will be no roundoff problems if this...
Definition: Geo_Box.h:506
Vec3P getVertexPos(int vx) const
Use bits in the vertex number to pick the signs, with 0=negative, 1=positive:
Definition: Geo_Box.h:268
This class, along with its sister class CoordinateDirection, provides convenient manipulation of the ...
Definition: CoordinateAxis.h:53
#define SimTK_ERRCHK3(cond, whereChecked, fmt, a1, a2, a3)
Definition: ExceptionMacros.h:330
OrientedBox_ & setTransform(const TransformP &newX_FB)
Change the pose of this box.
Definition: Geo_Box.h:550
TODO: A 3d box oriented and positioned with respect to an unspecified frame F.
Definition: Geo.h:63
A 3d box aligned with an unspecified frame F and centered at a given point measured from that frame&#39;s...
Definition: Geo.h:62
const CoordinateAxis::YCoordinateAxis YAxis
Constant representing the Y coordinate axis; will implicitly convert to the integer 1 when used in a ...
unsigned char square(unsigned char u)
Definition: Scalar.h:349
Vec3P & updCenter()
Definition: Geo_Box.h:558
const Vec3P & getHalfLengths() const
Definition: Geo_Box.h:563
void getFaceEdges(int fx, int e[4]) const
A face has four edges, ordered by the vertex ordering: v0-v1, v1-v2, v2-v3, v3-v1.
Definition: Geo_Box.h:407
static int getNumFaces()
Use bits in the vertex number to pick the signs, with 0=negative, 1=positive:
Definition: Geo_Box.h:257
bool intersectsSphere(const Geo::Sphere_< P > &sphere) const
Given a sphere with center measured and expressed in the box frame, return true if the box and sphere...
Definition: Geo_Box.h:208
const Box_< P > & getBox() const
Definition: Geo_Box.h:565
int getDirection() const
Returns 1 or -1 to indicate the direction along the coordinate axis returned by getAxis().
Definition: CoordinateAxis.h:277
int findSupportVertex(const Vec3P &d) const
Find the vertex (0-7) that is furthest in the direction d, which is given in the box frame...
Definition: Geo_Box.h:275
RealP findDistanceSqrToPoint(const Vec3P &pt) const
Return the square of the distance from this box to a given point whose location is measured from and ...
Definition: Geo_Box.h:155
static int getNumVertices()
Use bits in the vertex number to pick the signs, with 0=negative, 1=positive:
Definition: Geo_Box.h:255
Includes internal headers providing declarations for the basic SimTK Core classes, including Simmatrix.
AlignedBox_ & setCenter(const Vec3P &center)
Change the center location of this box.
Definition: Geo_Box.h:475
UnitVec3P getEdgeNormal(int ex) const
Edge normals point diagonally outwards from the edge.
Definition: Geo_Box.h:306
CoordinateDirection getEdgeCoordinateDirection(int ex) const
Return the direction of an edge, going from its first vertex to its second vertex, as a CoordinateDirection.
Definition: Geo_Box.h:323
Vec3P findClosestPointOfSolidBox(const Vec3P &pt, bool &ptWasInside) const
Given a point location in the box frame, return the closest point of the solid box, and a flag saying whether the given point was inside the box, using the same definition of "inside" as the containsPoint() method.
Definition: Geo_Box.h:121
Box_(const Vec3P &halfLengths)
Construct a Box with the given nonnegative half-dimensions.
Definition: Geo_Box.h:62
const RotationP & getOrientation() const
Definition: Geo_Box.h:559
This class represents the rotate-and-shift transform which gives the location and orientation of a ne...
Definition: Transform.h:43
RealP getOrderedHalfLength(int i) const
Get lengths in order shortest to longest; 0 is shortest, 2 is longest.
Definition: Geo_Box.h:91
ELEM max(const VectorBase< ELEM > &v)
Definition: VectorMath.h:251
Vec3P findClosestPointOnSurface(const Vec3P &pt, bool &ptWasInside) const
Given a point location in the box frame, return the closest point on the box surface, and a flag saying whether the given point was inside the box, using the same definition of "inside" as the containsPoint() method.
Definition: Geo_Box.h:135
Defines geometric primitive shapes and algorthms.
RealP findVolume() const
Calculate the volume of this box.
Definition: Geo_Box.h:103
void getEdgeFaces(int ex, int f[2], int w[2]) const
An edge is between two faces.
Definition: Geo_Box.h:415
#define SimTK_OOSQRT2
One over the square root of 2; also half the square root of 2 since 1/sqrt(2) == 2^(-1/2) == sqrt(2)/...
Definition: Constants.h:184
void getEdgeVertices(int ex, int v[2]) const
An edge connects two vertices.
Definition: Geo_Box.h:386
Vec3P & updCenter()
Return a writable reference to the center location of this box.
Definition: Geo_Box.h:485
RealP findDistanceSqrToSphere(const Geo::Sphere_< P > &sphere) const
Return the square of the distance from this box to a given sphere whose center location is measured f...
Definition: Geo_Box.h:180
Box_ & addToHalfLengths(const Vec3P &incr)
Change the half-dimensions of this box by adding the given vector.
Definition: Geo_Box.h:77
OrientedBox_ & setHalfLengths(const Vec3P &halfLengths)
Change the dimensions of this box.
Definition: Geo_Box.h:554
Box_< P > & updBox()
Definition: Geo_Box.h:491
UnitVec3P getFaceNormal(int fx) const
Return the outward normal for the given face as a unit vector in the box frame.
Definition: Geo_Box.h:361
RowVectorBase< typename CNT< ELEM >::TAbs > abs(const RowVectorBase< ELEM > &v)
Definition: VectorMath.h:120
const Vec3P & getHalfLengths() const
Return the half-lengths of this box as a Vec3 from the center to the first quadrant vertex...
Definition: Geo_Box.h:88
CoordinateDirection getFaceCoordinateDirection(int fx) const
Return the outward normal for the given face as a CoordinateDirection.
Definition: Geo_Box.h:341
const Box_< P > & getBox() const
Definition: Geo_Box.h:490
This is the header file that every Simmath compilation unit should include first. ...
const Vec3P & getCenter() const
Return the location of the center of this box (box frame origin).
Definition: Geo_Box.h:483
OrientedBox_ & stretchBoundary()
Stretch this box in place by a small amount to ensure that there will be no roundoff problems if this...
Definition: Geo_Box.h:582
const TransformP & getTransform() const
Definition: Geo_Box.h:561
UnitVec3P getVertexNormal(int vx) const
Vertex normals point diagonally outwards from the box corners.
Definition: Geo_Box.h:285
RotationP & updOrientation()
Definition: Geo_Box.h:560
const Vec3P & getCenter() const
Definition: Geo_Box.h:557
#define SimTK_INDEXCHECK(ix, ub, where)
Definition: ExceptionMacros.h:145
A geometric primitive representing a sphere by its radius and center point, and a collection of spher...
Definition: Geo.h:56
AlignedBox_(const Vec3P &center, const Vec3P &halfLengths)
Construct an AlignedBox with the given center location and half-dimensions.
Definition: Geo_Box.h:471
OrientedBox_()
Construct an uninitialized OrientedBox object; the dimensions and pose will be garbage.
Definition: Geo_Box.h:537
OrientedBox_(const TransformP &X_FB, const Geo::Box_< P > &box)
Construct an OrientedBox with the given box shape with positioned and oriented according to the given...
Definition: Geo_Box.h:541
Vec3P getEdgeCenter(int ex) const
Return the center point of the specified edge, in the box frame.
Definition: Geo_Box.h:292
Vec3P getFaceCenter(int fx) const
Return the center point position for the given face.
Definition: Geo_Box.h:351
OrientedBox_(const TransformP &X_FB, const Vec3P &halfLengths)
Construct an OrientedBox with the given location and half-dimensions.
Definition: Geo_Box.h:545
bool intersectsAlignedBox(const Geo::AlignedBox_< P > &aab) const
Given an aligned box with center measured and expressed in the from of this box, return true if the t...
Definition: Geo_Box.h:221
Box_()
Construct an uninitialized Box object; the dimensions will be garbage.
Definition: Geo_Box.h:59
UnitVec3P getEdgeDirection(int ex) const
Return a unit vector aligned with the selected edge, pointing in the direction from the first vertex ...
Definition: Geo_Box.h:335
#define SimTK_SIMMATH_EXPORT
Definition: SimTKmath/include/simmath/internal/common.h:64
#define SimTK_SQRT3
The square root of 3.
Definition: Constants.h:214
Vec3P findSupportPoint(const Vec3 &d) const
Find a supporting point on the surface of the box in the given direction, which must be expressed in ...
Definition: Geo_Box.h:171
RealP getRadius() const
Get the sphere&#39;s radius.
Definition: Geo_Sphere.h:111
const Vec3P & getCenter() const
Get the location of the sphere&#39;s center point.
Definition: Geo_Sphere.h:107
RealP findArea() const
Calculate the surface area of this box, ignoring degeneracy (meaning that all pairs of sides are coun...
Definition: Geo_Box.h:106
const CoordinateAxis::XCoordinateAxis XAxis
Constant representing the X coordinate axis; will implicitly convert to the integer 0 when used in a ...