SCM

SCM Repository

[matrix] Diff of /pkg/tests/dpoMatrix.R
ViewVC logotype

Diff of /pkg/tests/dpoMatrix.R

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 508, Fri Feb 4 18:22:28 2005 UTC revision 509, Fri Feb 4 23:34:01 2005 UTC
# Line 11  Line 11 
11  str(f9 <- as(chol(h9), "dtrMatrix"))  str(f9 <- as(chol(h9), "dtrMatrix"))
12  ## h9 now has factorization  ## h9 now has factorization
13  stopifnot(names(h9@factors) == "Cholesky")  stopifnot(names(h9@factors) == "Cholesky")
14  rcond(h9)  stopifnot(all.equal(rcond(h9), 9.0938e-13))
15  rcond(f9)  stopifnot(all.equal(rcond(f9), 9.1272e-7, tol = 1e-6))# more precision fails
16  str(h9)# has 'rcond' and 'factors'  str(h9)# has 'rcond' and 'factors'
17  options(digits=4)  options(digits=4)
18  crossprod(f9)# looks the same as  (cf9 <- crossprod(f9))# looks the same as  h9 :
19  h9 # but not internally!  stopifnot(all.equal(as(h9, "matrix"),
20  ##   i.e. this is all wrong : all.equal(h9, crossprod(f9))                      as(cf9,"matrix"), tol= 1e-15))
21    

Legend:
Removed from v.508  
changed lines
  Added in v.509

root@r-forge.r-project.org
ViewVC Help
Powered by ViewVC 1.0.0  
Thanks to:
Vienna University of Economics and Business Powered By FusionForge