Skip to main content

JAX-enabled continuous-time Riccati equation solver

Project description

logo

LQRax

LQRax is JAX-enabled continuous-time Riccati equation solver for LQR problems.

  • It accelerates numerical simulation through JAX's scan mechanism;
  • It allows rapid prototyping of iterative LQR (iLQR) for nonlinear control by supporting auto-differentiation on the loss function and dynamics;
  • It also enables batch-based large-scale optimal control using JAX's vmap mechanism.

This repo is currently under active development.

Install

Follow the instructions to install JAX before installing this package.

To install: pip install lqrax

Usage

There are two modules: LQR and iLQR,

The LQR module solves the following time-varying LQR problem:

$$ \arg\min_{u(t)} \int_0^T \Big[ (x(t)-x_{ref}(t))^\top Q (x(t)-x_{ref}(t)) + u(t)^\top R u(t) \Big] dt $$ $$ \text{s.t. } \dot{x}(t) = A(t) x(t) + B(t) u(t), \quad x(0) = x_0 $$

The iLQR module solves a different time-varying LQR problem:

$$ \arg\min_{v(t)} \int_0^T \Big[ z(t)^\top Q z(t) + v(t)^\top R v(t) + z(t)^\top a(t) + v(t)^\top b(t) \Big] dt $$ $$ \text{s.t. } \dot{z}(t) = A(t) z(t) + B(t) v(t), \quad z(0) = 0. $$

This formulation is often used as the sub-problem for iterative linear quadratic regulator (iLQR) to calculate the steepest descent direction on the control for a general nonlinear control problem:

$$ \arg\min_{u(t)} \int_0^T l(x(t), u(t)) dt, \text{ s.t. } \dot{x}(t) = f(x(t), u(t)), $$

where the $z(t)$ and $v(t)$ are perturbations on the system's state $x(t)$ and control $u(t)$, and $A(t)$ and $B(t)$ are the linearized system dynamics $f(x(t), u(t))$ on the current system trajectory with respect to the state and control.

Examples

Example Code Example Code
LQR for double integrator
drawing
[Jupyter Notebook]
[Google Colab]
Nonlinear control with iLQR
drawing
[Jupyter Notebook]
[Google Colab]
Multi-agent iLQGames
drawing
[Jupyter Notebook]
[Google Colab]
Ergodic control
drawing
[Jupyter Notebook]
[Google Colab]

Please also checkout Linear Quadratic Flow Matching that uses this package (the LQR module).

Copyright and License

The implementations contained herein are copyright (C) 2024 - 2025 by Max Muchen Sun, and are distributed under the terms of the GNU General Public License (GPL) version 3 (or later). Please see the LICENSE for more information.

If you use the package in your research, please cite this repository. You can see the citation information at the right side panel under "About". The BibTeX file is attached below:

@software{sun_lqrax_2025,
    author = {["Sun"], Max Muchen},
    license = {GPL-3.0},
    month = march,
    title = {{LQRax: JAX-enabled continuous-time LQR solver}},
    url = {https://github.com/MaxMSun/lqrax},
    version = {0.0.6},
    year = {2025}
}

Contact: msun@u.northwestern.edu

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

lqrax-0.0.6.tar.gz (16.9 kB view details)

Uploaded Source

Built Distribution

If you're not sure about the file name format, learn more about wheel file names.

lqrax-0.0.6-py3-none-any.whl (17.9 kB view details)

Uploaded Python 3

File details

Details for the file lqrax-0.0.6.tar.gz.

File metadata

  • Download URL: lqrax-0.0.6.tar.gz
  • Upload date:
  • Size: 16.9 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.1.0 CPython/3.10.14

File hashes

Hashes for lqrax-0.0.6.tar.gz
Algorithm Hash digest
SHA256 101e9492a8749d0122d1ebee4114137d9ace8af9b5b9771828c45d63bc204f29
MD5 5eaa21d013410acadc16c7d05be4366a
BLAKE2b-256 1d1c4f4eac1d9b8f08b6226633563efea0eba5252b4b9c2c5819a0cf6135fec8

See more details on using hashes here.

File details

Details for the file lqrax-0.0.6-py3-none-any.whl.

File metadata

  • Download URL: lqrax-0.0.6-py3-none-any.whl
  • Upload date:
  • Size: 17.9 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.1.0 CPython/3.10.14

File hashes

Hashes for lqrax-0.0.6-py3-none-any.whl
Algorithm Hash digest
SHA256 e01903414bc0669330d84031e2bbc7da714af1533dcc0cd2af7b57bb14e4b058
MD5 37e7b8818f6c3650cc103888edf0f6c2
BLAKE2b-256 72f9f15a86a278f0ae5e2e566512ac3c60384e4cf7c4fe0b02f15b28360e5d45

See more details on using hashes here.

Supported by

AWS Cloud computing and Security Sponsor Datadog Monitoring Depot Continuous Integration Fastly CDN Google Download Analytics Pingdom Monitoring Sentry Error logging StatusPage Status page