I'm trying to use the acb_calc_integrate function from the Arb C library to calculate the integral of a Fortran function through Fortran-C interoperability.
However, I'm encountering a segmentation fault when integrating a Fortran program with the polynomial function f(z) = 1 + z + z^2 + z^3 as the integrand. This issue doesn't occur when the polynomial function is f(z) = 1 + z + z^2, and the same Fortran program works fine. This issue is not present when integrating in equivalent C program.
I'm looking for assistance to understand why this segmentation fault is happening in Fortran and how to fix it. I have detailed the issue in an issue report on the Flint GitHub repository here:
https://github.com/flintlib/flint/issues/2036.
Any help would be appreciated.