MikaelSlevinsky / SincFun.jl

Sinc numerical methods in Julia
Other
3 stars 0 forks source link

Add new type ConformalMap #12

Closed MikaelSlevinsky closed 9 years ago

MikaelSlevinsky commented 9 years ago

Has data u0::T and u::Vector{T} to unify this and this code.

Have routines evaluate and differentiate. Derivatives could be implemented by updating another parameter in the ConformalMap, say n::Int, that alters the way it's evaluated.

MikaelSlevinsky commented 9 years ago

Closed by https://github.com/MikaelSlevinsky/SincFun.jl/commit/a3e4068592191fcb1d8cd8ecf673a33ba3283bd5.

Now this works!

julia> using SincFun

julia> h = ConformalMap(1.0,rand(5))
ConformalMap{Float64}(1.0,[0.794062,0.0445882,0.237025,0.27319,0.289326],[],0)

julia> x = linspace(-5,5,101);

julia> @time h[x];
elapsed time: 1.8547e-5 seconds (4264 bytes allocated)

julia> hp = differentiate(h)
ConformalMap{Float64}(1.0,[0.0445882,0.474051,0.81957,1.1573],[],1)

julia> @time hp[x];
elapsed time: 1.5722e-5 seconds (4264 bytes allocated)

julia> h3 = differentiate(h,3)
ConformalMap{Float64}(1.0,[1.63914,6.94382],[],3)

julia> @time h3[x];
elapsed time: 2.4084e-5 seconds (4264 bytes allocated)