SCM

SCM Repository

[matrix] Diff of /pkg/tests/matr-exp.R
ViewVC logotype

Diff of /pkg/tests/matr-exp.R

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

revision 519, Sat Feb 5 14:25:56 2005 UTC revision 520, Sat Feb 5 14:41:29 2005 UTC
# Line 13  Line 13 
13  assert.EQ.mat(e2, rbind(c(-0.73575875814475, 0.55181909965810),  assert.EQ.mat(e2, rbind(c(-0.73575875814475, 0.55181909965810),
14                          c(-1.47151759908826, 1.10363824071557)),                          c(-1.47151759908826, 1.10363824071557)),
15                tol = 1e-14)                tol = 1e-14)
16    ## The ``surprising identity''      det(exp(A)) == exp( tr(A) )
17    ## or                           log det(exp(A)) == tr(A) :
18    stopifnot(all.equal(determinant(e2)$modulus, sum(diag(m2))))
19    
20  m3 <- Matrix(cbind(0,rbind(6*diag(3),0)), nc = 4)  m3 <- Matrix(cbind(0,rbind(6*diag(3),0)), nc = 4)
21  e3 <- expm(m3)  e3 <- expm(m3)

Legend:
Removed from v.519  
changed lines
  Added in v.520

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