Mercurial > repos > public > sbplib_julia

comparison src/Grids/EquidistantGrid.jl @ 358:64ad8ec0eae0 feature/equidistant_grid/subgrid

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change name from subgrid to restrict
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 28 Sep 2020 22:00:04 +0200
parents 8257cc75ea6b
children dacbcba33d7d
comparison
equal deleted inserted replaced
353:8257cc75ea6b 358:64ad8ec0eae0
65 h = spacing(grid) 65 h = spacing(grid)
66 return broadcast(I -> grid.limit_lower .+ (I.-1).*h, indices) 66 return broadcast(I -> grid.limit_lower .+ (I.-1).*h, indices)
67 end 67 end
68 68
69 """ 69 """
70 subgrid(::EquidistantGrid, dim) 70 restrict(::EquidistantGrid, dim)
71 71
72 Pick out given dimensions from the grid and return a grid for them 72 Pick out given dimensions from the grid and return a grid for them
73 """ 73 """
74 function subgrid(grid::EquidistantGrid, dim) 74 function restrict(grid::EquidistantGrid, dim)
75 size = grid.size[dim] 75 size = grid.size[dim]
76 limit_lower = grid.limit_lower[dim] 76 limit_lower = grid.limit_lower[dim]
77 limit_upper = grid.limit_upper[dim] 77 limit_upper = grid.limit_upper[dim]
78 78
79 return EquidistantGrid(size, limit_lower, limit_upper) 79 return EquidistantGrid(size, limit_lower, limit_upper)
80 end 80 end
81 export subgrid 81 export restrict
82 82
83 function pointsalongdim(grid::EquidistantGrid, dim::Integer) 83 function pointsalongdim(grid::EquidistantGrid, dim::Integer)
84 @assert dim<=dimension(grid) 84 @assert dim<=dimension(grid)
85 @assert dim>0 85 @assert dim>0
86 points = collect(range(grid.limit_lower[dim],stop=grid.limit_upper[dim],length=grid.size[dim])) 86 points = collect(range(grid.limit_lower[dim],stop=grid.limit_upper[dim],length=grid.size[dim]))