| Anonymous | Login | Signup for a new account | 2024-11-22 09:21 MSK |
| Main | My View | View Issues | Change Log | Roadmap | Docs |
| Viewing Issue Simple Details [ Jump to Notes ] | [ View Advanced ] [ Issue History ] [ Print ] | ||||||
| ID | Category | Severity | Reproducibility | Date Submitted | Last Update | ||
| 0000010 | [ALGLIB] Spec.functions | major | always | 2009-03-10 21:02 | 2009-07-23 12:42 | ||
| Reporter | SergeyB | View Status | public | ||||
| Assigned To | SergeyB | ||||||
| Priority | normal | Resolution | fixed | ||||
| Status | resolved | Product Version | |||||
| Summary | 0000010: FIXED: modified Besel function returns incorrect results | ||||||
| Description |
From: Kors van der Werf < korsvanderwerf@yahoo.com > Hi, Very nice, ythis ALGLib. I use the VB version for now, and it seems the Modified BesselI funtion for orders of negative half integers do not give proper results, can you have a look? Already order = -0.5 fails, and can be compared to the well known sqrt(2/(PI x)) cosh(x) For x=1 it should be 1.2312, I get 0.293525. Could you ahve a look please. |
||||||
| Additional Information | |||||||
| Programming language | All | ||||||
| Attached Files | |||||||
|
|
|||||||
 Relationships |
|
|
Notes |
|
|
(0000006) SergeyB (administrator) 2009-03-11 10:21 |
Confirmed: should be 1.2312, returns 0.2935 |
|
(0000007) SergeyB (administrator) 2009-03-11 12:15 edited on: 2009-03-11 14:12 |
Fixed. There was an error in original Cephes code when calculating 1F1(v,2*v,2*x) with v=0. Straightforward formula for 1F1 may return 1.0 or NAN (depending on implementation, 1.0 in Cephes), more elaborate calculations gives us 0.5*(1+exp(2*x)). This error can only happen when v=-0.5, should be OK with other values of v. |
|
(0000008) SergeyB (administrator) 2009-03-11 23:36 |
Other negative half-integer values should be fixed also... |
|
(0000010) SergeyB (administrator) 2009-04-19 11:07 |
Looks like that this subroutine can't calculate modified Bessel functions of negative orders. It is not a bug - it is a limitation of formula used. I've modified documentation... |
|
Issue History |
|||
| Date Modified | Username | Field | Change |
| 2009-03-10 21:02 | SergeyB | New Issue | |
| 2009-03-10 21:02 | SergeyB | Programming language | => VBA |
| 2009-03-10 21:50 | SergeyB | Status | new => assigned |
| 2009-03-10 21:50 | SergeyB | Assigned To | => SergeyB |
| 2009-03-10 21:54 | SergeyB | Assigned To | SergeyB => |
| 2009-03-10 21:55 | SergeyB | Status | assigned => new |
| 2009-03-11 10:20 | SergeyB | Status | new => assigned |
| 2009-03-11 10:20 | SergeyB | Assigned To | => SergeyB |
| 2009-03-11 10:21 | SergeyB | Note Added: 0000006 | |
| 2009-03-11 10:21 | SergeyB | Status | assigned => confirmed |
| 2009-03-11 10:21 | SergeyB | Reproducibility | have not tried => always |
| 2009-03-11 12:04 | SergeyB | Programming language | VBA => All |
| 2009-03-11 12:15 | SergeyB | Note Added: 0000007 | |
| 2009-03-11 12:15 | SergeyB | Status | confirmed => resolved |
| 2009-03-11 12:15 | SergeyB | Fixed in Version | => Next minor release |
| 2009-03-11 12:15 | SergeyB | Resolution | open => fixed |
| 2009-03-11 12:17 | SergeyB | Note Edited: 0000007 | |
| 2009-03-11 14:12 | SergeyB | Note Edited: 0000007 | |
| 2009-03-11 23:36 | SergeyB | Note Added: 0000008 | |
| 2009-03-11 23:36 | SergeyB | Status | resolved => feedback |
| 2009-03-11 23:36 | SergeyB | Resolution | fixed => reopened |
| 2009-04-19 11:07 | SergeyB | Note Added: 0000010 | |
| 2009-04-19 11:07 | SergeyB | Status | feedback => resolved |
| 2009-04-19 11:07 | SergeyB | Resolution | reopened => fixed |
| 2009-04-20 09:42 | SergeyB | Summary | modified Besel function returns incorrect results => FIXED: modified Besel function returns incorrect results |
| 2009-04-20 09:47 | SergeyB | Target Version | => Upcoming release |
| 2009-07-23 12:42 | SergeyB | Category | => Special functions |
| Mantis 1.1.6[^] Copyright © 2000 - 2008 Mantis Group |  |
