Mercurial > octave
comparison liboctave/numeric/lo-specfun.cc @ 21066:258c787cd9ce
maint: Use "FIXME:" consistently in code base.
* file-editor-tab.h, mxarray.in.h, pr-output.cc, symtab.h, ov-bool-sparse.h,
ov-bool.h, ov-float.h, ov-scalar.h, ov-str-mat.cc, oct-parse.in.yy, f77-fcn.h,
lo-specfun.cc, randmtzig.c, randpoisson.c, lo-ieee.h, oct-rl-edit.c,
sparse-util.h: Use "FIXME:" consistently in code base.
| author | Rik <rik@octave.org> |
|---|---|
| date | Thu, 14 Jan 2016 07:02:41 -0800 |
| parents | f7084eae3318 |
| children | dfcb9d74b253 |
comparison
equal
deleted
inserted
replaced
| 21065:e1ee2203efe0 | 21066:258c787cd9ce |
|---|---|
| 3930 { | 3930 { |
| 3931 // Remember that psifn and dpsifn return scales values | 3931 // Remember that psifn and dpsifn return scales values |
| 3932 // When n is 1: do nothing since ((-1)**(n+1)/gamma(n+1)) == 1 | 3932 // When n is 1: do nothing since ((-1)**(n+1)/gamma(n+1)) == 1 |
| 3933 // When n is 0: change sign since ((-1)**(n+1)/gamma(n+1)) == -1 | 3933 // When n is 0: change sign since ((-1)**(n+1)/gamma(n+1)) == -1 |
| 3934 if (n > 1) | 3934 if (n > 1) |
| 3935 // FIXME xgamma here is a killer for our precision since it grows | 3935 // FIXME: xgamma here is a killer for our precision since it grows |
| 3936 // way too fast | 3936 // way too fast. |
| 3937 ans = ans / (pow (-1.0, n + 1) / xgamma (double (n+1))); | 3937 ans = ans / (pow (-1.0, n + 1) / xgamma (double (n+1))); |
| 3938 else if (n == 0) | 3938 else if (n == 0) |
| 3939 ans = -ans; | 3939 ans = -ans; |
| 3940 } | 3940 } |
| 3941 else if (ierr == 2) | 3941 else if (ierr == 2) |