The R Project SVN R

Rev

Rev 8431 | Rev 42302 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 8431 Rev 24049
Line 39... Line 39...
39
    /* Uses these GLOBALS to save time when many rv's are generated : */
39
    /* Uses these GLOBALS to save time when many rv's are generated : */
40
    static double beta, gamma, delta, k1, k2;
40
    static double beta, gamma, delta, k1, k2;
41
    static double olda = -1.0;
41
    static double olda = -1.0;
42
    static double oldb = -1.0;
42
    static double oldb = -1.0;
43
 
43
 
44
    if (aa <= 0. || bb <= 0.)
44
    if (aa <= 0. || bb <= 0. || (!R_FINITE(aa) && !R_FINITE(bb)))
45
	ML_ERR_return_NAN;
45
	ML_ERR_return_NAN;
46
 
46
 
-
 
47
    if (!R_FINITE(aa))
-
 
48
    	return 1.0;
-
 
49
 
-
 
50
    if (!R_FINITE(bb))
-
 
51
    	return 0.0;
-
 
52
 
47
    /* Test if we need new "initializing" */
53
    /* Test if we need new "initializing" */
48
    qsame = (olda == aa) && (oldb == bb);
54
    qsame = (olda == aa) && (oldb == bb);
49
    if (!qsame) { olda = aa; oldb = bb; }
55
    if (!qsame) { olda = aa; oldb = bb; }
50
 
56
 
51
    a = fmin2(aa, bb);
57
    a = fmin2(aa, bb);