forked from Kitware/VTK
-
Notifications
You must be signed in to change notification settings - Fork 0
/
vtkMinimalStandardRandomSequence.cxx
123 lines (104 loc) · 3.67 KB
/
vtkMinimalStandardRandomSequence.cxx
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
/*=========================================================================
Program: Visualization Toolkit
Module: vtkMinimalStandardRandomSequence.cxx
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
#include "vtkMinimalStandardRandomSequence.h"
#include "vtkObjectFactory.h"
#include <cassert>
vtkStandardNewMacro(vtkMinimalStandardRandomSequence);
const int VTK_K_A = 16807;
const int VTK_K_M = 2147483647; // Mersenne prime 2^(31)-1
const int VTK_K_Q = 127773; // M/A
const int VTK_K_R = 2836; // M%A
//------------------------------------------------------------------------------
vtkMinimalStandardRandomSequence::vtkMinimalStandardRandomSequence()
{
this->Seed = 1;
}
//------------------------------------------------------------------------------
vtkMinimalStandardRandomSequence::~vtkMinimalStandardRandomSequence() = default;
//------------------------------------------------------------------------------
void vtkMinimalStandardRandomSequence::SetSeedOnly(int value)
{
this->Seed = value;
// fit the seed to the valid range [1,2147483646]
if (this->Seed < 1)
{
this->Seed += 2147483646;
}
else
{
if (this->Seed == 2147483647)
{
this->Seed = 1;
}
}
}
//------------------------------------------------------------------------------
void vtkMinimalStandardRandomSequence::SetSeed(int value)
{
this->SetSeedOnly(value);
// the first random number after setting the seed is proportional to the
// seed value. To help solve this, call Next() a few times.
// This doesn't ruin the repeatability of Next().
this->Next();
this->Next();
this->Next();
}
//------------------------------------------------------------------------------
int vtkMinimalStandardRandomSequence::GetSeed()
{
return this->Seed;
}
//------------------------------------------------------------------------------
double vtkMinimalStandardRandomSequence::GetValue()
{
double result = static_cast<double>(this->Seed) / VTK_K_M;
assert("post: unit_range" && result >= 0.0 && result <= 1.0);
return result;
}
//------------------------------------------------------------------------------
void vtkMinimalStandardRandomSequence::Next()
{
int hi = this->Seed / VTK_K_Q;
int lo = this->Seed % VTK_K_Q;
this->Seed = VTK_K_A * lo - VTK_K_R * hi;
if (this->Seed <= 0)
{
this->Seed += VTK_K_M;
}
}
//------------------------------------------------------------------------------
double vtkMinimalStandardRandomSequence::GetRangeValue(double rangeMin, double rangeMax)
{
double result;
if (rangeMin == rangeMax)
{
result = rangeMin;
}
else
{
result = rangeMin + this->GetValue() * (rangeMax - rangeMin);
}
assert("post: valid_result" &&
((rangeMin <= rangeMax && result >= rangeMin && result <= rangeMax) ||
(rangeMax <= rangeMin && result >= rangeMax && result <= rangeMin)));
return result;
}
//------------------------------------------------------------------------------
double vtkMinimalStandardRandomSequence::GetNextRangeValue(double rangeMin, double rangeMax)
{
this->Next();
return this->GetRangeValue(rangeMin, rangeMax);
}
//------------------------------------------------------------------------------
void vtkMinimalStandardRandomSequence::PrintSelf(ostream& os, vtkIndent indent)
{
this->Superclass::PrintSelf(os, indent);
}