Attempting to split the code into several parts and avoid some repeating code blocks. Does not work yet.

Reogranisation into several submodules with separate functions. Some improvement in code readability by removing repeated code.

Better display of not converged points.
Add a few plot "other_kwargs" options.

Change default setting for optimization.
Add loop to check if it's really stable after optimization.

Typo in Eq. 19

Default behaviour is use_GaussLMCI=False
Add MC impedance to _PTBL_alves

Merge modif from mbtrack_v0.7

See merge request PA/collective-effects/robinsonmodes!1

Switch optimization to generic methods.
Save data plot by default.
Add many plot options using other_kwargs.
Fix coutours when there is not converged data.
Add kwargs for opti and add method to save name.

Fix scan_RoQ_Q0 and scan_psi_HC_Q0 as QL was a fixed value (not modfied when Q0 is modified).
Add save and load functions for data recording.
Add max_azi and max_rad as kwargs for __PTBL_alves

Remove set MC.theta phase in Avles for now (introduced bug ...)

Different ways to compute PTBL:
+ "Venturini" uses T.He threshold
+ "Alves" uses longitudinal mode coupling from pycolleff

Merge code for active and passive HC.
Add documentation.