Adapt ITPA backgrounds to new kinetic init
Bug description / feature request:
In struphy.kinetic_backgrounds.maxwellians
, the class Maxwellian6DITPA
can be eliminated. Instead, we can add the density
n = c3*np.exp(-c2/c1*np.tanh((eta1 - c0)/c2))
as a new class to the module struphy.initial.perturbations.py
.
Moreover, the class Maxwellian5DUniform
should be renamed to Maxwellian5D
and function in an analogous way to Maxwellian6D
. Thus the class Maxwellian5DITPA
can be eliminated as well.
Expected behavior:
5D kinetic backgrounds should function in analogy to 6D backgrounds.
Proposed solution:
See above.