Mercurial > repos > public > sbplib_julia

annotate src/SbpOperators/boundaryops/boundary_restriction.jl @ 1288:7de1df0aad6a refactor/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add component_type function to Grid
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 03 Mar 2023 15:42:05 +0100
parents 1cc45207817e
children e94ddef5e72f
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
510
db64cfe4d9de Start sketching on 1D boundary restriction operators.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
1 """
1281
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
2 boundary_restriction(g, closure_stencil::Stencil, boundary)
580
0779713f95a2 Add constructor that infers T and N
Jonatan Werpers <jonatan@werpers.com>
parents: 579
diff changeset
3
995
1ba8a398af9c Rename types
Jonatan Werpers <jonatan@werpers.com>
parents: 966
diff changeset
4 Creates boundary restriction operators `e` as `LazyTensor`s on `boundary`
510
db64cfe4d9de Start sketching on 1D boundary restriction operators.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
5
925
6b47a9ee1632 Add functionality for creating operators from a stencil set. Note: Tests are not updated yet!
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 871
diff changeset
6 `e` is the restriction of a grid function to `boundary` using a `Stencil` `closure_stencil`.
936
22c80fb36400 Fix docs and add some references
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 925
diff changeset
7 `e'` is the prolongation of a grid function on `boundary` to the whole grid using the same `closure_stencil`.
1281
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
8 On a one-dimensional grid, `e` is a `BoundaryOperator`. On a multi-dimensional grid, `e` is the inflation of
947
38d1752a9aff Reformat "See also:"
Jonatan Werpers <jonatan@werpers.com>
parents: 936
diff changeset
9 a `BoundaryOperator`.
38d1752a9aff Reformat "See also:"
Jonatan Werpers <jonatan@werpers.com>
parents: 936
diff changeset
10
1104
c0ab81e4c39c Update docs
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1100
diff changeset
11 See also: [`BoundaryOperator`](@ref), [`LazyTensors.inflate`](@ref).
610
e40e7439d1b4 Add a general boundary operator and make BoundaryRestriction a specialization of it.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 588
diff changeset
12 """
1281
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
13 #TODO: Check docstring
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
14 function boundary_restriction(g::TensorGrid, stencil_set::StencilSet, boundary::TensorGridBoundary)
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
15 op = boundary_restriction(g.grids[grid_id(boundary)], stencil_set, boundary_id(boundary))
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
16 return LazyTensors.inflate(op, size(g), grid_id(boundary))
827
beae513ef8b3 Fix bug in boundary_restriction when grids and stencil does not have matching type
Jonatan Werpers <jonatan@werpers.com>
parents: 777
diff changeset
17 end
925
6b47a9ee1632 Add functionality for creating operators from a stencil set. Note: Tests are not updated yet!
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 871
diff changeset
18
1281
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
19 function boundary_restriction(g::EquidistantGrid, stencil_set::StencilSet, boundary)
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
20 closure_stencil = parse_stencil(stencil_set["e"]["closure"])
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
21 converted_stencil = convert(Stencil{eltype(g)}, closure_stencil)
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
22 return BoundaryOperator(g, converted_stencil, boundary)
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1104
diff changeset
23 end