-
Notifications
You must be signed in to change notification settings - Fork 406
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Extend Support for Mimic Joint #2441
base: devel
Are you sure you want to change the base?
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
👋 Hi,
This is a reminder message to assign an extra build label to this Pull Request if needed.
By default, this PR will be build with minimal build options (URDF support and Python bindings)
The possible extra labels are:
- build_collision (build Pinocchio with coal support)
- build_casadi (build Pinocchio with CasADi support)
- build_autodiff (build Pinocchio with CppAD support)
- build_codegen (build Pinocchio with CppADCodeGen support)
- build_extra (build Pinocchio with extra algorithms)
- build_mpfr (build Pinocchio with Boost.Multiprecision support)
- build_sdf (build Pinocchio with SDF parser)
- build_accelerate (build Pinocchio with APPLE Accelerate framework support)
- build_all (build Pinocchio with ALL the options stated above)
Thanks.
The Pinocchio development team.
c43fad1
to
52358ed
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I will provide a detailed review later.
3f25772
to
5ea8402
Compare
Hi, I started using mimic joints, so I'm following this PR with interest ! How/Where can I report any bug I find ? |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Maybe you could add typedef
for mimic joints ?
typedef JointModelMimicTpl<Scalar, Options, JointCollectionTpl> JointModelMimic;
typedef JointDataMimicTpl<Scalar, Options, JointCollectionTpl> JointDataMimic;
res.setIndexes(id(), idx_q(), idx_v()); | ||
m_jmodel_ref.template cast<NewScalar>(), ScalarCast<NewScalar, Scalar>::cast(m_scaling), | ||
ScalarCast<NewScalar, Scalar>::cast(m_offset)); | ||
res.setIndexes(id(), idx_q(), idx_v(), idx_j()); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Isn't a setMimicIndices
missing ?
See abussy-aldebaran@5a8f205
and abussy-aldebaran@afef317
bd466b3
to
b90f6e0
Compare
Hi, Thanks a lot for your interest in this PR, ! In order to simplify the PR and the review process, if you want to add new unittests, please do it directly on this branch. If you don't feel confortable with this solution, we'll find another way. As for bugs, please add comments here, and we'll take a look at it. Thanks again ! |
What do you mean ? Opening a PR in your forked repository ? |
ea956c8
to
083d240
Compare
geom_data.oMg[i] = | ||
(static_cast<GeometryModel::SE3>(data.oMi[joint_id]) | ||
* geom_model.geometryObjects[i].placement); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is the static_cast mandatory ? This code build without it.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Did you try with CppAd ? I removed all the other cast, and ci with cppAd, does not work
} | ||
|
||
for (JointIndex j = 1; int(j) < model.njoints; ++j) | ||
{ | ||
if (boost::get<JointModelMimicTpl<Scalar, Options, JointCollectionTpl>>(&model.joints[j])) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Why you don't check this joint ?
Add a comment to explains or use the if to jump/change invalid check for mimic.
ColsBlock Ag_cols = jmodel.jointCols(data.Ag); | ||
motionSet::inertiaAction(data.oYcrb[i], J_cols, Ag_cols); | ||
ColsBlock Ag_cols = jmodel.jointVelCols(data.Ag); | ||
motionSet::inertiaAction<ADDTO>(data.oYcrb[i], J_cols, Ag_cols); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Did you fix a bug in this pass ?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
No, the default behavior is not changed
Before this backward pass, Ag is set to zero (line 177 of this file).
The columns of the DoF that does not involve any mimicking, are selected only once in this pass and thus "added to" only once. So for those Dof, the ADDTO act the same as a SETTO.
(But for the mimic DoF this allows for summing the results)
assert( | ||
(std::is_same<JointModel, JointModelMimicTpl<Scalar, Options, JointCollectionTpl>>::value | ||
== false) | ||
&& "Algorithm not supported for mimic joints"); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This assert is nice, but how to warn the user in Python binding ?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
All the asserts were changed to PINOCCHIO_THROW, so that it works in release mode (and in python)
J_cols.transpose() * data.dFda.middleCols(jmodel.idx_v(), data.nvSubtree[i]); | ||
data.oYcrb[parent] += data.oYcrb[i]; | ||
|
||
if (ContactMode) | ||
{ | ||
jmodel.jointVelocitySelector(data.nle).noalias() = | ||
jmodel.jointVelocityFromNvSelector(data.nle).noalias() = |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Did we support constrained dynamics ?
assert(secondary_id > primary_id && "Mimicking joint id is before the primary."); | ||
|
||
size_t ancestor_prim, ancestor_sec; | ||
findCommonAncestor(model, primary_id, secondary_id, ancestor_prim, ancestor_sec); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Instead of computing this data at each call, it's maybe possible to store it in JointDataMimicTpl
.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I agree with you that computing every time seems like a waste, but I am not sure that it is worth storing it.
When JointDataMimicTpl
is instantiated, it doesn't have any available information regarding the kinematic structure of the whole model. (The JointModelMimicTpl
does not contain this information either)
… and mimiced -> mimicked for naming consistency
Following our IRL meeting, here is the todo list:
The following point are Nice To Have and should be addressed in separate issues / PRs:
|
Overview
This PR enhances support for the mimic joint.$q_{mimic} = \alpha * q_{primary} + \beta$ (where $\alpha, \beta$ are defined by the user)
Mimic joint refer to a joint which configuration depends on the configuration of another joint (that we named the primary joint). The relation is
Key Updates:
URDF Parsing:
The mimic tag can now be parsed from urdf, and add the corresponding mimic joint to the model.
A
mimic
flag has been added to the parser to enable this feature or not (default tofalse
) in order to preserve backward compatibility.Robot and Joint Model Modifications:
nq
andnv
(resp.idx_q
,idx_v
) values have now been extended with anvExtended
(respidx_vExtended
) value, holding the dimension of the extended robot/joint tangeant space. (for all classical jointnv==nvExtened
but for mimicnv==0
andnvExtended==primary.nvExtended
). The suffix vExtended stand for Jacobian as those values are mainly used for computing jacobians rows and columns.Features of the joint mimic:
Joint
Not all the joints are able to be mimicked. Please see the following list for the joints that can be mimicked:
- Prismatics
- Revolutes
- Helicoidals
Please also note that the mimic joint must be after its primary (i.e. id_mimicking > id_mimicked) in the kinematic tree for it to be created in pinocchio. If that's not the case, the model will not be created.
Algorithms
The following algorithms support the new joint mimic :
- RNEA
- CRBA
- Forward Kinematics
- Jacobians and Frames
- Centroidal Algorithm
Testing & Validation:
Limitations:
Feedback Needed:
The naming of the new member variable and functions is pretty arbitrary, we are eager to change them to better suit pinocchio naming standards.
Non exhaustively, here are a list of the lesser quality one (that should most probably be renamed) :
idx_j
/nj
: representing the dimension and index in the "extended model"jointCols
,jointRows
, etc method has been split intojointVelCols
andjointJacCols
(respectively if they refer toidx_v
/nv
oridx_j
/nj
)SimilarlyjointConfigSelector
was splitted intojointConfigFromDofSelector
andjointConfigFromNqSelector
to meet the need for mimic joint to select config either from the parent joint or don't select anythingSimilarlyjointVelocitySelector
was splitted intojointVelocityFromDofSelector
andjointConfigFromNqSelector
to meet the need for mimic joint to select config either from the parent joint or don't select anythingSpecial Thanks:
A special thank you to @EtienneAr for his guidance and advice throughout this development process.