python-unsiotools 1.0.1

Python wrapper to unsiotools

ABOUT

UNSIOTOOLS contains a collection of unsio based programs and depends of the Nemo package. You have access to falcON algorithm to compute density and gravity.

Installing python wrapper

pip install python-unsiotools -U

To get some docstring help

# Help on falcON class
pydoc unsiotools.simulations.cfalcon

Usage

• In the example below, we load a RAMSES simulation, and compute density and gravity on HALO particles

import unsio.input as uns_in  # unsio reading module
import unsiotools.simulations.cfalcon as falcon

myfile="/home/jcl/output_00004" # input RAMSES simulation
# we instantiate a CUNS_IN object
my_in=uns_in.CUNS_IN(myfile,"halo") # We select components HALO
#
# Reading
#
if my_in.nextFrame(): # load snapshot
  # read halo positions
  status,pos=my_in.getData("halo","pos")
  # read halo mass
  status,mass=my_in.getData("halo","mass")
  # read time simulation
  status,timex=my_in.getData("time")

# compute density
cf=falcon.CFalcon()
ok,rho,hsml=cf.getDensity(pos,mass)
print("Rho=",rho)

# compute gravity
eps=0.05 # softening
ok,acc,phi=cf.getGravity(pos,mass,eps)
print (ok,acc,phi)

• In this more simple example, we compute density and gravity on random data. Note that data must be in float32 format

import unsiotools.simulations.cfalcon as falcon
import numpy as np

pos=np.float32(np.random.random_sample((300,)))
mass=np.float32(np.random.random_sample((100,)))

cf=falcon.CFalcon()
# density
ok,rho,hsml=cf.getDensity(pos,mass)
print(f"ok={ok}\nrho={rho}")
# gravity
eps=0.05
ok,acc,phi=cf.getGravity(pos,mass,eps)
print (f"ok={ok}\nacc={acc}\nphi={phi}")

Licence

UNSIOTOOLS is open source and released under the terms of the CeCILL2 Licence

Webpage

PLease visit :

• UNSIO project home page
• UNSIO Python reading manual
• UNSIO Python writing manual
• UNSIO Pypi page
• NEMO home page

Copyright

Copyright Jean-Charles LAMBERT
Jean-Charles.Lambert_at_lam.fr