Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 66:543b7a5ab831 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add type parameter for Laplace operand
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 17 Jan 2019 15:32:53 +0100
parents 7054230b639c
children e4fa13137d12
comparison
equal deleted inserted replaced
65:7054230b639c 66:543b7a5ab831
45 u = zeros(eltype(v), size(v)) 45 u = zeros(eltype(v), size(v))
46 apply!(D,v,u) 46 apply!(D,v,u)
47 return u 47 return u
48 end 48 end
49 49
50 struct Laplace{Dim} <: DiffOp 50 struct Laplace{Dim,T<:Real} <: DiffOp
51 grid::Grid.EquidistantGrid{Dim} 51 grid::Grid.EquidistantGrid{Dim,T}
52 a::Real 52 a::T
53 op::D2{Float64} 53 op::D2{Float64}
54 end 54 end
55 55
56 # u = L*v 56 # u = L*v
57 function apply(L::Laplace{1}, v::AbstractVector, i::Int) 57 function apply(L::Laplace{1}, v::AbstractVector, i::Int)