Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 48:be53b58581a6

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug with call by reference in Laplace1D apply
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 11 Jan 2019 11:16:22 +0100
parents ef060ab3b035
children 947f7579ba9c
comparison
equal deleted inserted replaced
47:8c4cfa680e5a 48:be53b58581a6
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 h = grid.spacings(L.grid)[1] 33 h = grid.spacings(L.grid)[1]
34 apply!(L.op, u, v, h) 34 apply!(L.op, u, v, h)
35 u .= L.a * u
35 return nothing 36 return nothing
36 end 37 end