markusbattarbee
commented
4 months ago I was now testing this, and the memory footprint is massive. The cause is, I think, that this project doesn't inherit the TriAxisSearch class, thus defaulting to the base project v-space initialization which in fact initializes all blocks (and because shrink_to_fit is bypassed, the footprint remains huge).
The solution is to make this inherit TriAxisSearch, but that requires that it implements a getV0() function which impacts initialization in general.
Fixes and updates to the existing 'KHB' Vlasiator project.
The updated version has a pressure balance across the simulation domain and options for a multichromatic initial perturbation with random phases.
Also added an example config file to samples/, which can be used as a quick test/demonstration.