Mercurial > repos > public > sbplib_julia
annotate src/Grids/manifolds.jl @ 1903:04c251bccbd4 feature/grids/manifolds
Merge feature/grids/parameter_spaces
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Sat, 01 Feb 2025 22:17:39 +0100 |
| parents | 20bd8887db0d |
| children | e7f8d11c4670 |
| 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 (c::Chart)(ξ) = c.mapping(ξ) |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
13 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
|
14 |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
15 """ |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
16 jacobian(c::Chart, ξ) |
| 1558 | 17 |
|
1581
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
18 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
|
19 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
|
20 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
|
21 chart itself. |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
22 ```julia |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
23 c = Chart(f, ps) |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
24 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
|
25 ``` |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
26 or |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
27 ```julia |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
28 c = Chart(f, ps) |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
29 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
|
30 ``` |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
31 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
|
32 """ |
|
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
33 jacobian(c::Chart, ξ) = jacobian(c.mapping, ξ) |
|
1644
e213bd857f3f
Add some todos and tbds
Jonatan Werpers <jonatan@werpers.com>
parents:
1625
diff
changeset
|
34 # TBD: Can we register a error hint for when jacobian is called with a function that doesn't have a registered jacobian? |
|
1581
f77c5309dd2b
Rename ConcreteChart to Chart and remove the abstarct chart type
Jonatan Werpers <jonatan@werpers.com>
parents:
1580
diff
changeset
|
35 |
| 1558 | 36 |
|
1867
de4b4f2aee4f
Add some tests for CartesianAtlas
Jonatan Werpers <jonatan@werpers.com>
parents:
1844
diff
changeset
|
37 # TBD: Should Charts, parameterspaces, Atlases, have boundary names? |
|
1644
e213bd857f3f
Add some todos and tbds
Jonatan Werpers <jonatan@werpers.com>
parents:
1625
diff
changeset
|
38 |
| 1558 | 39 """ |
| 40 Atlas | |
| 41 | |
| 42 A collection of charts and their connections. | |
|
1844
1987347752ef
Add stub for connections(::Atlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1782
diff
changeset
|
43 Should implement methods for `charts` and `connections`. |
| 1558 | 44 """ |
| 45 abstract type Atlas end | |
| 46 | |
| 47 """ | |
| 48 charts(::Atlas) | |
| 49 | |
| 50 The colloction of charts in the atlas. | |
| 51 """ | |
| 52 function charts end | |
| 53 | |
| 54 """ | |
|
1844
1987347752ef
Add stub for connections(::Atlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1782
diff
changeset
|
55 connections(::Atlas) |
| 1558 | 56 |
| 57 TBD: What exactly should this return? | |
| 58 """ | |
|
1844
1987347752ef
Add stub for connections(::Atlas)
Jonatan Werpers <jonatan@werpers.com>
parents:
1782
diff
changeset
|
59 function connections end |
| 1558 | 60 |
| 61 struct CartesianAtlas <: Atlas | |
| 62 charts::Matrix{Chart} | |
| 63 end | |
| 64 | |
| 65 charts(a::CartesianAtlas) = a.charts | |
|
1782
614f731af685
Add stubs for connections method
Jonatan Werpers <jonatan@werpers.com>
parents:
1781
diff
changeset
|
66 connections(a::CartesianAtlas) = nothing |
| 1558 | 67 |
| 68 struct UnstructuredAtlas <: Atlas | |
| 69 charts::Vector{Chart} | |
| 70 connections | |
| 71 end | |
| 72 | |
| 73 charts(a::UnstructuredAtlas) = a.charts | |
|
1782
614f731af685
Add stubs for connections method
Jonatan Werpers <jonatan@werpers.com>
parents:
1781
diff
changeset
|
74 connections(a::UnstructuredAtlas) = nothing |
| 1558 | 75 |
| 76 | |
| 77 ### | |
| 78 # Geometry | |
| 79 ### | |
| 80 | |
| 81 abstract type Curve end | |
| 82 abstract type Surface end | |
| 83 | |
| 84 | |
| 85 struct Line{PT} <: Curve | |
| 86 p::PT | |
| 87 tangent::PT | |
| 88 end | |
| 89 | |
| 90 (c::Line)(s) = c.p + s*c.tangent | |
| 91 | |
| 92 | |
| 93 struct LineSegment{PT} <: Curve | |
| 94 a::PT | |
| 95 b::PT | |
| 96 end | |
| 97 | |
| 98 (c::LineSegment)(s) = (1-s)*c.a + s*c.b | |
| 99 | |
| 100 | |
|
1625
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
101 function linesegments(ps...) |
|
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
102 return [LineSegment(ps[i], ps[i+1]) for i ∈ 1:length(ps)-1] |
|
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
103 end |
|
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
104 |
|
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
105 |
|
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
106 function polygon_edges(ps...) |
|
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
107 n = length(ps) |
|
1869
20bd8887db0d
Fix Aqua and JET fails
Jonatan Werpers <jonatan@werpers.com>
parents:
1867
diff
changeset
|
108 return [LineSegment(ps[i], ps[mod1(i+1,n)]) for i ∈ eachindex(ps)] |
|
1625
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
109 end |
|
84c3b9d71218
Add linesegments function
Jonatan Werpers <jonatan@werpers.com>
parents:
1581
diff
changeset
|
110 |
| 1558 | 111 struct Circle{T,PT} <: Curve |
| 112 c::PT | |
| 113 r::T | |
| 114 end | |
| 115 | |
|
1869
20bd8887db0d
Fix Aqua and JET fails
Jonatan Werpers <jonatan@werpers.com>
parents:
1867
diff
changeset
|
116 function (C::Circle)(θ) |
|
20bd8887db0d
Fix Aqua and JET fails
Jonatan Werpers <jonatan@werpers.com>
parents:
1867
diff
changeset
|
117 (;c, r) = C |
|
20bd8887db0d
Fix Aqua and JET fails
Jonatan Werpers <jonatan@werpers.com>
parents:
1867
diff
changeset
|
118 c + r*@SVector[cos(θ), sin(θ)] |
|
20bd8887db0d
Fix Aqua and JET fails
Jonatan Werpers <jonatan@werpers.com>
parents:
1867
diff
changeset
|
119 end |
| 1558 | 120 |
| 121 struct TransfiniteInterpolationSurface{T1,T2,T3,T4} <: Surface | |
| 122 c₁::T1 | |
| 123 c₂::T2 | |
| 124 c₃::T3 | |
| 125 c₄::T4 | |
| 126 end | |
| 127 | |
| 128 function (s::TransfiniteInterpolationSurface)(u,v) | |
| 129 c₁, c₂, c₃, c₄ = s.c₁, s.c₂, s.c₃, s.c₄ | |
| 130 P₀₀ = c₁(0) | |
| 131 P₁₀ = c₂(0) | |
| 132 P₁₁ = c₃(0) | |
| 133 P₀₁ = c₄(0) | |
| 134 return (1-v)*c₁(u) + u*c₂(v) + v*c₃(1-u) + (1-u)*c₄(1-v) - ( | |
| 135 (1-u)*(1-v)*P₀₀ + u*(1-v)*P₁₀ + u*v*P₁₁ + (1-u)*v*P₀₁ | |
| 136 ) | |
| 137 end | |
| 138 | |
| 139 function (s::TransfiniteInterpolationSurface)(ξ̄::AbstractArray) | |
| 140 s(ξ̄...) | |
| 141 end | |
| 142 | |
|
1644
e213bd857f3f
Add some todos and tbds
Jonatan Werpers <jonatan@werpers.com>
parents:
1625
diff
changeset
|
143 # TODO: Implement jacobian() for the different mapping helpers |
| 1558 | 144 |