From 99100bc76417c1ab40115939198181724d13618d Mon Sep 17 00:00:00 2001
From: Emmanuel Christophe <emmanuel.christophe@orfeo-toolbox.org>
Date: Sat, 7 May 2011 11:16:31 -0700
Subject: [PATCH] STYLE: otbGamma
---
Code/Common/otbGamma.h | 141 +++++++++++++++++++++--------------------
1 file changed, 74 insertions(+), 67 deletions(-)
diff --git a/Code/Common/otbGamma.h b/Code/Common/otbGamma.h
index 25f8419ea8..893d681a4d 100644
--- a/Code/Common/otbGamma.h
+++ b/Code/Common/otbGamma.h
@@ -34,75 +34,82 @@ namespace otb
*/
class Gamma {
public:
- double gamma(double x)
- {
- int i, k, m;
- double ga, gr, r, z;
+ double gamma(double x)
+ {
+ int m;
+ double ga, gr, r, z;
- static double g[] = {
- 1.0,
- 0.5772156649015329,
- -0.6558780715202538,
- -0.420026350340952e-1,
- 0.1665386113822915,
- -0.421977345555443e-1,
- -0.9621971527877e-2,
- 0.7218943246663e-2,
- -0.11651675918591e-2,
- -0.2152416741149e-3,
- 0.1280502823882e-3,
- -0.201348547807e-4,
- -0.12504934821e-5,
- 0.1133027232e-5,
- -0.2056338417e-6,
- 0.6116095e-8,
- 0.50020075e-8,
- -0.11812746e-8,
- 0.1043427e-9,
- 0.77823e-11,
- -0.36968e-11,
- 0.51e-12,
- -0.206e-13,
- -0.54e-14,
- 0.14e-14};
+ static double g[] = {
+ 1.0,
+ 0.5772156649015329,
+ -0.6558780715202538,
+ -0.420026350340952e-1,
+ 0.1665386113822915,
+ -0.421977345555443e-1,
+ -0.9621971527877e-2,
+ 0.7218943246663e-2,
+ -0.11651675918591e-2,
+ -0.2152416741149e-3,
+ 0.1280502823882e-3,
+ -0.201348547807e-4,
+ -0.12504934821e-5,
+ 0.1133027232e-5,
+ -0.2056338417e-6,
+ 0.6116095e-8,
+ 0.50020075e-8,
+ -0.11812746e-8,
+ 0.1043427e-9,
+ 0.77823e-11,
+ -0.36968e-11,
+ 0.51e-12,
+ -0.206e-13,
+ -0.54e-14,
+ 0.14e-14};
- if (x > 171.0) return 1e308; // This value is an overflow flag.
- if (x == (int)x) {
- if (x > 0.0) {
- ga = 1.0; // use factorial
- for (i=2; i<x; i++) {
- ga *= i;
- }
- }
- else
- ga = 1e308;
- }
- else {
- if (vcl_abs(x) > 1.0) {
- z = vcl_abs(x);
- m = (int)z;
- r = 1.0;
- for (k=1; k<=m; k++) {
- r *= (z-k);
- }
- z -= m;
- }
- else
- z = x;
- gr = g[24];
- for (k=23; k>=0; k--) {
- gr = gr*z+g[k];
- }
- ga = 1.0/(gr*z);
- if (vcl_abs(x) > 1.0) {
- ga *= r;
- if (x < 0.0) {
- ga = -otb::CONST_PI/(x*ga*vcl_sin(otb::CONST_PI*x));
- }
- }
- }
- return ga;
- }
+ if (x > 171.0) return 1e308; // This value is an overflow flag.
+ if (x == (int)x)
+ {
+ if (x > 0.0)
+ {
+ ga = 1.0; // use factorial
+ for (int i = 2; i < x; i++)
+ {
+ ga *= i;
+ }
+ }
+ else ga = 1e308;
+ }
+ else
+ {
+ if (vcl_abs(x) > 1.0)
+ {
+ z = vcl_abs(x);
+ m = (int)z;
+ r = 1.0;
+ for (int k = 1; k <= m; k++)
+ {
+ r *= (z-k);
+ }
+ z -= m;
+ }
+ else z = x;
+ gr = g[24];
+ for (int k = 23; k >= 0; k--)
+ {
+ gr = gr*z+g[k];
+ }
+ ga = 1.0/(gr*z);
+ if (vcl_abs(x) > 1.0)
+ {
+ ga *= r;
+ if (x < 0.0)
+ {
+ ga = -otb::CONST_PI/(x*ga*vcl_sin(otb::CONST_PI*x));
+ }
+ }
+ }
+ return ga;
+ }
};
} // end namespace otb
--
GitLab