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

added things

Browse source
This commit is contained in:
EiffL 2022-10-20 21:22:39 -07:00
parent ebab27b4c3
commit 70600fb7f7
2 changed files with 8 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

16
jaxpm/distributed_ops.py
View file

@ -147,24 +147,24 @@ def cic_paint(pos, mesh_shape, halo_size=0):
# Perform halo exchange
# Halo exchange along x
left = lax.pshuffle(mesh[-halo_size:],
left = lax.pshuffle(mesh[-2*halo_size:],
perm=range(mesh_size['nx'])[::-1],
axis_name='x')
right = lax.pshuffle(mesh[:halo_size],
right = lax.pshuffle(mesh[:2*halo_size],
perm=range(mesh_size['nx'])[::-1],
axis_name='x')
mesh = mesh.at[:halo_size].add(left)
mesh = mesh.at[-halo_size:].add(right)
mesh = mesh.at[:2*halo_size].add(left)
mesh = mesh.at[-2*halo_size:].add(right)
# Halo exchange along y
left = lax.pshuffle(mesh[:, -halo_size:],
left = lax.pshuffle(mesh[:, -2*halo_size:],
perm=range(mesh_size['ny'])[::-1],
axis_name='y')
right = lax.pshuffle(mesh[:, :halo_size],
right = lax.pshuffle(mesh[:, :2*halo_size],
perm=range(mesh_size['ny'])[::-1],
axis_name='y')
mesh = mesh.at[:, :halo_size].add(left)
mesh = mesh.at[:, -halo_size:].add(right)
mesh = mesh.at[:, :2*halo_size].add(left)
mesh = mesh.at[:, -2*halo_size:].add(right)
# removing halo and returning mesh
return mesh[halo_size:-halo_size, halo_size:-halo_size]

1
jaxpm/distributed_pm.py
View file

@ -64,7 +64,6 @@ def lpt(cosmo, initial_conditions, positions, a):
Computes first order LPT displacement
"""
initial_force = pm_forces(positions, delta_k=initial_conditions)
print(initial_force.shape)
a = jnp.atleast_1d(a)
dx = dops.scalar_multiply(initial_force, growth_factor(cosmo, a))
p = dops.scalar_multiply(dx, a**2 * growth_rate(cosmo, a) *