Deprecated since version 0.1.
paint_old(pos, mesh, weights=1.0, mode='raise', period=None, transform=None)¶
CIC approximation (trilinear), painting points to Nmesh, each point has a weight given by weights. This does not give density. pos is supposed to be row vectors. aka for 3d input pos.shape is (?, 3).
pos[:, i] should have been normalized in the range of [ 0, mesh.shape[i] )
thus z is the fast moving index
- mode can be :
- “raise” : raise exceptions if a particle is painted
- outside the mesh
“ignore”: ignore particle contribution outside of the mesh
period can be a scalar or of length len(mesh.shape). if period is given the particles are wrapped by the period.
transform is a function that transforms pos to mesh units: transform(pos[:, 3]) -> meshpos[:, 3]
readout_old(mesh, pos, mode='raise', period=None, transform=None, out=None)¶
CIC approximation, reading out mesh values at pos, see document of paint.