LCOV - code coverage report
Current view: top level - defines - random.hpp (source / functions) Hit Total Coverage
Test: clean.info Lines: 18 18 100.0 %
Date: 2013-05-24 Functions: 5 5 100.0 %
Branches: 5 6 83.3 %

           Branch data     Line data    Source code
       1                 :            : /* This software is distributed under BSD 3-clause license (see LICENSE file).
       2                 :            :  *
       3                 :            :  * Copyright (c) 2012-2013 Sergey Lisitsyn
       4                 :            :  */
       5                 :            : 
       6                 :            : #ifndef TAPKEE_DEFINES_RANDOM_H_
       7                 :            : #define TAPKEE_DEFINES_RANDOM_H_
       8                 :            : 
       9                 :            : #include <cstdlib>
      10                 :            : #include <algorithm>
      11                 :            : #include <limits>
      12                 :            : 
      13                 :            : namespace tapkee
      14                 :            : {
      15                 :            : 
      16                 :       4998 : inline IndexType uniform_random_index() 
      17                 :            : {
      18                 :            : #ifdef CUSTOM_UNIFORM_RANDOM_INDEX_FUNCTION
      19                 :            :         return CUSTOM_UNIFORM_RANDOM_INDEX_FUNCTION % std::numeric_limits<IndexType>::max();
      20                 :            : #else
      21                 :       4998 :         return std::rand();
      22                 :            : #endif
      23                 :            : }
      24                 :            : 
      25                 :       4998 : inline IndexType uniform_random_index_bounded(IndexType upper)
      26                 :            : {
      27                 :       4998 :         return uniform_random_index() % upper;
      28                 :            : }
      29                 :            : 
      30                 :         33 : inline ScalarType uniform_random()
      31                 :            : {
      32                 :            : #ifdef CUSTOM_UNIFORM_RANDOM_FUNCTION
      33                 :            :         return CUSTOM_UNIFORM_RANDOM_FUNCTION;
      34                 :            : #else
      35                 :         33 :         return std::rand()/((double)RAND_MAX+1);
      36                 :            : #endif
      37                 :            : }
      38                 :            : 
      39                 :        383 : inline ScalarType gaussian_random()
      40                 :            : {
      41                 :            : #ifdef CUSTOM_GAUSSIAN_RANDOM_FUNCTION
      42                 :            :         return CUSTOM_GAUSSIAN_RANDOM_FUNCTION;
      43                 :            : #else
      44                 :            :         ScalarType x, y, radius;
      45 [ +  + ][ -  + ]:        383 :         do {
                 [ +  + ]
      46                 :        383 :                 x = 2*(std::rand()/((double)RAND_MAX+1)) - 1;
      47                 :        383 :                 y = 2*(std::rand()/((double)RAND_MAX+1)) - 1;
      48                 :        383 :                 radius = (x * x) + (y * y);
      49                 :            :         } while ((radius >= 1.0) || (radius == 0.0));
      50                 :        308 :         radius = std::sqrt(-2 * std::log(radius) / radius);
      51                 :        308 :         x *= radius;
      52                 :        308 :         y *= radius;
      53                 :        308 :         return x;
      54                 :            : #endif
      55                 :            : }
      56                 :            : 
      57                 :            : template <class RAI>
      58                 :        102 : inline void random_shuffle(RAI first, RAI last)
      59                 :            : {
      60                 :        102 :         std::random_shuffle(first,last,uniform_random_index_bounded);
      61                 :        102 : }
      62                 :            : 
      63                 :            : }
      64                 :            : 
      65                 :            : #endif
      66                 :            : 
      67                 :            : 

Generated by: LCOV version 1.9