Commit Graph

2 Commits

Author SHA1 Message Date
Michele Guerini Rocco
3a2dd3697f
src/gray_plasma.f90: use mollifier to interpolate density
This methods ensures n_e(ψ) is everywhere positive and C²-continuous.

The density boundary is now dependent on the width w of the mollifier,
which also controls the smoothness of the curve, specifically:
  ψ_bnd = ψ_last + w
where ψ_last = 1.01 currently is an extra point added to improve the
convergence of the mollified density to the original data as w→0.
2025-02-21 21:37:30 +01:00
Michele Guerini Rocco
38a8edd439
doc: reorganize example input files 2024-02-16 00:38:43 +01:00