Mercurial > repos > public > sbplib_julia

comparison src/Grids/curvilinear_grid.jl @ 1504:976f5784d7b9 feature/grids/curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Spelling
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 16 Feb 2024 14:31:45 +0100
parents 704a84eef8b6
children 63101a1cd0e6
comparison
equal deleted inserted replaced
1503:704a84eef8b6 1504:976f5784d7b9
30 boundary_indices(g::CurvilinearGrid, id::TensorGridBoundary) = boundary_indices(g.logicalgrid, id) 30 boundary_indices(g::CurvilinearGrid, id::TensorGridBoundary) = boundary_indices(g.logicalgrid, id)
31 31
32 function boundary_grid(g::CurvilinearGrid, id::TensorGridBoundary) 32 function boundary_grid(g::CurvilinearGrid, id::TensorGridBoundary)
33 b_indices = boundary_indices(g.logicalgrid, id) 33 b_indices = boundary_indices(g.logicalgrid, id)
34 34
35 # Calculate indices of needed jacobian combonents 35 # Calculate indices of needed jacobian components
36 D = ndims(g) 36 D = ndims(g)
37 all_indices = SVector{D}(1:D) 37 all_indices = SVector{D}(1:D)
38 free_variable_indices = deleteat(all_indices, grid_id(id)) 38 free_variable_indices = deleteat(all_indices, grid_id(id))
39 jacobian_components = (:, free_variable_indices) 39 jacobian_components = (:, free_variable_indices)
40 40