docs: def _multiply(X, coef):
      """ Multiple X by coef element-wise, preserving sparsity. """
      if sp.issparse(X):
          return X.multiply(sp.csr_matrix(coef))
      else:
          return np.multiply(X, coef)