Function assertValidPinv

test/unit-tests/function/matrix/pinv.test.js:15–32

* Tests whether `A_` is Moore–Penrose inverse of `A`. * * @param {Matrix | Array} A A two dimensional matrix or array. * @param {Matrix | Array} A_ A two dimensional matrix or array. * * @return {Boolean} Returns true if `A_` is a valid Moore–Penrose inverse of `A`

(A, A_)

 * @return {Boolean} Returns true if `A_` is a valid Moore–Penrose inverse of `A`
 */
function assertValidPinv (A, A_) {
  const Asize = math.size(A).valueOf()

  const rows = Asize[0]
  const cols = Asize[1]

  // sizes match
  assert.deepStrictEqual(math.size(A_).valueOf(), [cols, rows])

  // A A_ A = A
  approxDeepEqual(math.add(math.multiply(A, math.multiply(A_, A)), math.Complex(1, 1)).valueOf(), math.add(A, math.Complex(1, 1)).valueOf())
  // A_ A A_ = A_
  approxDeepEqual(math.add(math.multiply(A_, math.multiply(A, A_)), math.Complex(1, 1)).valueOf(), math.add(A_, math.Complex(1, 1)).valueOf())
  // (A A_)* = A A_
  approxDeepEqual(math.add(math.ctranspose(math.multiply(A, A_)), math.Complex(1, 1)).valueOf(), math.add(math.multiply(A, A_), math.Complex(1, 1)).valueOf())
  // (A_ A)* = A_ A
  approxDeepEqual(math.add(math.ctranspose(math.multiply(A_, A)), math.Complex(1, 1)).valueOf(), math.add(math.multiply(A_, A), math.Complex(1, 1)).valueOf())
}

describe('pinv', function () {
  function check (A, A_, strict = false) {