1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
|
// Begin License:
// Copyright (C) 2006-2014 Tobias Sargeant (tobias.sargeant@gmail.com).
// All rights reserved.
//
// This file is part of the Carve CSG Library (http://carve-csg.com/)
//
// This file may be used under the terms of either the GNU General
// Public License version 2 or 3 (at your option) as published by the
// Free Software Foundation and appearing in the files LICENSE.GPL2
// and LICENSE.GPL3 included in the packaging of this file.
//
// This file is provided "AS IS" with NO WARRANTY OF ANY KIND,
// INCLUDING THE WARRANTIES OF DESIGN, MERCHANTABILITY AND FITNESS FOR
// A PARTICULAR PURPOSE.
// End:
#pragma once
#include <carve/carve.hpp>
#include <carve/math_constants.hpp>
#include <math.h>
namespace carve {
namespace geom {
template<unsigned ndim> struct vector;
}
}
namespace carve {
namespace math {
struct Matrix3;
int cubic_roots(double c3, double c2, double c1, double c0, double *roots);
void eigSolveSymmetric(const Matrix3 &m,
double &l1, carve::geom::vector<3> &e1,
double &l2, carve::geom::vector<3> &e2,
double &l3, carve::geom::vector<3> &e3);
void eigSolve(const Matrix3 &m, double &l1, double &l2, double &l3);
static inline bool ZERO(double x) { return fabs(x) < carve::EPSILON; }
static inline double radians(double deg) { return deg * M_PI / 180.0; }
static inline double degrees(double rad) { return rad * 180.0 / M_PI; }
static inline double ANG(double x) {
return (x < 0) ? x + M_TWOPI : x;
}
template<typename T>
static inline const T &clamp(const T &val, const T &min, const T &max) {
if (val < min) return min;
if (val > max) return max;
return val;
}
}
}
|