We recently [encountered](https://github.com/cp2k/cp2k/pull/3577#issuecomment-2261386023) a floating point exception in libgrpp. I believe the problem occurs when [this branch](https://github.com/aoleynichenko/libgrpp/blob/6e63e88f75385b811837efbbb143b2dcd83b00e8/libgrpp/specfunc_scaled_mod_sph_bessel.c#L83) returns [zero](https://github.com/aoleynichenko/libgrpp/blob/6e63e88f75385b811837efbbb143b2dcd83b00e8/libgrpp/specfunc_scaled_mod_sph_bessel.c#L92), which then leads to a [division by zero](https://github.com/aoleynichenko/libgrpp/blob/6e63e88f75385b811837efbbb143b2dcd83b00e8/libgrpp/screening.c#L356).
We recently encountered a floating point exception in libgrpp.
I believe the problem occurs when this branch returns zero, which then leads to a division by zero.