31 #ifndef OPENVDB_MATH_VEC4_HAS_BEEN_INCLUDED
32 #define OPENVDB_MATH_VEC4_HAS_BEEN_INCLUDED
35 #include <openvdb/Exceptions.h>
46 template<
typename T>
class Mat3;
49 class Vec4:
public Tuple<4, T>
59 explicit Vec4(T val) { this->mm[0] = this->mm[1] = this->mm[2] = this->mm[3] = val; }
71 template <
typename Source>
81 template<
typename Source>
84 this->mm[0] =
static_cast<T
>(v[0]);
85 this->mm[1] =
static_cast<T
>(v[1]);
86 this->mm[2] =
static_cast<T
>(v[2]);
87 this->mm[3] =
static_cast<T
>(v[3]);
91 T&
x() {
return this->mm[0]; }
92 T&
y() {
return this->mm[1]; }
93 T&
z() {
return this->mm[2]; }
94 T&
w() {
return this->mm[3]; }
97 T
x()
const {
return this->mm[0]; }
98 T
y()
const {
return this->mm[1]; }
99 T
z()
const {
return this->mm[2]; }
100 T
w()
const {
return this->mm[3]; }
118 this->mm[0] = x; this->mm[1] = y; this->mm[2] = z; this->mm[3] = w;
125 this->mm[0] = 0; this->mm[1] = 0; this->mm[2] = 0; this->mm[3] = 0;
130 template<
typename Source>
164 template <
typename T0,
typename T1>
167 this->mm[0] = v1[0] + v2[0];
168 this->mm[1] = v1[1] + v2[1];
169 this->mm[2] = v1[2] + v2[2];
170 this->mm[3] = v1[3] + v2[3];
178 template <
typename T0,
typename T1>
181 this->mm[0] = v1[0] - v2[0];
182 this->mm[1] = v1[1] - v2[1];
183 this->mm[2] = v1[2] - v2[2];
184 this->mm[3] = v1[3] - v2[3];
191 template <
typename T0,
typename T1>
194 this->mm[0] = scale * v[0];
195 this->mm[1] = scale * v[1];
196 this->mm[2] = scale * v[2];
197 this->mm[3] = scale * v[3];
202 template <
typename T0,
typename T1>
205 this->mm[0] = v[0] / scalar;
206 this->mm[1] = v[1] / scalar;
207 this->mm[2] = v[2] / scalar;
208 this->mm[3] = v[3] / scalar;
216 return (this->mm[0]*v.mm[0] + this->mm[1]*v.mm[1]
217 + this->mm[2]*v.mm[2] + this->mm[3]*v.mm[3]);
224 this->mm[0]*this->mm[0] +
225 this->mm[1]*this->mm[1] +
226 this->mm[2]*this->mm[2] +
227 this->mm[3]*this->mm[3]);
235 return (this->mm[0]*this->mm[0] + this->mm[1]*this->mm[1]
236 + this->mm[2]*this->mm[2] + this->mm[3]*this->mm[3]);
268 template <
typename S>
271 this->mm[0] *= scalar;
272 this->mm[1] *= scalar;
273 this->mm[2] *= scalar;
274 this->mm[3] *= scalar;
279 template <
typename S>
282 this->mm[0] *= v1[0];
283 this->mm[1] *= v1[1];
284 this->mm[2] *= v1[2];
285 this->mm[3] *= v1[3];
291 template <
typename S>
294 this->mm[0] /= scalar;
295 this->mm[1] /= scalar;
296 this->mm[2] /= scalar;
297 this->mm[3] /= scalar;
302 template <
typename S>
305 this->mm[0] /= v1[0];
306 this->mm[1] /= v1[1];
307 this->mm[2] /= v1[2];
308 this->mm[3] /= v1[3];
313 template <
typename S>
316 this->mm[0] += scalar;
317 this->mm[1] += scalar;
318 this->mm[2] += scalar;
319 this->mm[3] += scalar;
324 template <
typename S>
327 this->mm[0] += v1[0];
328 this->mm[1] += v1[1];
329 this->mm[2] += v1[2];
330 this->mm[3] += v1[3];
335 template <
typename S>
338 this->mm[0] -= scalar;
339 this->mm[1] -= scalar;
340 this->mm[2] -= scalar;
341 this->mm[3] -= scalar;
346 template <
typename S>
349 this->mm[0] -= v1[0];
350 this->mm[1] -= v1[1];
351 this->mm[2] -= v1[2];
352 this->mm[3] -= v1[3];
364 return isnan(this->mm[0]) || isnan(this->mm[1])
365 || isnan(this->mm[2]) || isnan(this->mm[3]);
371 return isinf(this->mm[0]) || isinf(this->mm[1])
372 || isinf(this->mm[2]) || isinf(this->mm[3]);
378 return finite(this->mm[0]) && finite(this->mm[1])
379 && finite(this->mm[2]) && finite(this->mm[3]);
388 template <
typename T0,
typename T1>
399 template <
typename T0,
typename T1>
403 template <
typename S,
typename T>
408 template <
typename S,
typename T>
417 template <
typename T0,
typename T1>
429 template <
typename S,
typename T>
439 template <
typename S,
typename T>
448 template <
typename T0,
typename T1>
453 result(v0[0]/v1[0], v0[1]/v1[1], v0[2]/v1[2], v0[3]/v1[3]);
458 template <
typename T0,
typename T1>
467 template <
typename S,
typename T>
476 template <
typename T0,
typename T1>
485 template <
typename S,
typename T>
499 template <
typename T>
510 template <
typename T>
526 #if DWREAL_IS_DOUBLE == 1
530 #endif // DWREAL_IS_DOUBLE
536 #endif // OPENVDB_MATH_VEC4_HAS_BEEN_INCLUDED