| ... | ... | @@ -18,8 +18,8 @@ The aim is to retain the intact region where a physics of interest happens and r |
|
|
|
# Testing the stabilizator
|
|
|
|
In the plot below we present the total energy evolution as we run the simulation. The curves correspond to different settings of the stabilization parameters. It is clearly seen that the energy is conserved once the parameters are set correctly, in this case:
|
|
|
|
```c
|
|
|
|
#define ASLDA_STABILIZATION_RETAIN_ABOVE_DENSITY 5.0e-4
|
|
|
|
#define ASLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY 5.0e-6
|
|
|
|
#define SLDA_STABILIZATION_RETAIN_ABOVE_DENSITY 5.0e-4
|
|
|
|
#define SLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY 5.0e-6
|
|
|
|
```
|
|
|
|
Otherwise, we may get unstable integration after some time, which typically is manifested as non-conservation of the energy or particle number.
|
|
|
|
 |
