| Directory: | ./ |
|---|---|
| File: | unittest/center-of-mass-derivatives.cpp |
| Date: | 2025-02-12 21:03:38 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 32 | 32 | 100.0% |
| Branches: | 84 | 166 | 50.6% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | // | ||
| 2 | // Copyright (c) 2019-2021 CNRS INRIA | ||
| 3 | // | ||
| 4 | |||
| 5 | #include "pinocchio/algorithm/center-of-mass.hpp" | ||
| 6 | #include "pinocchio/algorithm/joint-configuration.hpp" | ||
| 7 | #include "pinocchio/algorithm/compute-all-terms.hpp" | ||
| 8 | #include "pinocchio/algorithm/center-of-mass-derivatives.hpp" | ||
| 9 | #include "pinocchio/multibody/sample-models.hpp" | ||
| 10 | |||
| 11 | #include <boost/test/unit_test.hpp> | ||
| 12 | #include <boost/utility/binary.hpp> | ||
| 13 | |||
| 14 | BOOST_AUTO_TEST_SUITE(BOOST_TEST_MODULE) | ||
| 15 | |||
| 16 |
33/66✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 8 taken 1 times.
✗ Branch 9 not taken.
✓ Branch 11 taken 1 times.
✗ Branch 12 not taken.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 19 taken 1 times.
✗ Branch 20 not taken.
✓ Branch 23 taken 1 times.
✗ Branch 24 not taken.
✓ Branch 27 taken 1 times.
✗ Branch 28 not taken.
✓ Branch 30 taken 1 times.
✗ Branch 31 not taken.
✓ Branch 33 taken 1 times.
✗ Branch 34 not taken.
✓ Branch 36 taken 1 times.
✗ Branch 37 not taken.
✓ Branch 41 taken 1 times.
✗ Branch 42 not taken.
✓ Branch 45 taken 1 times.
✗ Branch 46 not taken.
✓ Branch 48 taken 1 times.
✗ Branch 49 not taken.
✓ Branch 52 taken 1 times.
✗ Branch 53 not taken.
✓ Branch 55 taken 1 times.
✗ Branch 56 not taken.
✓ Branch 58 taken 1 times.
✗ Branch 59 not taken.
✓ Branch 61 taken 1 times.
✗ Branch 62 not taken.
✓ Branch 66 taken 1 times.
✗ Branch 67 not taken.
✓ Branch 70 taken 1 times.
✗ Branch 71 not taken.
✓ Branch 73 taken 1 times.
✗ Branch 74 not taken.
✓ Branch 77 taken 1 times.
✗ Branch 78 not taken.
✓ Branch 80 taken 1 times.
✗ Branch 81 not taken.
✓ Branch 83 taken 1 times.
✗ Branch 84 not taken.
✓ Branch 86 taken 1 times.
✗ Branch 87 not taken.
✓ Branch 91 taken 1 times.
✗ Branch 92 not taken.
✓ Branch 95 taken 1 times.
✗ Branch 96 not taken.
✓ Branch 98 taken 1 times.
✗ Branch 99 not taken.
✓ Branch 102 taken 1 times.
✗ Branch 103 not taken.
✓ Branch 105 taken 1 times.
✗ Branch 106 not taken.
✓ Branch 108 taken 1 times.
✗ Branch 109 not taken.
✓ Branch 111 taken 1 times.
✗ Branch 112 not taken.
✓ Branch 116 taken 1 times.
✗ Branch 117 not taken.
|
4 | BOOST_AUTO_TEST_CASE(test_kinematics_derivatives_vcom) |
| 17 | { | ||
| 18 | using namespace Eigen; | ||
| 19 | using namespace pinocchio; | ||
| 20 | |||
| 21 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | Model model; |
| 22 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | buildModels::humanoid(model); |
| 23 | |||
| 24 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | Data data_ref(model); |
| 25 | |||
| 26 |
2/4✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
2 | model.lowerPositionLimit.head<3>().fill(-1.); |
| 27 |
2/4✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
2 | model.upperPositionLimit.head<3>().fill(1.); |
| 28 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | VectorXd q = randomConfiguration(model); |
| 29 |
2/4✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
2 | VectorXd vq(VectorXd::Random(model.nv)); |
| 30 |
2/4✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
2 | VectorXd aq(VectorXd::Random(model.nv)); |
| 31 | |||
| 32 | // Approximate dvcom_dq by finite diff. | ||
| 33 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | centerOfMass(model, data_ref, q, vq); |
| 34 |
1/2✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
|
2 | const Eigen::Vector3d vcom0 = data_ref.vcom[0]; |
| 35 | 2 | const double alpha = 1e-8; | |
| 36 |
2/4✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
2 | Eigen::VectorXd dq = VectorXd::Zero(model.nv); |
| 37 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | Data::Matrix3x dvcom_dqn(3, model.nv); |
| 38 | |||
| 39 |
2/2✓ Branch 0 taken 34 times.
✓ Branch 1 taken 1 times.
|
70 | for (int k = 0; k < model.nv; ++k) |
| 40 | { | ||
| 41 |
1/2✓ Branch 1 taken 34 times.
✗ Branch 2 not taken.
|
68 | dq[k] = alpha; |
| 42 |
2/4✓ Branch 1 taken 34 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 34 times.
✗ Branch 5 not taken.
|
68 | centerOfMass(model, data_ref, integrate(model, q, dq), vq); |
| 43 |
4/8✓ Branch 2 taken 34 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 34 times.
✗ Branch 6 not taken.
✓ Branch 8 taken 34 times.
✗ Branch 9 not taken.
✓ Branch 11 taken 34 times.
✗ Branch 12 not taken.
|
68 | dvcom_dqn.col(k) = (data_ref.vcom[0] - vcom0) / alpha; |
| 44 |
1/2✓ Branch 1 taken 34 times.
✗ Branch 2 not taken.
|
68 | dq[k] = 0; |
| 45 | } | ||
| 46 | |||
| 47 | { | ||
| 48 | // Compute dvcom_dq using the algorithm | ||
| 49 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | Data data(model); |
| 50 |
2/4✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
2 | Data::Matrix3x dvcom_dq = Data::Matrix3x::Zero(3, model.nv); |
| 51 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | centerOfMass(model, data, q, vq); |
| 52 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | getCenterOfMassVelocityDerivatives(model, data, dvcom_dq); |
| 53 | |||
| 54 | // Check that algo result and finite-diff approx are similar. | ||
| 55 |
7/14✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 11 taken 1 times.
✗ Branch 12 not taken.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 17 taken 1 times.
✗ Branch 18 not taken.
✓ Branch 20 taken 1 times.
✗ Branch 21 not taken.
✗ Branch 25 not taken.
✓ Branch 26 taken 1 times.
|
2 | BOOST_CHECK(dvcom_dq.isApprox(dvcom_dqn, sqrt(alpha))); |
| 56 | 2 | } | |
| 57 | |||
| 58 | { | ||
| 59 | // Compute dvcom_dq using the algorithm | ||
| 60 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | Data data(model); |
| 61 |
2/4✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
|
2 | Data::Matrix3x dvcom_dq = Data::Matrix3x::Zero(3, model.nv); |
| 62 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | computeAllTerms(model, data, q, vq); |
| 63 |
1/2✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
|
2 | getCenterOfMassVelocityDerivatives(model, data, dvcom_dq); |
| 64 | |||
| 65 | // Check that algo result and finite-diff approx are similar. | ||
| 66 |
7/14✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 11 taken 1 times.
✗ Branch 12 not taken.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 17 taken 1 times.
✗ Branch 18 not taken.
✓ Branch 20 taken 1 times.
✗ Branch 21 not taken.
✗ Branch 25 not taken.
✓ Branch 26 taken 1 times.
|
2 | BOOST_CHECK(dvcom_dq.isApprox(dvcom_dqn, sqrt(alpha))); |
| 67 | 2 | } | |
| 68 | 2 | } | |
| 69 | |||
| 70 | BOOST_AUTO_TEST_SUITE_END() | ||
| 71 |