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.
Note to myself: don't fix signs while listening to talks!!
Other mods/fixes:
* 2LPT is now supported also with distinct baryon/DM
power spectra. This creates now an overhead when
using Eisenstein&Hu or BBKS. But who cares.
* Reverted to determinant formulation of LLA rather than
second order expansion. This is much more accurate.
* Added baryon density filter. Baryon density field is now
consistent with CIC DM density. This is achieved by
convolving with an additional grid cell filter (and
no deconvolution).
* Tested LLA, skewness is now exact to a few per cent
but tends to be systematically higher.
* 2LPT needs a bit less memory now in multigrid mode
(FFT based 2LPT still needs humongous amounts).
* Cleaned up code output.
* The DM transfer function is now deconvolved with the cell assignment
function in order to restore power near the Nyquist wavenumber.
* Added an option to change the random number sample cube size
[random]/cubesize. This should be set >128 to generate un-
correlated numbers. True fix pending.
* Removed some unused code
* Increased Buffer size for Gadget-2 output plug-in
Fixed a bug with top grid kernel normalization, this reduces errors a bit.
A non-cubic grid is now extended to a cubic grid in the density calculation. This should reduce artefacts due to TF truncation.
Added a new (undocumented and internal) parameter that forces grid centering on the specified refinement region also in unigrid mode. This should in the future also obey that the shift comes in multiples of the top grid cells (which are however unknown).
Fixed a bug in the generic output plugin that caused non-cubic grids to be cut to cubic grids.
Fixed a bug causing wrong shift-backs of particle positions for Gadget.
Disabled optimized unigrid mode as it did not work in all circumstances (issue posted).
