Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/boundary_restriction.jl @ 588:0c411d865d66 feature/boundary_ops

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 02 Dec 2020 13:33:54 +0100
parents 0e1a95b35999
children e40e7439d1b4
comparison
equal deleted inserted replaced
587:6e02c497ccd5 588:0c411d865d66
14 r = region(boundary) 14 r = region(boundary)
15 d = dim(boundary) 15 d = dim(boundary)
16 e = BoundaryRestriction(restrict(grid, d), closureStencil, r) 16 e = BoundaryRestriction(restrict(grid, d), closureStencil, r)
17 17
18 # Create 1D IdentityMappings for each coordinate direction 18 # Create 1D IdentityMappings for each coordinate direction
19 one_d_grids = restrict.(Ref(grid), Tuple(1:D)) 19 one_d_grids = restrict.(Ref(grid), Tuple(1:Dim))
20 Is = IdentityMapping{T}.(size.(one_d_grids)) 20 Is = IdentityMapping{T}.(size.(one_d_grids))
21 21
22 # Formulate the correct outer product sequence of the identity mappings and 22 # Formulate the correct outer product sequence of the identity mappings and
23 # the boundary restriction operator 23 # the boundary restriction operator
24 parts = Base.setindex(Is, e, d) 24 parts = Base.setindex(Is, e, d)