Tools for reconstructing the velocities in a LCDM universe using the Jonker-Volgenant algorithm
Project description
void-reconstructor
------------------
This is a python program that provides the tools for reconstructing the velocities in a LCDM universe using the Jonker-Volgenant algorithm. The input are the positions of the particles in the final configuration. The LAPJV algorithm will then reconstruct the initial positions. Once the initial positions are known, we can calculate the Lagrangian and Eulerian velocities.
To use this program:
from void-reconstructor import reconstruction as rec
from void_reconstructor.volume_division import *
Then, a Reconstructor object can be created: R = rec.Reconstructor(X_final, ID, masses)
Where X_final is an array of the final positions. For example, if 1000 particles are reconstructed, the shape of X_final is (1000,3) where X_final[:,0] are all x-coordinates for the 1000 particles.
ID has shape (1000,) and contains a unique number for each particle. This can for example simply be np.range(0,1000). This makes sure that positions, velocities etc. are calculated for the correct particle.
masses contains the mass for each particle.
The grid on which the particles are reconstructed can be determined using R.get_grid
------------
Modifications
1.2.0
- Added README
- Fixed an error with mirror_positions()
1.2.2
- Fixed return of box_sizes
1.2.4
- Added a boolean to choose if the reconstruction should be done considering periodic boundary conditions, default is True. (R.periodic_condition)
1.2.6
- Added the possibility to add fake particles in full_reconstruction_nested
1.2.8
- Different selection method for homogeneous grid
- Fixed a bug with get_grid (only applicable if not a square grid)
1.2.10
- fake_particle is now N_fake_particle: N_part/N_fake_particle are added
------------------
This is a python program that provides the tools for reconstructing the velocities in a LCDM universe using the Jonker-Volgenant algorithm. The input are the positions of the particles in the final configuration. The LAPJV algorithm will then reconstruct the initial positions. Once the initial positions are known, we can calculate the Lagrangian and Eulerian velocities.
To use this program:
from void-reconstructor import reconstruction as rec
from void_reconstructor.volume_division import *
Then, a Reconstructor object can be created: R = rec.Reconstructor(X_final, ID, masses)
Where X_final is an array of the final positions. For example, if 1000 particles are reconstructed, the shape of X_final is (1000,3) where X_final[:,0] are all x-coordinates for the 1000 particles.
ID has shape (1000,) and contains a unique number for each particle. This can for example simply be np.range(0,1000). This makes sure that positions, velocities etc. are calculated for the correct particle.
masses contains the mass for each particle.
The grid on which the particles are reconstructed can be determined using R.get_grid
------------
Modifications
1.2.0
- Added README
- Fixed an error with mirror_positions()
1.2.2
- Fixed return of box_sizes
1.2.4
- Added a boolean to choose if the reconstruction should be done considering periodic boundary conditions, default is True. (R.periodic_condition)
1.2.6
- Added the possibility to add fake particles in full_reconstruction_nested
1.2.8
- Different selection method for homogeneous grid
- Fixed a bug with get_grid (only applicable if not a square grid)
1.2.10
- fake_particle is now N_fake_particle: N_part/N_fake_particle are added
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Close
Hashes for void_reconstructor-1.2.10.tar.gz
Algorithm | Hash digest | |
---|---|---|
SHA256 | ae6ac1f7d0d738b1813b32399bf664a103a323e1668f273da98247e85b4b9127 |
|
MD5 | 12212aa597832ee252d0c7ce50b96385 |
|
BLAKE2b-256 | 7b9511934797915a0ddd50c3bafb4942a69c78a90b5c9121c16c0644479cd1b1 |