Two squares falling onto the ground under gravity, contacting with each other, and then compressed by a ceiling is simulated with an inversion-free hyperelastic potential and IPC with implicit Euler time integration. Each time step is solved by minimizing the Incremental Potential with the projected Newton method.

pip install numpy scipy pygame

python simulator.py