ArbNumerics.jl
ArbNumerics.jl copied to clipboard
extended precision math, accurate and performant
You might want to support the `base` keyword as in JuliaLang/julia#42428 1. For `precision`, the simplest thing is to overload `Base._precision` instead of `Base.precision`, if `isdefined(Base, :_precision)`; that way you...
`setinterval` and `setball` should accept Float64, Float32.
by doing `promote` stuff simple arithmetic allocates a lot: ```julia julia> z = ArbComplex(2) 2.0 + 0im julia> @time z+1; 0.000023 seconds (11 allocations: 600 bytes) ``` by exploiting `arb`s...
I guess arbs `slong` translates to `Clong` on the julia side: ```julia julia> w = typemax(Cint) +1 2147483648 julia> ArbComplex(w) ERROR: InexactError: trunc(Int32, 2147483648) (...) julia> z = ArbComplex{128}() 0...
integration
Arb has [acb_calc](http://arblib.org/acb_calc.html), where `acb_calc_integrate` is provided. I tried to wrap the function as follows ```julia #src/float/integrate.jl struct acb_calc_integrate_opt_struct deg_limit::Cint #slong eval_limit::Cint #slong depth_limit::Cint #slong use_heap::Cint #int verbose::Cint #int function...
Add test cases for the `ArbComplex`type. A few bug and simplifications were made to the source code, as revealed by test failures, and to simplify test coverage. Some test cases...
to this package (or ArbFloat or ArbReal)?
Are the examples in the README.md correct? With ArbNumerics-1.3.3, I get errors like ```julia julia> ArbFloat(2, 32) ERROR: MethodError: no method matching ArbFloat(::Int64, ::Int64) Closest candidates are: ArbFloat(::T) where T
Maybe I'm doing it wrong, but it looks like there is a problem with airyaiprime and BigFloats: ``` _ _ _(_)_ | Documentation: https://docs.julialang.org (_) | (_) (_) | _...
Several methods are added for `promote_type` in [`src/libarb/promote.jl`](https://github.com/JeffreySarnoff/ArbNumerics.jl/blob/8ad37c2eb61e9006eb09be7b2c35961078245a4a/src/libarb/promote.jl). Xref JuliaLang/julia#54138