comparison src/Grids/mapped_grid.jl @ 1796:4d55ac2892a2 feature/grids/curvilinear

Improve plotting of grids so that it can handle embedded grids
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 01 Oct 2024 11:24:40 +0200
parents 08710ebe940d
children 25dbb1c02ddc
comparison
equal deleted inserted replaced
1795:a38911d39682 1796:4d55ac2892a2
70 Base.size(g::MappedGrid, d) = size(g.logical_grid, d) 70 Base.size(g::MappedGrid, d) = size(g.logical_grid, d)
71 71
72 boundary_identifiers(g::MappedGrid) = boundary_identifiers(g.logical_grid) 72 boundary_identifiers(g::MappedGrid) = boundary_identifiers(g.logical_grid)
73 boundary_indices(g::MappedGrid, id::TensorGridBoundary) = boundary_indices(g.logical_grid, id) 73 boundary_indices(g::MappedGrid, id::TensorGridBoundary) = boundary_indices(g.logical_grid, id)
74 74
75 # Review: Error when calling plot(boundary_grid(g, id))
76 # Currently need to collect first, i.e., plot(collect(boundary_grid(g, id)))
77 function boundary_grid(g::MappedGrid, id::TensorGridBoundary) 75 function boundary_grid(g::MappedGrid, id::TensorGridBoundary)
78 b_indices = boundary_indices(g.logical_grid, id) 76 b_indices = boundary_indices(g.logical_grid, id)
79 77
80 # Calculate indices of needed jacobian components 78 # Calculate indices of needed jacobian components
81 D = ndims(g) 79 D = ndims(g)