1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
|
//
// Creator: http://www.dicelocksecurity.com
// Version: vers.5.0.0.1
//
// Copyright © 2008-2011 DiceLock Security, LLC. All rights reserved.
//
// DISCLAIMER
//
// THIS SOFTWARE IS PROVIDED "AS IS" AND ANY EXPRESSED OR IMPLIED WARRANTIES,
// INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY
// AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
// REGENTS OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
// OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
// WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
// OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
//
// DICELOCK IS A REGISTERED TRADEMARK OR TRADEMARK OF THE OWNERS.
//
// DICELOCK IS PROTECTED BY US PATENT 7508945 AND EUROPEAN PATENT 1182777 WHERE APPLICABLE.
//
// LICENSE INFORMATION CAN BE OBTAINED AT CORPORATE WEB SITE
//
#ifndef MATHEMATICALFUNCTIONS_HPP
#define MATHEMATICALFUNCTIONS_HPP
#ifdef DICELOCKCIPHER_EXPORTS
#define CLASS_DECLSPEC __declspec(dllexport)
#else
#define CLASS_DECLSPEC __declspec(dllimport)
#endif
#define MAX(x,y) ((x) < (y) ? (y) : (x))
#define MIN(x,y) ((x) > (y) ? (y) : (x))
#define isNonPositive(x) ((x) <= 0.e0 ? 1 : 0)
#define isPositive(x) ((x) > 0.e0 ? 1 : 0)
#define isNegative(x) ((x) < 0.e0 ? 1 : 0)
#define isGreaterThanOne(x) ((x) > 1.e0 ? 1 : 0)
#define isZero(x) ((x) == 0.e0 ? 1 : 0)
#define isOne(x) ((x) == 1.e0 ? 1 : 0)
namespace DiceLockSecurity {
namespace RandomTest {
enum MathematicalErrors {
None,
Overflow,
Underflow,
MathematicianNan
};
class MathematicalFunctions {
private:
static const double MACHEP; // 2**-53
static const double MAXLOG; // log(MAXNUM)
static const double MAXNUM; // 2**1024*(1-MACHEP)
static const double PI; // pi
static const double LOGPI; // log(pi)
static const double LS2PI; //( log( sqrt( 2*pi ) )
static const double big;
static const double biginv;
// A[]: Stirling's formula expansion of log gamma
// B[], C[]: log gamma function between 2 and 3
static double A_lgam[];
static double B_lgam[];
static double C_lgam[];
static const double MAXLGM;
protected:
bool Error;
MathematicalErrors MathError;
public:
// Constructor, default
CLASS_DECLSPEC MathematicalFunctions();
// Destructor
CLASS_DECLSPEC ~MathematicalFunctions();
// Logarithm of gamma function
CLASS_DECLSPEC double LGamma(double);
// incomplete gamma function
CLASS_DECLSPEC double IGamma(double, double);
// Complemented incomplete gamma integral
CLASS_DECLSPEC double IGammaC(double, double);
// Evaluate polynomial of degree N
CLASS_DECLSPEC double PolEvl(double,double *, int);
// N
// Evaluate polynomial when coefficient of x is 1.0.
CLASS_DECLSPEC double P1Evl(double, double *, int);
// Error function in double precision
CLASS_DECLSPEC double ErF(double);
// Error function in double precision
CLASS_DECLSPEC double ErFc(double);
// Statistical Normal function
CLASS_DECLSPEC double Normal(double);
// Class common error handling member
CLASS_DECLSPEC int GetError();
// Class common error handling member
CLASS_DECLSPEC MathematicalErrors GetMathError();
};
}
}
#endif
|
