Skip to content

JAX-FEM in the Mosaic benchmark suite #92

Description

@dionhaefner

We just released Mosaic (github.com/pasteurlabs/mosaic), a benchmark for differentiable physics solvers. Each solver runs in a container behind a uniform apply / vjp interface, so one harness compares gradient accuracy and cost across JAX, PyTorch, Julia, and C++-adjoint codes on shared tasks.

JAX-FEM is one of the benchmarked solvers on the structural-mesh and thermal-mesh domains (results if you're interested).

We did our best to do it justice, but we don't know your code as well as you do and have surely made suboptimal calls. For example, we solve the (linear) elasticity and heat problems with the default direct UMFPACK solver and take the implicit-function-theorem adjoint JAX-FEM exposes via ad_wrapper. If a different linear solver, or a different way of invoking the adjoint, is what you'd recommend, we'd love to know.

The setup is self-contained under mosaic/tesseracts/structural-mesh/jax-fem/ and mosaic/tesseracts/thermal-mesh/jax-fem/. If something looks odd, corrections or a PR are very welcome, here on GitHub or on the Tesseract Forum.

Mostly we wanted you to know your work is represented, and to get it right.

Cheers,
Dion on behalf of the Mosaic authors

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions