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 Split this topic September 17, 2021, 12:23am 8

A post was split to a new topic: Problem upgrading GR package

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

NaN returned by gauss()

Related topics