The R Project SVN R-packages

Rev

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

Rev 1577 Rev 1631
Line 237... Line 237...
237
 
237
 
238
void make_array_triangular(double *x, SEXP from);
238
void make_array_triangular(double *x, SEXP from);
239
 
239
 
240
SEXP Matrix_expand_pointers(SEXP pP);
240
SEXP Matrix_expand_pointers(SEXP pP);
241
 
241
 
-
 
242
 
-
 
243
/** 
-
 
244
 * Elementwise increment dest by src
-
 
245
 * 
-
 
246
 * @param dest vector to be incremented
-
 
247
 * @param src vector to be added to dest
-
 
248
 * @param n length of vectors
-
 
249
 * 
-
 
250
 * @return dest
-
 
251
 */
-
 
252
static R_INLINE double*
-
 
253
vecIncrement(double dest[], const double src[], int n) {
-
 
254
    int i;
-
 
255
    for (i = 0; i < n; i++) dest[i] += src[i];
-
 
256
    return dest;
-
 
257
}
-
 
258
 
-
 
259
/** 
-
 
260
 * Elementwise sum of src1 and src2 into dest
-
 
261
 * 
-
 
262
 * @param dest vector to be incremented
-
 
263
 * @param src1 vector to be added
-
 
264
 * @param src1 second vector to be added
-
 
265
 * @param n length of vectors
-
 
266
 * 
-
 
267
 * @return dest
-
 
268
 */
-
 
269
static R_INLINE double*
-
 
270
vecSum(double dest[], const double src1[], const double src2[],
-
 
271
       int n) {
-
 
272
    int i;
-
 
273
    for (i = 0; i < n; i++) dest[i] = src1[i] + src2[i];
-
 
274
    return dest;
-
 
275
}
-
 
276
 
242
#ifdef __cplusplus
277
#ifdef __cplusplus
243
}
278
}
244
#endif
279
#endif
245
 
280
 
246
#endif /* MATRIX_MUTILS_H_ */
281
#endif /* MATRIX_MUTILS_H_ */