Planning along Differentiable Charts of Constraint Manifolds with General-Purpose IK Solvers

MIT Computer Science and Artificial Intelligence Laboratory
Under review

*Indicates Equal Contribution

Abstract

Planning trajectories for robot manipulators under kinematic equality constraints restricts feasible motions to a measure-zero submanifold of the configuration space, requiring special algorithmic treatment. A promising strategy is parametrizing the set of feasible configurations using analytic inverse kinematics (IK). Bespoke analytic IK functions can be written to be differentiable, a necessary property for gradient-based trajectory optimization. But the vast majority of IK functions are computed by automated meta-solvers like IKFast, and are difficult to modify for differentiability. We present a new approach for computing gradients of analytic IK parameterizations: we leverage the inverse function theorem to recover the desired gradients from the ordinary forward kinematic Jacobian. Furthermore, we present a least-squares domain extension and an optimization-amenable description of the reachability constraint, which preserves gradient signal outside the reachable workspace. We demonstrate the efficacy of our approach through numerical experiments and downstream tasks, including a hardware demonstration of an RB-Y1 picking up a box and placing it on a table.

Method

Kinematic equality constraints — carrying a box with two hands, opening a door, keeping feet planted, or any policy that reasons about the robot as a floating hand — confine a robot to a measure-zero constraint manifold. Rather than carry that constraint through the optimization, we parameterize the manifold with analytic inverse kinematics: optimize over end-effector pose and self-motion, and let IK map back to joint angles. The constraint then holds by construction, leaving a positive-measure feasible set bounded only by inequalities.

This needs derivatives of the IK map, and almost every manipulator gets its IK from a meta-solver like IKFast, whose machine-generated code is hopeless to differentiate through. We never touch it. Analytic IK is the local inverse of the augmented forward kinematics \(\mathrm{FK}_A\), which is trivially differentiable, so the inverse function theorem recovers the chart's derivative from a linear solve:

\[ D\mathrm{FK}_A(q)\, \frac{\partial q}{\partial y} = \left( \frac{\partial X}{\partial y},\ \frac{\partial \psi}{\partial y} \right)^{\!\top} \]

This is the same linear algebra that differential IK (DiffIK) uses, but put to a different purpose. DiffIK resolves redundancy locally, one step at a time, and can get stuck tracking a difficult end-effector path. Because the augmented Jacobian pins the self-motion down explicitly, our charts stay globally consistent, so the optimizer can plan motions that avoid joint limits, workspace boundaries, and kinematic singularities outright.

Analytic IK is also undefined outside the reachable workspace, where optimizers routinely wander. We extend it as a least-squares projection whose exact derivative collapses to the inverse function theorem result whenever the target is reachable, and replace the usual reachability constraint — active at every feasible point — with one that only activates near the workspace boundary:

\[ -\log\det\!\left( J J^\top + \epsilon I \right) \le \tau \]

The Domain Extension maps non-reachable end-effector targets to the boundary of the workspace. The Boundary Reachability constraint stays slack in the interior and becomes active at the boundary.

How much do we lose versus bespoke autodiff?

Log-log plot of gradient error versus partial derivative vector size.
Log-log plot of gradient computation runtime versus partial derivative vector size.

Error (left) and runtime (right) versus autodiff, over 10,000 sampled reachable configurations. Axes are log-scale, shading is the 5–95% range, and outliers with higher error are marked in red.

The bimanual KUKA iiwa of Cohn et al. (ICRA 2024) has a handwritten analytic IK that was instrumented for autodiff, giving us a ground truth. The two agree to numerical precision, and our method is faster once the partial size passes \(2^6\). Downstream, our gradients paired with the old direct reachability constraint cost about 4× on the IrisNp2 convex region growth that feeds graph-of-convex-sets planning, but paired with boundary reachability and residual damping they land at 1.15 s against the bespoke baseline's 1.03 s. Boundary reachability is what makes the difference: it handles unreachable optimizer iterates robustly, while direct reachability stays active even in the interior of the feasible set.

Table I — Downstream runtimes (seconds) by gradient strategy.

