comparison SbpOperators/src/readoperator.jl @ 272:51e7de109c25 boundary_conditions

Remove redundant property inverseQuadratureClosure
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 05 Dec 2019 11:57:39 +0100
parents f89718833620
children
comparison
equal deleted inserted replaced
271:ecd49ffe0bc8 272:51e7de109c25
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
25 eClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["e"][1]), boundarySize)) 24 eClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["e"][1]), boundarySize))
26 dClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["d1"][1]), boundarySize)) 25 dClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["d1"][1]), boundarySize))
27 26
28 d2 = D2( 27 d2 = D2(
29 quadratureClosure, 28 quadratureClosure,
30 inverseQuadratureClosure,
31 innerStencil, 29 innerStencil,
32 closureStencils, 30 closureStencils,
33 eClosure, 31 eClosure,
34 dClosure, 32 dClosure,
35 even 33 even