The R Project SVN R

Rev

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

Rev 77685 Rev 90186
Line 30... Line 30...
30
{
30
{
31
#ifdef IEEE_754
31
#ifdef IEEE_754
32
    if (ISNAN(p) || ISNAN(shape) || ISNAN(scale))
32
    if (ISNAN(p) || ISNAN(shape) || ISNAN(scale))
33
	return p + shape + scale;
33
	return p + shape + scale;
34
#endif
34
#endif
35
    if (shape <= 0 || scale <= 0) ML_WARN_return_NAN;
35
    if (shape < 0 || scale <= 0) ML_WARN_return_NAN;
36
 
36
 
37
    R_Q_P01_boundaries(p, 0, ML_POSINF);
37
    R_Q_P01_boundaries(p, 0, ML_POSINF);
38
 
38
 
39
    return scale * pow(- R_DT_Clog(p), 1./shape) ;
39
    return scale * pow(- R_DT_Clog(p), 1./shape) ;
40
}
40
}