docs: public DoubleMatrix invert()
      {
          if (A == null)
          {
              throw new IllegalArgumentException("decompose() not called");
          }
          DoubleMatrix IA = getInstance(rows(), columns());
          lupInvert(A, P, IA);
          return IA;
      }