Not needed any more due the SynchrotronRadiation change.

Fix bug on Beam.distance_between_bunches for the last bunch when the first bunches are empty.
Change Beam.distance_between_bunches to compute the array only at initalization with Beam.update_distance_between_bunches
Beam.mpi_share_distributions can now share all bunch dimensions.

Correct time decay !
Now ok versus for loop (decay then beam loading)

Typo in the type of data to share for bunch profile. Not sure if there was an effect.

In track method of cavity resonator, there was an additional phasor decay of bin length.
To avoid that the good boundaries are set bins[0] instead of center[0].
For now the bins are also shared in Beam.mpi_share_distributions

Add single core tracking for CavityResonator
Correct missing ref_frame is CavityResonator.track
Change phase of generator voltage to be the bin phase
Correct generator phase for non integer m
Change init_phasor track time to 10 times
Change phasor_evolution k to change order of steps

Add phasor initialization (on the "rf" reference)
Change beam phasor evolution to allow for "beam" or "rf" reference

Correct case when where is no mp in the bin
Remove useless test

Add tracking to CavityResonator based on the mbtrack implementation but using a analytic formula for speed-up
Add method to Beam to share longitudinal distributions
Add methods to Mpi to look for next and previous bunches
Add start of a CavityMonitor, to improve to deal with multi-cores

Main change
- Variable local_beta, local_alpha, and local_gamma are made property objects.

Minor change
- Add an option to plot phase space in different style.

Remove part which was causing "damping".
We have to re-process the wake function at each track call and not to reuse the previously processed wake function

WakePotential -> add method to plot wake potential for a perfect gaussian bunch
utilities -> new function gaussian_bunch which gives the gaussian bunch profile
utilities -> read_IW2D now also reads wake function from IW2D
utilities -> new function read_IW2D_folder which reads the full results folder of IW2D into a WakeField element.
resistive_wall -> change the sign of the wake functions to agree with the code convention

Bugfix due to switch to numpy arrays.
Allow to save less mp than the full number in the bunch.

Bugfix2 : PhaseSpaceMonitor

Fix dipole_moment due to switch to numpy arrays for Bunch.binning

Update docstrings
Use pd.DataFrame representation for bunch __repr__
Add the track_alive option for Bunch and Beam.init_beam

Change track method to allow for Beam objects
If pynaff can not find the tune, the tune value is nan.
Nanmean and nanstd are used to avoid to bias the tune value.
! Untested if some particles are lost !

WakePotentialMonitor now use interpolation on a fixed grid to be able to save the wake potential data (which size change from turn to turn).

1st test, not complete

Now based on two time arrays!
The bunch profile is computed from its own time array and interpolated into the wake function time array.
The input wake functions must be uniformly sampled!
Reference computation now gives loss factor and kick factor.
Tested with Wlong and Wxdip/Wydip -> OK
Untested for Wxquad/Wyquad.

Add CavityResonator class
Add BeamLoadingVlasov class, to do tested and developped further.
From BeamLoadingVlasov, only the beam_equilibrium method is well tested and benchmarked.
The cannonical_transform and solveB methods have to be tested and benchmarked further.

Compared to mbtrack, the head/tail convention is reversed which leads to this sign change.
In this python code, the head is for negative times, t<0.
The sign change allows successful tracking benchmark with mbtrack for pure inductive and resonator impedance

Modify the docstrings and add some precisions.

A method to compare energy loss from a wake potential to a loss factor.
Minor changes :
- sigma method in Ring object has been edited.
- add WakePotential in __init__ in tracking module.

Correct scale in plot_wakedata
Avoid circular import by moving import in function

Remove mean of tau array during wake function interpolation to avoid wrong interpolation if tau is not centered on zero

Add transverse wake potential to WakePotential and move class to dedicated file
Merge Resonator and TransverseResonator into Resonator with dedicated file
wakefunction_from_imp -> change trimming to allow for a float value
wakefunction_from_imp -> transverse impedance divided by 1j, to be checked with impedance convention but ok with benchmark
Wakefield -> add wake_components

Fix mbi_threshold to have sigma in [s}
Add plot_critical_mass to instabilities (ions)
Add sigma(position) to synchrotron to compute the beam size at any ring location

- Deconvolution methods in WakeFunction and Impedance class were improved.
- Analytical impedance calculation method was added to Resonator class.
- Wake potential calculation method in WakePotential2 class was corrected.

- New Wakefunction class to store a wake function object.
- New deconvolve_wp and deconvolve_imp methods were added in Wakefunction class.
- New  get-impedance method was added in Impedance class.
- New Resonator class to calculate the analytical resonator wake function.
- New Wakepotential2 class to take in wake function objects, then compute the wake potential and do the tracking.

Add ABI threshold calculation
Bugfix synchrotron_tune computation
Add zoom option to ImpedanceModel.plot_area