Mercurial > repos > public > sbplib_julia
changeset 627:9f27f451d0a0 feature/volume_and_boundary_operators
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 | 8581610da4f1 |
| children | 316dbfd31d35 |
| files | src/SbpOperators/boundaryops/boundary_operator.jl src/SbpOperators/volumeops/volume_operator.jl |
| diffstat | 2 files changed, 4 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- a/src/SbpOperators/boundaryops/boundary_operator.jl Thu Dec 31 08:11:57 2020 +0100 +++ b/src/SbpOperators/boundaryops/boundary_operator.jl Thu Dec 31 08:13:04 2020 +0100 @@ -10,6 +10,8 @@ the boundary restriction operator in the y-direction direction is `Ix⊗op⊗Iz`. """ function boundary_operator(grid::EquidistantGrid{Dim,T}, closure_stencil::Stencil{T}, boundary::CartesianBoundary) where {Dim,T} + #TODO:Check that dim(boundary) <= Dim? + # Create 1D boundary operator r = region(boundary) d = dim(boundary)
--- a/src/SbpOperators/volumeops/volume_operator.jl Thu Dec 31 08:11:57 2020 +0100 +++ b/src/SbpOperators/volumeops/volume_operator.jl Thu Dec 31 08:13:04 2020 +0100 @@ -9,6 +9,8 @@ the boundary restriction operator in the y-direction direction is `Ix⊗op⊗Iz`. """ function volume_operator(grid::EquidistantGrid{Dim,T}, inner_stencil::Stencil{T}, closure_stencils::NTuple{M,Stencil{T}}, parity, direction) where {Dim,T,M} + #TODO: Check that direction <= Dim? + # Create 1D volume operator in along coordinate direction op = VolumeOperator(restrict(grid, direction), inner_stencil, closure_stencils, parity) # Create 1D IdentityMappings for each coordinate direction