diff options
author | Michele Calgaro <michele.calgaro@yahoo.it> | 2021-05-23 20:48:35 +0900 |
---|---|---|
committer | Michele Calgaro <michele.calgaro@yahoo.it> | 2021-05-29 15:17:38 +0900 |
commit | d63c9d696eb6e2539528b99afc21f4086c9defe3 (patch) | |
tree | b3bfc97a66431a12cdd8f9379c0072673ede43df /karbon/core/vglobal.cpp | |
parent | 5363fe3c36504c37bdc6dcfafd5f71daeae251e8 (diff) | |
download | koffice-d63c9d69.tar.gz koffice-d63c9d69.zip |
Renaming of files in preparation for code style tools.
Signed-off-by: Michele Calgaro <michele.calgaro@yahoo.it>
(cherry picked from commit 8b78a8791bc539bcffe7159f9d9714d577cb3d7d)
Diffstat (limited to 'karbon/core/vglobal.cpp')
-rw-r--r-- | karbon/core/vglobal.cpp | 90 |
1 files changed, 90 insertions, 0 deletions
diff --git a/karbon/core/vglobal.cpp b/karbon/core/vglobal.cpp new file mode 100644 index 000000000..b118e5e26 --- /dev/null +++ b/karbon/core/vglobal.cpp @@ -0,0 +1,90 @@ +/* This file is part of the KDE project + Copyright (C) 2001, The Karbon Developers + Copyright (C) 2002, The Karbon Developers + + This library is free software; you can redistribute it and/or + modify it under the terms of the GNU Library General Public + License as published by the Free Software Foundation; either + version 2 of the License, or (at your option) any later version. + + This library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Library General Public License for more details. + + You should have received a copy of the GNU Library General Public License + along with this library; see the file COPYING.LIB. If not, write to + the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, + * Boston, MA 02110-1301, USA. +*/ + + +#include <math.h> + +#include "vglobal.h" + + +int +VGlobal::binomialCoeff( unsigned n, unsigned k ) +{ + return + static_cast<int>( + 0.5 + + exp( + factorialLn( n ) - + factorialLn( k ) - + factorialLn( n - k ) ) ); +} + +double +VGlobal::factorialLn( unsigned n ) +{ + const unsigned cacheSize = 100; + + // A static array is initalized to zero. + static double cache[ cacheSize ]; + + + if( n <= 1 ) + return 0.0; + + if( n <= cacheSize - 1 ) + { + return cache[ n ] + ? cache[ n ] + : ( cache[ n ] = gammaLn( n + 1.0 ) ); + } + else + { + return gammaLn( n + 1.0 ); + } +} + +double +VGlobal::gammaLn( double x ) +{ + static double coeff[ 6 ] = + { + 76.18009172947146, + -86.50532032941677, + 24.01409824083091, + -1.231739572450155, + 0.1208650973866179e-2, + -0.5395239384953e-5 + }; + + double y = x; + + double tmp = x + 5.5; + tmp -= ( x + 0.5 ) * log( tmp ); + + double ser = 1.000000000190015; + + for( int i = 0; i < 5; ++i ) + { + ser += coeff[ i ] / ++y; + } + + return -tmp + log( 2.5066282746310005 * ser / x ); +} + |