Displaying 1 result from an estimated 1 matches for "copymatrixdata".
2008 Apr 05
2
Adding a Matrix Exponentiation Operator
...ies for this implementation of
exponentiation by squaring. Is there a way to cut down on the number of
copies I am making here (I am assuming that the lhs and rhs of matprod()
must be different instances).
Any feedback appreciated !
Thanks
Rory
<snip>
/* Convenience function */
static void copyMatrixData(SEXP a, SEXP b, int nrows, int ncols, int mode) {
for (int i=0; i < ncols; ++i)
for (int j=0; j < nrows; ++j)
REAL(b)[i * nrows + j] = REAL(a)[i * nrows + j];
}
SEXP do_matexp(SEXP call, SEXP op, SEXP args, SEXP rho)
{
int nrows, ncols;
SEXP matrix, tmp, dims, dims2;...