On Tue, Aug 9, 2022 at 10:23 AM Patrick Dupre <pdu...@gmx.com> wrote:
> Hello, > > I get an error > gsl: qag.c:247: ERROR: roundoff error prevents tolerance from being achieve > Then, I tried to see where is this file > qag.c > > Unfortunately > dnf repoquery --whatprovides 'qag.c' > > does not provide the information > qag.c is from a C translation of QUADPACK, a Fortran library for numerical integration. The error is very unlikely to be a bug -- numerical quadrature routines often allow you to specify the desired error tolerance of the result, and will issue an error if they are not able to provide the required tolerance. https://www.gnu.org/software/gsl/doc/html/integration.html has: int gsl_integration_qag(const gsl_function <https://www.gnu.org/software/gsl/doc/html/roots.html#c.gsl_function> *f, double a, double b, double epsabs, double epsrel, size_t limit, int key, gsl_integration_workspace <https://www.gnu.org/software/gsl/doc/html/integration.html#c.gsl_integration_workspace> *workspace, double *result, double *abserr) <https://www.gnu.org/software/gsl/doc/html/integration.html#c.gsl_integration_qag> This function applies an integration rule adaptively until an estimate of the integral of [image: f] over [image: (a,b)] is achieved within the desired absolute and relative error limits, epsabs <https://www.gnu.org/software/gsl/doc/html/integration.html#c.gsl_integration_qag> and epsrel <https://www.gnu.org/software/gsl/doc/html/integration.html#c.gsl_integration_qag>. The function returns the final approximation, result <https://www.gnu.org/software/gsl/doc/html/integration.html#c.gsl_integration_qag>, and an estimate of the absolute error, abserr <https://www.gnu.org/software/gsl/doc/html/integration.html#c.gsl_integration_qag>. The integration rule is determined by the value of key <https://www.gnu.org/software/gsl/doc/html/integration.html#c.gsl_integration_qag>, which should be chosen from the following symbolic names, -- George N. White III
_______________________________________________ users mailing list -- users@lists.fedoraproject.org To unsubscribe send an email to users-le...@lists.fedoraproject.org Fedora Code of Conduct: https://docs.fedoraproject.org/en-US/project/code-of-conduct/ List Guidelines: https://fedoraproject.org/wiki/Mailing_list_guidelines List Archives: https://lists.fedoraproject.org/archives/list/users@lists.fedoraproject.org Do not reply to spam, report it: https://pagure.io/fedora-infrastructure/new_issue
