Compute the covariance matrix cov = inv (J^T J) by QRP^T decomposition of J
Returns (Value, idx).
(Value, idx)