on of the high resolution region.
added output of comoving size of high-res region to screen
added output of densest point also when using the generic output plugin
via option [random]/restart = yes/no
if set to yes, it will try to load the data, but all files have to exist,
and the refinement hierarchy cannot be changed.
also the 'LINGER++' one was present (they used the same class name, so which one
was actually used depended on the compile order rather than the user choice).
* Added new transfer function type 'music' which reads in a file with
density and velocity transfer functions for baryons and CDM.
* Powerspectrum evolution in two-fluid sim appears to scale as
predicted now at the per cent level
* Tool to compute velocity transfer functions is not included as of
now. Let's see..
* Velocity potentials are not yet supported with 2LPT.
* Fixed a bug in the coarse grid density correction that lead to a shift in the correction in some rare cases
* Fixed a bug in restricting fine grid noise when the double size fine grid (for FFT) is actually larger than the domain.
* Added new log file routines that use printf syntax
* Added log file output to most parts of the code
* Removed the stats file output, it is now appended to the log file
* Fixed a normalisation bug in the hybrid Poisson solver
* Added calculation of the exact T(r=0) by summing modes rather than by analytical integration
* Fixed a bug in the deconvolution that caused the Nyquist mode to be excluded in the deconv.
* Fixed a bug in the random number generator when using single thereaded FFTW.
New conservative transfer function scheme.
New random number generator that preserves coarse Fourier modes instead of applying Hoffman-Ribak
Some minor bugs fixed.
* Moved hybrid poisson solver from convolution_kernel.cc to poisson.cc
* Fine-tuned deconvolution schemes for various running modes
* Added lots of doxygen documentation
improves the solution.
Also fixed another bug related to the previous grid placement correction with align_top,
size corrected grids have to be enforced to have an even offset.
LLA baryon density might deviate from the mean density by up to a few percent.
This is now corrected explcitly by subtracting the mean after calculation
Grid sizes have to be a multiple of 4 during density convolution when using
isolated boundary conditions. Otherwise the padding region is not coextensive
with coarse grid cells. This is now enforced during density generation but
is readjusted afterwards.