Mercurial > repos > public > sbplib_julia
diff src/Grids/EquidistantGrid.jl @ 682:3ed922e95a35 feature/boundary_quads
Make boundary_grid throw if invalid boundary identifier is passed.
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Mon, 08 Feb 2021 20:12:09 +0100 |
| parents | 1d3e29ffc6c6 |
| children | 984f6fe42708 |
line wrap: on
line diff
--- a/src/Grids/EquidistantGrid.jl Mon Feb 08 18:44:44 2021 +0100 +++ b/src/Grids/EquidistantGrid.jl Mon Feb 08 20:12:09 2021 +0100 @@ -122,10 +122,14 @@ grid is a zero-dimensional grid. """ function boundary_grid(grid::EquidistantGrid,id::CartesianBoundary) + dims = collect(1:dimension(grid)) + orth_dims = dims[dims .!= dim(id)] + if orth_dims == dims + throw(DomainError("boundary identifier not matching grid")) + end dims = collect(1:dimension(grid)) orth_dims = dims[dims .!= dim(id)] return restrict(grid,orth_dims) end export boundary_grid boundary_grid(::EquidistantGrid{1},::CartesianBoundary{1}) = EquidistantGrid((),(),()) -boundary_grid(::EquidistantGrid{1},::CartesianBoundary) = throw(DimensionMismatch("dimension of Grid and BoundaryIdentifier not matching"))