RobinHankin/onion

bug in matrix division

Closed this issue · 0 comments

Dividing by a matrix should be elementwise but uses matrix inverses:

> Him/matrix(1,2,2)
Error in solve.default(e2) : 
  Lapack routine dgesv: system is exactly singular: U[2,2] = 0
>