Aquila-Consortium/JaxPM_highres
0
0
Fork 1
forked from guilhem_lavaux/JaxPM
Code Pull requests Activity

Update jaxpm/pm.py

Browse source 
Co-authored-by: Francois Lanusse <EiffL@users.noreply.github.com>
This commit is contained in:
Denise Lanzieri 2022-05-18 09:59:59 +02:00 • committed by GitHub
parent dd71d359e0
commit 3df0c05d29
1 changed files with 3 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
jaxpm/pm.py
View file

@ -75,6 +75,9 @@ def make_ode_fn(mesh_shape):
def pgd_correction(pos, params):
"""
improve the short-range interactions of PM-Nbody simulations with potential gradient descent method, based on https://arxiv.org/abs/1804.00671
args:
pos: particle positions [npart, 3]
params: [alpha, kl, ks] pgd parameters
"""
kvec = fftk(mesh_shape)