r102:81563e019fa4
Thu, 20 Mar 2008 17:15:35
[Not Reviewed]
0 comments (0 inline)
0
2
0
| ... | ... |
@@ -566,24 +566,44 @@ |
| 566 | 566 |
/// |
| 567 | 567 |
/// Assign operator. The generated sequence will be identical to |
| 568 | 568 |
/// the other sequence. It can be used to save the current state |
| 569 | 569 |
/// of the generator and later use it to generate the same |
| 570 | 570 |
/// sequence. |
| 571 | 571 |
Random& operator=(const Random& other) {
|
| 572 | 572 |
if (&other != this) {
|
| 573 | 573 |
core.copyState(other.core); |
| 574 | 574 |
} |
| 575 | 575 |
return *this; |
| 576 | 576 |
} |
| 577 | 577 |
|
| 578 |
/// \brief Seeding random sequence |
|
| 579 |
/// |
|
| 580 |
/// Seeding the random sequence. The current number type will be |
|
| 581 |
/// converted to the architecture word type. |
|
| 582 |
template <typename Number> |
|
| 583 |
void seed(Number seed) {
|
|
| 584 |
_random_bits::Initializer<Number, Word>::init(core, seed); |
|
| 585 |
} |
|
| 586 |
|
|
| 587 |
/// \brief Seeding random sequence |
|
| 588 |
/// |
|
| 589 |
/// Seeding the random sequence. The given range should contain |
|
| 590 |
/// any number type and the numbers will be converted to the |
|
| 591 |
/// architecture word type. |
|
| 592 |
template <typename Iterator> |
|
| 593 |
void seed(Iterator begin, Iterator end) {
|
|
| 594 |
typedef typename std::iterator_traits<Iterator>::value_type Number; |
|
| 595 |
_random_bits::Initializer<Number, Word>::init(core, begin, end); |
|
| 596 |
} |
|
| 597 |
|
|
| 578 | 598 |
/// \brief Returns a random real number from the range [0, 1) |
| 579 | 599 |
/// |
| 580 | 600 |
/// It returns a random real number from the range [0, 1). The |
| 581 | 601 |
/// default Number type is \c double. |
| 582 | 602 |
template <typename Number> |
| 583 | 603 |
Number real() {
|
| 584 | 604 |
return _random_bits::RealConversion<Number, Word>::convert(core); |
| 585 | 605 |
} |
| 586 | 606 |
|
| 587 | 607 |
double real() {
|
| 588 | 608 |
return real<double>(); |
| 589 | 609 |
} |
| ... | ... |
@@ -15,23 +15,31 @@ |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#include <lemon/random.h> |
| 20 | 20 |
#include "test_tools.h" |
| 21 | 21 |
|
| 22 | 22 |
///\file \brief Test cases for random.h |
| 23 | 23 |
/// |
| 24 | 24 |
///\todo To be extended |
| 25 | 25 |
/// |
| 26 | 26 |
|
| 27 |
int seed_array[] = {1, 2};
|
|
| 28 |
|
|
| 27 | 29 |
int main() |
| 28 | 30 |
{
|
| 29 | 31 |
double a=lemon::rnd(); |
| 30 | 32 |
check(a<1.0&&a>0.0,"This should be in [0,1)"); |
| 31 | 33 |
a=lemon::rnd.gauss(); |
| 32 | 34 |
a=lemon::rnd.gamma(3.45,0); |
| 33 | 35 |
a=lemon::rnd.gamma(4); |
| 34 | 36 |
//Does gamma work with integer k? |
| 35 | 37 |
a=lemon::rnd.gamma(4.0,0); |
| 36 | 38 |
a=lemon::rnd.poisson(.5); |
| 39 |
|
|
| 40 |
lemon::rnd.seed(100); |
|
| 41 |
lemon::rnd.seed(seed_array, seed_array + |
|
| 42 |
(sizeof(seed_array) / sizeof(seed_array[0]))); |
|
| 43 |
|
|
| 44 |
return 0; |
|
| 37 | 45 |
} |
0 comments (0 inline)