| Directory: | ./ |
|---|---|
| File: | include/pinocchio/multibody/joint/joint-revolute-unbounded-unaligned.hpp |
| Date: | 2025-02-12 21:03:38 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 68 | 73 | 93.2% |
| Branches: | 73 | 158 | 46.2% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | // | ||
| 2 | // Copyright (c) 2019-2020 INRIA | ||
| 3 | // | ||
| 4 | |||
| 5 | #ifndef __pinocchio_multibody_joint_revolute_unbounded_unaligned_hpp__ | ||
| 6 | #define __pinocchio_multibody_joint_revolute_unbounded_unaligned_hpp__ | ||
| 7 | |||
| 8 | #include "pinocchio/fwd.hpp" | ||
| 9 | #include "pinocchio/spatial/inertia.hpp" | ||
| 10 | #include "pinocchio/math/rotation.hpp" | ||
| 11 | #include "pinocchio/math/matrix.hpp" | ||
| 12 | |||
| 13 | #include "pinocchio/multibody/joint/joint-revolute-unaligned.hpp" | ||
| 14 | |||
| 15 | namespace pinocchio | ||
| 16 | { | ||
| 17 | |||
| 18 | template<typename Scalar, int Options = 0> | ||
| 19 | struct JointRevoluteUnboundedUnalignedTpl; | ||
| 20 | |||
| 21 | template<typename _Scalar, int _Options> | ||
| 22 | struct traits<JointRevoluteUnboundedUnalignedTpl<_Scalar, _Options>> | ||
| 23 | { | ||
| 24 | enum | ||
| 25 | { | ||
| 26 | NQ = 2, | ||
| 27 | NV = 1 | ||
| 28 | }; | ||
| 29 | typedef _Scalar Scalar; | ||
| 30 | enum | ||
| 31 | { | ||
| 32 | Options = _Options | ||
| 33 | }; | ||
| 34 | |||
| 35 | typedef Eigen::Matrix<Scalar, NQ, 1, Options> ConfigVector_t; | ||
| 36 | typedef Eigen::Matrix<Scalar, NV, 1, Options> TangentVector_t; | ||
| 37 | |||
| 38 | typedef JointDataRevoluteUnboundedUnalignedTpl<Scalar, Options> JointDataDerived; | ||
| 39 | typedef JointModelRevoluteUnboundedUnalignedTpl<Scalar, Options> JointModelDerived; | ||
| 40 | typedef JointMotionSubspaceRevoluteUnalignedTpl<Scalar, Options> Constraint_t; | ||
| 41 | typedef SE3Tpl<Scalar, Options> Transformation_t; | ||
| 42 | typedef MotionRevoluteUnalignedTpl<Scalar, Options> Motion_t; | ||
| 43 | typedef MotionZeroTpl<Scalar, Options> Bias_t; | ||
| 44 | typedef Eigen::Matrix<Scalar, 6, NV, Options> F_t; | ||
| 45 | |||
| 46 | // [ABA] | ||
| 47 | typedef Eigen::Matrix<Scalar, 6, NV, Options> U_t; | ||
| 48 | typedef Eigen::Matrix<Scalar, NV, NV, Options> D_t; | ||
| 49 | typedef Eigen::Matrix<Scalar, 6, NV, Options> UD_t; | ||
| 50 | |||
| 51 | PINOCCHIO_JOINT_DATA_BASE_ACCESSOR_DEFAULT_RETURN_TYPE | ||
| 52 | }; | ||
| 53 | |||
| 54 | template<typename _Scalar, int _Options> | ||
| 55 | struct traits<JointDataRevoluteUnboundedUnalignedTpl<_Scalar, _Options>> | ||
| 56 | { | ||
| 57 | typedef JointRevoluteUnboundedUnalignedTpl<_Scalar, _Options> JointDerived; | ||
| 58 | typedef _Scalar Scalar; | ||
| 59 | }; | ||
| 60 | |||
| 61 | template<typename _Scalar, int _Options> | ||
| 62 | struct traits<JointModelRevoluteUnboundedUnalignedTpl<_Scalar, _Options>> | ||
| 63 | { | ||
| 64 | typedef JointRevoluteUnboundedUnalignedTpl<_Scalar, _Options> JointDerived; | ||
| 65 | typedef _Scalar Scalar; | ||
| 66 | }; | ||
| 67 | |||
| 68 | template<typename _Scalar, int _Options> | ||
| 69 | struct JointDataRevoluteUnboundedUnalignedTpl | ||
| 70 | : public JointDataBase<JointDataRevoluteUnboundedUnalignedTpl<_Scalar, _Options>> | ||
| 71 | { | ||
| 72 | EIGEN_MAKE_ALIGNED_OPERATOR_NEW | ||
| 73 | typedef JointRevoluteUnboundedUnalignedTpl<_Scalar, _Options> JointDerived; | ||
| 74 | PINOCCHIO_JOINT_DATA_TYPEDEF_TEMPLATE(JointDerived); | ||
| 75 | 837 | PINOCCHIO_JOINT_DATA_BASE_DEFAULT_ACCESSOR | |
| 76 | |||
| 77 | ConfigVector_t joint_q; | ||
| 78 | TangentVector_t joint_v; | ||
| 79 | |||
| 80 | Transformation_t M; | ||
| 81 | Constraint_t S; | ||
| 82 | Motion_t v; | ||
| 83 | Bias_t c; | ||
| 84 | |||
| 85 | // [ABA] specific data | ||
| 86 | U_t U; | ||
| 87 | D_t Dinv; | ||
| 88 | UD_t UDinv; | ||
| 89 | D_t StU; | ||
| 90 | |||
| 91 | 76 | JointDataRevoluteUnboundedUnalignedTpl() | |
| 92 |
3/5✓ Branch 1 taken 72 times.
✓ Branch 2 taken 4 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 4 times.
✗ Branch 6 not taken.
|
76 | : joint_q(Scalar(1), Scalar(0)) |
| 93 |
3/5✓ Branch 1 taken 4 times.
✓ Branch 2 taken 72 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
|
76 | , joint_v(TangentVector_t::Zero()) |
| 94 |
1/2✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
|
76 | , M(Transformation_t::Identity()) |
| 95 |
3/5✓ Branch 1 taken 4 times.
✓ Branch 2 taken 72 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
|
76 | , S(Constraint_t::Vector3::Zero()) |
| 96 |
3/6✓ Branch 1 taken 76 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 76 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
|
76 | , v(Constraint_t::Vector3::Zero(), (Scalar)0) |
| 97 |
3/5✓ Branch 1 taken 4 times.
✓ Branch 2 taken 72 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
|
76 | , U(U_t::Zero()) |
| 98 |
3/5✓ Branch 1 taken 4 times.
✓ Branch 2 taken 72 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
|
76 | , Dinv(D_t::Zero()) |
| 99 |
3/5✓ Branch 1 taken 4 times.
✓ Branch 2 taken 72 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
|
76 | , UDinv(UD_t::Zero()) |
| 100 |
3/5✓ Branch 2 taken 4 times.
✓ Branch 3 taken 72 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 4 times.
✗ Branch 6 not taken.
|
152 | , StU(D_t::Zero()) |
| 101 | { | ||
| 102 | 76 | } | |
| 103 | |||
| 104 | template<typename Vector3Like> | ||
| 105 | 16 | JointDataRevoluteUnboundedUnalignedTpl(const Eigen::MatrixBase<Vector3Like> & axis) | |
| 106 |
3/5✓ Branch 1 taken 13 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
|
16 | : joint_q(Scalar(1), Scalar(0)) |
| 107 |
3/5✓ Branch 1 taken 1 times.
✓ Branch 2 taken 13 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
16 | , joint_v(TangentVector_t::Zero()) |
| 108 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
16 | , M(Transformation_t::Identity()) |
| 109 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
16 | , S(axis) |
| 110 |
2/4✓ Branch 1 taken 14 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
16 | , v(axis, (Scalar)0) |
| 111 |
3/5✓ Branch 1 taken 1 times.
✓ Branch 2 taken 13 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
16 | , U(U_t::Zero()) |
| 112 |
3/5✓ Branch 1 taken 1 times.
✓ Branch 2 taken 13 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
16 | , Dinv(D_t::Zero()) |
| 113 |
3/5✓ Branch 1 taken 1 times.
✓ Branch 2 taken 13 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
16 | , UDinv(UD_t::Zero()) |
| 114 |
3/5✓ Branch 2 taken 1 times.
✓ Branch 3 taken 13 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
|
32 | , StU(D_t::Zero()) |
| 115 | { | ||
| 116 | 16 | } | |
| 117 | |||
| 118 | 144 | static std::string classname() | |
| 119 | { | ||
| 120 |
1/2✓ Branch 2 taken 144 times.
✗ Branch 3 not taken.
|
144 | return std::string("JointDataRevoluteUnboundedUnalignedTpl"); |
| 121 | } | ||
| 122 | 3 | std::string shortname() const | |
| 123 | { | ||
| 124 | 3 | return classname(); | |
| 125 | } | ||
| 126 | |||
| 127 | }; // struct JointDataRevoluteUnboundedUnalignedTpl | ||
| 128 | |||
| 129 | PINOCCHIO_JOINT_CAST_TYPE_SPECIALIZATION(JointModelRevoluteUnboundedUnalignedTpl); | ||
| 130 | |||
| 131 | template<typename _Scalar, int _Options> | ||
| 132 | struct JointModelRevoluteUnboundedUnalignedTpl | ||
| 133 | : public JointModelBase<JointModelRevoluteUnboundedUnalignedTpl<_Scalar, _Options>> | ||
| 134 | { | ||
| 135 | EIGEN_MAKE_ALIGNED_OPERATOR_NEW | ||
| 136 | typedef JointRevoluteUnboundedUnalignedTpl<_Scalar, _Options> JointDerived; | ||
| 137 | PINOCCHIO_JOINT_TYPEDEF_TEMPLATE(JointDerived); | ||
| 138 | typedef Eigen::Matrix<Scalar, 3, 1, Options> Vector3; | ||
| 139 | |||
| 140 | typedef JointModelBase<JointModelRevoluteUnboundedUnalignedTpl> Base; | ||
| 141 | using Base::id; | ||
| 142 | using Base::idx_q; | ||
| 143 | using Base::idx_v; | ||
| 144 | using Base::setIndexes; | ||
| 145 | |||
| 146 | 92 | JointModelRevoluteUnboundedUnalignedTpl() | |
| 147 |
1/2✓ Branch 3 taken 92 times.
✗ Branch 4 not taken.
|
92 | : axis(Vector3::UnitX()) |
| 148 | { | ||
| 149 | 92 | } | |
| 150 | |||
| 151 | JointModelRevoluteUnboundedUnalignedTpl(const Scalar & x, const Scalar & y, const Scalar & z) | ||
| 152 | : axis(x, y, z) | ||
| 153 | { | ||
| 154 | normalize(axis); | ||
| 155 | assert(isUnitary(axis) && "Rotation axis is not unitary"); | ||
| 156 | } | ||
| 157 | |||
| 158 | template<typename Vector3Like> | ||
| 159 | 33 | JointModelRevoluteUnboundedUnalignedTpl(const Eigen::MatrixBase<Vector3Like> & axis) | |
| 160 | 33 | : axis(axis) | |
| 161 | { | ||
| 162 | EIGEN_STATIC_ASSERT_VECTOR_ONLY(Vector3Like); | ||
| 163 |
5/11✓ Branch 1 taken 1 times.
✓ Branch 2 taken 20 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 21 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✓ Branch 9 taken 1 times.
✗ Branch 10 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
|
33 | assert(isUnitary(axis) && "Rotation axis is not unitary"); |
| 164 | 33 | } | |
| 165 | |||
| 166 | 16 | JointDataDerived createData() const | |
| 167 | { | ||
| 168 | 16 | return JointDataDerived(axis); | |
| 169 | } | ||
| 170 | |||
| 171 | ✗ | const std::vector<bool> hasConfigurationLimit() const | |
| 172 | { | ||
| 173 | ✗ | return {false, false}; | |
| 174 | } | ||
| 175 | |||
| 176 | ✗ | const std::vector<bool> hasConfigurationLimitInTangent() const | |
| 177 | { | ||
| 178 | ✗ | return {false}; | |
| 179 | } | ||
| 180 | |||
| 181 | using Base::isEqual; | ||
| 182 | 15 | bool isEqual(const JointModelRevoluteUnboundedUnalignedTpl & other) const | |
| 183 | { | ||
| 184 |
3/4✓ Branch 1 taken 13 times.
✓ Branch 2 taken 2 times.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
|
15 | return Base::isEqual(other) && internal::comparison_eq(axis, other.axis); |
| 185 | } | ||
| 186 | |||
| 187 | template<typename ConfigVector> | ||
| 188 | 32 | void calc(JointDataDerived & data, const typename Eigen::MatrixBase<ConfigVector> & qs) const | |
| 189 | { | ||
| 190 |
1/2✓ Branch 3 taken 22 times.
✗ Branch 4 not taken.
|
32 | data.joint_q = qs.template segment<NQ>(idx_q()); |
| 191 | |||
| 192 | 32 | const Scalar & ca = data.joint_q(0); | |
| 193 | 32 | const Scalar & sa = data.joint_q(1); | |
| 194 | |||
| 195 | 32 | toRotationMatrix(axis, ca, sa, data.M.rotation()); | |
| 196 | 32 | } | |
| 197 | |||
| 198 | template<typename TangentVector> | ||
| 199 | void | ||
| 200 | 1 | calc(JointDataDerived & data, const Blank, const typename Eigen::MatrixBase<TangentVector> & vs) | |
| 201 | const | ||
| 202 | { | ||
| 203 | 1 | data.joint_v[0] = vs[idx_v()]; | |
| 204 | 1 | data.v.angularRate() = data.joint_v[0]; | |
| 205 | 1 | } | |
| 206 | |||
| 207 | template<typename ConfigVector, typename TangentVector> | ||
| 208 | 20 | void calc( | |
| 209 | JointDataDerived & data, | ||
| 210 | const typename Eigen::MatrixBase<ConfigVector> & qs, | ||
| 211 | const typename Eigen::MatrixBase<TangentVector> & vs) const | ||
| 212 | { | ||
| 213 | 20 | calc(data, qs.derived()); | |
| 214 |
0/4✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
20 | data.joint_v[0] = vs[idx_v()]; |
| 215 | 20 | data.v.angularRate() = data.joint_v[0]; | |
| 216 | 20 | } | |
| 217 | |||
| 218 | template<typename VectorLike, typename Matrix6Like> | ||
| 219 | 5 | void calc_aba( | |
| 220 | JointDataDerived & data, | ||
| 221 | const Eigen::MatrixBase<VectorLike> & armature, | ||
| 222 | const Eigen::MatrixBase<Matrix6Like> & I, | ||
| 223 | const bool update_I) const | ||
| 224 | { | ||
| 225 |
3/6✓ Branch 2 taken 5 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 5 times.
✗ Branch 6 not taken.
✓ Branch 8 taken 5 times.
✗ Branch 9 not taken.
|
5 | data.U.noalias() = I.template middleCols<3>(Motion::ANGULAR) * axis; |
| 226 |
0/10✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
|
10 | data.Dinv[0] = |
| 227 |
3/8✗ Branch 1 not taken.
✓ Branch 2 taken 5 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 5 times.
✗ Branch 6 not taken.
✓ Branch 8 taken 5 times.
✗ Branch 9 not taken.
|
5 | Scalar(1) / (axis.dot(data.U.template segment<3>(Motion::ANGULAR)) + armature[0]); |
| 228 |
2/4✓ Branch 2 taken 5 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 5 times.
✗ Branch 6 not taken.
|
5 | data.UDinv.noalias() = data.U * data.Dinv; |
| 229 | |||
| 230 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 5 times.
|
5 | if (update_I) |
| 231 | ✗ | PINOCCHIO_EIGEN_CONST_CAST(Matrix6Like, I).noalias() -= data.UDinv * data.U.transpose(); | |
| 232 | 5 | } | |
| 233 | |||
| 234 | 166 | static std::string classname() | |
| 235 | { | ||
| 236 |
1/2✓ Branch 2 taken 166 times.
✗ Branch 3 not taken.
|
166 | return std::string("JointModelRevoluteUnboundedUnaligned"); |
| 237 | } | ||
| 238 | 24 | std::string shortname() const | |
| 239 | { | ||
| 240 | 24 | return classname(); | |
| 241 | } | ||
| 242 | |||
| 243 | /// \returns An expression of *this with the Scalar type casted to NewScalar. | ||
| 244 | template<typename NewScalar> | ||
| 245 | 5 | JointModelRevoluteUnboundedUnalignedTpl<NewScalar, Options> cast() const | |
| 246 | { | ||
| 247 | typedef JointModelRevoluteUnboundedUnalignedTpl<NewScalar, Options> ReturnType; | ||
| 248 |
1/2✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
|
5 | ReturnType res(axis.template cast<NewScalar>()); |
| 249 |
4/8✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
5 | res.setIndexes(id(), idx_q(), idx_v()); |
| 250 | 5 | return res; | |
| 251 | } | ||
| 252 | |||
| 253 | // data | ||
| 254 | /// | ||
| 255 | /// \brief axis of rotation of the joint. | ||
| 256 | /// | ||
| 257 | Vector3 axis; | ||
| 258 | }; // struct JointModelRevoluteUnboundedUnalignedTpl | ||
| 259 | |||
| 260 | } // namespace pinocchio | ||
| 261 | |||
| 262 | #include <boost/type_traits.hpp> | ||
| 263 | |||
| 264 | namespace boost | ||
| 265 | { | ||
| 266 | template<typename Scalar, int Options> | ||
| 267 | struct has_nothrow_constructor< | ||
| 268 | ::pinocchio::JointModelRevoluteUnboundedUnalignedTpl<Scalar, Options>> | ||
| 269 | : public integral_constant<bool, true> | ||
| 270 | { | ||
| 271 | }; | ||
| 272 | |||
| 273 | template<typename Scalar, int Options> | ||
| 274 | struct has_nothrow_copy<::pinocchio::JointModelRevoluteUnboundedUnalignedTpl<Scalar, Options>> | ||
| 275 | : public integral_constant<bool, true> | ||
| 276 | { | ||
| 277 | }; | ||
| 278 | |||
| 279 | template<typename Scalar, int Options> | ||
| 280 | struct has_nothrow_constructor< | ||
| 281 | ::pinocchio::JointDataRevoluteUnboundedUnalignedTpl<Scalar, Options>> | ||
| 282 | : public integral_constant<bool, true> | ||
| 283 | { | ||
| 284 | }; | ||
| 285 | |||
| 286 | template<typename Scalar, int Options> | ||
| 287 | struct has_nothrow_copy<::pinocchio::JointDataRevoluteUnboundedUnalignedTpl<Scalar, Options>> | ||
| 288 | : public integral_constant<bool, true> | ||
| 289 | { | ||
| 290 | }; | ||
| 291 | } // namespace boost | ||
| 292 | |||
| 293 | #endif // ifndef __pinocchio_multibody_joint_revolute_unbounded_unaligned_hpp__ | ||
| 294 |