|
|
|
For reporting basic quantities (like particle number, energy, etc) W-SLDA code used `wlog` file. It is txt file produced by function `int logger( )` from [logger.h](https://gitlab.fizyka.pw.edu.pl/gabrielw/wslda/-/blob/public/st-project-template/logger.h) file. This file can be modified by the user depending on a specifics of a given problem.
|
|
|
|
|
|
|
|
### Example
|
|
|
|
In addition to defined by default variables, the user needs to print out variable:
|
|
|
|
```math
|
|
|
|
\Delta_{\textrm{avg}}(t) = \dfrac{\int |\Delta(\bm{r},t)|(n_a(\bm{r},t)+n_b(\bm{r},t))\,d\bm{r}}{\int (n_a(\bm{r},t)+n_b(\bm{r},t))\,d\bm{r}}
|
|
|
|
```
|
|
|
|
which is (weighted) average value of delta. This can be done by adding to `logger(..)`:
|
|
|
|
```c
|
|
|
|
int ixyz;
|
|
|
|
int lNXYZ = h_potentials.blocklength;
|
|
|
|
double delta_avg=0.0, Ntot = 0.0;
|
|
|
|
for(ixyz=0; ixyz<lNXYZ; ixyz++) delta_avg+=cabs(h_potentials.delta[ixyz])*(h_densities.rho_a[ixyz]+h_densities.rho_b[ixyz]);
|
|
|
|
for(ixyz=0; ixyz<lNXYZ; ixyz++) Ntot+=h_densities.rho_a[ixyz]+h_densities.rho_b[ixyz];
|
|
|
|
delta_avg /= Ntot;
|
|
|
|
```
|
|
|
|
and supplementing `printf` statement appropriate entry, like:
|
|
|
|
```c
|
|
|
|
// add entry
|
|
|
|
fprintf(log, "%6d %18.10g %18.10g
|
|
|
|
%18.10g %18.10g %18.10g
|
|
|
|
%18.10g %18.10g %18.10g
|
|
|
|
%18.10g %18.10g %18.10g
|
|
|
|
%18.10g %18.10g %18.10g
|
|
|
|
%18.10g %18.10g %18.10g %10.2f %20s\n",
|
|
|
|
<-- HERE
|
|
|
|
it, // 1
|
|
|
|
npart[SPINA], // 2
|
|
|
|
npart[SPINB], // 3
|
|
|
|
npart[SPINA]+npart[SPINB], // 4
|
|
|
|
E_tot/Effg, // 5
|
|
|
|
energy[EKIN]/Effg, // 6
|
|
|
|
energy[EPOT]/Effg, // 7
|
|
|
|
energy[EPAIR]/Effg, // 8
|
|
|
|
energy[ECURRENT]/Effg, // 9
|
|
|
|
energy[EPOTEXT]/Effg, //10
|
|
|
|
energy[EPAIREXT]/Effg, //11
|
|
|
|
energy[EVELEXT]/Effg, //12
|
|
|
|
mu[SPINA]/eF, //13
|
|
|
|
mu[SPINB]/eF, //14
|
|
|
|
kF, // 15
|
|
|
|
eF, // 16
|
|
|
|
Effg, // 17
|
|
|
|
delta_avg, // 18 <-- HERE
|
|
|
|
logger_get_time_from_last_entry(), //19
|
|
|
|
buffer
|
|
|
|
);
|
|
|
|
|
|
|
|
``` |
|
|
\ No newline at end of file |
