docs: public void solve(DoubleMatrix b, DoubleMatrix x)
      {
          if (A == null)
          {
              throw new IllegalArgumentException("decompose() not called");
          }
          lupSolve(A, P, b, x);
      }