Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 43:ef060ab3b035

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
remove stride and remove some bugs
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 10 Jan 2019 16:30:54 +0100
parents bb841977d198
children be53b58581a6
comparison
equal deleted inserted replaced
40:8b04efde1a46 43:ef060ab3b035
28 op 28 op
29 end 29 end
30 30
31 # u = L*v 31 # u = L*v
32 function apply!(L::Laplace1D, u::AbstractVector, v::AbstractVector) 32 function apply!(L::Laplace1D, u::AbstractVector, v::AbstractVector)
33 N = closureSize(L.op) 33 h = grid.spacings(L.grid)[1]
34 M = length(v) 34 apply!(L.op, u, v, h)
35
36 h = scaling(L.grid)
37
38 apply!(L.op, u, v, grid.spacings(L.grid)[1], 1, L.grid.numberOfPointsPerDim, stride=1)
39 return nothing 35 return nothing
40 end 36 end