Displaying 1 result from an estimated 1 matches for "tr_inv".
2006 Oct 13
3
Rmpi performance
...ned successfully. 0 failed.
and we want to compute the function op_mat (see below) first with the
master and then with the cluster using system.time for checking the
computational performance.
op_mat = function(mat) {
+ inv = solve(mat)
+ det_inv = det(inversa)
+ tr_inv = sum(diag(inversa))
+ return(list(c(det=det_inv,tr=tr_inv)))
+ }
>nn = 3000
>XX = matrix(rnorm(nn*nn),nn,nn)
# with the master
> system.time(op_matrici(XX))
[1] 42.283 1.883 44.168 0.000 0.000
# with the cluster
> system.time(clusterCall(cl,op_matrici,XX))
[1] 11.523 12....