I ran some tests and it turns out the asymptotic expansion can be truncated a bit more
while retaining machine accuracy. You may (to reduce code and table size in ducc) or may not (to keep similarity with fastgl) simplify
your functions. The following is enough to maintain machine accuracy according to my tests.
NB: The "fix" commit should be squashed with the previous one.
I ran some tests and it turns out the asymptotic expansion can be truncated a bit more while retaining machine accuracy. You may (to reduce code and table size in ducc) or may not (to keep similarity with fastgl) simplify your functions. The following is enough to maintain machine accuracy according to my tests.
NB: The "fix" commit should be squashed with the previous one.