changeset 30941:272614f05636

besselj.cc: Add more BISTs for the function airy (bug #62321) besselj.cc: Also tighten tolerance on tests
author Arun Giridhar <arungiridhar@gmail.com>
date Sun, 17 Apr 2022 17:58:14 -0400
parents 84944164799e
children 84e7222b6b5c
files libinterp/corefcn/besselj.cc
diffstat 1 files changed, 33 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
--- 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 <K must be an integer value> (airy ("foo", 2, false))