Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/volume_operator.jl @ 627:9f27f451d0a0 feature/volume_and_boundary_operators

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add todo for checking valid inputs to boundary_operator and volume_operators
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 31 Dec 2020 08:13:04 +0100
parents 332f65c1abf3
children 86bb3606b215
comparison
equal deleted inserted replaced
626:8581610da4f1 627:9f27f451d0a0
7 When `Dim>1`, the `VolumeOperator` `op` is inflated by the outer product 7 When `Dim>1`, the `VolumeOperator` `op` is inflated by the outer product
8 of `IdentityMappings` in orthogonal coordinate directions, e.g for `Dim=3`, 8 of `IdentityMappings` in orthogonal coordinate directions, e.g for `Dim=3`,
9 the boundary restriction operator in the y-direction direction is `Ix⊗op⊗Iz`. 9 the boundary restriction operator in the y-direction direction is `Ix⊗op⊗Iz`.
10 """ 10 """
11 function volume_operator(grid::EquidistantGrid{Dim,T}, inner_stencil::Stencil{T}, closure_stencils::NTuple{M,Stencil{T}}, parity, direction) where {Dim,T,M} 11 function volume_operator(grid::EquidistantGrid{Dim,T}, inner_stencil::Stencil{T}, closure_stencils::NTuple{M,Stencil{T}}, parity, direction) where {Dim,T,M}
12 #TODO: Check that direction <= Dim?
13
12 # Create 1D volume operator in along coordinate direction 14 # Create 1D volume operator in along coordinate direction
13 op = VolumeOperator(restrict(grid, direction), inner_stencil, closure_stencils, parity) 15 op = VolumeOperator(restrict(grid, direction), inner_stencil, closure_stencils, parity)
14 # Create 1D IdentityMappings for each coordinate direction 16 # Create 1D IdentityMappings for each coordinate direction
15 one_d_grids = restrict.(Ref(grid), Tuple(1:Dim)) 17 one_d_grids = restrict.(Ref(grid), Tuple(1:Dim))
16 Is = IdentityMapping{T}.(size.(one_d_grids)) 18 Is = IdentityMapping{T}.(size.(one_d_grids))