joint_acceleration_backward_difference.cpp
Go to the documentation of this file.
1 //
2 // Copyright (c) 2018, University of Edinburgh
3 // All rights reserved.
4 //
5 // Redistribution and use in source and binary forms, with or without
6 // modification, are permitted provided that the following conditions are met:
7 //
8 // * Redistributions of source code must retain the above copyright notice,
9 // this list of conditions and the following disclaimer.
10 // * Redistributions in binary form must reproduce the above copyright
11 // notice, this list of conditions and the following disclaimer in the
12 // documentation and/or other materials provided with the distribution.
13 // * Neither the name of nor the names of its contributors may be used to
14 // endorse or promote products derived from this software without specific
15 // prior written permission.
16 //
17 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
18 // AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
19 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
20 // ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
21 // LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
22 // CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
23 // SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
24 // INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
25 // CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
26 // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
27 // POSSIBILITY OF SUCH DAMAGE.
28 //
29 
31 
32 REGISTER_TASKMAP_TYPE("JointAccelerationBackwardDifference", exotica::JointAccelerationBackwardDifference);
33 
34 namespace exotica
35 {
37 {
38  scene_ = scene;
39 
40  // Get ndof
41  N_ = scene_->GetKinematicTree().GetNumControlledJoints();
42 
43  // Set binomial coefficient parameters
45 
46  // Init each col of q_ with start state
47  q_.resize(N_, 2);
48  if (parameters_.StartState.rows() == 0)
49  {
50  q_.setZero(N_, 2);
51  }
52  else if (parameters_.StartState.rows() == N_)
53  {
54  for (int i = 0; i < 2; ++i)
55  q_.col(i) = parameters_.StartState;
56  }
57  else
58  {
59  ThrowPretty("Wrong size for StartState!");
60  }
61 
62  // Init qbd_
64 
65  // Init identity matrix
66  I_ = Eigen::MatrixXd::Identity(N_, N_);
67 }
68 
70 {
71  // Input check
72  if (joint_state.rows() != N_) ThrowNamed("Wrong size for joint_state!");
73 
74  // Push back previous joint states
75  q_.col(1) = q_.col(0);
76  q_.col(0) = joint_state;
77 
78  // Compute new qbd_
80 }
81 
83 {
84  // Input check
85  if (phi.rows() != N_) ThrowNamed("Wrong size of phi!");
86 
87  // Estimate second time derivative
88  phi = x + qbd_;
89 }
90 
92 {
93  // Input check
94  if (jacobian.rows() != N_ || jacobian.cols() != N_) ThrowNamed("Wrong size of jacobian! " << N_);
95 
96  // Estimate second time derivative and set Jacobian to identity matrix
97  Update(x, phi);
98  jacobian = I_;
99 }
100 
102 {
103  return N_;
104 }
105 } // namespace exotica
Eigen::MatrixXd q_
Log of previous two joint states.
Eigen::Ref< Eigen::VectorXd > VectorXdRef
Eigen::Vector2d backward_difference_params_
Binomial cooeficient parameters.
#define ThrowPretty(m)
std::shared_ptr< Scene > ScenePtr
#define ThrowNamed(m)
JointAccelerationBackwardDifferenceInitializer parameters_
const Eigen::Ref< const Eigen::VectorXd > & VectorXdRefConst
REGISTER_TASKMAP_TYPE("JointAccelerationBackwardDifference", exotica::JointAccelerationBackwardDifference)
void Update(Eigen::VectorXdRefConst x, Eigen::VectorXdRef phi) override
void SetPreviousJointState(Eigen::VectorXdRefConst joint_state)
Logs previous joint state. SetPreviousJointState must be called after solve is called in a Python/C++...
Eigen::Ref< Eigen::MatrixXd > MatrixXdRef
Eigen::VectorXd qbd_
x+qbd_ is a simplifed estimate of the second time derivative.
Time-derivative estimation by backward differencing. JointAccelerationBackwardDifference uses backwar...


exotica_core_task_maps
Author(s): Yiming Yang, Michael Camilleri
autogenerated on Sat Apr 10 2021 02:36:09