I’m using the function gauss() from QuadGK in order to calculate a non-adaptive quadrature and am irritated about the appearance of NaN beyond a certain N. Specifically, gauss(1093,38.68598671891722,150000.0) returns an array that is partly populated with NaN, whereas gauss(1092,38.68598671891722,15…

[image] trg818: I’m using the function gauss() from QuadGK in order to calculate a non-adaptive quadrature and am irritated about the appearance of NaN beyond a certain N . The bug is fixed in the latest version (2.4.2) of QuadGK: julia> all.(isfinite, gauss(1093,38.68598671891722,150000.0))…

NaN returned by gauss()

Specific Domains Numerics

stevengj September 16, 2021, 7:33pm 7

The bug is fixed in the latest version (2.4.2) of QuadGK:

julia> all.(isfinite, gauss(1093,38.68598671891722,150000.0))
(true, true)

1 Like

Topic		Replies	Views
QuadGK returns NaNs Numerics question , quadgk , integral	4	380	October 24, 2023
Handling an overflow example in QuadGK Numerics integral	5	478	July 16, 2023
LAPACK error from quadgk integration Numerics question	5	535	December 29, 2020
Unexpected behavior of QuadGK with extended precision floats Numerics integral , type-piracy	4	590	December 7, 2021
Integrand produced a NaN. Domain error using QuadGK General Usage quadgk , integral	3	1702	April 30, 2021

NaN returned by gauss()

Related topics