Skip to content
DraconiticResearch
Search Sources API
186,546 studies31 sourcesUpdated 3 min ago
PreprintCC BY-NC-ND 4.0arXiv

On the Compatibility of DEIM Subspaces for Nonlinear Forces Reduction in Projective Dynamics

Abstract

In physics-based simulations, the evaluation of nonlinear constraint forces represents a primary source of computational complexity. Within the framework of projective dynamics, the solver alternates between estimating these forces and updating vertex positions for the subsequent frame. This work investigates the potential of snapshot-based subspace methods to accelerate the computation of constraint terms while maintaining the fidelity of nonlinear dynamics. We focus specifically on the Discrete Empirical Interpolation Method (DEIM), examining whether it can identify a reduced set of constrained mesh elements that enable accurate interpolation and approximation of forces across the full system. Our experiments demonstrate that DEIM is neither intuitive nor consistently effective for constraint projection reduction. In particular, its performance degrades when constraints involve multiple vertices. To validate these observations, we present comparisons between full-order and reduced-order simulations on both surface meshes embedded in three dimensions and volumetric meshes. Furthermore, we report reconstruction errors for various classes of constraints, providing quantitative evidence of the method limitations in this context.

Related

  1. Covering Projective Height Balls by Subspaces in Rigid Adelic Spaces
    Ruida Di, Runjie Hu
    arXivPreprintarXiv licenceFull textAlgebra & number theoryarXiv
  2. Fermat's Last Theorem in star-invariant subspaces II: small boundary spectrum, projective plane curves, and the cubic Fermat equation
    Stephan Ramon Garcia
    arXivPreprintarXiv licenceFull textAnalysis & PDEsarXiv
  3. Equidistant Affine Subspaces: Explicit Bounds, Codimension-Two Reduction, and Transverse Products
    Jia Li
    arXivPreprintarXiv licenceFull textCombinatorics & discrete mathematicsGeometry & topologyarXiv
  4. Nonlinear Lebesgue spaces: Dense subspaces, completeness and separability
    Guillaume Sérieys, Alain Trouvé
    HALPreprintCC BY-NC-SA 4.0Full textAnalysis & PDEsHAL
  5. Admissible Subspaces for the Collective Haar Twirl: Schur-Weyl Classification and Heredity under Tensor-Factor Reduction
    Sayantan Ganguly
    arXivPreprintCC BY 4.0Full textMathematical physicsQuantum physics & quantum informationarXiv
BibTeX
@misc{abdelhafez2026compatibility,
  title = {{On the Compatibility of DEIM Subspaces for Nonlinear Forces Reduction in Projective Dynamics}},
  author = {Shaimaa Monem Abdelhafez and Igor Pontes Duff and Peter Benner},
  howpublished = {arXiv},
  year = {2026},
  month = sep,
  doi = {10.48550/arxiv.2609.40206},
  eprint = {2609.40206},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2609.40206},
  note = {CC BY-NC-ND 4.0},
}
RIS
TY  - UNPB
TI  - On the Compatibility of DEIM Subspaces for Nonlinear Forces Reduction in Projective Dynamics
AU  - Shaimaa Monem Abdelhafez
AU  - Igor Pontes Duff
AU  - Peter Benner
PY  - 2026
DA  - 2026/09/30
PB  - arXiv
DO  - 10.48550/arxiv.2609.40206
UR  - https://arxiv.org/abs/2609.40206
AB  - In physics-based simulations, the evaluation of nonlinear constraint forces represents a primary source of computational complexity. Within the framework of projective dynamics, the solver alternates between estimating these forces and updating vertex positions for the subsequent frame. This work investigates the potential of snapshot-based subspace methods to accelerate the computation of constraint terms while maintaining the fidelity of nonlinear dynamics. We focus specifically on the Discrete Empirical Interpolation Method (DEIM), examining whether it can identify a reduced set of constrained mesh elements that enable accurate interpolation and approximation of forces across the full system. Our experiments demonstrate that DEIM is neither intuitive nor consistently effective for constraint projection reduction. In particular, its performance degrades when constraints involve multiple vertices. To validate these observations, we present comparisons between full-order and reduced-order simulations on both surface meshes embedded in three dimensions and volumetric meshes. Furthermore, we report reconstruction errors for various classes of constraints, providing quantitative evidence of the method limitations in this context.
ER  -