# HG changeset patch # User Arun Giridhar # Date 1650232694 14400 # Node ID 272614f05636c236d3f92d18d5bc30593a6cc6e4 # Parent 84944164799ed15f47fcb8d59d8bc0a931e37cda besselj.cc: Add more BISTs for the function airy (bug #62321) besselj.cc: Also tighten tolerance on tests diff -r 84944164799e -r 272614f05636 libinterp/corefcn/besselj.cc --- a/libinterp/corefcn/besselj.cc Sun Apr 17 16:57:39 2022 -0400 +++ b/libinterp/corefcn/besselj.cc Sun Apr 17 17:58:14 2022 -0400 @@ -719,15 +719,40 @@ /* %!test <*62321> -%! assert (airy (0, 1, false), 0.1352924163128814, 1e-14) -%! assert (airy (0, -1, false), 0.5355608832923521, 1e-14) -%! assert (airy (0, i, false), 0.3314933054321411 - 0.3174498589684437i, 1e-14) -%! assert (airy (0, -i, false), 0.3314933054321411 + 0.3174498589684437i, 1e-14) -%! assert (airy (0, 1, true), 0.2635136447491401, 1e-14) -%! assert (airy (0, -1, true), 0.4208904755499093 - 0.3311746779333462i, 1e-14) -%! assert (airy (0, i, true), 0.2743053542644657 - 0.08256069414005915i, 1e-14) -%! assert (airy (0, -i, true), 0.2743053542644657 + 0.08256069414005915i, 1e-14) +%! assert (airy (0, +1, false), 0.1352924163128814, 1e-15) +%! assert (airy (0, -1, false), 0.5355608832923521, 1e-15) +%! assert (airy (0, +i, false), 0.3314933054321411 - 0.3174498589684437i, 1e-15) +%! assert (airy (0, -i, false), 0.3314933054321411 + 0.3174498589684437i, 1e-15) +%! assert (airy (0, +1, true), 0.2635136447491401, 1e-15) +%! assert (airy (0, -1, true), 0.4208904755499093 - 0.3311746779333462i, 1e-15) +%! assert (airy (0, +i, true), 0.2743053542644657 - 0.0825606941400591i, 1e-15) +%! assert (airy (0, -i, true), 0.2743053542644657 + 0.0825606941400591i, 1e-15) +%! assert (airy (1, +1, false), -0.1591474412967932, 1e-15) +%! assert (airy (1, -1, false), -0.0101605671166451, 1e-15) +%! assert (airy (1, +i, false), -0.4324926598418070 + 0.0980478562292432i, 1e-15) +%! assert (airy (1, -i, false), -0.4324926598418070 - 0.0980478562292432i, 1e-15) +%! assert (airy (1, +1, true), -0.3099768889605148, 1e-15) +%! assert (airy (1, -1, true), -0.0079850602592406 + 0.0062829878869967i, 1e-15) +%! assert (airy (1, +i, true), -0.2682792317254506 - 0.0680655016581369i, 1e-15) +%! assert (airy (1, -i, true), -0.2682792317254506 + 0.0680655016581369i, 1e-15) +%! assert (airy (2, +1, false), 1.2074235949528715, 1e-15) +%! assert (airy (2, -1, false), 0.1039973894969446, 1e-15) +%! assert (airy (2, +i, false), 0.6488582083303949 + 0.3449586347680484i, 1e-15) +%! assert (airy (2, -i, false), 0.6488582083303949 - 0.3449586347680484i, 1e-15) +%! assert (airy (2, +1, true), 0.6199119435726786, 1e-15) +%! assert (airy (2, -1, true), 0.1039973894969446, 1e-15) +%! assert (airy (2, +i, true), 0.4049686654663902 + 0.2152973271658010i, 1e-15) +%! assert (airy (2, -i, true), 0.4049686654663902 - 0.2152973271658010i, 1e-15) +%! assert (airy (3, +1, false), 0.9324359333927754, 1e-15) +%! assert (airy (3, -1, false), 0.5923756264227923, 1e-15) +%! assert (airy (3, +i, false), 0.1350266467108190 - 0.1288373867812549i, 1e-15) +%! assert (airy (3, -i, false), 0.1350266467108190 + 0.1288373867812549i, 1e-15) +%! assert (airy (3, +1, true), 0.4787285706049846, 1e-15) +%! assert (airy (3, -1, true), 0.5923756264227923, 1e-15) +%! assert (airy (3, +i, true), 0.0842735134099415 - 0.0804106412111761i, 1e-15) +%! assert (airy (3, -i, true), 0.0842735134099415 + 0.0804106412111761i, 1e-15) +Input validation tests %!error (airy ()) %!error (airy (0, 1, 2, 3)) %!error (airy ("foo", 2, false))