utils.test_linalg fails when use_symmetry_reduced_integrals is set True

[woodsp-ibm]

In looking to remedy the amount of deprecation messages in the unit test logs I had also changed the settings over to the defaults I would expect us to move to in the next version. So I had updated the settings file, e.g.

        self._use_symmetry_reduced_integrals: bool = True

This however causes utils.test_linalg to fail (works fine when set False), as follows:

  File "test/utils/test_linalg.py", line 104, in test_double_factorized_error_threshold_max_vecs
    diag_coulomb_mats, orbital_rotations = double_factorized(
  File "qiskit_nature/utils/linalg.py", line 312, in double_factorized
    reshaped_tensor = np.reshape(two_body_tensor, (n_modes**2, n_modes**2))
  File "<__array_function__ internals>", line 5, in reshape
  File "qiskit_nature/second_q/operators/symmetric_two_body.py", line 159, in __array_function__
    return super().__array_function__(func, types, args, kwargs)
  File "qiskit_nature/second_q/operators/tensor.py", line 179, in __array_function__
    ret = func(*new_args, **kwargs)
  File "<__array_function__ internals>", line 5, in reshape
  File "python3.9/site-packages/numpy/core/fromnumeric.py", line 298, in reshape
    return _wrapfunc(a, 'reshape', newshape, order=order)
  File "python3.9/site-packages/numpy/core/fromnumeric.py", line 57, in _wrapfunc
    return bound(*args, **kwds)
ValueError: cannot reshape array of size 231 into shape (36,36)