perazz / fortran-bessels

Fortran port of the Bessels.jl repository
MIT License
10 stars 2 forks source link

resolving the potential loss of precision #1

Closed shahmoradi closed 2 years ago

shahmoradi commented 2 years ago

Hi @perazz , Thanks for starting this. Looking at the code, it came to my attention that [real(BK :: ...] may not necessarily respect the requested precision (in fact, compilers often do not). This pull request converts all such declarations to explicit suffixes real numbers to ensure the desired precision in the constants. I look forward to seeing and helping with the full implementation of procedures.