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

rafael.guerra September 16, 2021, 2:12pm 3

Is using BigFloats an option?
gauss(BigFloat, 1093, 38.68598671891722, 150000.0)

Topic		Replies	Views
QuadGK returns NaNs Numerics question , quadgk , integral	4	349	October 24, 2023
Handling an overflow example in QuadGK Numerics integral	5	417	July 16, 2023
Integrand produced a NaN. Domain error using QuadGK General Usage quadgk , integral	3	1653	April 30, 2021
Advice on approximating Complicated integral General Usage	4	1176	April 27, 2020
Unexpected behavior of QuadGK with extended precision floats Numerics integral , type-piracy	4	566	December 7, 2021

NaN returned by gauss()

Related topics