The R Project SVN R

Rev

Rev 68947 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 68947 Rev 77685
Line 31... Line 31...
31
#include "nmath.h"
31
#include "nmath.h"
32
 
32
 
33
double rnorm(double mu, double sigma)
33
double rnorm(double mu, double sigma)
34
{
34
{
35
    if (ISNAN(mu) || !R_FINITE(sigma) || sigma < 0.)
35
    if (ISNAN(mu) || !R_FINITE(sigma) || sigma < 0.)
36
	ML_ERR_return_NAN;
36
	ML_WARN_return_NAN;
37
    if (sigma == 0. || !R_FINITE(mu))
37
    if (sigma == 0. || !R_FINITE(mu))
38
	return mu; /* includes mu = +/- Inf with finite sigma */
38
	return mu; /* includes mu = +/- Inf with finite sigma */
39
    else
39
    else
40
	return mu + sigma * norm_rand();
40
	return mu + sigma * norm_rand();
41
}
41
}