Pyoptimalestimation

Minor fix to ensure pyOptimalEstimation works with numpy >= 1.24

Contribution by deweatherman for computing the Jacobian in a more efficient way. The code applies automatically the selected scheme with the following priority:

`userJacobian` (if userJacobian function is provided) -> multipleProfiles Jacobian (if `multipleForwardKwArgs` is provided) -> normal way to estimate the Jacobian iteratively (if none of the new inputs is provided)

The two new capabilities are:

`multipleForwardKwArgs`: Use the capability of most forward models to simulate several input profiles at the same time: so instead of looping through the perturbed profiles, the forward model is called once for all the perturbed profiles; this simply uses any under the hood optimisations of the forward solver. Requires arguments to set up the forward model for multiple profiles (we call them multipleForwardKwArgs)

`userJacobian`: If the user has its own function to compute Jacobians in a better way, they can provide it in the same way as for the forward model. This exploits the fact that most forward models include efficient ways to compute Jacobians using the same internal data structures as the forward model. Requires the function definition (e.g. similar to the forward function definition, we call it `userJacobian`)

* Added new option `convergenceTest` to determine whether convergence test should be done in x or y space. Options `x`, `y`, or `auto` (use space of smaller dimension). Default `x`, because experience shows that the retrieval converges faster without a notable change in the quality of the solution.
* Added tests ensuring that covariance matrices are symmetric
* Bugfix in the convergence test in x-space

Lot's of minor bugfixes

No major changes.