... | @@ -5,8 +5,8 @@ The ASLDA functional contains terms that may be a source of numerical instabilit |
... | @@ -5,8 +5,8 @@ The ASLDA functional contains terms that may be a source of numerical instabilit |
|
|
|
|
|
In order to obtain stable integration, one needs to set correct in `predefines.h` values of
|
|
In order to obtain stable integration, one needs to set correct in `predefines.h` values of
|
|
```c
|
|
```c
|
|
#define ASLDA_STABILIZATION_RETAIN_ABOVE_DENSITY 1.0e-5
|
|
#define SLDA_STABILIZATION_RETAIN_ABOVE_DENSITY 1.0e-5
|
|
#define ASLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY 1.0e-7
|
|
#define SLDA_STABILIZATION_EXCLUDE_BELOW_DENISTY 1.0e-7
|
|
```
|
|
```
|
|
|
|
|
|
One can do this by inspecting the initial state. In the examples below, we study the dynamics of two vortices confined in a tube. The color map presents a region where the regularization is not applied, by white we plot the region that is excluded when computing the problematic terms, while by black we plot the transition region.
|
|
One can do this by inspecting the initial state. In the examples below, we study the dynamics of two vortices confined in a tube. The color map presents a region where the regularization is not applied, by white we plot the region that is excluded when computing the problematic terms, while by black we plot the transition region.
|
... | | ... | |