31 #ifndef OPENVDB_MATH_VEC3_HAS_BEEN_INCLUDED
32 #define OPENVDB_MATH_VEC3_HAS_BEEN_INCLUDED
35 #include <openvdb/Exceptions.h>
44 template<
typename T>
class Mat3;
47 class Vec3:
public Tuple<3, T>
57 explicit Vec3(T val) { this->mm[0] = this->mm[1] = this->mm[2] = val; }
68 template <
typename Source>
77 template<
typename Source>
80 this->mm[0] =
static_cast<T
>(v[0]);
81 this->mm[1] =
static_cast<T
>(v[1]);
82 this->mm[2] =
static_cast<T
>(v[2]);
85 template<
typename Other>
88 this->mm[0] =
static_cast<T
>(v[0]);
89 this->mm[1] =
static_cast<T
>(v[1]);
90 this->mm[2] =
static_cast<T
>(v[2]);
94 T&
x() {
return this->mm[0]; }
95 T&
y() {
return this->mm[1]; }
96 T&
z() {
return this->mm[2]; }
99 T
x()
const {
return this->mm[0]; }
100 T
y()
const {
return this->mm[1]; }
101 T
z()
const {
return this->mm[2]; }
116 this->mm[0] = x; this->mm[1] = y; this->mm[2] = z;
124 this->mm[0] = 0; this->mm[1] = 0; this->mm[2] = 0;
129 template<
typename Source>
141 bool eq(
const Vec3<T> &v, T eps = static_cast<T>(1.0e-7))
const
154 template <
typename T0,
typename T1>
157 this->mm[0] = v1[0] + v2[0];
158 this->mm[1] = v1[1] + v2[1];
159 this->mm[2] = v1[2] + v2[2];
166 template <
typename T0,
typename T1>
169 this->mm[0] = v1[0] - v2[0];
170 this->mm[1] = v1[1] - v2[1];
171 this->mm[2] = v1[2] - v2[2];
178 template <
typename T0,
typename T1>
181 this->mm[0] = scale * v[0];
182 this->mm[1] = scale * v[1];
183 this->mm[2] = scale * v[2];
188 template <
typename T0,
typename T1>
191 this->mm[0] = v[0] /
scale;
192 this->mm[1] = v[1] /
scale;
193 this->mm[2] = v[2] /
scale;
202 this->mm[0]*v.mm[0] +
203 this->mm[1]*v.mm[1] +
210 return static_cast<T
>(sqrt(
double(
211 this->mm[0]*this->mm[0] +
212 this->mm[1]*this->mm[1] +
213 this->mm[2]*this->mm[2])));
222 this->mm[0]*this->mm[0] +
223 this->mm[1]*this->mm[1] +
224 this->mm[2]*this->mm[2];
230 return Vec3<T>(this->mm[1]*v.mm[2] - this->mm[2]*v.mm[1],
231 this->mm[2]*v.mm[0] - this->mm[0]*v.mm[2],
232 this->mm[0]*v.mm[1] - this->mm[1]*v.mm[0]);
241 this->mm[0] = v1.mm[1]*v2.mm[2] - v1.mm[2]*v2.mm[1];
242 this->mm[1] = v1.mm[2]*v2.mm[0] - v1.mm[0]*v2.mm[2];
243 this->mm[2] = v1.mm[0]*v2.mm[1] - v1.mm[1]*v2.mm[0];
248 template <
typename S>
251 this->mm[0] *= scalar;
252 this->mm[1] *= scalar;
253 this->mm[2] *= scalar;
258 template <
typename S>
261 this->mm[0] *= v1[0];
262 this->mm[1] *= v1[1];
263 this->mm[2] *= v1[2];
268 template <
typename S>
271 this->mm[0] /= scalar;
272 this->mm[1] /= scalar;
273 this->mm[2] /= scalar;
278 template <
typename S>
281 this->mm[0] /= v1[0];
282 this->mm[1] /= v1[1];
283 this->mm[2] /= v1[2];
288 template <
typename S>
291 this->mm[0] += scalar;
292 this->mm[1] += scalar;
293 this->mm[2] += scalar;
298 template <
typename S>
301 this->mm[0] += v1[0];
302 this->mm[1] += v1[1];
303 this->mm[2] += v1[2];
308 template <
typename S>
311 this->mm[0] -= scalar;
312 this->mm[1] -= scalar;
313 this->mm[2] -= scalar;
318 template <
typename S>
321 this->mm[0] -= v1[0];
322 this->mm[1] -= v1[1];
323 this->mm[2] -= v1[2];
363 T component(
const Vec3<T> &onto, T eps=1.0e-7)
const
368 return dot(onto)*(T(1)/l);
375 T l = onto.lengthSqr();
378 return onto*(dot(onto)*(T(1)/l));
389 if ( fabs(this->mm[0]) >= fabs(this->mm[1]) ) {
391 l = this->mm[0]*this->mm[0] + this->mm[2]*this->mm[2];
392 l =
static_cast<T
>(T(1)/sqrt(
double(l)));
393 u.mm[0] = -this->mm[2]*l;
395 u.mm[2] = +this->mm[0]*l;
398 l = this->mm[1]*this->mm[1] + this->mm[2]*this->mm[2];
399 l =
static_cast<T
>(T(1)/sqrt(
double(l)));
401 u.mm[1] = +this->mm[2]*l;
402 u.mm[2] = -this->mm[1]*l;
409 bool isNan()
const {
return isnan(this->mm[0]) || isnan(this->mm[1]) || isnan(this->mm[2]); }
412 bool isInfinite()
const
414 return isinf(this->mm[0]) || isinf(this->mm[1]) || isinf(this->mm[2]);
418 bool isFinite()
const
420 return finite(this->mm[0]) && finite(this->mm[1]) && finite(this->mm[2]);
429 template <
typename T0,
typename T1>
437 template <
typename T0,
typename T1>
441 template <
typename S,
typename T>
445 template <
typename S,
typename T>
454 template <
typename T0,
typename T1>
463 template <
typename S,
typename T>
470 template <
typename S,
typename T>
479 template <
typename T0,
typename T1>
487 template <
typename T0,
typename T1>
496 template <
typename S,
typename T>
505 template <
typename T0,
typename T1>
514 template <
typename S,
typename T>
524 template <
typename T>
528 return static_cast<T
>(atan2(c.length(), v1.dot(v2)));
531 template <
typename T>
537 template <
typename T>
548 template <
typename T>
582 template <
typename T>
592 template <
typename T>
606 #if DWREAL_IS_DOUBLE == 1
610 #endif // DWREAL_IS_DOUBLE
616 #endif // OPENVDB_MATH_VEC3_HAS_BEEN_INCLUDED