Sheehan Olver

Fix stietljes on Line() with square root singularities

1

`1/sqrt(x+im)`

Use AppleAccelerate for plotting?

11

This suggests about a 2x speed increase: julia> r=rand(10000000);@time sqrt(r); 0.049254 seconds (6 allocations: 76.294 MB, 11.57% gc time) julia> r=rand(10000000);@time AppleAccelerate.sqrt(r); 0.026610 seconds (8 allocations: 76.294 MB, 21.70% gc...

Speedup Cauchy on complicated contours

2

One approach is to embed the domain between an inner and outer ellipse (or in the simplest case, an annulus). Then (mapped) Laurent series can be used instead of the...

OffHilbert -> Stieltjes?

2

Should OffHilbert be multiplied by -π and renamed Stieltjes? This matches the function call stieltjes which is now overwritten.

Add dct and r2r

1

It's not clear why `fft` is here but not `dct` or `r2r`.

This avoids the following stack overflow: ```julia [1] circcopy!(::Array{Float64,2}, ::Array{Float64,2}) at ./multidimensional.jl:1071 [2] copy1 at /Users/solver/.julia/packages/AbstractFFTs/JAxy0/src/definitions.jl:38 [inlined] [3] *(::SphericalHarmonics.SphericalHarmonicTransform{Float64}, ::Array{Float64,2}) at /Users/solver/.julia/packages/AbstractFFTs/JAxy0/src/definitions.jl:208 (repeats 51959 times) ``` Haven't had time to...

Add copy

11

This PR adds a copy method. The easiest way to do this was to call `deepcopy`, @ararslan @timholy is there any hidden issue in doing this?

Conversion to complex should throw an error

1

The following behaviour is problematic: ```julia julia> d = convert(AbstractInterval{ComplexF64}, 0..1) 0.0 + 0.0im..1.0 + 0.0im julia> 0.5 in d ERROR: MethodError: no method matching isless(::Complex{Float64}, ::Float64) Closest candidates are:...

Support broadcasting for `+`, `-`, `*`, `/` in 0.7?

2

While it makes sense to _not_ override `+`, etc. to mean interval arithmetic, overriding broadcast notation seems less confusion. I'd propose that broadcasting is meant in the sense of set...