radial.jl 717 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
"""
    Radial( nodes::NodesType, epsilon, gamma )

Radial inteprolation

"""
mutable struct Radial <: InterpolationType

    nx         :: Int64
    epsilon    :: Float64
    xk         :: Array{Float64,1}
    colloc_mat :: Array{Float64,2}

Pierre Navaro's avatar
Pierre Navaro committed
14
    function Radial( nodes::NodesType, epsilon::Real )
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
        
        xk         = nodes.xk
        nx         = nodes.nx
        colloc_mat = evaluate_radial(xk, xk, nx, epsilon)
        new( nx, epsilon, xk, colloc_mat )
        
    end

end

"""
    Radial( xe )

returns evaluation matrix
"""
function (interp::Radial)( xe::Array{Float64,1} )

    nx      = interp.nx
    epsilon = interp.epsilon
    xk      = interp.xk
    evaluate_radial(xk, xe, nx, epsilon)

end