inst/tinytest/test_pinv.R

In lintools: Manipulation of Linear Systems of (in)Equalities

## Pseudoinverse

## Computing pseudoinverse
 # example from appendix of Schaum's outline (2009) on linear algebra
  A <- matrix(c(
     1,  1, -1,  2,
     2,  2, -1,  3,
    -1, -1,  2, -3
  ),byrow=TRUE,nrow=3)
  Aplus55 <- matrix(c(
     1, 18,  15,
     1, 18,  15,
    -2, 19,  25,
     3, -1, -10
    ),byrow=TRUE, nrow=4)

  expect_equal(pinv(A),Aplus55/55)

## bugfixes
  # this crashed because there's only one s.v. and 'diag' 
  # reacts differently when presented a single number.
  pinv(matrix(c(1,-1,0,0,0)))