tests/internals/linalg.r

In wrathematics/fmlr: Fused Matrix Library for R

# add
test = linalg_add(x=x, y=y)$to_robj()
truth = xr + yr
cmp(test, truth)

test = linalg_add(transx=TRUE, transy=TRUE, x=x, y=y)$to_robj()
truth = t(xr) + t(yr)
cmp(test, truth)



# matmult
test = linalg_matmult(transx=FALSE, transy=TRUE, x=x, y=y)$to_robj()
truth = xr %*% t(yr)
cmp(test, truth)

test = linalg_matmult(transx=TRUE, transy=FALSE, x=x, y=y)$to_robj()
truth = t(xr) %*% yr
cmp(test, truth)

test = linalg_matmult(transx=FALSE, transy=FALSE, x=x, y=v)$to_robj()
truth = xr %*% vr
cmp(test, truth)



# crossprod
test = linalg_crossprod(x=x)$to_robj()
truth = crossprod(xr)
cmp(test[lower.tri(test, diag=TRUE)], truth[lower.tri(truth, diag=TRUE)])

test = linalg_tcrossprod(x=x)$to_robj()
truth = tcrossprod(xr)
cmp(test[lower.tri(test, diag=TRUE)], truth[lower.tri(truth, diag=TRUE)])



# xpose
test = linalg_xpose(x=x)$to_robj()
truth = t(xr)
cmp(test, truth)