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
126
127
128
129
130
131
132
133
134
135

//
// Creator:    http://www.dicelocksecurity.com
// Version:    vers.4.0.0.1
//
// Copyright © 2008-2010 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
// ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
// 
// DICELOCK IS A REGISTERED TRADEMARK OR TRADEMARK OF THE OWNERS.
// 
 
#include <stdexcept>
#include <stdlib.h>
#include <math.h>
#include "frequencyTest.h"
 
 
using namespace std;
 
 
namespace DiceLockSecurity {
  
  namespace RandomTest {
 
  // Random Test Class enumerator name
  const RandomTests FrequencyTest::test = Frequency;
  // Random Test Class minimum stream length
  const unsigned int  FrequencyTest::minimumLength = 100;
 
  // Constructor, default 
  FrequencyTest::FrequencyTest() {
 
    sum = 0;
    sumDiv_n = 0.0;
  }
 
 
  // Constructor with a MathematicalFunctions object instantiated 
  FrequencyTest::FrequencyTest(MathematicalFunctions* mathFuncObj) {
 
    sum = 0;
    sumDiv_n = 0.0;
  }
 
  // Destructor
  FrequencyTest::~FrequencyTest() {
 
    sum = 0;
    sumDiv_n = 0.0;
  }
  
  // Gets the BaseRandomTest random state of the last executed BaseCryptoRandomStream
  bool FrequencyTest::IsRandom(void) {
 
    return BaseRandomTest::IsRandom();
  }
 
  // Tests randomness of the BaseCryptoRandomStream and returns the random value
  bool FrequencyTest::IsRandom(BaseCryptoRandomStream* bitStream) {
    unsigned long i;
    double f, s_obs, sum;
    double sqrt2 = 1.41421356237309504880;
    unsigned short int bitTemp;
 
    if (bitStream->GetBitLength() < this->GetMinimumLength()) {
      this->error = InsufficientNumberOfBits;
      this->random = false;
      return this->random;
    }
    bitStream->SetBitPosition(0);
    this->error = NoError;
    sum = 0.0;
    for(i = 0; i < bitStream->GetBitLength(); i++) {
      bitTemp = (unsigned short int)bitStream->GetBitPosition(i);
      sum += ((2 * bitTemp) - 1);
    }
    s_obs = fabs(sum)/sqrt((double)bitStream->GetBitLength());
    f = s_obs/sqrt2;
    this->pValue = this->mathFuncs->ErFc(f);
    if (this->pValue < this->GetAlpha()) { 
      this->random = false;
    }
    else {
      this->random = true;
    }
    this->sum = (int)sum;
    this->sumDiv_n = sum/bitStream->GetBitLength();
    return this->random;
  }
 
  // Initializes the object
  void FrequencyTest::Initialize(void) {
 
    BaseRandomTest::Initialize();
    sum = 0;
    sumDiv_n = 0.0;
  }
 
  // Gets the type of the object
  RandomTests FrequencyTest::GetType(void) {
 
    return this->test;
  }
 
  // Gets the minimum random stream length
  unsigned int FrequencyTest::GetMinimumLength(void) {
 
    return this->minimumLength;
  }
 
  // Gets the "sum" result
  int FrequencyTest::GetSum(void) {
 
    return this->sum;
  }
 
  // Gets the "sumDiv_n" result
  double FrequencyTest::GetSumDiv_n(void) {
 
    return this->sumDiv_n;
  }
  }
}