Mercurial > repos > public > sbplib_julia

comparison SbpOperators/src/readoperator.jl @ 260:f89718833620 boundary_conditions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Store the inverse quadrature closure for D2. Implement the stencil application for the inverse quadrature
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 26 Nov 2019 08:18:23 -0800
parents 7cb4492ccd60
children 51e7de109c25
comparison
equal deleted inserted replaced
259:5571d2c5bf0f 260:f89718833620
19 r = (1-i,width-i) 19 r = (1-i,width-i)
20 closureStencils = (closureStencils..., Stencil(r, stencilWeights)) 20 closureStencils = (closureStencils..., Stencil(r, stencilWeights))
21 end 21 end
22 22
23 quadratureClosure = pad_tuple(stringToTuple(Float64, h["closure"][1]), boundarySize) 23 quadratureClosure = pad_tuple(stringToTuple(Float64, h["closure"][1]), boundarySize)
24 inverseQuadratureClosure = 1.0 ./ quadratureClosure
24 eClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["e"][1]), boundarySize)) 25 eClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["e"][1]), boundarySize))
25 dClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["d1"][1]), boundarySize)) 26 dClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["d1"][1]), boundarySize))
26 27
27 d2 = D2( 28 d2 = D2(
28 quadratureClosure, 29 quadratureClosure,
30 inverseQuadratureClosure,
29 innerStencil, 31 innerStencil,
30 closureStencils, 32 closureStencils,
31 eClosure, 33 eClosure,
32 dClosure, 34 dClosure,
33 even 35 even