Rcpp Version 1.0.9
rt.h
Go to the documentation of this file.
1
// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; indent-tabs-mode: nil; -*-
2
//
3
// rt.h: Rcpp R/C++ interface class library --
4
//
5
// Copyright (C) 2010 - 2016 Douglas Bates, Dirk Eddelbuettel and Romain Francois
6
//
7
// This file is part of Rcpp.
8
//
9
// Rcpp is free software: you can redistribute it and/or modify it
10
// under the terms of the GNU General Public License as published by
11
// the Free Software Foundation, either version 2 of the License, or
12
// (at your option) any later version.
13
//
14
// Rcpp is distributed in the hope that it will be useful, but
15
// WITHOUT ANY WARRANTY; without even the implied warranty of
16
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17
// GNU General Public License for more details.
18
//
19
// You should have received a copy of the GNU General Public License
20
// along with Rcpp. If not, see <http://www.gnu.org/licenses/>.
21
22
#ifndef Rcpp__stats__random_rt_h
23
#define Rcpp__stats__random_rt_h
24
25
namespace
Rcpp
{
26
namespace
stats {
27
28
class
TGenerator
:
public
::Rcpp::Generator
<double> {
29
public
:
30
31
TGenerator
(
double
df_ ) :
df
(df_),
df_2
(df_/2.0) {}
32
33
inline
double
operator()
()
const
{
34
/* Some compilers (including MW6) evaluated this from right to left
35
return norm_rand() / sqrt(rchisq(df) / df); */
36
double
num =
norm_rand
();
37
38
// return num / sqrt(rchisq(df) / df);
39
// replaced by the followoing line to skip the test in
40
// rchisq because we already know
41
return
num /
::sqrt
( ::Rf_rgamma(
df_2
, 2.0) /
df
);
42
}
43
44
private
:
45
double
df
,
df_2
;
46
};
47
48
}
// stats
49
}
// Rcpp
50
51
#endif
Rcpp::Generator
Definition:
random.h:28
Rcpp::stats::TGenerator
Definition:
rt.h:28
Rcpp::stats::TGenerator::df
double df
Definition:
rt.h:45
Rcpp::stats::TGenerator::df_2
double df_2
Definition:
rt.h:45
Rcpp::stats::TGenerator::TGenerator
TGenerator(double df_)
Definition:
rt.h:31
Rcpp::stats::TGenerator::operator()
double operator()() const
Definition:
rt.h:33
R::norm_rand
double norm_rand(void)
Definition:
Rmath.h:30
Rcpp::algorithm::sqrt
void sqrt(InputIterator begin, InputIterator end, OutputIterator out)
Definition:
algorithm.h:479
Rcpp
Rcpp API.
Definition:
algo.h:28
inst
include
Rcpp
stats
random
rt.h
Generated on Sat Jul 9 2022 09:14:52 for Rcpp Version 1.0.9 by
1.9.1