Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in
  • SWIFTsim SWIFTsim
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 57
    • Issues 57
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 21
    • Merge requests 21
  • Deployments
    • Deployments
    • Releases
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • SWIFT
  • SWIFTsimSWIFTsim
  • Issues
  • #642
Closed
Open
Issue created Dec 27, 2019 by Matthieu Schaller@matthieuOwner

Crash in planetary SESAME equation of state

The function SESAME_pressure_from_internal_energy() can return infinite values which then propagate in the code and lead to a crash (eventually by having particles with NaN positions as a result of infinite accelerations). In some cases the denominator in the equation for intp_u_2 turns out to be zero, meaning there are two consecutive values in table_log_u_rho_T which are equal. For instance the entries 59241 and 59242.

If the tables are correct then that function should be made robust to this. Or something needs to be done to prevent crashes down the line.

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information
Assignee
Assign to
Time tracking