Matrix-Free Delassus Operations: Scalable and Memory-Efficient Algorithms
Preprint
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 resourcelimited 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 2x to over 400x, in contact-rich scenarios.Keywords
All Related Versions
This publication has 18 references indexed in Scilit:
- Constrained Articulated Body Algorithms for Closed-Loop MechanismsIEEE Transactions on Robotics, 2026
- Constrained Articulated Body Dynamics AlgorithmsIEEE Transactions on Robotics, 2024
- PV-OSIMr: A Lowest Order Complexity Algorithm for Computing the Delassus MatrixIEEE Robotics and Automation Letters, 2024
- Efficient Constrained Dynamics Algorithms Based on an Equivalent LQR Formulation Using Gauss' Principle of Least ConstraintIEEE Transactions on Robotics, 2023
- 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
- Realistic haptic rendering of interacting deformable objects in virtual environmentsIEEE Transactions on Visualization and Computer Graphics, 2005
- Efficient recursive algorithm for the operational space inertia matrix of branching mechanismsAdvanced Robotics, 2001
- A Very Efficient Algorithm for the Simulation of Robots and Similar Multibody Systems without Inversion of the Mass MatrixIFAC Proceedings Volumes, 1986
- On-Line Computational Scheme for Mechanical ManipulatorsJournal of Dynamic Systems, Measurement, and Control, 1980