rename cond -> rcond

There is a slight mess with the naming condition number in scipy and
numpy. scipy uses word `cond` and numpy `rcond`. rcond makes actually
more sense because it is a lower bound on the reciprocal condition
number.

Author: agoscinski

src/equisolve/numpy/models/linear_model.py

@@ -106,7 +106,7 @@ def _solve(
         y: TensorBlock,
         alpha: TensorBlock,
         sample_weight: TensorBlock,
-        cond: Optional[float] = None,
+        rcond: Optional[float] = None,
     ) -> TensorBlock:
         """A regularized solver using ``np.linalg.lstsq``."""
         self._used_auto_solver = None
@@ -197,9 +197,9 @@ def _solve(
             # and b is [y*sqrt(w), 0]
             X_eff = np.vstack([sqrt_sw_arr * X_arr, np.diag(np.sqrt(alpha_arr))])
             y_eff = np.hstack([y_arr * sqrt_sw_arr.flatten(), np.zeros(num_properties)])
-            if cond is None:
-                cond = max(X_arr.shape) * np.finfo(X_arr.dtype.char.lower()).eps
-            w = scipy.linalg.lstsq(X_eff, y_eff, cond=cond, overwrite_a=True)[0].ravel()
+            if rcond is None:
+                rcond = max(X_arr.shape) * np.finfo(X_arr.dtype.char.lower()).eps
+            w = scipy.linalg.lstsq(X_eff, y_eff, cond=rcond, overwrite_a=True)[0].ravel()
         else:
             raise ValueError(
                 f"Unknown solver {self._solver!r} only 'auto', 'cholesky',"
@@ -227,7 +227,7 @@ def fit(
         alpha: Union[float, TensorMap] = 1.0,
         sample_weight: Union[float, TensorMap] = None,
         solver="auto",
-        cond: Optional[float] = None,
+        rcond: Optional[float] = None,
     ) -> None:
         """Fit a regression model to each block in `X`.
 
@@ -263,7 +263,7 @@ def fit(
               on the dual problem (X@X.T) w_dual = y,
               the primal weights are obtained by w = X.T @ w_dual
             - **"lstsq"**: using :func:`scipy.linalg.lstsq` on the linear system X w = y
-        :param cond:
+        :param rcond:
             Cut-off ratio for small singular values during the fit. For the purposes of
             rank determination, singular values are treated as zero if they are smaller
             than `cond` times the largest singular value in "weights" matrix. Only
@@ -306,7 +306,7 @@ def fit(
             alpha_block = alpha.block(key)
             sw_block = sample_weight.block(key)
 
-            weight_block = self._solve(X_block, y_block, alpha_block, sw_block, cond)
+            weight_block = self._solve(X_block, y_block, alpha_block, sw_block, rcond)
 
             weights_blocks.append(weight_block)