comparison src/Grids/equidistant_grid.jl @ 1842:f2b32da29b73 feature/grids/manifolds

Introduce mapped_grid for ParameterSpace and simplify implementation of equidistant_grid(::Chart)
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 22 Oct 2024 15:32:33 +0200
parents 2a8a2b52a112
children 81559cb7b11c
comparison
equal deleted inserted replaced
1841:1c58005429fd 1842:f2b32da29b73
153 153
154 equidistant_grid(d::Interval, size::Int) = equidistant_grid(limits(d)..., size) 154 equidistant_grid(d::Interval, size::Int) = equidistant_grid(limits(d)..., size)
155 equidistant_grid(hb::HyperBox, dims::Vararg{Int}) = equidistant_grid(limits(hb)..., dims...) 155 equidistant_grid(hb::HyperBox, dims::Vararg{Int}) = equidistant_grid(limits(hb)..., dims...)
156 156
157 function equidistant_grid(c::Chart, dims::Vararg{Int}) 157 function equidistant_grid(c::Chart, dims::Vararg{Int})
158 lg = equidistant_grid(parameterspace(c), dims...) 158 mapped_grid(c, ξ->jacobian(c,ξ), parameterspace(c), dims...)
159 return MappedGrid(
160 lg,
161 map(c,lg),
162 map(ξ->jacobian(c, ξ), lg),
163 )
164 end 159 end
165 160
166 161
167 CartesianBoundary{D,BID} = TensorGridBoundary{D,BID} # TBD: What should we do about the naming of this boundary? 162 CartesianBoundary{D,BID} = TensorGridBoundary{D,BID} # TBD: What should we do about the naming of this boundary?
168 163