ttricco / sarracen

A Python library for smoothed particle hydrodynamics (SPH) analysis and visualization.
https://sarracen.readthedocs.io
GNU General Public License v3.0
15 stars 18 forks source link

Minimum smoothing length for fast interpolation #55

Closed ttricco closed 1 year ago

ttricco commented 1 year ago

Fast interpolation treats a pixel as just the point at the centre. Particles that don't overlap with any pixel centre can fall in the cracks between pixels, potentially leading to artefacts.

See issue #49 for discussion on this point. (Thanks @JBorrow.)

JBorrow commented 1 year ago

Thanks! FWIW, my choice for dealing with this is to treat particles below the resolution limit as if you are just doing PIC (i.e. just pixel += mass for the nearest pixel; this is super fast).