GitHub

Phase-Space DTFE

Documentation for PhaseSpaceDTFE.

Density field

Installation

The PhaseSpaceDTFE package can be installed with the Julia package manager. From the Julia REPL, type ] to enter the Pkg REPL mode and run:

pkg> add PhaseSpaceDTFE

Or, equivalently, via the Pkg API:

julia> import Pkg; Pkg.add("PhaseSpaceDTFE")

Usage

Given the initial (coords_q) and final (coords_x) particle positions and velocities vels of an $N$-body simulation, we estimate the density, velocity and number of streams fields as follows:

using PhaseSpaceDTFE

m       = 1.  # particle mass
depth   = 5
sim_box = SimBox(L, Ni)

ps_dtfe_sb = ps_dtfe_subbox(coords_q, coords_x, vels, m, depth, sim_box)

Range      = 0.:0.2:100.
coords_arr = [[L/2., y, z] for y in Range, z in Range]
density_field  = density_subbox(coords_arr, ps_dtfe_sb)
nstreams_field = numberOfStreams_subbox(coords_arr, ps_dtfe_sb)
velocity_field = velocitySum_subbox(coords_arr, ps_dtfe_sb)

Please have a look at the Tutorial page for more details.

Contributors

This code was written by:

We thank:

  • Bram Alferink