Grids of `AttractorViaRecurrences` unecessarily limit the number type to `Float64`
The grids used in AttractorsViaRecurrences convert all numbers to Float64. E.g., we have:
function RegularGrid(grid::NTuple)
D = length(grid)
grid_steps = SVector{D,Float64}(step.(grid))
grid_maxima = SVector{D,Float64}(maximum.(grid))
grid_minima = SVector{D,Float64}(minimum.(grid))
return RegularGrid(grid_steps, grid_minima, grid_maxima, grid)
end
This is unecessary. We should simply extract the number type T and use vectors of type T. This can impact performance if different number types are used.
Are there a large number of places in DynamicalSystems.jl where this lower precision type specification is hard-coded? I'm asking because I'm experimenting with root finding algorithm basins and I need to use very high precision for the standards of that field. I don't mind changing a few functions by hand but if it's the whole package .... Thank you.
No other place. See also the function 'fixedpoints' from ChaosTools.jl that uses interval arethmetic for max precision