Skip to content

GitLab

  • Menu
Projects Groups Snippets
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in
  • wslda wslda
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 0
    • Issues 0
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 0
    • Merge requests 0
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Monitor
    • Monitor
    • Incidents
  • Packages & Registries
    • Packages & Registries
    • Package Registry
    • Container Registry
    • Infrastructure Registry
  • Analytics
    • Analytics
    • CI/CD
    • Repository
    • Value stream
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • wtools
  • wsldawslda
  • Wiki
  • Monitoring of conservation laws

Monitoring of conservation laws · Changes

Page history
Update Monitoring of conservation laws authored Feb 20, 2026 by Gabriel Wlazłowski's avatar Gabriel Wlazłowski
Hide whitespace changes
Inline Side-by-side
Monitoring-of-conservation-laws.md
View page @ 86995442
**VERSION>=2022.04.05**
The W-SLDA Toolkit uses conserved quantities to monitor the stability of time integration. The quantities of special interest are:
* _Particle number_ ($`N_a+N_b`$): this quantity is conserved almost for all cases. The special case when the total particle number may not be conserved is the situation when an external pairing field is activated.
* _Particle number_ ($`N_a+N_b`$): this quantity is conserved for almost all cases. The special case in which the total particle number may not be conserved arises when an external pairing field is applied.
* _Total energy_: if there is no time-dependent external potential, the energy is conserved.
To avoid situations where the code uses computing resources for generating trajectories that are incorrect (violation of conservation laws) there is build in mechanism for monitoring conservation laws. This option is controlled by the following tags in the input file:
......@@ -17,7 +17,7 @@ To avoid situations where the code uses computing resources for generating traje
# Nconservation_stop 10000 # Stop to monitor total particle number conservation at time*eF, default=1e12 (infinity)
# Nconservation_tol 0.05 # if |[N(t)-N(start_t)]/N(start_t)|>tol then the code will terminate, default=0.05
```
Note that the energy conservation checking is disabled by default while particle number conservation monitoring is activated. If the violation (measured as deviation from the initial value) is detected, the code terminates. The information will be displays in `stdout` file:
Note that the energy conservation checking is disabled by default, while particle number conservation monitoring is activated. If the violation (measured as deviation from the initial value) is detected, the code terminates. The information will be displayed in the `stdout` file:
```
# SIMULATION INSTABILITY CRITERIA MET FOR PARTICLE NUMBER!!!
# INITIAL VALUE=19.000000, PRESENT VALUE=20.351277, RELATIVE CHANGE=0.07112 [> 0.050000]
......@@ -29,7 +29,7 @@ If the checkpoint tag is selected
```bash
checkpoint 1 # do checkpoint at the end of the calculation?
```
then immediately after terminating the code, the checkpoint files will be created
then, immediately after terminating the code, the checkpoint files will be created
```
# SIMULATION INSTABILITY CRITERIA MET FOR PARTICLE NUMBER!!!
# INITIAL VALUE=19.000000, PRESENT VALUE=20.351277, RELATIVE CHANGE=0.07112 [> 0.050000]
......
Clone repository

Content of Documentation
Official webpage
W-BSK Toolkit