0.6

Bug correction

• solve issue with optimized (scipy bisect) radial solver (see PR #8)

Changes

• switch radial solvers: radial corresponds now to the corrected optimized radial solver using scipy for bisect search, radial_legacy is the old one
• Store radial and azimuthal functions of the modes in the radial solver in modes0.data[<ind_mode>]['radial_func'] and modes0.data[<ind_mode>]['azimuthal_func'], can be used to apply to your mesh, e.g.:

modes = solver.solve(mode='radial_test', ...)
X, Y = np.meshgrid(...)
TH = np.arctan2(Y, X)
R = np.sqrt(X**2 + Y**2)
ind_mode = 0
psi_r = modes.data[ind_mode]['radial_func'](R)
psi_theta = modes.data[ind_mode]['azimuthal_func'](TH)

plt.figure()
plt.imshow(np.real(R*TH))

• in the radial solver, argument min_radius_bc is now in units of wavelength, defaults to 4.