{"action":"created","issue":{"url":"https://api.github.com/repos/scipy/scipy/issues/5227","labels_url":"https://api.github.com/repos/scipy/scipy/issues/5227/labels{/name}","comments_url":"https://api.github.com/repos/scipy/scipy/issues/5227/comments","events_url":"https://api.github.com/repos/scipy/scipy/issues/5227/events","html_url":"https://github.com/scipy/scipy/issues/5227","id":104917266,"number":5227,"title":"Unexpected behavior of pro_rad1 for x < 1.0","user":{"login":"josePhoenix","id":172886,"avatar_url":"https://avatars.githubusercontent.com/u/172886?v=3","gravatar_id":"","url":"https://api.github.com/users/josePhoenix","html_url":"https://github.com/josePhoenix","followers_url":"https://api.github.com/users/josePhoenix/followers","following_url":"https://api.github.com/users/josePhoenix/following{/other_user}","gists_url":"https://api.github.com/users/josePhoenix/gists{/gist_id}","starred_url":"https://api.github.com/users/josePhoenix/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/josePhoenix/subscriptions","organizations_url":"https://api.github.com/users/josePhoenix/orgs","repos_url":"https://api.github.com/users/josePhoenix/repos","events_url":"https://api.github.com/users/josePhoenix/events{/privacy}","received_events_url":"https://api.github.com/users/josePhoenix/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":3,"created_at":"2015-09-04T15:50:27Z","updated_at":"2015-09-04T16:47:19Z","closed_at":null,"body":"The docstring for `scipy.special.pro_rad1` says \"Computes the prolate sheroidal radial function of the first kind and its derivative (with respect to x) for mode parameters m>=0 and n>=m, spheroidal parameter c and |x| < 1.0.\" However, when trying to compute the function with m = n = 0 and c = 1 on the interval (-1.0, 1.0), I get an array full of `nan` values.\r\n\r\nLooking at the wrapper code in https://github.com/scipy/scipy/blob/d1203c280a8661536d64f17e6a90d7ef000302d5/scipy/special/specfun_wrappers.c#L953, is it possible that the comparison is wrong? That is, it returns `nan` when `(x <= 1.0)`, when it should be testing for abs(x) > 1.\r\n\r\n```\r\nIn [5]: scipy.special.pro_rad1(0, 0, 10, 200)\r\nOut[5]: (0.00046907225742714526, -0.0017336833438237117)\r\n\r\nIn [6]: scipy.special.pro_rad1(0, 0, 1, 0.9)\r\nOut[6]: (nan, nan)\r\n```"},"comment":{"url":"https://api.github.com/repos/scipy/scipy/issues/comments/137788305","html_url":"https://github.com/scipy/scipy/issues/5227#issuecomment-137788305","issue_url":"https://api.github.com/repos/scipy/scipy/issues/5227","id":137788305,"user":{"login":"josePhoenix","id":172886,"avatar_url":"https://avatars.githubusercontent.com/u/172886?v=3","gravatar_id":"","url":"https://api.github.com/users/josePhoenix","html_url":"https://github.com/josePhoenix","followers_url":"https://api.github.com/users/josePhoenix/followers","following_url":"https://api.github.com/users/josePhoenix/following{/other_user}","gists_url":"https://api.github.com/users/josePhoenix/gists{/gist_id}","starred_url":"https://api.github.com/users/josePhoenix/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/josePhoenix/subscriptions","organizations_url":"https://api.github.com/users/josePhoenix/orgs","repos_url":"https://api.github.com/users/josePhoenix/repos","events_url":"https://api.github.com/users/josePhoenix/events{/privacy}","received_events_url":"https://api.github.com/users/josePhoenix/received_events","type":"User","site_admin":false},"created_at":"2015-09-04T16:47:19Z","updated_at":"2015-09-04T16:47:19Z","body":"For example, http://reference.wolfram.com/language/ref/SpheroidalS1.html"}}