Rm-tools

-minor fix to uncertainties generated by RMpeakfit3D

-changed internal handling of Stokes I models; fixed problem when using user-defined model.

-updates to do_fitIcube tool: removed deprecated cutoff keywords, fixed tqdm
progress bars, some minor speed and efficiency improvements.
-relaxed numpy version restriction, since the numpy-astropy version conflict
has been resolved.
-added tqdm as install dependency.
-removed internal copies of external packages (corner, emcee, nestle) used by QU-fitting
-added pre-commit auto-formatting system
-escaped all LaTeX characters
-implemented auto-formatting, using black, and updated all existing code to that standard.
-removed RMS noise estimate for FDF -- it was inacurrate and misleading
-fixed corMAD noise estimate for FDF -- it was estimating MAD of PI rather than Q and U,
and thus was biased low. Now uses Q and U of FDF, and reproduces correct value in
signal-free and RM-cleaned simulations. Still biased high in dirty spectra; no fixing that.
-old versions of core rmsynth and RMSF calculation functions are now officially deprecated
-fixed crashing when trying to generate RMsynth1d plots for source with negative Stokes I model
-improved FITS headers of 3D tools output
-1D: removed RMS noise estimate (for being misleading); removed _chan outputs (for being low value and confusing)
-fixed estimation of empirical FDF noise. (1D RMclean)
-QU-fitting: fixed problems with de-rotated angle wrapping
-Integrated nufft tool for the Fourier transform math (makes rmsynth faster)
-added options to rmsynth3d to decrease RAM use
-new QU-fitting models: m5 (single Burn slab), m6 (double Burn slab), m7 (single mixed emission slab)

-progress bars updated to use tdqm, which should make them faster.
-Integrated super-resolution method of Cotton & Rudnick 2023, which creates
narrower RMSF shapes by setting lambda^2_0 = 0.
-Added options to dofitIcube tool to save more outputs (coefficient maps,
errors, etc).
-Made RMSF theoretical width calculations consistent: Dickey et al. 2019
now used throughout.

-merged in bwdepol branch. This does not change any core functionality in the
existing RM-Tools programs, but adds two new command-line tools: rmtools_bwpredict
and rmtools_bwdepol, which predict and correct for bandwidth depolarization
respectively.
-fixed a bug with JSON outputs not handling numpy variables. This was partially
due to a problem with asscalar() in the utils using a depcreciated numpy function,
but partially from other causes.
-an improved version of the Stokes I cube fitting routine has been developed
by Lerato Sebokolody, which includes parallelization.

-changed Stokes I fitting to always use 64 bit values (even if data is 32 bit).
This is because the fitter was found to be vulnerable to numerical issues
when fitting log-polynomials in 32-bits (thanks to Shannon Vanderwoude for
discovering the problem).
-added parameter errors of Stokes I fit as output to Stokes I fitDict, to
RMsynth1D outputs, and QU-fitting outputs. These are the 1-sigma errors
computed by the mpfit routine.
-the renormalize_StokesI_model function is now depreciated, since it cannot
transform the errors. It is left in for the time being in case it is needed
for any legacy reasons.
-the make_freq_file.py/rmtools_freqfile tool was found to sometime produce
truncated output (not all channels included). This was fixed and changed
to a more rigourous calculation (using WCS, as suggested by Alec Thomson).
-modified do_fitIcube.py/rmtools_fitIcube to output maps of the fit coefficients
and their errors. Thanks to Biny Sebastian for providing this code.
-merged pull request from Alec Thomson that changes the fractional polarization
reported by rmsynth1d to use the bias-corrected polarized intensity.
!!!Previous versions used the non-corrected polarized intensity, from this
version onwards the corrected polarized intensity will be used!!!
-merged pull request from Alec Thomson that implements windowed cleaning
and adds a command line flag to set the windowed clean threshold.
-merged pull request from Alec Thomson that changed the QU-fitting routine
to use the bilby package for the nested sampling. This package allows the
user to specify which sampler they wish to use -- the default is dynesty,
but users can choose to use pymultinest or other supported samplers.
-removed support for Python 3.6, as a consequence of the change to bilby.
From this point forward, Python 3.7 is the minimum version supported
(at least for the QU-fitting portion of the package).