DualNumbers.jl
DualNumbers.jl copied to clipboard
Julia package for representing dual numbers and for performing dual algebra
now `div()` and `rem()` are not supported: ``` julia> x = Dual(1.1, 2.2); julia> y = Dual(3.3, 4.4); julia> div(y, x) ERROR: MethodError: no method matching div(::Dual128, ::Dual128, ::RoundingMode{:ToZero}) julia>...
I'd really like to be able to differentiate complex functions with ForwardDiff. For almost all analytic functions, the derivative can be taken in the same way as in the real...
Well, I looked at https://github.com/JuliaDiff/DualNumbers.jl/issues/45, so I thought this PR might help. The code is basically the same, but I exported `Dual`, `value` and `partials`. I also added pretty printing,...
https://github.com/JuliaDiff/DualNumbers.jl/blob/4603cc1e07dc6363d394875046bca74cd7ce6fc6/src/dual.jl#L178
This matches the docstring of `complex`, so now ```julia julia> typeof(complex(zero(Dual128))) == complex(Dual128) true ```
This package works a treat for evaluating the derivative functions at complex locations: ``` julia> using DualNumbers julia> f(x) = (1, -2 *x) .* exp(-x^2) # value and derivative f...
According to the Julia Base docs, `==` for floating-point types should fail for `NaN`s but not for zeros with opposite signs and `isequal` should have the opposite behaviour for these...
I noticed, in Julia 1.0.0, the following: ```julia julia> @show imε; imε = imɛ ``` The output produced by the `@show` macro uses 2 different characters for the epsilon: first,...