Mercurial > octave

comparison liboctave/numeric/randpoisson.c @ 21066:258c787cd9ce

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 4197fc428c7d
children 26f85aa072de
comparison
equal deleted inserted replaced
21065:e1ee2203efe0 21066:258c787cd9ce
533 /* numerical recipes */ 533 /* numerical recipes */
534 poisson_rejection (L, &ret, 1); 534 poisson_rejection (L, &ret, 1);
535 } 535 }
536 else if (INFINITE(L)) 536 else if (INFINITE(L))
537 { 537 {
538 /* FIXME R uses NaN, but the normal approx. suggests that as 538 /* FIXME: R uses NaN, but the normal approximation suggests that
539 * limit should be inf. Which is correct? */ 539 * limit should be Inf. Which is correct? */
540 ret = NAN; 540 ret = NAN;
541 } 541 }
542 else 542 else
543 { 543 {
544 /* normal approximation: from Phys. Rev. D (1994) v50 p1284 */ 544 /* normal approximation: from Phys. Rev. D (1994) v50 p1284 */
606 /* numerical recipes */ 606 /* numerical recipes */
607 poisson_rejection_float (L, &ret, 1); 607 poisson_rejection_float (L, &ret, 1);
608 } 608 }
609 else if (INFINITE(L)) 609 else if (INFINITE(L))
610 { 610 {
611 /* FIXME R uses NaN, but the normal approx. suggests that as 611 /* FIXME: R uses NaN, but the normal approximation suggests that
612 * limit should be inf. Which is correct? */ 612 * limit should be Inf. Which is correct? */
613 ret = NAN; 613 ret = NAN;
614 } 614 }
615 else 615 else
616 { 616 {
617 /* normal approximation: from Phys. Rev. D (1994) v50 p1284 */ 617 /* normal approximation: from Phys. Rev. D (1994) v50 p1284 */