Mercurial > repos > public > sbplib_julia

comparison src/Grids/mapped_grid.jl @ 1570:063a2bfb03da feature/grids/curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bugs in calls to equidistant_grid after merge with default
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 25 Apr 2024 13:57:42 +0200
parents 5d32ecb98db8
children c2620118f110 64452a678e7a
comparison
equal deleted inserted replaced
1560:5d32ecb98db8 1570:063a2bfb03da
51 # TBD: refine and coarsen could be implemented once we have a simple manifold implementation. 51 # TBD: refine and coarsen could be implemented once we have a simple manifold implementation.
52 # Before we do, we should consider the overhead of including such a field in the mapped grid struct. 52 # Before we do, we should consider the overhead of including such a field in the mapped grid struct.
53 53
54 function mapped_grid(x, J, size...) 54 function mapped_grid(x, J, size...)
55 D = length(size) 55 D = length(size)
56 lg = equidistant_grid(size, ntuple(i->0., D), ntuple(i->1., D)) 56 lg = equidistant_grid(ntuple(i->0., D), ntuple(i->1., D), size...)
57 return MappedGrid( 57 return MappedGrid(
58 lg, 58 lg,
59 map(x,lg), 59 map(x,lg),
60 map(J,lg), 60 map(J,lg),
61 ) 61 )