Modify

Ticket #5560 (closed Bugs: fixed)

Opened 3 years ago

Last modified 3 years ago

cyl_bessel_j return NaN for low numbers

Reported by: agribov_no_s_p_a_m@… Owned by: johnmaddock
Milestone: To Be Determined Component: math
Version: Boost 1.45.0 Severity: Problem
Keywords: Cc:

Description

Calculate cyl_bessel_j<double, double> for v = 10.0 and x = 1.0e-100. Return value is NaN instead of 0.0. Problem happens in file: \math\special_functions\detail\bessel_jn.hpp Line: 92. ratio become NaN.

Attachments

Change History

comment:1 Changed 3 years ago by johnmaddock

Confirmed, investigating.

comment:2 Changed 3 years ago by johnmaddock

  • Status changed from new to closed
  • Resolution set to fixed

(In [72221]) Update Bessel functions to use series representations when x is small compared to v and also to prevent the generation of NaN's in such cases. Update docs to match. Fixes #5560.

comment:3 Changed 3 years ago by johnmaddock

  • Status changed from closed to reopened
  • Resolution fixed deleted

Reopened: there's still an issue for v < 0 and |v| >> x. Will be harder to fix though :-(

comment:4 Changed 3 years ago by johnmaddock

  • Status changed from reopened to closed
  • Resolution set to fixed

(In [72415]) Add new series evaluation forms for the Bessel J and Y functions. Also take extra care when evaluating to avoid spurious overflow/underflow. Changed check_series_iterations to know the real-number type so messages are properly formatted. Updated docs to match and regenerated. Fixes #5560.

View

Add a comment

Modify Ticket

Change Properties
<Author field>
Action
as closed
The resolution will be deleted. Next status will be 'reopened'
Author


E-mail address and user name can be saved in the Preferences.

 
Note: See TracTickets for help on using tickets.