Mercurial > repos > public > sbplib_julia

comparison sbpD2.jl @ 170:6ba2238a9687 boundary_conditions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Avoid division in apply_d
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 21 May 2019 14:40:07 +0200
parents 6b6d921e8f05
children
comparison
equal deleted inserted replaced
169:24ee4def7ffb 170:6ba2238a9687
101 function apply_e(op::D2, v::AbstractVector, ::Type{Upper}) 101 function apply_e(op::D2, v::AbstractVector, ::Type{Upper})
102 apply(flip(op.eClosure),v,length(v)) 102 apply(flip(op.eClosure),v,length(v))
103 end 103 end
104 104
105 105
106 function apply_d(op::D2, h::Real, v::AbstractVector, ::Type{Lower}) 106 function apply_d(op::D2, h_inv::Real, v::AbstractVector, ::Type{Lower})
107 -apply(op.dClosure,v,1)/h 107 -h_inv*apply(op.dClosure,v,1)
108 end 108 end
109 109
110 function apply_d(op::D2, h::Real, v::AbstractVector, ::Type{Upper}) 110 function apply_d(op::D2, h_inv::Real, v::AbstractVector, ::Type{Upper})
111 -apply(flip(op.dClosure),v,length(v))/h 111 -h_inv*apply(flip(op.dClosure),v,length(v))
112 end 112 end
113 113
114 function readSectionedFile(filename)::Dict{String, Vector{String}} 114 function readSectionedFile(filename)::Dict{String, Vector{String}}
115 f = open(filename) 115 f = open(filename)
116 sections = Dict{String, Vector{String}}() 116 sections = Dict{String, Vector{String}}()