Add a dict_dtype optional parameter to Monitor.monitor_init to allow to choose the dtype to use to save the data.
Improve CavityMonitor to allow for saving beam and cavity phasor for all bunches
Add a plot_cavitydata function
Fix a bug on CavityResonator.track -> nturn += 1 at correct time
Change CavityResonator.track to allow for monitoring the cavity/beam phasor at each bunch

Change back invariant name to cs_invariant
Simplify SkewQuadrupole track method

Add comments.
Remove usage of bins, use only center.

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

Remove for-loops inside the function.

- cs_invariant method in Bunch is renamed to 'action'.
- Action is now saved together with BunchMonitor and can be plotted in plot_bunchdata

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.