zizhou

Reordered local vertices in load_mesh / flipped tets in mesh

4

When a mesh in msh22 format is loaded, if I call `mesh->cell_vertex(c, v)` to get the cell list, the local vertices are sorted. Here're the first few cells in my...

Differentiable refactor

1

Main changes: 1. Wrap `std::vector` into a class and hide for loops over the vector inside. 2. Remove `first_derivative_unweighted` in `AdjointForm`, it may cause NAN if the weight is zero...

"Augment Lagrangian fails to solve & Wrong lumped mass matrix" for higher-order elements

1

**Describe the bug** The BCLagrangianForm uses the lumped mass matrix as the weights for Augmented Lagrangian. For linear elements, the lumped mass matrix is diagonally dominant, so the weights are...

Can run physics simulations through JSON input, and perform differentiable physics and solve inverse optimizations (https://github.com/Huangzizhou/DiffIPC-data).