Mercurial > repos > public > sbplib_julia
changeset 583:6db96ef13a29 feature/boundary_ops
Fix bug in boundary_restriction
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Wed, 02 Dec 2020 08:18:09 +0100 |
| parents | aa44edea36fc |
| children | 0e1a95b35999 |
| files | src/SbpOperators/boundaryops/boundary_restriction.jl |
| diffstat | 1 files changed, 2 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
diff -r aa44edea36fc -r 6db96ef13a29 src/SbpOperators/boundaryops/boundary_restriction.jl --- a/src/SbpOperators/boundaryops/boundary_restriction.jl Tue Dec 01 19:41:16 2020 +0100 +++ b/src/SbpOperators/boundaryops/boundary_restriction.jl Wed Dec 02 08:18:09 2020 +0100 @@ -8,7 +8,8 @@ d = dim(boundary) e = BoundaryRestriction(restrict(grid, d), closureStencil, r) - one_d_grids = restrict.(Ref(grid), tuple(1:D)) + one_d_grids = restrict.(Ref(grid), Tuple(1:D)) + Is = IdentityMapping{T}.(size.(one_d_grids)) parts = Base.setindex(Is, e, d) return foldl(⊗, parts)