LUP decomposition seems to violate API

[vincechoo750]

The API for LUP decomposition seems to indicate that

(defn LU-contract [A]
 (let [{:keys [L U P]} (lu A)]
  (= (mmul P A) (mmul L U))))

should return true.

When the implementation is vectorz, such as the case of

    [[ 1 -2  3]]
A =  [ 2 -5  12]
     [ 0  2 -10]]

it is not the case that (LU-contract A) returns true.

Instead, it appears to be the case that (= (mmul A) (mmul P L U)).

Not sure if this is a typo in the API docs or in the code itself - the tests in the tests folder seem to indicate that (= (mmul A) (mmul P L U)) is the intended behavior.