Matrix-Free Delassus Operations: Scalable and Memory-Efficient Algorithms
- 8 June 2026
- journal article
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Robotics and Automation Letters
- Vol. 11 (8) , 9375-9382
Abstract
The Delassus matrix, closely related to the operational-space inertia matrix, is a fundamental quantity in robotics with applications in simulation, system identification, and control. Traditional approaches compute and store this matrix explicitly, either in sparse or dense form. In this work, we depart from this convention by treating the Delassus matrix as a matrix-free operator. We derive efficient algorithms with low computational complexity that multiply the Delassus matrix or its damped inverse by an input vector or matrix. Unlike approaches based on explicit matrix construction, our method achieves a linear memory footprint, making it scalable to problems with thousands of constraints and suitable for execution on resource-limited hardware. We implement these matrix-free operations on top of the open-source Pinocchio library and evaluate their performance against state-of-the-art methods that rely on explicit matrix computation. Our benchmarks demonstrate substantial speedups, ranging from approximately $2\times$ to $4\times$ for single-robot scenarios, up to $10\times$ when joint friction is enabled, and over $400\times$ in the Pyramid+Robots contact-rich multi-body scene.Keywords
All Related Versions
Funding Information
- European Union’s Horizon Europe research and innovation programme
- Marie Skłodowska-Curie Postdoctoral Fellowship (101211945)
- ARTIFACT project (101165695)
- ARTIFACT project (101070165)
- French Government under the management of Agence Nationale de la Recherche (ANR-22-CE33-0007-01)
- PR[AI]RIE-PSAI” AI Cluster (ANR-23-IACL-0008)
- France 2030 Program with the references Organic Robotics Program
- PIQ Program
- Management of Agence de Programme du Numérique
This publication has 28 references indexed in Scilit:
- Convex and analytically-invertible dynamics with contacts and constraints: Theory and implementation in MuJoCo2014
- MuJoCo: A physics engine for model-based control2012
- A reduced-order recursive algorithm for the computation of the operational-space inertia matrix2012
- A matrix-free cone complementarity approach for solving large-scale, nonsmooth, rigid body dynamicsComputer Methods in Applied Mechanics and Engineering, 2011
- Rigid Body Dynamics Algorithms2008
- Nonsmooth Mechanics1999
- On the limited memory BFGS method for large scale optimizationMathematical Programming, 1989
- A unified approach for motion and force control of robot manipulators: The operational space formulationIEEE Journal on Robotics and Automation, 1987
- A Very Efficient Algorithm for the Simulation of Robots and Similar Multibody Systems without Inversion of the Mass MatrixIFAC Proceedings Volumes, 1986
- The conjugate-gradient method for solving linear systems1956