author | Timo Bingmann <timo.bingmann@student.kit.edu> |
Wed, 04 Mar 2009 18:55:32 +0100 | |
changeset 4256 | cc06d903ca09 |
parent 4249 | c2d4a91bcc8a |
child 4317 | 0a250f44e0ed |
permissions | -rw-r--r-- |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1 |
/* -*- Mode:C++; c-file-style:"gnu"; indent-tabs-mode:nil; -*- */ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
2 |
// |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
3 |
// Copyright (c) 2006 Georgia Tech Research Corporation |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
4 |
// |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
5 |
// This program is free software; you can redistribute it and/or modify |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
6 |
// it under the terms of the GNU General Public License version 2 as |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
7 |
// published by the Free Software Foundation; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
8 |
// |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
9 |
// This program is distributed in the hope that it will be useful, |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
10 |
// but WITHOUT ANY WARRANTY; without even the implied warranty of |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
11 |
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
12 |
// GNU General Public License for more details. |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
13 |
// |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
14 |
// You should have received a copy of the GNU General Public License |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
15 |
// along with this program; if not, write to the Free Software |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
16 |
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
17 |
// |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
18 |
// Author: Rajib Bhattacharjea<raj.b@gatech.edu> |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
19 |
// Author: Hadi Arbabi<marbabi@cs.odu.edu> |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
20 |
// |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
21 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
22 |
#include <iostream> |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
23 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
24 |
#include <math.h> |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
25 |
#include <stdlib.h> |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
26 |
#include <sys/time.h> // for gettimeofday |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
27 |
#include <unistd.h> |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
28 |
#include <iostream> |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
29 |
#include <sys/types.h> |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
30 |
#include <sys/stat.h> |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
31 |
#include <fcntl.h> |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
32 |
#include <sstream> |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
33 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
34 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
35 |
#include "assert.h" |
4235
c070d2fca46d
Fix previous changeset
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4234
diff
changeset
|
36 |
#include "config.h" |
4234
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
37 |
#include "integer.h" |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
38 |
#include "random-variable.h" |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
39 |
#include "rng-stream.h" |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
40 |
#include "fatal-error.h" |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
41 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
42 |
using namespace std; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
43 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
44 |
namespace ns3{ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
45 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
46 |
//----------------------------------------------------------------------------- |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
47 |
// Seed Manager |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
48 |
//----------------------------------------------------------------------------- |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
49 |
|
4234
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
50 |
uint32_t SeedManager::GetSeed() |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
51 |
{ |
4234
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
52 |
uint32_t s[6]; |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
53 |
RngStream::GetPackageSeed (s); |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
54 |
NS_ASSERT( |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
55 |
s[0] == s[1] && |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
56 |
s[0] == s[2] && |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
57 |
s[0] == s[3] && |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
58 |
s[0] == s[4] && |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
59 |
s[0] == s[5] |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
60 |
); |
7ec503ed040d
Almost have env variable support for RNG right
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4227
diff
changeset
|
61 |
return s[0]; |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
62 |
} |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
63 |
|
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
64 |
void SeedManager::SetSeed(uint32_t seed) |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
65 |
{ |
4235
c070d2fca46d
Fix previous changeset
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4234
diff
changeset
|
66 |
Config::SetGlobal("RngSeed", IntegerValue(seed)); |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
67 |
} |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
68 |
|
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
69 |
void SeedManager::SetRun(uint32_t run) |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
70 |
{ |
4235
c070d2fca46d
Fix previous changeset
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4234
diff
changeset
|
71 |
Config::SetGlobal("RngRun", IntegerValue(run)); |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
72 |
} |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
73 |
|
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
74 |
uint32_t SeedManager::GetRun() |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
75 |
{ |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
76 |
return RngStream::GetPackageRun (); |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
77 |
} |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
78 |
|
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
79 |
bool SeedManager::CheckSeed (uint32_t seed) |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
80 |
{ |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
81 |
return RngStream::CheckSeed(seed); |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
82 |
} |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
83 |
|
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
84 |
//----------------------------------------------------------------------------- |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
85 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
86 |
// RandomVariableBase methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
87 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
88 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
89 |
class RandomVariableBase |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
90 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
91 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
92 |
RandomVariableBase (); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
93 |
RandomVariableBase (const RandomVariableBase &o); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
94 |
virtual ~RandomVariableBase(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
95 |
virtual double GetValue() = 0; |
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
96 |
virtual uint32_t GetInteger(); |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
97 |
virtual RandomVariableBase* Copy(void) const = 0; |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
98 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
99 |
protected: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
100 |
RngStream* m_generator; //underlying generator being wrapped |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
101 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
102 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
103 |
RandomVariableBase::RandomVariableBase() |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
104 |
: m_generator(NULL) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
105 |
{ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
106 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
107 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
108 |
RandomVariableBase::RandomVariableBase(const RandomVariableBase& r) |
1811
85b79acd72ad
Fixed randomvariable memory bug that caused a segfault
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1806
diff
changeset
|
109 |
:m_generator(0) |
360
7bffd987426c
Fixed major memory bug and docs in RandomVariable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
346
diff
changeset
|
110 |
{ |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
111 |
if (r.m_generator) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
112 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
113 |
m_generator = new RngStream(*r.m_generator); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
114 |
} |
360
7bffd987426c
Fixed major memory bug and docs in RandomVariable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
346
diff
changeset
|
115 |
} |
7bffd987426c
Fixed major memory bug and docs in RandomVariable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
346
diff
changeset
|
116 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
117 |
RandomVariableBase::~RandomVariableBase() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
118 |
{ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
119 |
delete m_generator; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
120 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
121 |
|
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
122 |
uint32_t RandomVariableBase::GetInteger() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
123 |
{ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
124 |
return (uint32_t)GetValue(); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
125 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
126 |
|
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
127 |
//------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
128 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
129 |
RandomVariable::RandomVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
130 |
: m_variable (0) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
131 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
132 |
RandomVariable::RandomVariable(const RandomVariable&o) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
133 |
: m_variable (o.m_variable->Copy ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
134 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
135 |
RandomVariable::RandomVariable (const RandomVariableBase &variable) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
136 |
: m_variable (variable.Copy ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
137 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
138 |
RandomVariable & |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
139 |
RandomVariable::operator = (const RandomVariable &o) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
140 |
{ |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
141 |
if (&o == this) |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
142 |
{ |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
143 |
return *this; |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
144 |
} |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
145 |
delete m_variable; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
146 |
m_variable = o.m_variable->Copy (); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
147 |
return *this; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
148 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
149 |
RandomVariable::~RandomVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
150 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
151 |
delete m_variable; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
152 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
153 |
double |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
154 |
RandomVariable::GetValue (void) const |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
155 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
156 |
return m_variable->GetValue (); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
157 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
158 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
159 |
uint32_t |
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
160 |
RandomVariable::GetInteger (void) const |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
161 |
{ |
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
162 |
return m_variable->GetInteger (); |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
163 |
} |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
164 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
165 |
RandomVariableBase * |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
166 |
RandomVariable::Peek (void) const |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
167 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
168 |
return m_variable; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
169 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
170 |
|
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
171 |
|
2445
e1b65471f159
convert RandomVariable to new helpers.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2439
diff
changeset
|
172 |
ATTRIBUTE_VALUE_IMPLEMENT (RandomVariable); |
e1b65471f159
convert RandomVariable to new helpers.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2439
diff
changeset
|
173 |
ATTRIBUTE_CHECKER_IMPLEMENT (RandomVariable); |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
174 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
175 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
176 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
177 |
// UniformVariableImpl |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
178 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
179 |
class UniformVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
180 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
181 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
182 |
* Creates a uniform random number generator in the |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
183 |
* range [0.0 .. 1.0). |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
184 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
185 |
UniformVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
186 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
187 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
188 |
* Creates a uniform random number generator with the specified range |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
189 |
* \param s Low end of the range |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
190 |
* \param l High end of the range |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
191 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
192 |
UniformVariableImpl(double s, double l); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
193 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
194 |
UniformVariableImpl(const UniformVariableImpl& c); |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
195 |
|
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
196 |
double GetMin (void) const; |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
197 |
double GetMax (void) const; |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
198 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
199 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
200 |
* \return A value between low and high values specified by the constructor |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
201 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
202 |
virtual double GetValue(); |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
203 |
|
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
204 |
/** |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
205 |
* \return A value between low and high values specified by parameters |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
206 |
*/ |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
207 |
virtual double GetValue(double s, double l); |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
208 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
209 |
virtual RandomVariableBase* Copy(void) const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
210 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
211 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
212 |
double m_min; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
213 |
double m_max; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
214 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
215 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
216 |
UniformVariableImpl::UniformVariableImpl() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
217 |
: m_min(0), m_max(1.0) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
218 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
219 |
UniformVariableImpl::UniformVariableImpl(double s, double l) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
220 |
: m_min(s), m_max(l) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
221 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
222 |
UniformVariableImpl::UniformVariableImpl(const UniformVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
223 |
: RandomVariableBase(c), m_min(c.m_min), m_max(c.m_max) { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
224 |
|
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
225 |
double |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
226 |
UniformVariableImpl::GetMin (void) const |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
227 |
{ |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
228 |
return m_min; |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
229 |
} |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
230 |
double |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
231 |
UniformVariableImpl::GetMax (void) const |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
232 |
{ |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
233 |
return m_max; |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
234 |
} |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
235 |
|
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
236 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
237 |
double UniformVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
238 |
{ |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
239 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
240 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
241 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
242 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
243 |
return m_min + m_generator->RandU01() * (m_max - m_min); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
244 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
245 |
|
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
246 |
double UniformVariableImpl::GetValue(double s, double l) |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
247 |
{ |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
248 |
if(!m_generator) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
249 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
250 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
251 |
} |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
252 |
return s + m_generator->RandU01() * (l-s); |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
253 |
} |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
254 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
255 |
RandomVariableBase* UniformVariableImpl::Copy() const |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
256 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
257 |
return new UniformVariableImpl(*this); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
258 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
259 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
260 |
UniformVariable::UniformVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
261 |
: RandomVariable (UniformVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
262 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
263 |
UniformVariable::UniformVariable(double s, double l) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
264 |
: RandomVariable (UniformVariableImpl (s, l)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
265 |
{} |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
266 |
|
4242
3d4788bd1ddd
improve coherency between UniformVariable and RandomVariable, explain why the extra overload is necessary.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4240
diff
changeset
|
267 |
double UniformVariable::GetValue(void) const |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
268 |
{ |
4242
3d4788bd1ddd
improve coherency between UniformVariable and RandomVariable, explain why the extra overload is necessary.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4240
diff
changeset
|
269 |
return this->RandomVariable::GetValue (); |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
270 |
} |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
271 |
|
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
272 |
double UniformVariable::GetValue(double s, double l) |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
273 |
{ |
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
274 |
return ((UniformVariableImpl*)Peek())->GetValue(s,l); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
275 |
} |
442
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
276 |
|
4256
cc06d903ca09
Adding UniformVariable::GetInteger and using it to fix bug in RealRandomStream::GetNext.
Timo Bingmann <timo.bingmann@student.kit.edu>
parents:
4249
diff
changeset
|
277 |
uint32_t UniformVariable::GetInteger (uint32_t s, uint32_t l) |
cc06d903ca09
Adding UniformVariable::GetInteger and using it to fix bug in RealRandomStream::GetNext.
Timo Bingmann <timo.bingmann@student.kit.edu>
parents:
4249
diff
changeset
|
278 |
{ |
cc06d903ca09
Adding UniformVariable::GetInteger and using it to fix bug in RealRandomStream::GetNext.
Timo Bingmann <timo.bingmann@student.kit.edu>
parents:
4249
diff
changeset
|
279 |
NS_ASSERT(s <= l); |
cc06d903ca09
Adding UniformVariable::GetInteger and using it to fix bug in RealRandomStream::GetNext.
Timo Bingmann <timo.bingmann@student.kit.edu>
parents:
4249
diff
changeset
|
280 |
return static_cast<uint32_t>( GetValue(s, l+1) ); |
cc06d903ca09
Adding UniformVariable::GetInteger and using it to fix bug in RealRandomStream::GetNext.
Timo Bingmann <timo.bingmann@student.kit.edu>
parents:
4249
diff
changeset
|
281 |
} |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
282 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
283 |
//----------------------------------------------------------------------------- |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
284 |
//----------------------------------------------------------------------------- |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
285 |
// ConstantVariableImpl methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
286 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
287 |
class ConstantVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
288 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
289 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
290 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
291 |
* Construct a ConstantVariableImpl RNG that returns zero every sample |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
292 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
293 |
ConstantVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
294 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
295 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
296 |
* Construct a ConstantVariableImpl RNG that returns the specified value |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
297 |
* every sample. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
298 |
* \param c Unchanging value for this RNG. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
299 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
300 |
ConstantVariableImpl(double c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
301 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
302 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
303 |
ConstantVariableImpl(const ConstantVariableImpl& c) ; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
304 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
305 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
306 |
* \brief Specify a new constant RNG for this generator. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
307 |
* \param c New constant value for this RNG. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
308 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
309 |
void NewConstant(double c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
310 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
311 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
312 |
* \return The constant value specified |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
313 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
314 |
virtual double GetValue(); |
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
315 |
virtual uint32_t GetInteger(); |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
316 |
virtual RandomVariableBase* Copy(void) const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
317 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
318 |
double m_const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
319 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
320 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
321 |
ConstantVariableImpl::ConstantVariableImpl() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
322 |
: m_const(0) { } |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
323 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
324 |
ConstantVariableImpl::ConstantVariableImpl(double c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
325 |
: m_const(c) { }; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
326 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
327 |
ConstantVariableImpl::ConstantVariableImpl(const ConstantVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
328 |
: RandomVariableBase(c), m_const(c.m_const) { } |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
329 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
330 |
void ConstantVariableImpl::NewConstant(double c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
331 |
{ m_const = c;} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
332 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
333 |
double ConstantVariableImpl::GetValue() |
442
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
334 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
335 |
return m_const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
336 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
337 |
|
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
338 |
uint32_t ConstantVariableImpl::GetInteger() |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
339 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
340 |
return (uint32_t)m_const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
341 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
342 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
343 |
RandomVariableBase* ConstantVariableImpl::Copy() const |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
344 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
345 |
return new ConstantVariableImpl(*this); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
346 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
347 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
348 |
ConstantVariable::ConstantVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
349 |
: RandomVariable (ConstantVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
350 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
351 |
ConstantVariable::ConstantVariable(double c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
352 |
: RandomVariable (ConstantVariableImpl (c)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
353 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
354 |
void |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
355 |
ConstantVariable::SetConstant(double c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
356 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
357 |
*this = ConstantVariable (c); |
442
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
358 |
} |
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
359 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
360 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
361 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
362 |
// SequentialVariableImpl methods |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
363 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
364 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
365 |
class SequentialVariableImpl : public RandomVariableBase { |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
366 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
367 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
368 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
369 |
* \brief Constructor for the SequentialVariableImpl RNG. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
370 |
* |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
371 |
* The four parameters define the sequence. For example |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
372 |
* SequentialVariableImpl(0,5,1,2) creates a RNG that has the sequence |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
373 |
* 0, 0, 1, 1, 2, 2, 3, 3, 4, 4, 0, 0 ... |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
374 |
* \param f First value of the sequence. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
375 |
* \param l One more than the last value of the sequence. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
376 |
* \param i Increment between sequence values |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
377 |
* \param c Number of times each member of the sequence is repeated |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
378 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
379 |
SequentialVariableImpl(double f, double l, double i = 1, uint32_t c = 1); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
380 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
381 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
382 |
* \brief Constructor for the SequentialVariableImpl RNG. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
383 |
* |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
384 |
* Differs from the first only in that the increment parameter is a |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
385 |
* random variable |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
386 |
* \param f First value of the sequence. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
387 |
* \param l One more than the last value of the sequence. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
388 |
* \param i Reference to a RandomVariableBase for the sequence increment |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
389 |
* \param c Number of times each member of the sequence is repeated |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
390 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
391 |
SequentialVariableImpl(double f, double l, const RandomVariable& i, uint32_t c = 1); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
392 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
393 |
SequentialVariableImpl(const SequentialVariableImpl& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
394 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
395 |
~SequentialVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
396 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
397 |
* \return The next value in the Sequence |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
398 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
399 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
400 |
virtual RandomVariableBase* Copy(void) const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
401 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
402 |
double m_min; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
403 |
double m_max; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
404 |
RandomVariable m_increment; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
405 |
uint32_t m_consecutive; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
406 |
double m_current; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
407 |
uint32_t m_currentConsecutive; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
408 |
}; |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
409 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
410 |
SequentialVariableImpl::SequentialVariableImpl(double f, double l, double i, uint32_t c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
411 |
: m_min(f), m_max(l), m_increment(ConstantVariable(i)), m_consecutive(c), |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
412 |
m_current(f), m_currentConsecutive(0) |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
413 |
{} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
414 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
415 |
SequentialVariableImpl::SequentialVariableImpl(double f, double l, const RandomVariable& i, uint32_t c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
416 |
: m_min(f), m_max(l), m_increment(i), m_consecutive(c), |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
417 |
m_current(f), m_currentConsecutive(0) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
418 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
419 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
420 |
SequentialVariableImpl::SequentialVariableImpl(const SequentialVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
421 |
: RandomVariableBase(c), m_min(c.m_min), m_max(c.m_max), |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
422 |
m_increment(c.m_increment), m_consecutive(c.m_consecutive), |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
423 |
m_current(c.m_current), m_currentConsecutive(c.m_currentConsecutive) |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
424 |
{} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
425 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
426 |
SequentialVariableImpl::~SequentialVariableImpl() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
427 |
{} |
360
7bffd987426c
Fixed major memory bug and docs in RandomVariable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
346
diff
changeset
|
428 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
429 |
double SequentialVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
430 |
{ // Return a sequential series of values |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
431 |
double r = m_current; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
432 |
if (++m_currentConsecutive == m_consecutive) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
433 |
{ // Time to advance to next |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
434 |
m_currentConsecutive = 0; |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
435 |
m_current += m_increment.GetValue(); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
436 |
if (m_current >= m_max) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
437 |
m_current = m_min + (m_current - m_max); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
438 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
439 |
return r; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
440 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
441 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
442 |
RandomVariableBase* SequentialVariableImpl::Copy() const |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
443 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
444 |
return new SequentialVariableImpl(*this); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
445 |
} |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
446 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
447 |
SequentialVariable::SequentialVariable(double f, double l, double i, uint32_t c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
448 |
: RandomVariable (SequentialVariableImpl (f, l, i, c)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
449 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
450 |
SequentialVariable::SequentialVariable(double f, double l, const RandomVariable& i, uint32_t c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
451 |
: RandomVariable (SequentialVariableImpl (f, l, i, c)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
452 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
453 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
454 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
455 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
456 |
// ExponentialVariableImpl methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
457 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
458 |
class ExponentialVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
459 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
460 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
461 |
* Constructs an exponential random variable with a mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
462 |
* value of 1.0. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
463 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
464 |
ExponentialVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
465 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
466 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
467 |
* \brief Constructs an exponential random variable with a specified mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
468 |
* \param m Mean value for the random variable |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
469 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
470 |
explicit ExponentialVariableImpl(double m); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
471 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
472 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
473 |
* \brief Constructs an exponential random variable with spefified |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
474 |
* \brief mean and upper limit. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
475 |
* |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
476 |
* Since exponential distributions can theoretically return unbounded values, |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
477 |
* it is sometimes useful to specify a fixed upper limit. Note however when |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
478 |
* the upper limit is specified, the true mean of the distribution is |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
479 |
* slightly smaller than the mean value specified. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
480 |
* \param m Mean value of the random variable |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
481 |
* \param b Upper bound on returned values |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
482 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
483 |
ExponentialVariableImpl(double m, double b); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
484 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
485 |
ExponentialVariableImpl(const ExponentialVariableImpl& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
486 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
487 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
488 |
* \return A random value from this exponential distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
489 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
490 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
491 |
virtual RandomVariableBase* Copy(void) const; |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
492 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
493 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
494 |
double m_mean; // Mean value of RV |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
495 |
double m_bound; // Upper bound on value (if non-zero) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
496 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
497 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
498 |
ExponentialVariableImpl::ExponentialVariableImpl() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
499 |
: m_mean(1.0), m_bound(0) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
500 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
501 |
ExponentialVariableImpl::ExponentialVariableImpl(double m) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
502 |
: m_mean(m), m_bound(0) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
503 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
504 |
ExponentialVariableImpl::ExponentialVariableImpl(double m, double b) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
505 |
: m_mean(m), m_bound(b) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
506 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
507 |
ExponentialVariableImpl::ExponentialVariableImpl(const ExponentialVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
508 |
: RandomVariableBase(c), m_mean(c.m_mean), m_bound(c.m_bound) { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
509 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
510 |
double ExponentialVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
511 |
{ |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
512 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
513 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
514 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
515 |
} |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
516 |
while(1) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
517 |
{ |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
518 |
double r = -m_mean*log(m_generator->RandU01()); |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
519 |
if (m_bound == 0 || r <= m_bound) return r; |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
520 |
//otherwise, try again |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
521 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
522 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
523 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
524 |
RandomVariableBase* ExponentialVariableImpl::Copy() const |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
525 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
526 |
return new ExponentialVariableImpl(*this); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
527 |
} |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
528 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
529 |
ExponentialVariable::ExponentialVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
530 |
: RandomVariable (ExponentialVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
531 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
532 |
ExponentialVariable::ExponentialVariable(double m) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
533 |
: RandomVariable (ExponentialVariableImpl (m)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
534 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
535 |
ExponentialVariable::ExponentialVariable(double m, double b) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
536 |
: RandomVariable (ExponentialVariableImpl (m, b)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
537 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
538 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
539 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
540 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
541 |
// ParetoVariableImpl methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
542 |
class ParetoVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
543 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
544 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
545 |
* Constructs a pareto random variable with a mean of 1 and a shape |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
546 |
* parameter of 1.5 |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
547 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
548 |
ParetoVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
549 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
550 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
551 |
* Constructs a pareto random variable with specified mean and shape |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
552 |
* parameter of 1.5 |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
553 |
* \param m Mean value of the distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
554 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
555 |
explicit ParetoVariableImpl(double m); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
556 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
557 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
558 |
* Constructs a pareto random variable with the specified mean value and |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
559 |
* shape parameter. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
560 |
* \param m Mean value of the distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
561 |
* \param s Shape parameter for the distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
562 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
563 |
ParetoVariableImpl(double m, double s); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
564 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
565 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
566 |
* \brief Constructs a pareto random variable with the specified mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
567 |
* \brief value, shape (alpha), and upper bound. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
568 |
* |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
569 |
* Since pareto distributions can theoretically return unbounded values, |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
570 |
* it is sometimes useful to specify a fixed upper limit. Note however |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
571 |
* when the upper limit is specified, the true mean of the distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
572 |
* is slightly smaller than the mean value specified. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
573 |
* \param m Mean value |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
574 |
* \param s Shape parameter |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
575 |
* \param b Upper limit on returned values |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
576 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
577 |
ParetoVariableImpl(double m, double s, double b); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
578 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
579 |
ParetoVariableImpl(const ParetoVariableImpl& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
580 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
581 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
582 |
* \return A random value from this Pareto distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
583 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
584 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
585 |
virtual RandomVariableBase* Copy() const; |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
586 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
587 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
588 |
double m_mean; // Mean value of RV |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
589 |
double m_shape; // Shape parameter |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
590 |
double m_bound; // Upper bound on value (if non-zero) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
591 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
592 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
593 |
ParetoVariableImpl::ParetoVariableImpl() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
594 |
: m_mean(1.0), m_shape(1.5), m_bound(0) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
595 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
596 |
ParetoVariableImpl::ParetoVariableImpl(double m) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
597 |
: m_mean(m), m_shape(1.5), m_bound(0) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
598 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
599 |
ParetoVariableImpl::ParetoVariableImpl(double m, double s) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
600 |
: m_mean(m), m_shape(s), m_bound(0) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
601 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
602 |
ParetoVariableImpl::ParetoVariableImpl(double m, double s, double b) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
603 |
: m_mean(m), m_shape(s), m_bound(b) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
604 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
605 |
ParetoVariableImpl::ParetoVariableImpl(const ParetoVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
606 |
: RandomVariableBase(c), m_mean(c.m_mean), m_shape(c.m_shape), |
360
7bffd987426c
Fixed major memory bug and docs in RandomVariable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
346
diff
changeset
|
607 |
m_bound(c.m_bound) { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
608 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
609 |
double ParetoVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
610 |
{ |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
611 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
612 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
613 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
614 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
615 |
double scale = m_mean * ( m_shape - 1.0) / m_shape; |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
616 |
while(1) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
617 |
{ |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
618 |
double r = (scale * ( 1.0 / pow(m_generator->RandU01(), 1.0 / m_shape))); |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
619 |
if (m_bound == 0 || r <= m_bound) return r; |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
620 |
//otherwise, try again |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
621 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
622 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
623 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
624 |
RandomVariableBase* ParetoVariableImpl::Copy() const |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
625 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
626 |
return new ParetoVariableImpl(*this); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
627 |
} |
442
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
628 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
629 |
ParetoVariable::ParetoVariable () |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
630 |
: RandomVariable (ParetoVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
631 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
632 |
ParetoVariable::ParetoVariable(double m) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
633 |
: RandomVariable (ParetoVariableImpl (m)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
634 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
635 |
ParetoVariable::ParetoVariable(double m, double s) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
636 |
: RandomVariable (ParetoVariableImpl (m, s)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
637 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
638 |
ParetoVariable::ParetoVariable(double m, double s, double b) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
639 |
: RandomVariable (ParetoVariableImpl (m, s, b)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
640 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
641 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
642 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
643 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
644 |
// WeibullVariableImpl methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
645 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
646 |
class WeibullVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
647 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
648 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
649 |
* Constructs a weibull random variable with a mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
650 |
* value of 1.0 and a shape (alpha) parameter of 1 |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
651 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
652 |
WeibullVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
653 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
654 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
655 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
656 |
* Constructs a weibull random variable with the specified mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
657 |
* value and a shape (alpha) parameter of 1.5. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
658 |
* \param m mean value of the distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
659 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
660 |
WeibullVariableImpl(double m) ; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
661 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
662 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
663 |
* Constructs a weibull random variable with the specified mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
664 |
* value and a shape (alpha). |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
665 |
* \param m Mean value for the distribution. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
666 |
* \param s Shape (alpha) parameter for the distribution. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
667 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
668 |
WeibullVariableImpl(double m, double s); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
669 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
670 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
671 |
* \brief Constructs a weibull random variable with the specified mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
672 |
* \brief value, shape (alpha), and upper bound. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
673 |
* Since WeibullVariableImpl distributions can theoretically return unbounded values, |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
674 |
* it is sometimes usefull to specify a fixed upper limit. Note however |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
675 |
* that when the upper limit is specified, the true mean of the distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
676 |
* is slightly smaller than the mean value specified. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
677 |
* \param m Mean value for the distribution. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
678 |
* \param s Shape (alpha) parameter for the distribution. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
679 |
* \param b Upper limit on returned values |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
680 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
681 |
WeibullVariableImpl(double m, double s, double b); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
682 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
683 |
WeibullVariableImpl(const WeibullVariableImpl& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
684 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
685 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
686 |
* \return A random value from this Weibull distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
687 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
688 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
689 |
virtual RandomVariableBase* Copy(void) const; |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
690 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
691 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
692 |
double m_mean; // Mean value of RV |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
693 |
double m_alpha; // Shape parameter |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
694 |
double m_bound; // Upper bound on value (if non-zero) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
695 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
696 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
697 |
WeibullVariableImpl::WeibullVariableImpl() : m_mean(1.0), m_alpha(1), m_bound(0) { } |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
698 |
WeibullVariableImpl::WeibullVariableImpl(double m) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
699 |
: m_mean(m), m_alpha(1), m_bound(0) { } |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
700 |
WeibullVariableImpl::WeibullVariableImpl(double m, double s) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
701 |
: m_mean(m), m_alpha(s), m_bound(0) { } |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
702 |
WeibullVariableImpl::WeibullVariableImpl(double m, double s, double b) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
703 |
: m_mean(m), m_alpha(s), m_bound(b) { }; |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
704 |
WeibullVariableImpl::WeibullVariableImpl(const WeibullVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
705 |
: RandomVariableBase(c), m_mean(c.m_mean), m_alpha(c.m_alpha), |
360
7bffd987426c
Fixed major memory bug and docs in RandomVariable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
346
diff
changeset
|
706 |
m_bound(c.m_bound) { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
707 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
708 |
double WeibullVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
709 |
{ |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
710 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
711 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
712 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
713 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
714 |
double exponent = 1.0 / m_alpha; |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
715 |
while(1) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
716 |
{ |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
717 |
double r = m_mean * pow( -log(m_generator->RandU01()), exponent); |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
718 |
if (m_bound == 0 || r <= m_bound) return r; |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
719 |
//otherwise, try again |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
720 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
721 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
722 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
723 |
RandomVariableBase* WeibullVariableImpl::Copy() const |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
724 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
725 |
return new WeibullVariableImpl(*this); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
726 |
} |
442
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
727 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
728 |
WeibullVariable::WeibullVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
729 |
: RandomVariable (WeibullVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
730 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
731 |
WeibullVariable::WeibullVariable(double m) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
732 |
: RandomVariable (WeibullVariableImpl (m)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
733 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
734 |
WeibullVariable::WeibullVariable(double m, double s) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
735 |
: RandomVariable (WeibullVariableImpl (m, s)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
736 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
737 |
WeibullVariable::WeibullVariable(double m, double s, double b) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
738 |
: RandomVariable (WeibullVariableImpl (m, s, b)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
739 |
{} |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
740 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
741 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
742 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
743 |
// NormalVariableImpl methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
744 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
745 |
class NormalVariableImpl : public RandomVariableBase { // Normally Distributed random var |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
746 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
747 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
748 |
static const double INFINITE_VALUE; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
749 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
750 |
* Constructs an normal random variable with a mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
751 |
* value of 0 and variance of 1. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
752 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
753 |
NormalVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
754 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
755 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
756 |
* \brief Construct a normal random variable with specified mean and variance |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
757 |
* \param m Mean value |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
758 |
* \param v Variance |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
759 |
* \param b Bound. The NormalVariableImpl is bounded within +-bound of the mean. |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
760 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
761 |
NormalVariableImpl(double m, double v, double b = INFINITE_VALUE); |
442
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
762 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
763 |
NormalVariableImpl(const NormalVariableImpl& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
764 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
765 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
766 |
* \return A value from this normal distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
767 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
768 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
769 |
virtual RandomVariableBase* Copy(void) const; |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
770 |
|
4249
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
771 |
double GetMean (void) const; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
772 |
double GetVariance (void) const; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
773 |
double GetBound (void) const; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
774 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
775 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
776 |
double m_mean; // Mean value of RV |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
777 |
double m_variance; // Mean value of RV |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
778 |
double m_bound; // Bound on value's difference from the mean (absolute value) |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
779 |
bool m_nextValid; // True if next valid |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
780 |
double m_next; // The algorithm produces two values at a time |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
781 |
static bool m_static_nextValid; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
782 |
static double m_static_next; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
783 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
784 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
785 |
bool NormalVariableImpl::m_static_nextValid = false; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
786 |
double NormalVariableImpl::m_static_next; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
787 |
const double NormalVariableImpl::INFINITE_VALUE = 1e307; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
788 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
789 |
NormalVariableImpl::NormalVariableImpl() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
790 |
: m_mean(0.0), m_variance(1.0), m_bound(INFINITE_VALUE), m_nextValid(false){} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
791 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
792 |
NormalVariableImpl::NormalVariableImpl(double m, double v, double b/*=INFINITE_VALUE*/) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
793 |
: m_mean(m), m_variance(v), m_bound(b), m_nextValid(false) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
794 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
795 |
NormalVariableImpl::NormalVariableImpl(const NormalVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
796 |
: RandomVariableBase(c), m_mean(c.m_mean), m_variance(c.m_variance), |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
797 |
m_bound(c.m_bound) { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
798 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
799 |
double NormalVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
800 |
{ |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
801 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
802 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
803 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
804 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
805 |
if (m_nextValid) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
806 |
{ // use previously generated |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
807 |
m_nextValid = false; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
808 |
return m_next; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
809 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
810 |
while(1) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
811 |
{ // See Simulation Modeling and Analysis p. 466 (Averill Law) |
3292 | 812 |
// for algorithm; basically a Box-Muller transform: |
813 |
// http://en.wikipedia.org/wiki/Box-Muller_transform |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
814 |
double u1 = m_generator->RandU01(); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
815 |
double u2 = m_generator->RandU01();; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
816 |
double v1 = 2 * u1 - 1; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
817 |
double v2 = 2 * u2 - 1; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
818 |
double w = v1 * v1 + v2 * v2; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
819 |
if (w <= 1.0) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
820 |
{ // Got good pair |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
821 |
double y = sqrt((-2 * log(w))/w); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
822 |
m_next = m_mean + v2 * y * sqrt(m_variance); |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
823 |
//if next is in bounds, it is valid |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
824 |
m_nextValid = fabs(m_next-m_mean) <= m_bound; |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
825 |
double x1 = m_mean + v1 * y * sqrt(m_variance); |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
826 |
//if x1 is in bounds, return it |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
827 |
if (fabs(x1-m_mean) <= m_bound) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
828 |
{ |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
829 |
return x1; |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
830 |
} |
4227
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
831 |
//otherwise try and return m_next if it is valid |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
832 |
else if (m_nextValid) |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
833 |
{ |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
834 |
m_nextValid = false; |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
835 |
return m_next; |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
836 |
} |
b30027eeb387
Michelle's bounds patch
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4224
diff
changeset
|
837 |
//otherwise, just run this loop again |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
838 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
839 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
840 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
841 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
842 |
RandomVariableBase* NormalVariableImpl::Copy() const |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
843 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
844 |
return new NormalVariableImpl(*this); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
845 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
846 |
|
4249
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
847 |
double |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
848 |
NormalVariableImpl::GetMean (void) const |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
849 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
850 |
return m_mean; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
851 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
852 |
|
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
853 |
double |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
854 |
NormalVariableImpl::GetVariance (void) const |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
855 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
856 |
return m_variance; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
857 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
858 |
|
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
859 |
double |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
860 |
NormalVariableImpl::GetBound (void) const |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
861 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
862 |
return m_bound; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
863 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
864 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
865 |
NormalVariable::NormalVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
866 |
: RandomVariable (NormalVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
867 |
{} |
4223
86a97665dcb0
Fix old code that crept in
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4218
diff
changeset
|
868 |
NormalVariable::NormalVariable(double m, double v) |
86a97665dcb0
Fix old code that crept in
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4218
diff
changeset
|
869 |
: RandomVariable (NormalVariableImpl (m, v)) |
86a97665dcb0
Fix old code that crept in
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4218
diff
changeset
|
870 |
{} |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
871 |
NormalVariable::NormalVariable(double m, double v, double b) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
872 |
: RandomVariable (NormalVariableImpl (m, v, b)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
873 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
874 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
875 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
876 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
877 |
class EmpiricalVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
878 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
879 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
880 |
* Constructor for the EmpiricalVariableImpl random variables. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
881 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
882 |
explicit EmpiricalVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
883 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
884 |
virtual ~EmpiricalVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
885 |
EmpiricalVariableImpl(const EmpiricalVariableImpl& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
886 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
887 |
* \return A value from this empirical distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
888 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
889 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
890 |
virtual RandomVariableBase* Copy(void) const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
891 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
892 |
* \brief Specifies a point in the empirical distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
893 |
* \param v The function value for this point |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
894 |
* \param c Probability that the function is less than or equal to v |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
895 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
896 |
virtual void CDF(double v, double c); // Value, prob <= Value |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
897 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
898 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
899 |
class ValueCDF { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
900 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
901 |
ValueCDF(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
902 |
ValueCDF(double v, double c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
903 |
ValueCDF(const ValueCDF& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
904 |
double value; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
905 |
double cdf; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
906 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
907 |
virtual void Validate(); // Insure non-decreasing emiprical values |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
908 |
virtual double Interpolate(double, double, double, double, double); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
909 |
bool validated; // True if non-decreasing validated |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
910 |
std::vector<ValueCDF> emp; // Empicical CDF |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
911 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
912 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
913 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
914 |
// ValueCDF methods |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
915 |
EmpiricalVariableImpl::ValueCDF::ValueCDF() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
916 |
: value(0.0), cdf(0.0){ } |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
917 |
EmpiricalVariableImpl::ValueCDF::ValueCDF(double v, double c) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
918 |
: value(v), cdf(c) { } |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
919 |
EmpiricalVariableImpl::ValueCDF::ValueCDF(const ValueCDF& c) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
920 |
: value(c.value), cdf(c.cdf) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
921 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
922 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
923 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
924 |
// EmpiricalVariableImpl methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
925 |
EmpiricalVariableImpl::EmpiricalVariableImpl() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
926 |
: validated(false) { } |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
927 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
928 |
EmpiricalVariableImpl::EmpiricalVariableImpl(const EmpiricalVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
929 |
: RandomVariableBase(c), validated(c.validated), emp(c.emp) { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
930 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
931 |
EmpiricalVariableImpl::~EmpiricalVariableImpl() { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
932 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
933 |
double EmpiricalVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
934 |
{ // Return a value from the empirical distribution |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
935 |
// This code based (loosely) on code by Bruce Mah (Thanks Bruce!) |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
936 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
937 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
938 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
939 |
} |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
940 |
if (emp.size() == 0) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
941 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
942 |
return 0.0; // HuH? No empirical data |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
943 |
} |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
944 |
if (!validated) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
945 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
946 |
Validate(); // Insure in non-decreasing |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
947 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
948 |
double r = m_generator->RandU01(); |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
949 |
if (r <= emp.front().cdf) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
950 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
951 |
return emp.front().value; // Less than first |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
952 |
} |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
953 |
if (r >= emp.back().cdf) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
954 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
955 |
return emp.back().value; // Greater than last |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
956 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
957 |
// Binary search |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
958 |
std::vector<ValueCDF>::size_type bottom = 0; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
959 |
std::vector<ValueCDF>::size_type top = emp.size() - 1; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
960 |
while(1) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
961 |
{ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
962 |
std::vector<ValueCDF>::size_type c = (top + bottom) / 2; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
963 |
if (r >= emp[c].cdf && r < emp[c+1].cdf) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
964 |
{ // Found it |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
965 |
return Interpolate(emp[c].cdf, emp[c+1].cdf, |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
966 |
emp[c].value, emp[c+1].value, |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
967 |
r); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
968 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
969 |
// Not here, adjust bounds |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
970 |
if (r < emp[c].cdf) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
971 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
972 |
top = c - 1; |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
973 |
} |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
974 |
else |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
975 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
976 |
bottom = c + 1; |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
977 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
978 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
979 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
980 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
981 |
RandomVariableBase* EmpiricalVariableImpl::Copy() const |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
982 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
983 |
return new EmpiricalVariableImpl(*this); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
984 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
985 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
986 |
void EmpiricalVariableImpl::CDF(double v, double c) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
987 |
{ // Add a new empirical datapoint to the empirical cdf |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
988 |
// NOTE. These MUST be inserted in non-decreasing order |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
989 |
emp.push_back(ValueCDF(v, c)); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
990 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
991 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
992 |
void EmpiricalVariableImpl::Validate() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
993 |
{ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
994 |
ValueCDF prior; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
995 |
for (std::vector<ValueCDF>::size_type i = 0; i < emp.size(); ++i) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
996 |
{ |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
997 |
ValueCDF& current = emp[i]; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
998 |
if (current.value < prior.value || current.cdf < prior.cdf) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
999 |
{ // Error |
1020
ebcb38b7ff7d
In some places, print debug messages / errors to stderr, not stdout.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
942
diff
changeset
|
1000 |
cerr << "Empirical Dist error," |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1001 |
<< " current value " << current.value |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1002 |
<< " prior value " << prior.value |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1003 |
<< " current cdf " << current.cdf |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1004 |
<< " prior cdf " << prior.cdf << endl; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1005 |
NS_FATAL_ERROR("Empirical Dist error"); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1006 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1007 |
prior = current; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1008 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1009 |
validated = true; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1010 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1011 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1012 |
double EmpiricalVariableImpl::Interpolate(double c1, double c2, |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1013 |
double v1, double v2, double r) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1014 |
{ // Interpolate random value in range [v1..v2) based on [c1 .. r .. c2) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1015 |
return (v1 + ((v2 - v1) / (c2 - c1)) * (r - c1)); |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1016 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1017 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1018 |
EmpiricalVariable::EmpiricalVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1019 |
: RandomVariable (EmpiricalVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1020 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1021 |
EmpiricalVariable::EmpiricalVariable (const RandomVariableBase &variable) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1022 |
: RandomVariable (variable) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1023 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1024 |
void |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1025 |
EmpiricalVariable::CDF(double v, double c) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1026 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1027 |
EmpiricalVariableImpl *impl = dynamic_cast<EmpiricalVariableImpl *> (Peek ()); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1028 |
NS_ASSERT (impl); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1029 |
impl->CDF (v, c); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1030 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1031 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1032 |
|
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1033 |
//----------------------------------------------------------------------------- |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1034 |
//----------------------------------------------------------------------------- |
2965
4b28e9740e3b
get rid of Attribute class. Use AttributeValue subclasses directly.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2947
diff
changeset
|
1035 |
// IntegerValue EmpiricalVariableImpl methods |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1036 |
class IntEmpiricalVariableImpl : public EmpiricalVariableImpl { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1037 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1038 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1039 |
IntEmpiricalVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1040 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1041 |
virtual RandomVariableBase* Copy(void) const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1042 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1043 |
* \return An integer value from this empirical distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1044 |
*/ |
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
1045 |
virtual uint32_t GetInteger(); |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1046 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1047 |
virtual double Interpolate(double, double, double, double, double); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1048 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1049 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1050 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1051 |
IntEmpiricalVariableImpl::IntEmpiricalVariableImpl() { } |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1052 |
|
2439
4a0b22a3f5fa
UintValue -> Uinteger, IntValue -> Integer
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2436
diff
changeset
|
1053 |
uint32_t IntEmpiricalVariableImpl::GetInteger() |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1054 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1055 |
return (uint32_t)GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1056 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1057 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1058 |
RandomVariableBase* IntEmpiricalVariableImpl::Copy() const |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1059 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1060 |
return new IntEmpiricalVariableImpl(*this); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1061 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1062 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1063 |
double IntEmpiricalVariableImpl::Interpolate(double c1, double c2, |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1064 |
double v1, double v2, double r) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1065 |
{ // Interpolate random value in range [v1..v2) based on [c1 .. r .. c2) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1066 |
return ceil(v1 + ((v2 - v1) / (c2 - c1)) * (r - c1)); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1067 |
} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1068 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1069 |
IntEmpiricalVariable::IntEmpiricalVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1070 |
: EmpiricalVariable (IntEmpiricalVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1071 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1072 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1073 |
//----------------------------------------------------------------------------- |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1074 |
//----------------------------------------------------------------------------- |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1075 |
// DeterministicVariableImpl |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1076 |
class DeterministicVariableImpl : public RandomVariableBase |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1077 |
{ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1078 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1079 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1080 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1081 |
* \brief Constructor |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1082 |
* |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1083 |
* Creates a generator that returns successive elements of the d array |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1084 |
* on successive calls to ::Value(). Note that the d pointer is copied |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1085 |
* for use by the generator (shallow-copy), not its contents, so the |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1086 |
* contents of the array d points to have to remain unchanged for the use |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1087 |
* of DeterministicVariableImpl to be meaningful. |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1088 |
* \param d Pointer to array of random values to return in sequence |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1089 |
* \param c Number of values in the array |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1090 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1091 |
explicit DeterministicVariableImpl(double* d, uint32_t c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1092 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1093 |
virtual ~DeterministicVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1094 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1095 |
* \return The next value in the deterministic sequence |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1096 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1097 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1098 |
virtual RandomVariableBase* Copy(void) const; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1099 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1100 |
uint32_t count; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1101 |
uint32_t next; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1102 |
double* data; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1103 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1104 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1105 |
DeterministicVariableImpl::DeterministicVariableImpl(double* d, uint32_t c) |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1106 |
: count(c), next(c), data(d) |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1107 |
{ // Nothing else needed |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1108 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1109 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1110 |
DeterministicVariableImpl::~DeterministicVariableImpl() { } |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1111 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1112 |
double DeterministicVariableImpl::GetValue() |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1113 |
{ |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1114 |
if (next == count) |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1115 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1116 |
next = 0; |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1117 |
} |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1118 |
return data[next++]; |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1119 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1120 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1121 |
RandomVariableBase* DeterministicVariableImpl::Copy() const |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1122 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1123 |
return new DeterministicVariableImpl(*this); |
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1124 |
} |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1125 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1126 |
DeterministicVariable::DeterministicVariable(double* d, uint32_t c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1127 |
: RandomVariable (DeterministicVariableImpl (d, c)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1128 |
{} |
438
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1129 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1130 |
//----------------------------------------------------------------------------- |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1131 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1132 |
// LogNormalVariableImpl |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1133 |
class LogNormalVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1134 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1135 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1136 |
* \param mu mu parameter of the lognormal distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1137 |
* \param sigma sigma parameter of the lognormal distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1138 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1139 |
LogNormalVariableImpl (double mu, double sigma); |
438
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1140 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1141 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1142 |
* \return A random value from this distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1143 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1144 |
virtual double GetValue (); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1145 |
virtual RandomVariableBase* Copy(void) const; |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
1146 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1147 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1148 |
double m_mu; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1149 |
double m_sigma; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1150 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1151 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1152 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1153 |
RandomVariableBase* LogNormalVariableImpl::Copy () const |
438
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1154 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1155 |
return new LogNormalVariableImpl (m_mu, m_sigma); |
438
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1156 |
} |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1157 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1158 |
LogNormalVariableImpl::LogNormalVariableImpl (double mu, double sigma) |
442
96d3e7dc8bb2
Added static RNGs, like ExponentialVariable::GetSingleValue(mean)
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
438
diff
changeset
|
1159 |
:m_mu(mu), m_sigma(sigma) |
438
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1160 |
{ |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1161 |
} |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1162 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1163 |
// The code from this function was adapted from the GNU Scientific |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1164 |
// Library 1.8: |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1165 |
/* randist/lognormal.c |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1166 |
* |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1167 |
* Copyright (C) 1996, 1997, 1998, 1999, 2000 James Theiler, Brian Gough |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1168 |
* |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1169 |
* This program is free software; you can redistribute it and/or modify |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1170 |
* it under the terms of the GNU General Public License as published by |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1171 |
* the Free Software Foundation; either version 2 of the License, or (at |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1172 |
* your option) any later version. |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1173 |
* |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1174 |
* This program is distributed in the hope that it will be useful, but |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1175 |
* WITHOUT ANY WARRANTY; without even the implied warranty of |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1176 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1177 |
* General Public License for more details. |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1178 |
* |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1179 |
* You should have received a copy of the GNU General Public License |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1180 |
* along with this program; if not, write to the Free Software |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1181 |
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1182 |
*/ |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1183 |
/* The lognormal distribution has the form |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1184 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1185 |
p(x) dx = 1/(x * sqrt(2 pi sigma^2)) exp(-(ln(x) - zeta)^2/2 sigma^2) dx |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1186 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1187 |
for x > 0. Lognormal random numbers are the exponentials of |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1188 |
gaussian random numbers */ |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1189 |
double |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1190 |
LogNormalVariableImpl::GetValue () |
438
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1191 |
{ |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
1192 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1193 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1194 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1195 |
} |
438
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1196 |
double u, v, r2, normal, z; |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1197 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1198 |
do |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1199 |
{ |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1200 |
/* choose x,y in uniform square (-1,-1) to (+1,+1) */ |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1201 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1202 |
u = -1 + 2 * m_generator->RandU01 (); |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1203 |
v = -1 + 2 * m_generator->RandU01 (); |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1204 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1205 |
/* see if it is in the unit circle */ |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1206 |
r2 = u * u + v * v; |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1207 |
} |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1208 |
while (r2 > 1.0 || r2 == 0); |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1209 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1210 |
normal = u * sqrt (-2.0 * log (r2) / r2); |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1211 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1212 |
z = exp (m_sigma * normal + m_mu); |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1213 |
|
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1214 |
return z; |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1215 |
} |
515bb5663cf3
Randomvariable Lognormal added
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
423
diff
changeset
|
1216 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1217 |
LogNormalVariable::LogNormalVariable (double mu, double sigma) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1218 |
: RandomVariable (LogNormalVariableImpl (mu, sigma)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1219 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1220 |
|
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1221 |
//----------------------------------------------------------------------------- |
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1222 |
//----------------------------------------------------------------------------- |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1223 |
// TriangularVariableImpl methods |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1224 |
class TriangularVariableImpl : public RandomVariableBase { |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1225 |
public: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1226 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1227 |
* Creates a triangle distribution random number generator in the |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1228 |
* range [0.0 .. 1.0), with mean of 0.5 |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1229 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1230 |
TriangularVariableImpl(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1231 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1232 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1233 |
* Creates a triangle distribution random number generator with the specified |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1234 |
* range |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1235 |
* \param s Low end of the range |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1236 |
* \param l High end of the range |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1237 |
* \param mean mean of the distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1238 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1239 |
TriangularVariableImpl(double s, double l, double mean); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1240 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1241 |
TriangularVariableImpl(const TriangularVariableImpl& c); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1242 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1243 |
/** |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1244 |
* \return A value from this distribution |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1245 |
*/ |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1246 |
virtual double GetValue(); |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1247 |
virtual RandomVariableBase* Copy(void) const; |
4218
debf1a8a96d3
Proposed RNG API changes
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
3679
diff
changeset
|
1248 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1249 |
private: |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1250 |
double m_min; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1251 |
double m_max; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1252 |
double m_mode; //easier to work with the mode internally instead of the mean |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1253 |
//they are related by the simple: mean = (min+max+mode)/3 |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1254 |
}; |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1255 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1256 |
TriangularVariableImpl::TriangularVariableImpl() |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1257 |
: m_min(0), m_max(1), m_mode(0.5) { } |
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1258 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1259 |
TriangularVariableImpl::TriangularVariableImpl(double s, double l, double mean) |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1260 |
: m_min(s), m_max(l), m_mode(3.0*mean-s-l) { } |
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1261 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1262 |
TriangularVariableImpl::TriangularVariableImpl(const TriangularVariableImpl& c) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1263 |
: RandomVariableBase(c), m_min(c.m_min), m_max(c.m_max), m_mode(c.m_mode) { } |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1264 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1265 |
double TriangularVariableImpl::GetValue() |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1266 |
{ |
1806
34e115600ef2
Fixed randomvariable bugs and modified unit tests so that there is no non-deterministic failure on the tests
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
1020
diff
changeset
|
1267 |
if(!m_generator) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1268 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1269 |
m_generator = new RngStream(); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1270 |
} |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1271 |
double u = m_generator->RandU01(); |
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1272 |
if(u <= (m_mode - m_min) / (m_max - m_min) ) |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1273 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1274 |
return m_min + sqrt(u * (m_max - m_min) * (m_mode - m_min) ); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1275 |
} |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1276 |
else |
4224
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1277 |
{ |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1278 |
return m_max - sqrt( (1-u) * (m_max - m_min) * (m_max - m_mode) ); |
40e5d5af3c71
Commit Mathieu's coding style cleanup of random-variable
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
4223
diff
changeset
|
1279 |
} |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1280 |
} |
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1281 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1282 |
RandomVariableBase* TriangularVariableImpl::Copy() const |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1283 |
{ |
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1284 |
return new TriangularVariableImpl(*this); |
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1285 |
} |
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1286 |
|
2336
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1287 |
TriangularVariable::TriangularVariable() |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1288 |
: RandomVariable (TriangularVariableImpl ()) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1289 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1290 |
TriangularVariable::TriangularVariable(double s, double l, double mean) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1291 |
: RandomVariable (TriangularVariableImpl (s,l,mean)) |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1292 |
{} |
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1293 |
|
28ce210b91bb
bug 133: automate RandomVariable memory management.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
1811
diff
changeset
|
1294 |
|
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1295 |
std::ostream &operator << (std::ostream &os, const RandomVariable &var) |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1296 |
{ |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1297 |
RandomVariableBase *base = var.Peek (); |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1298 |
ConstantVariableImpl *constant = dynamic_cast<ConstantVariableImpl *> (base); |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1299 |
if (constant != 0) |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1300 |
{ |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1301 |
os << "Constant:" << constant->GetValue (); |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1302 |
return os; |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1303 |
} |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1304 |
UniformVariableImpl *uniform = dynamic_cast<UniformVariableImpl *> (base); |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1305 |
if (uniform != 0) |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1306 |
{ |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1307 |
os << "Uniform:" << uniform->GetMin () << ":" << uniform->GetMax (); |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1308 |
return os; |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1309 |
} |
4249
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1310 |
NormalVariableImpl *normal = dynamic_cast<NormalVariableImpl *> (base); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1311 |
if (normal != 0) |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1312 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1313 |
os << "Normal:" << normal->GetMean () << ":" << normal->GetVariance (); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1314 |
double bound = normal->GetBound (); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1315 |
if (bound != NormalVariableImpl::INFINITE_VALUE) |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1316 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1317 |
os << ":" << bound; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1318 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1319 |
return os; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1320 |
} |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1321 |
// XXX: support other distributions |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1322 |
os.setstate (std::ios_base::badbit); |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1323 |
return os; |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1324 |
} |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1325 |
std::istream &operator >> (std::istream &is, RandomVariable &var) |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1326 |
{ |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1327 |
std::string value; |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1328 |
is >> value; |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1329 |
std::string::size_type tmp; |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1330 |
tmp = value.find (":"); |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1331 |
if (tmp == std::string::npos) |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1332 |
{ |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1333 |
is.setstate (std::ios_base::badbit); |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1334 |
return is; |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1335 |
} |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1336 |
std::string type = value.substr (0, tmp); |
2628
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1337 |
value = value.substr (tmp + 1, value.npos); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1338 |
if (type == "Constant") |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1339 |
{ |
2628
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1340 |
istringstream iss (value); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1341 |
double constant; |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1342 |
iss >> constant; |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1343 |
var = ConstantVariable (constant); |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1344 |
} |
2628
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1345 |
else if (type == "Uniform") |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1346 |
{ |
2628
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1347 |
if (value.size () == 0) |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1348 |
{ |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1349 |
var = UniformVariable (); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1350 |
} |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1351 |
else |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1352 |
{ |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1353 |
tmp = value.find (":"); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1354 |
if (tmp == value.npos) |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1355 |
{ |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1356 |
NS_FATAL_ERROR ("bad Uniform value: " << value); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1357 |
} |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1358 |
istringstream issA (value.substr (0, tmp)); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1359 |
istringstream issB (value.substr (tmp + 1, value.npos)); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1360 |
double a, b; |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1361 |
issA >> a; |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1362 |
issB >> b; |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1363 |
var = UniformVariable (a, b); |
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1364 |
} |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1365 |
} |
4249
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1366 |
else if (type == "Normal") |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1367 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1368 |
if (value.size () == 0) |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1369 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1370 |
var = NormalVariable (); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1371 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1372 |
else |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1373 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1374 |
tmp = value.find (":"); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1375 |
if (tmp == value.npos) |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1376 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1377 |
NS_FATAL_ERROR ("bad Normal value: " << value); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1378 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1379 |
std::string::size_type tmp2; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1380 |
std::string sub = value.substr (tmp + 1, value.npos); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1381 |
tmp2 = sub.find (":"); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1382 |
if (tmp2 == value.npos) |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1383 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1384 |
istringstream issA (value.substr (0, tmp)); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1385 |
istringstream issB (sub); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1386 |
double a, b; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1387 |
issA >> a; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1388 |
issB >> b; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1389 |
var = NormalVariable (a, b); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1390 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1391 |
else |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1392 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1393 |
istringstream issA (value.substr (0, tmp)); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1394 |
istringstream issB (sub.substr (0, tmp2)); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1395 |
istringstream issC (sub.substr (tmp2 + 1, value.npos)); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1396 |
double a, b, c; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1397 |
issA >> a; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1398 |
issB >> b; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1399 |
issC >> c; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1400 |
var = NormalVariable (a, b, c); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1401 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1402 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1403 |
} |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1404 |
else |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1405 |
{ |
2628
88cdcd394491
Fix deserialization for Constant and Uniform random variables.
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
2445
diff
changeset
|
1406 |
NS_FATAL_ERROR ("RandomVariable deserialization not implemented for " << type); |
2409
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1407 |
// XXX: support other distributions. |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1408 |
} |
8744ec60e8e6
convert RandomVariable, Rectangle, Vector, and Time to the new Class Helper
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2384
diff
changeset
|
1409 |
return is; |
2384
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1410 |
} |
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1411 |
|
500ada6a4874
add Value support to RandomVariable.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
2336
diff
changeset
|
1412 |
|
942
b65368bf40f9
Added Triangularly distributed random variables
Raj Bhattacharjea <raj.b@gatech.edu>
parents:
665
diff
changeset
|
1413 |
|
346
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1414 |
}//namespace ns3 |
4a76f247e7dc
Add random number files, base class Application
Tom Henderson <tomh@tomh.org>
parents:
diff
changeset
|
1415 |
|
4243
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1416 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1417 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1418 |
#ifdef RUN_SELF_TESTS |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1419 |
#include "test.h" |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1420 |
#include <vector> |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1421 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1422 |
namespace ns3 { |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1423 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1424 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1425 |
class RandomVariableTest : public Test |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1426 |
{ |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1427 |
public: |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1428 |
RandomVariableTest () : Test ("RandomVariable") {} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1429 |
virtual bool RunTests (void) |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1430 |
{ |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1431 |
bool result = true; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1432 |
const double desired_mean = 1.0; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1433 |
const double desired_stddev = 1.0; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1434 |
double tmp = log (1 + (desired_stddev/desired_mean)*(desired_stddev/desired_mean)); |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1435 |
double sigma = sqrt (tmp); |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1436 |
double mu = log (desired_mean) - 0.5*tmp; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1437 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1438 |
// Test a custom lognormal instance |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1439 |
{ |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1440 |
LogNormalVariable lognormal (mu, sigma); |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1441 |
vector<double> samples; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1442 |
const int NSAMPLES = 10000; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1443 |
double sum = 0; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1444 |
for (int n = NSAMPLES; n; --n) |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1445 |
{ |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1446 |
double value = lognormal.GetValue (); |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1447 |
sum += value; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1448 |
samples.push_back (value); |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1449 |
} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1450 |
double obtained_mean = sum / NSAMPLES; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1451 |
sum = 0; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1452 |
for (vector<double>::iterator iter = samples.begin (); iter != samples.end (); iter++) |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1453 |
{ |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1454 |
double tmp = (*iter - obtained_mean); |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1455 |
sum += tmp*tmp; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1456 |
} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1457 |
double obtained_stddev = sqrt (sum / (NSAMPLES - 1)); |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1458 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1459 |
if (not (obtained_mean/desired_mean > 0.90 and obtained_mean/desired_mean < 1.10)) |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1460 |
{ |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1461 |
result = false; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1462 |
Failure () << "Obtained lognormal mean value " << obtained_mean << ", expected " << desired_mean << std::endl; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1463 |
} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1464 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1465 |
if (not (obtained_stddev/desired_stddev > 0.90 and obtained_stddev/desired_stddev < 1.10)) |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1466 |
{ |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1467 |
result = false; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1468 |
Failure () << "Obtained lognormal stddev value " << obtained_stddev << |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1469 |
", expected " << desired_stddev << std::endl; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1470 |
} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1471 |
} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1472 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1473 |
// Test attribute serialization |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1474 |
{ |
4249
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1475 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1476 |
RandomVariableValue val; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1477 |
val.DeserializeFromString ("Uniform:0.1:0.2", MakeRandomVariableChecker ()); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1478 |
RandomVariable rng = val.Get (); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1479 |
NS_TEST_ASSERT_EQUAL (val.SerializeToString (MakeRandomVariableChecker ()), "Uniform:0.1:0.2"); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1480 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1481 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1482 |
RandomVariableValue val; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1483 |
val.DeserializeFromString ("Normal:0.1:0.2", MakeRandomVariableChecker ()); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1484 |
RandomVariable rng = val.Get (); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1485 |
NS_TEST_ASSERT_EQUAL (val.SerializeToString (MakeRandomVariableChecker ()), "Normal:0.1:0.2"); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1486 |
} |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1487 |
{ |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1488 |
RandomVariableValue val; |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1489 |
val.DeserializeFromString ("Normal:0.1:0.2:0.15", MakeRandomVariableChecker ()); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1490 |
RandomVariable rng = val.Get (); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1491 |
NS_TEST_ASSERT_EQUAL (val.SerializeToString (MakeRandomVariableChecker ()), "Normal:0.1:0.2:0.15"); |
c2d4a91bcc8a
Bug 485: implement deserialization of normal/gaussian random variables
Gustavo J. A. M. Carneiro <gjc@inescporto.pt>
parents:
4243
diff
changeset
|
1492 |
} |
4243
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1493 |
} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1494 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1495 |
return result; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1496 |
} |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1497 |
}; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1498 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1499 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1500 |
static RandomVariableTest g_random_variable_tests; |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1501 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1502 |
}//namespace ns3 |
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1503 |
|
d84e644c9692
add regression tests back
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
4242
diff
changeset
|
1504 |
#endif /* RUN_SELF_TESTS */ |