Antkeeper  0.0.1
euler-ik-constraint.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2023 Christopher J. Howard
3  *
4  * This file is part of Antkeeper source code.
5  *
6  * Antkeeper source code is free software: you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation, either version 3 of the License, or
9  * (at your option) any later version.
10  *
11  * Antkeeper source code is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with Antkeeper source code. If not, see <http://www.gnu.org/licenses/>.
18  */
19 
21 #include <algorithm>
22 #include <cmath>
23 
25 {
26  // Store w-component of quaternion
27  float old_w = q.w();
28 
29  // Derive Euler angles from quaternion
30  auto angles = math::euler_from_quat(m_rotation_sequence, q);
31 
32  // Constrain Euler angles
33  angles = math::clamp(angles, m_min_angles, m_max_angles);
34 
35  // Rebuild quaternion from constrained Euler angles
36  q = math::euler_to_quat(m_rotation_sequence, angles);
37 
38  // Restore quaternion sign
39  q.w() = std::copysign(q.w(), old_w);
40 }
void solve(math::fquat &q) override
Solves the constraint.
vec3< T > euler_from_quat(const quat< T > &q, T tolerance=T{1e-6})
Derives Euler angles from a unit quaternion.
quat< T > euler_to_quat(const vec3< T > &angles)
Constructs a quaternion from Euler angles.
constexpr vector< T, N > clamp(const vector< T, N > &x, const vector< T, N > &min, const vector< T, N > &max)
Clamps the values of a vector's elements.
Definition: vector.hpp:1069
Quaternion composed of a real scalar part and imaginary vector part.
Definition: quaternion.hpp:39