31 #ifndef OPENVDB_MATH_COORD_HAS_BEEN_INCLUDED 32 #define OPENVDB_MATH_COORD_HAS_BEEN_INCLUDED 39 namespace tbb {
class split; }
57 using Limits = std::numeric_limits<ValueType>;
62 explicit Coord(
const Vec3i& v): mVec{{v[0], v[1], v[2]}} {}
64 explicit Coord(
const Int32* v): mVec{{v[0], v[1], v[2]}} {}
117 return Coord(mVec[0] + dx, mVec[1] + dy, mVec[2] + dz);
137 return Coord(mVec[0] + rhs[0], mVec[1] + rhs[1], mVec[2] + rhs[2]);
141 return Coord(mVec[0] - rhs[0], mVec[1] - rhs[1], mVec[2] - rhs[2]);
145 Coord operator>> (
size_t n)
const {
return Coord(mVec[0]>>n, mVec[1]>>n, mVec[2]>>n); }
149 Coord operator& (
Int32 n)
const {
return Coord(mVec[0] & n, mVec[1] & n, mVec[2] & n); }
150 Coord operator| (
Int32 n)
const {
return Coord(mVec[0] | n, mVec[1] | n, mVec[2] | n); }
151 Coord& operator&= (
Int32 n) { mVec[0]&=n; mVec[1]&=n; mVec[2]&=n;
return *
this; }
152 Coord& operator|= (
Int32 n) { mVec[0]|=n; mVec[1]|=n; mVec[2]|=n;
return *
this; }
175 return (mVec[0] == rhs.mVec[0] && mVec[1] == rhs.mVec[1] && mVec[2] == rhs.mVec[2]);
182 return this->x() < rhs.
x() ? true : this->x() > rhs.
x() ? false
183 : this->y() < rhs.
y() ? true : this->y() > rhs.
y() ? false
184 : this->z() < rhs.
z() ? true :
false;
189 return this->x() < rhs.
x() ? true : this->x() > rhs.
x() ? false
190 : this->y() < rhs.
y() ? true : this->y() > rhs.
y() ? false
191 : this->z() <=rhs.
z() ? true :
false;
201 mVec[0] =
std::min(mVec[0], other.mVec[0]);
202 mVec[1] =
std::min(mVec[1], other.mVec[1]);
203 mVec[2] =
std::min(mVec[2], other.mVec[2]);
209 mVec[0] =
std::max(mVec[0], other.mVec[0]);
210 mVec[1] =
std::max(mVec[1], other.mVec[1]);
211 mVec[2] =
std::max(mVec[2], other.mVec[2]);
234 return (a[0] < b[0] || a[1] < b[1] || a[2] < b[2]);
243 void read(std::istream& is) { is.read(reinterpret_cast<char*>(mVec.data()),
sizeof(mVec)); }
246 os.write(reinterpret_cast<const char*>(mVec.data()),
sizeof(mVec));
250 std::array<Int32, 3> mVec;
270 template<
bool ZYXOrder>
281 operator bool()
const {
return ZYXOrder ? (mPos[0] <= mMax[0]) : (mPos[2] <= mMax[2]); }
287 return ((mPos == other.mPos) && (mMin == other.mMin) && (mMax == other.mMax));
292 template<
size_t a,
size_t b,
size_t c>
295 if (mPos[a] < mMax[a]) { ++mPos[a]; }
296 else if (mPos[b] < mMax[b]) { mPos[a] = mMin[a]; ++mPos[b]; }
297 else if (mPos[c] <= mMax[c]) { mPos[a] = mMin[a]; mPos[b] = mMin[b]; ++mPos[c]; }
299 Coord mPos, mMin, mMax;
313 : mMin(xMin, yMin, zMin), mMax(xMax, yMax, zMax)
320 assert(this->is_divisible());
321 const size_t n = this->maxExtent();
322 mMax[n] = (mMin[n] + mMax[n]) >> 1;
323 other.mMin[n] = mMax[n] + 1;
369 bool empty()
const {
return (mMin[0] > mMax[0] || mMin[1] > mMax[1] || mMin[2] > mMax[2]); }
371 operator bool()
const {
return !this->empty(); }
388 const Coord d = this->dim();
392 bool is_divisible()
const {
return mMin[0]<mMax[0] && mMin[1]<mMax[1] && mMin[2]<mMax[2]; }
395 size_t minExtent()
const {
return this->dim().minIndex(); }
398 size_t maxExtent()
const {
return this->dim().maxIndex(); }
403 return !(Coord::lessThan(xyz,mMin) || Coord::lessThan(mMax,xyz));
409 return !(Coord::lessThan(b.mMin,mMin) || Coord::lessThan(mMax,b.mMax));
415 return !(Coord::lessThan(mMax,b.mMin) || Coord::lessThan(b.mMax,mMin));
421 mMin.offset(-padding);
422 mMax.offset( padding);
428 return CoordBBox(mMin.offsetBy(-padding),mMax.offsetBy(padding));
434 mMin.minComponent(xyz);
435 mMax.maxComponent(xyz);
441 mMin.minComponent(bbox.
min());
442 mMax.maxComponent(bbox.
max());
447 mMin.maxComponent(bbox.
min());
448 mMax.minComponent(bbox.
max());
454 mMin.minComponent(min);
455 mMax.maxComponent(min.
offsetBy(dim-1));
466 assert(p !=
nullptr);
467 p->
reset(mMin.x(), mMin.y(), mMin.z()); ++p;
468 p->
reset(mMin.x(), mMin.y(), mMax.z()); ++p;
469 p->
reset(mMin.x(), mMax.y(), mMin.z()); ++p;
470 p->
reset(mMin.x(), mMax.y(), mMax.z()); ++p;
471 p->
reset(mMax.x(), mMin.y(), mMin.z()); ++p;
472 p->
reset(mMax.x(), mMin.y(), mMax.z()); ++p;
473 p->
reset(mMax.x(), mMax.y(), mMin.z()); ++p;
474 p->
reset(mMax.x(), mMax.y(), mMax.z());
490 void read(std::istream& is) { mMin.read(is); mMax.read(is); }
492 void write(std::ostream& os)
const { mMin.write(os); mMax.write(os); }
504 os << xyz.
asVec3i();
return os;
542 template <
typename T>
554 template <
typename T>
569 os << b.
min() <<
" -> " << b.
max();
577 #endif // OPENVDB_MATH_COORD_HAS_BEEN_INCLUDED Int32 operator[](size_t i) const
Definition: Coord.h:157
Coord offsetBy(Int32 dx, Int32 dy, Int32 dz) const
Definition: Coord.h:115
math::Vec3< Index32 > Vec3I
Definition: Types.h:76
bool operator!=(const Coord &rhs) const
Definition: Coord.h:177
static Coord floor(const Vec3< T > &xyz)
Return the largest integer coordinates that are not greater than xyz (node centered conversion)...
Definition: Coord.h:80
void read(std::istream &is)
Unserialize this bounding box from the given stream.
Definition: Coord.h:490
Coord operator+(const Coord &rhs) const
Definition: Coord.h:135
Coord & min()
Definition: Coord.h:337
bool operator>(const Coord &rhs) const
Lexicographic greater than.
Definition: Coord.h:194
Int32 z() const
Definition: Coord.h:156
General-purpose arithmetic and comparison routines, most of which accept arbitrary value types (or at...
Iterator & operator++()
Increment the iterator to point to the next coordinate.
Definition: Coord.h:279
std::numeric_limits< ValueType > Limits
Definition: Coord.h:57
Coord & setZ(Int32 z)
Definition: Coord.h:105
XYZIterator beginXYZ() const
Return an XYZ-order iterator that points to the minimum coordinate.
Definition: Coord.h:356
Coord()
Definition: Coord.h:59
Coord(Int32 xyz)
Definition: Coord.h:60
bool hasVolume() const
Return true if this bounding box is nonempty (i.e., encloses at least one coordinate).
Definition: Coord.h:373
void write(std::ostream &os) const
Definition: Coord.h:244
bool operator==(const CoordBBox &rhs) const
Definition: Coord.h:365
CoordBBox(ValueType xMin, ValueType yMin, ValueType zMin, ValueType xMax, ValueType yMax, ValueType zMax)
Construct from individual components of the min and max bounds.
Definition: Coord.h:311
bool is_divisible() const
Return true if this bounding box can be subdivided [mainly for use by TBB].
Definition: Coord.h:392
Vec3< typename promote< T, typename Coord::ValueType >::type > operator+(const Coord &v1, const Vec3< T > &v0)
Allow a Coord to be added to or subtracted from a Vec3.
Definition: Coord.h:530
size_t MinIndex(const Vec3T &v)
Return the index [0,1,2] of the smallest value in a 3D vector.
Definition: Math.h:886
void read(std::istream &is)
Definition: Coord.h:243
static Coord min()
Return the smallest possible coordinate.
Definition: Coord.h:67
Index64 volume() const
Return the integer volume of coordinates spanned by this bounding box.
Definition: Coord.h:386
Coord & operator-=(const Coord &rhs)
Definition: Coord.h:128
Int32 y() const
Definition: Coord.h:155
Coord & max()
Definition: Coord.h:338
uint32_t Index32
Definition: Types.h:55
ZYXIterator beginZYX() const
Return a ZYX-order iterator that points to the minimum coordinate.
Definition: Coord.h:354
Coord & operator>>=(size_t n)
Definition: Coord.h:148
ZYXIterator end() const
Return a ZYX-order iterator that points past the maximum coordinate.
Definition: Coord.h:359
CoordBBox(CoordBBox &other, const tbb::split &)
Splitting constructor for use in TBB ranges.
Definition: Coord.h:318
tbb::atomic< Index32 > i
Definition: LeafBuffer.h:71
int Ceil(float x)
Return the ceiling of x.
Definition: Math.h:806
Vec3i asVec3i() const
Definition: Coord.h:169
const Int32 * data() const
Definition: Coord.h:163
XYZIterator endXYZ() const
Return an XYZ-order iterator that points past the maximum coordinate.
Definition: Coord.h:363
void resetToCube(const Coord &min, ValueType dim)
Definition: Coord.h:342
Axis-aligned bounding box of signed integer coordinates.
Definition: Coord.h:261
static Coord maxComponent(const Coord &lhs, const Coord &rhs)
Return the component-wise maximum of the two Coords.
Definition: Coord.h:223
Coord(const Int32 *v)
Definition: Coord.h:64
Coord offsetBy(Int32 n) const
Definition: Coord.h:119
void expand(ValueType padding)
Pad this bounding box with the specified padding.
Definition: Coord.h:419
Iterator over the Coord domain covered by a CoordBBox.
Definition: Coord.h:271
Vec3< int32_t > Vec3i
Definition: Vec3.h:675
CoordBBox expandBy(ValueType padding) const
Return a new instance that is expanded by the specified padding.
Definition: Coord.h:426
Vec3< float > Vec3s
Definition: Vec3.h:677
Coord(const Vec3I &v)
Definition: Coord.h:63
CoordBBox & operator>>=(size_t n)
Bit-wise operations performed on both the min and max members.
Definition: Coord.h:482
Coord & operator<<=(size_t n)
Definition: Coord.h:147
Coord getEnd() const
Return the maximum coordinate plus one.
Definition: Coord.h:349
Vec3d getCenter() const
Return the floating-point position of the center of this bounding box.
Definition: Coord.h:376
bool isInside(const CoordBBox &b) const
Return true if the given bounding box is inside this bounding box.
Definition: Coord.h:407
Int32 * data()
Definition: Coord.h:164
uint64_t Index64
Definition: Types.h:56
Int32 & operator[](size_t i)
Definition: Coord.h:161
void expand(const Coord &xyz)
Expand this bounding box to enclose point (x, y, z).
Definition: Coord.h:432
bool operator<=(const Coord &rhs) const
Lexicographic less than or equal to.
Definition: Coord.h:187
const Int32 * asPointer() const
Definition: Coord.h:165
Int32 & x()
Definition: Coord.h:158
#define OPENVDB_VERSION_NAME
Definition: version.h:43
Coord extents() const
Definition: Coord.h:383
void translate(const Coord &t)
Translate this bounding box by .
Definition: Coord.h:458
Coord getStart() const
Return the minimum coordinate.
Definition: Coord.h:346
CoordBBox()
The default constructor produces an empty bounding box.
Definition: Coord.h:307
bool operator==(const Iterator &other) const
Return true if this iterator and the given iterator point to the same coordinate. ...
Definition: Coord.h:285
size_t minIndex() const
Return the index (0, 1 or 2) with the smallest value.
Definition: Coord.h:238
Coord Abs(const Coord &xyz)
Definition: Coord.h:509
Coord & operator+=(const Coord &rhs)
Definition: Coord.h:121
static CoordBBox inf()
Return an "infinite" bounding box, as defined by the Coord value range.
Definition: Coord.h:332
static Coord max()
Return the largest possible coordinate.
Definition: Coord.h:70
Coord & reset(Int32 xyz)
Reset all three coordinates with the same specified argument.
Definition: Coord.h:101
void asXYZ(Int32 &x, Int32 &y, Int32 &z) const
Definition: Coord.h:171
bool isInside(const Coord &xyz) const
Return true if point (x, y, z) is inside this bounding box.
Definition: Coord.h:401
Definition: Exceptions.h:39
Coord & setX(Int32 x)
Definition: Coord.h:103
void minComponent(const Coord &other)
Perform a component-wise minimum with the other Coord.
Definition: Coord.h:199
Coord & offset(Int32 dx, Int32 dy, Int32 dz)
Definition: Coord.h:107
Int32 & z()
Definition: Coord.h:160
const Coord & max() const
Definition: Coord.h:335
void intersect(const CoordBBox &bbox)
Intersect this bounding box with the given bounding box.
Definition: Coord.h:445
Vec3< double > Vec3d
Definition: Vec3.h:678
Int32 & y()
Definition: Coord.h:159
bool operator!=(const Iterator &other) const
Return true if this iterator and the given iterator point to different coordinates.
Definition: Coord.h:290
size_t maxIndex() const
Return the index (0, 1 or 2) with the largest value.
Definition: Coord.h:241
bool operator>=(const Coord &rhs) const
Lexicographic greater than or equal to.
Definition: Coord.h:196
static CoordBBox createCube(const Coord &min, ValueType dim)
Definition: Coord.h:326
void getCornerPoints(Coord *p) const
Populates an array with the eight corner points of this bounding box.
Definition: Coord.h:464
size_t maxExtent() const
Return the index (0, 1 or 2) of the longest axis.
Definition: Coord.h:398
bool empty() const
Return true if this bounding box is empty (i.e., encloses no coordinates).
Definition: Coord.h:369
Int32 x() const
Definition: Coord.h:154
void expand(const CoordBBox &bbox)
Union this bounding box with the given bounding box.
Definition: Coord.h:439
Coord operator-() const
Definition: Coord.h:143
Coord & reset(Int32 x, Int32 y, Int32 z)
Reset all three coordinates with the specified arguments.
Definition: Coord.h:93
size_t minExtent() const
Return the index (0, 1 or 2) of the shortest axis.
Definition: Coord.h:395
Int32 * asPointer()
Definition: Coord.h:166
Coord operator-(const Coord &rhs) const
Definition: Coord.h:139
static Coord round(const Vec3< T > &xyz)
Return xyz rounded to the closest integer coordinates (cell centered conversion). ...
Definition: Coord.h:74
Vec3I asVec3I() const
Definition: Coord.h:170
bool operator<(const Coord &rhs) const
Lexicographic less than.
Definition: Coord.h:180
void reset(const Coord &min, const Coord &max)
Definition: Coord.h:341
Vec3< typename promote< T, Coord::ValueType >::type > operator-(const Coord &v1, const Vec3< T > &v0)
Allow a Coord to be subtracted from a Vec3.
Definition: Coord.h:556
Signed (x, y, z) 32-bit integer coordinates.
Definition: Coord.h:48
bool operator==(const Coord &rhs) const
Definition: Coord.h:173
Vec3s asVec3s() const
Definition: Coord.h:168
float Round(float x)
Return x rounded to the nearest integer.
Definition: Math.h:769
int32_t Int32
Definition: Types.h:59
std::ostream & operator<<(std::ostream &os, const CoordBBox &b)
Definition: Coord.h:567
static Coord minComponent(const Coord &lhs, const Coord &rhs)
Return the component-wise minimum of the two Coords.
Definition: Coord.h:215
Iterator(const CoordBBox &b)
C-tor from a bounding box.
Definition: Coord.h:275
Coord dim() const
Return the dimensions of the coordinates spanned by this bounding box.
Definition: Coord.h:381
ZYXIterator begin() const
Return a ZYX-order iterator that points to the minimum coordinate.
Definition: Coord.h:352
void expand(const Coord &min, Coord::ValueType dim)
Union this bounding box with the cubical bounding box of the given size and with the given minimum co...
Definition: Coord.h:452
CoordBBox(const Coord &min, const Coord &max)
Construct a bounding box with the given min and max bounds.
Definition: Coord.h:309
#define OPENVDB_USE_VERSION_NAMESPACE
Definition: version.h:71
int32_t Int32
Definition: Coord.h:51
Int32 ValueType
Definition: Coord.h:56
bool operator!=(const CoordBBox &rhs) const
Definition: Coord.h:366
Coord & setY(Int32 y)
Definition: Coord.h:104
void maxComponent(const Coord &other)
Perform a component-wise maximum with the other Coord.
Definition: Coord.h:207
static Coord ceil(const Vec3< T > &xyz)
Return the largest integer coordinates that are not greater than xyz+1 (node centered conversion)...
Definition: Coord.h:87
static bool lessThan(const Coord &a, const Coord &b)
Definition: Coord.h:232
CoordBBox & operator<<=(size_t n)
Bit-wise operations performed on both the min and max members.
Definition: Coord.h:481
const Coord & min() const
Definition: Coord.h:334
int Floor(float x)
Return the floor of x.
Definition: Math.h:798
void reset()
Definition: Coord.h:340
uint32_t Index32
Definition: Coord.h:52
Vec3d asVec3d() const
Definition: Coord.h:167
bool hasOverlap(const CoordBBox &b) const
Return true if the given bounding box overlaps with this bounding box.
Definition: Coord.h:413
const Coord & operator*() const
Return a const reference to the coordinate currently pointed to.
Definition: Coord.h:283
Coord::ValueType ValueType
Definition: Coord.h:265
Coord(Int32 x, Int32 y, Int32 z)
Definition: Coord.h:61
Coord & offset(Int32 n)
Definition: Coord.h:114
uint64_t Index64
Definition: Coord.h:264
size_t MaxIndex(const Vec3T &v)
Return the index [0,1,2] of the largest value in a 3D vector.
Definition: Math.h:907
ZYXIterator endZYX() const
Return a ZYX-order iterator that points past the maximum coordinate.
Definition: Coord.h:361
Coord(const Vec3i &v)
Definition: Coord.h:62
void write(std::ostream &os) const
Serialize this bounding box to the given stream.
Definition: Coord.h:492