The R Project SVN R

Rev

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

Rev 22959 Rev 24049
Line 40... Line 40...
40 
    double ai;
 40 
    double ai;
41 
    int i;
 41 
    int i;
42 
    errno = 0;
 42 
    errno = 0;
43 
    for (i = 0; i < n; i++) {
 43 
    for (i = 0; i < n; i++) {
44 
	ai = a[i % na];
 44 
	ai = a[i % na];
45 
	if (R_FINITE(ai)) {
 - 
 
46 
	    x[i] = MATH_CHECK(f(ai));
 45 
	x[i] = MATH_CHECK(f(ai));
47 
	    if (!R_FINITE(x[i])) naflag = 1;
 46 
	if (!R_FINITE(x[i])) naflag = 1;
48 
	}
 - 
 
49 
	else x[i] = NA_REAL;
 - 
 
50 
    }
 47 
    }
51 
    return(naflag);
 48 
    return(naflag);
52 
}
 49 
}
53 
 
 50 
 
54 
#define RAND1(num,name) \
 51 
#define RAND1(num,name) \
Line 115... Line 112...
115 
    Rboolean naflag = FALSE;
 112 
    Rboolean naflag = FALSE;
116 
    errno = 0;
 113 
    errno = 0;
117 
    for (i = 0; i < n; i++) {
 114 
    for (i = 0; i < n; i++) {
118 
	ai = a[i % na];
 115 
	ai = a[i % na];
119 
	bi = b[i % nb];
 116 
	bi = b[i % nb];
120 
	if (R_FINITE(ai) && R_FINITE(bi)) {
 - 
 
121 
	    x[i] = MATH_CHECK(f(ai, bi));
 117 
	x[i] = MATH_CHECK(f(ai, bi));
122 
	    if (!R_FINITE(x[i])) naflag = 1;
 118 
	if (!R_FINITE(x[i])) naflag = 1;
123 
	}
 - 
 
124 
	else x[i] = NA_REAL;
 - 
 
125 
    }
 119 
    }
126 
    return(naflag);
 120 
    return(naflag);
127 
}
 121 
}
128 
 
 122 
 
129 
#define RAND2(num,name) \
 123 
#define RAND2(num,name) \
Line 202... Line 196...
202 
    errno = 0;
 196 
    errno = 0;
203 
    for (i = 0; i < n; i++) {
 197 
    for (i = 0; i < n; i++) {
204 
	ai = a[i % na];
 198 
	ai = a[i % na];
205 
	bi = b[i % nb];
 199 
	bi = b[i % nb];
206 
	ci = c[i % nc];
 200 
	ci = c[i % nc];
207 
	if (R_FINITE(ai) && R_FINITE(bi) && R_FINITE(ci)) {
 - 
 
208 
	    x[i] = MATH_CHECK(f(ai, bi, ci));
 201 
	x[i] = MATH_CHECK(f(ai, bi, ci));
209 
	    if (!R_FINITE(x[i])) naflag = TRUE;
 202 
	if (!R_FINITE(x[i])) naflag = TRUE;
210 
	}
 - 
 
211 
	else x[i] = NA_REAL;
 - 
 
212 
    }
 203 
    }
213 
    return(naflag);
 204 
    return(naflag);
214 
}
 205 
}
215 
 
 206 
 
216 
#define RAND3(num,name) \
 207 
#define RAND3(num,name) \