diff options
Diffstat (limited to 'src/test_rtnorm.cpp')
| -rw-r--r-- | src/test_rtnorm.cpp | 45 |
1 files changed, 0 insertions, 45 deletions
diff --git a/src/test_rtnorm.cpp b/src/test_rtnorm.cpp deleted file mode 100644 index 48ce933..0000000 --- a/src/test_rtnorm.cpp +++ /dev/null @@ -1,45 +0,0 @@ -// Example for using rtnorm -// -// Copyright (C) 2012 Guillaume Dollé, Vincent Mazet (LSIIT, CNRS/Université de Strasbourg) -// Licence: GNU General Public License Version 2 -// see http://www.gnu.org/licenses/old-licenses/gpl-2.0.txt -// -// Depends: LibGSL -// OS: Unix based system - - -#include <iostream> -#include <gsl/gsl_rng.h> - -#include "rtnorm.hpp" - - -int main() -{ - double a = 1; // Left bound - double b = 9; // Right bound - double mu = 2; // Mean - double sigma = 3; // Standard deviation - std::pair<double, double> s; // Output argument of rtnorm - int K = 1e5; // Number of random variables to generate - - //--- GSL random init --- - gsl_rng_env_setup(); // Read variable environnement - const gsl_rng_type* type = gsl_rng_default; // Default algorithm 'twister' - gsl_rng *gen = gsl_rng_alloc (type); // Rand generator allocation - - //--- generate and display the random numbers --- - //std::cout<<"# x p(x)"<<std::endl; - std::cout<<a<<" "<<b<<std::endl; - std::cout<<mu<<" "<<sigma<<std::endl; - for(int k=0; k<K; k++) - { - s = rtnorm(gen,a,b,mu,sigma); - std::cout<<s.first<<" "<<s.second<<std::endl; - } - - gsl_rng_free(gen); // GSL rand generator deallocation - - return 0; -} - |
