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	373	October 24, 2023
LAPACK error from quadgk integration Numerics question	5	532	December 29, 2020
Unexpected behavior of QuadGK with extended precision floats Numerics integral , type-piracy	4	590	December 7, 2021
QuadGK don't work in CUDA kernel General Usage question , package , cudanative , cuda	6	962	July 30, 2019
Integrand produced a NaN. Domain error using QuadGK General Usage quadgk , integral	3	1690	April 30, 2021

NaN returned by gauss()

Related topics