Gradient MethodIrisNp2TrajoptTOPPRA
Autodiff Baselines
Autodiff, Direct Reachability1.141.540.59
Autodiff, Probing Reachability1.030.980.60
Autodiff, Boundary Reachability1.081.240.58
Proposed: IFT, Direct Reachability
IFT, Zero Gradients4.111.350.75
IFT, Pseudoinverse4.121.310.81
IFT, LM Constant3.121.420.80
IFT, LM SVT2.221.490.80
IFT, Residual Damping1.581.760.81
IFT, Anisotropic Damping1.722.520.81
IFT, Full Newton1.711.430.82
Proposed: IFT, Boundary Reachability
IFT, Zero Gradients1.641.200.75
IFT, Pseudoinverse1.751.230.80
IFT, LM Constant2.511.310.80
IFT, LM SVT1.291.210.80
IFT, Residual Damping1.151.240.81
IFT, Anisotropic Damping1.151.230.81
IFT, Full Newton1.171.200.82

Grasp selection on a UR5e with EAIK

The full grasp selection scene with a UR5e, a mug, tables, and shelf obstacles.

The grasp selection scene, with shelf obstacles.

46.5%C-space baseline
61.4%IFT, direct
63.2%IFT, boundary

Success rate from a single initialization.

A UR5e must grasp a mug among shelf obstacles, with IK derived automatically by EAIK. We generate 100 problem instances by rejection sampling collision-free goal states, and solve each from 10 randomly sampled initializations shared across all three formulations. Reformulating the problem in minimal coordinates raises SNOPT's single-initialization success rate by nearly 17 points over the configuration-space baseline; with three-initialization multistart, the parameterized formulations reach 90–93% against the baseline's 86%. The IFT solutions cost slightly more, likely because the change of variables distorts the joint-centering objective.

Table II — UR5e grasp selection IK experimental results. The number in parentheses is restricted to three-way mutual successes.

Formulation Success Rate (Multistart) Mean CostMedian CostMean RuntimeMedian Runtime
1 init3 init
C-Space Baseline46.5%86%12.68 (12.81)10.94 (11.11)0.357 (0.365)0.105 (0.096)
IFT, Direct61.4%93%13.30 (12.57)12.65 (11.95)0.719 (0.720)0.100 (0.083)
IFT, Boundary63.2%90%13.53 (13.03)12.83 (12.66)0.776 (0.825)0.110 (0.095)

Whole-body pick-and-place on RB-Y1 hardware

20/20hardware trials
56.1 smean planning time
0.555 mmmean measured violation
2.235 mmmax measured violation

A Rainbow Robotics RB-Y1 picks a box off the floor with both hands and places it on a table, using IKFast-generated IK, while keeping its center of mass over its support polytope. Grasps come from an optimization IK formulation and are ranked by a cheap heuristic plus a plannability check; the motions are planned with bidirectional RRT and shortcutting, refined by trajectory optimization, and retimed with TOPPRA. Constrained legs of the motion use the parameterization; unconstrained legs stay in C-space. We ran it from 20 box placements on a 4×5 grid spaced 3 cm apart. All 20 planned and executed. The plans satisfy the constraint to floating point tolerance, so joint position control alone was enough — what the robot actually shows is controller tracking error, averaging just over half a millimeter.

Table III — RB-Y1 box pickup experimental results.

MetricMeanMax
Planning Statistics
Total Runtime (s)56.1135.0
Optimization IK Runtime (s)19.484.0
Trajectory Optimization Runtime (s)5.629.6
Hardware Statistics
Trajectory Duration (s)54.163.6
Measured Constraint Violation (mm)0.5552.235
Measured Constraint Violation (mrad)1.5976.809

BibTeX

@article{cohn2026planning,
  title={Planning along Differentiable Charts of Constraint Manifolds with General-Purpose IK Solvers},
  author={Cohn, Thomas and Shaw, Seiji and Biggie, Harel and Manderson, Travis and Roy, Nicholas and Tedrake, Russ},
  journal={arXiv preprint arXiv:2609.10905},
  year={2026}
}