Mercurial > repos > public > sbplib_julia
annotate src/Grids/manifolds.jl @ 1998:6dd00ea0511a feature/grids/manifolds
Add check if the logical coordinates are in the parameter space when calling a chart
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Fri, 25 Apr 2025 08:28:34 +0200 |
parents | ebe883807b62 |
children | a1b2453c02c9 |
rev | line source |
---|---|
1558 | 1 """ |
1581
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
2 Chart{D} |
1558 | 3 |
4 A parametrized description of a manifold or part of a manifold. | |
5 """ | |
1581
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
6 struct Chart{D, PST<:ParameterSpace{D}, MT} |
1564
35fe4375b35f
Export things and fix ConcreteChart and Simplex
Jonatan Werpers <jonatan@werpers.com>
parents:
1558
diff
changeset
|
7 mapping::MT |
1558 | 8 parameterspace::PST |
9 end | |
10 | |
1781
a73838c9ef94
Let Chart implement Base.ndims instead of domain_dim
Jonatan Werpers <jonatan@werpers.com>
parents:
1780
diff
changeset
|
11 Base.ndims(::Chart{D}) where D = D |
1581
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
12 parameterspace(c::Chart) = c.parameterspace |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
13 |
1998
6dd00ea0511a
Add check if the logical coordinates are in the parameter space when calling a chart
Jonatan Werpers <jonatan@werpers.com>
parents:
1987
diff
changeset
|
14 function (c::Chart)(ξ) |
6dd00ea0511a
Add check if the logical coordinates are in the parameter space when calling a chart
Jonatan Werpers <jonatan@werpers.com>
parents:
1987
diff
changeset
|
15 if ξ ∉ parameterspace(c) |
6dd00ea0511a
Add check if the logical coordinates are in the parameter space when calling a chart
Jonatan Werpers <jonatan@werpers.com>
parents:
1987
diff
changeset
|
16 throw(DomainError(ξ, "chart was called logical coordinates outside the parameterspace. If this was inteded, use the `mapping` field from the Chart struct instead.")) |
6dd00ea0511a
Add check if the logical coordinates are in the parameter space when calling a chart
Jonatan Werpers <jonatan@werpers.com>
parents:
1987
diff
changeset
|
17 end |
6dd00ea0511a
Add check if the logical coordinates are in the parameter space when calling a chart
Jonatan Werpers <jonatan@werpers.com>
parents:
1987
diff
changeset
|
18 return c.mapping(ξ) |
6dd00ea0511a
Add check if the logical coordinates are in the parameter space when calling a chart
Jonatan Werpers <jonatan@werpers.com>
parents:
1987
diff
changeset
|
19 end |
6dd00ea0511a
Add check if the logical coordinates are in the parameter space when calling a chart
Jonatan Werpers <jonatan@werpers.com>
parents:
1987
diff
changeset
|
20 |
1581
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
21 """ |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
22 jacobian(c::Chart, ξ) |
1558 | 23 |
1581
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
24 The jacobian of the mapping evaluated at `ξ`. This defers to the |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
25 implementation of `jacobian` for the mapping itself. If no implementation is |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
26 available one can easily be specified for either the mapping function or the |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
27 chart itself. |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
28 ```julia |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
29 c = Chart(f, ps) |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
30 jacobian(f::typeof(f), ξ) = f′(ξ) |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
31 ``` |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
32 or |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
33 ```julia |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
34 c = Chart(f, ps) |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
35 jacobian(c::typeof(c),ξ) = f′(ξ) |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
36 ``` |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
37 which will both allow calling `jacobian(c,ξ)`. |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
38 """ |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
39 jacobian(c::Chart, ξ) = jacobian(c.mapping, ξ) |
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
40 |
1946
c63116e2ec8e
Make boundaries() a method of boundary_identifiers
Jonatan Werpers <jonatan@werpers.com>
parents:
1945
diff
changeset
|
41 boundary_identifiers(c::Chart) = boundary_identifiers(parameterspace(c)) |
1558 | 42 |
1644
e213bd857f3f
Add some todos and tbds
Jonatan Werpers <jonatan@werpers.com>
parents:
1625
diff
changeset
|
43 |
1558 | 44 """ |
45 Atlas | |
46 | |
47 A collection of charts and their connections. | |
1844
1987347752ef
Add stub for connections(::Atlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1782
diff
changeset
|
48 Should implement methods for `charts` and `connections`. |
1558 | 49 """ |
50 abstract type Atlas end | |
51 | |
52 """ | |
53 charts(::Atlas) | |
54 | |
55 The colloction of charts in the atlas. | |
56 """ | |
57 function charts end | |
58 | |
59 """ | |
1844
1987347752ef
Add stub for connections(::Atlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1782
diff
changeset
|
60 connections(::Atlas) |
1558 | 61 |
1987
ebe883807b62
Make doc string more clear
Jonatan Werpers <jonatan@werpers.com>
parents:
1964
diff
changeset
|
62 Collection of 2-tuples of multiblock boundary identifiers. |
1558 | 63 """ |
1844
1987347752ef
Add stub for connections(::Atlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1782
diff
changeset
|
64 function connections end |
1558 | 65 |
1948 | 66 |
67 """ | |
68 CartesianAtlas{D,C<:Chart,AT<:AbstractArray{C,D}} <: Atlas | |
69 | |
70 An atlas where the charts are arranged and connected like an array. | |
71 """ | |
1932
ac63820222f8
Make field types of CartesianAtlas concrete
Jonatan Werpers <jonatan@werpers.com>
parents:
1931
diff
changeset
|
72 struct CartesianAtlas{D,C<:Chart,AT<:AbstractArray{C,D}} <: Atlas |
ac63820222f8
Make field types of CartesianAtlas concrete
Jonatan Werpers <jonatan@werpers.com>
parents:
1931
diff
changeset
|
73 charts::AT |
1558 | 74 end |
75 | |
76 charts(a::CartesianAtlas) = a.charts | |
1919
71d218593cac
Add implementation for size to CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1914
diff
changeset
|
77 |
1920
ae83c91286a2
Implement connections for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1919
diff
changeset
|
78 function connections(a::CartesianAtlas) |
ae83c91286a2
Implement connections for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1919
diff
changeset
|
79 c = Tuple{MultiBlockBoundary, MultiBlockBoundary}[] |
ae83c91286a2
Implement connections for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1919
diff
changeset
|
80 |
1928
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
81 for d ∈ 1:ndims(charts(a)) |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
82 Is = eachslice(CartesianIndices(charts(a)); dims=d) |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
83 for i ∈ 1:length(Is)-1 # For each interface between slices |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
84 for jk ∈ eachindex(Is[i]) # For each block in slice |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
85 Iᵢⱼₖ = Tuple(Is[i][jk]) |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
86 Iᵢ₊₁ⱼₖ = Tuple(Is[i+1][jk]) |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
87 push!(c, |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
88 ( |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
89 MultiBlockBoundary{Iᵢⱼₖ, CartesianBoundary{d,UpperBoundary}}(), |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
90 MultiBlockBoundary{Iᵢ₊₁ⱼₖ, CartesianBoundary{d,LowerBoundary}}(), |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
91 ) |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
92 ) |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
93 end |
96b3e9c516fe
Make cartesian atlas work for any dimension
Jonatan Werpers <jonatan@werpers.com>
parents:
1927
diff
changeset
|
94 end |
1920
ae83c91286a2
Implement connections for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1919
diff
changeset
|
95 end |
ae83c91286a2
Implement connections for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1919
diff
changeset
|
96 |
ae83c91286a2
Implement connections for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1919
diff
changeset
|
97 return c |
ae83c91286a2
Implement connections for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1919
diff
changeset
|
98 end |
1558 | 99 |
1948 | 100 """ |
101 boundary_identifiers(a::CartesianAtlas) | |
102 | |
103 All non-connected boundaries of the charts of `a`. | |
104 """ | |
1946
c63116e2ec8e
Make boundaries() a method of boundary_identifiers
Jonatan Werpers <jonatan@werpers.com>
parents:
1945
diff
changeset
|
105 function boundary_identifiers(a::CartesianAtlas) |
1940
1bd0e23499da
Change variable name in boundaries(::CartesianAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1937
diff
changeset
|
106 bs = MultiBlockBoundary[] |
1934
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
107 |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
108 for d ∈ 1:ndims(charts(a)) |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
109 Is = eachslice(CartesianIndices(charts(a)); dims=d) |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
110 |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
111 for (i,b) ∈ ((1,LowerBoundary),(length(Is),UpperBoundary)) # For first and last slice |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
112 for jk ∈ eachindex(Is[i]) # For each block in slice |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
113 Iᵢⱼₖ = Tuple(Is[i][jk]) |
1940
1bd0e23499da
Change variable name in boundaries(::CartesianAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1937
diff
changeset
|
114 push!(bs, |
1934
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
115 MultiBlockBoundary{Iᵢⱼₖ, CartesianBoundary{d,b}}(), |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
116 ) |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
117 end |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
118 end |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
119 end |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
120 |
1940
1bd0e23499da
Change variable name in boundaries(::CartesianAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1937
diff
changeset
|
121 return bs |
1934
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
122 end |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
123 |
1b012e2d5db1
Add function for getting boundaries from CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1932
diff
changeset
|
124 |
1948 | 125 """ |
126 UnstructuredAtlas{C<:Chart, CN<:Tuple{MultiBlockBoundary,MultiBlockBoundary}, ...} <: Atlas | |
127 | |
128 An atlas with connections determined by a vector `MultiBlockBoundary` pairs. | |
129 """ | |
1937
755fc0907e99
Add tetss for construcors and getters for UnstructuredAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1936
diff
changeset
|
130 struct UnstructuredAtlas{C<:Chart, CN<:Tuple{MultiBlockBoundary,MultiBlockBoundary}, CV<:AbstractVector{C}, CNV<:AbstractVector{CN}} <: Atlas |
755fc0907e99
Add tetss for construcors and getters for UnstructuredAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1936
diff
changeset
|
131 charts::CV |
755fc0907e99
Add tetss for construcors and getters for UnstructuredAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1936
diff
changeset
|
132 connections::CNV |
1558 | 133 end |
134 | |
135 charts(a::UnstructuredAtlas) = a.charts | |
1931
bb2c708d724d
Implement connections for unstructured atlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1929
diff
changeset
|
136 connections(a::UnstructuredAtlas) = a.connections |
1936
a7c9d04b57e0
Stub out boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1934
diff
changeset
|
137 |
1948 | 138 """ |
139 boundary_identifiers(a::UnstructuredAtlas) | |
140 | |
141 All non-connected boundaries of the charts of `a`. | |
142 """ | |
1946
c63116e2ec8e
Make boundaries() a method of boundary_identifiers
Jonatan Werpers <jonatan@werpers.com>
parents:
1945
diff
changeset
|
143 function boundary_identifiers(a::UnstructuredAtlas) |
1941
1e3089197397
Add tests for boundaries(::UnstructuredAtlas). Attempt implementation but more Charts functionality is needed
Jonatan Werpers <jonatan@werpers.com>
parents:
1940
diff
changeset
|
144 bs = MultiBlockBoundary[] |
1e3089197397
Add tests for boundaries(::UnstructuredAtlas). Attempt implementation but more Charts functionality is needed
Jonatan Werpers <jonatan@werpers.com>
parents:
1940
diff
changeset
|
145 |
1945
19efb9dacfff
Fix boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1944
diff
changeset
|
146 for (i,c) ∈ enumerate(charts(a)) |
1946
c63116e2ec8e
Make boundaries() a method of boundary_identifiers
Jonatan Werpers <jonatan@werpers.com>
parents:
1945
diff
changeset
|
147 for b ∈ boundary_identifiers(c) |
1945
19efb9dacfff
Fix boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1944
diff
changeset
|
148 mbb = MultiBlockBoundary{i,typeof(b)}() |
19efb9dacfff
Fix boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1944
diff
changeset
|
149 |
19efb9dacfff
Fix boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1944
diff
changeset
|
150 if !any(cn->mbb∈cn, connections(a)) |
19efb9dacfff
Fix boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1944
diff
changeset
|
151 push!(bs, mbb) |
1941
1e3089197397
Add tests for boundaries(::UnstructuredAtlas). Attempt implementation but more Charts functionality is needed
Jonatan Werpers <jonatan@werpers.com>
parents:
1940
diff
changeset
|
152 end |
1e3089197397
Add tests for boundaries(::UnstructuredAtlas). Attempt implementation but more Charts functionality is needed
Jonatan Werpers <jonatan@werpers.com>
parents:
1940
diff
changeset
|
153 end |
1e3089197397
Add tests for boundaries(::UnstructuredAtlas). Attempt implementation but more Charts functionality is needed
Jonatan Werpers <jonatan@werpers.com>
parents:
1940
diff
changeset
|
154 end |
1945
19efb9dacfff
Fix boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1944
diff
changeset
|
155 |
19efb9dacfff
Fix boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1944
diff
changeset
|
156 return bs |
1936
a7c9d04b57e0
Stub out boundaries(::UnstructuredAtlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1934
diff
changeset
|
157 end |