Displaying 2 results from an estimated 2 matches for "trprod2".
Did you mean:
trprod
2006 Mar 09
0
When calling external C-function repeatedly I get different results; can't figure out why..
...int nrx, int ncx,
double *y, int nry, int ncy, double *z)
{
char *transa = "N", *transb = "N";
double one = 1.0, zero = 0.0;
F77_CALL(dgemm)(transa, transb, &nrx, &ncy, &ncx, &one,
x, &nrx, y, &nry, &zero, z, &nrx);
}
SEXP trProd2(SEXP x, SEXP y)
{
int nrx, ncx, nry, ncy, mode, i;
SEXP xdims, ydims, ans, ans2, tr;
xdims = getAttrib(x, R_DimSymbol);
ydims = getAttrib(y, R_DimSymbol);
mode = REALSXP;
nrx = INTEGER(xdims)[0];
ncx = INTEGER(xdims)[1];
nry = INTEGER(ydims)[0];
ncy = INTEGER(ydims)[1];
PROTECT(...
2006 Mar 09
0
When calling external C-function repeatedly I get differentresults; can't figure out why..
...*y, int nry, int ncy, double *z) {
> char *transa = "N", *transb = "N";
> double one = 1.0, zero = 0.0;
> F77_CALL(dgemm)(transa, transb, &nrx, &ncy, &ncx, &one,
> x, &nrx, y, &nry, &zero, z, &nrx); }
>
> SEXP trProd2(SEXP x, SEXP y)
> {
> int nrx, ncx, nry, ncy, mode, i;
> SEXP xdims, ydims, ans, ans2, tr;
> xdims = getAttrib(x, R_DimSymbol);
> ydims = getAttrib(y, R_DimSymbol);
> mode = REALSXP;
> nrx = INTEGER(xdims)[0];
> ncx = INTEGER(xdims)[1];
> nry = INTEGER(ydi...