search

jacobwilliams / quadpack

Modern Fortran QUADPACK Library for 1D numerical quadrature
Other
62 stars 8 forks source link

Bug in dqawo #1

Closed jacobwilliams closed 2 years ago

jacobwilliams commented 2 years ago

I believe this line:

        if (Ier == 6) lvl = 0

should be :

        if (Ier == 6) lvl = 1

Which matches all the other routines. This is a bug that has been in there since 1983.