Mercurial > repos > public > sbplib_julia

diff 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
line wrap: on
line diff
--- a/SbpOperators/src/readoperator.jl	Fri Jun 28 14:17:13 2019 +0200
+++ b/SbpOperators/src/readoperator.jl	Tue Nov 26 08:18:23 2019 -0800
@@ -21,11 +21,13 @@
     end
 
     quadratureClosure = pad_tuple(stringToTuple(Float64, h["closure"][1]), boundarySize)
+    inverseQuadratureClosure = 1.0 ./ quadratureClosure
     eClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["e"][1]), boundarySize))
     dClosure = Stencil((0,boundarySize-1), pad_tuple(stringToTuple(Float64, d["d1"][1]), boundarySize))
 
     d2 = D2(
         quadratureClosure,
+        inverseQuadratureClosure,
         innerStencil,
         closureStencils,
         eClosure,