bugGNU Scientific Library - Bugs: bug #39292, possible error in...


You are not allowed to post comments on this tracker with your current authentication level.

bug #39292: possible error in gsl_sf_coulomb_wave_FG_e

Submitter:  Patrick Alken <psa>
Submitted:  Wed 19 Jun 2013 04:37:38 PM UTC
Category:  Runtime error Severity:  3 - Normal
Operating System:  Status:  None
Assigned to:  None Open/Closed:  Open

Wed 19 Jun 2013 04:37:38 PM UTC, original submission:  

From Alexey A. Illarionov:

Accidentally found a bug in specfun/coulomb.c . Trivial error in
solution of quadratic equation. None of the current test cases are
influenced by this bug since the value of C (actually N) determines
only the condition of applicability of "continued fraction+recursion"

=== modified file 'specfunc/coulomb.c'
- --- specfunc/coulomb.c 2007-07-02 18:34:24 +0000
+++ specfunc/coulomb.c 2012-10-26 04:27:18 +0000
@@ -1121,7 +1121,7 @@
      *           we must go at least as low as lam_G
     const double SMALL = GSL_SQRT_DBL_EPSILON;
- -    const double C = sqrt(1.0 + 4.0*x*(x-2.0*eta));
+    const double C = 0.5 sqrt(1.0 + 4.0*x(x-2.0*eta));
     const int N = ceil(lam_F - C + 0.5);
     const double lam_0   = lam_F - GSL_MAX(N, 0);
     const double lam_min = GSL_MIN(lam_0, lam_G);

Patrick Alken <psa>
Group administrator


(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -email is unavailable- added by psa (Submitted the item)

    No changes have been made to this item

    Back to the top

    Powered by Savane 3.13-a854.
    Corresponding source code