Mercurial > repos > public > sbplib_julia

comparison src/Grids/EquidistantGrid.jl @ 352:a18bd337a280 feature/equidistant_grid/subgrid

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add function for getting a one dimensional grid for a given dimension from a equidistant grid
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 27 Sep 2020 14:38:05 +0200
parents 01b851161018
children 8257cc75ea6b
comparison
equal deleted inserted replaced
350:28e71a861531 352:a18bd337a280
64 indices = Tuple.(CartesianIndices(grid.size)) 64 indices = Tuple.(CartesianIndices(grid.size))
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 function subgrid(grid::EquidistantGrid, dim::Integer)
70 size = grid.size[dim]
71 limit_lower = grid.limit_lower[dim]
72 limit_upper = grid.limit_upper[dim]
73
74 return EquidistantGrid(size, limit_lower, limit_upper)
75 end
76 export subgrid
77
69 function pointsalongdim(grid::EquidistantGrid, dim::Integer) 78 function pointsalongdim(grid::EquidistantGrid, dim::Integer)
70 @assert dim<=dimension(grid) 79 @assert dim<=dimension(grid)
71 @assert dim>0 80 @assert dim>0
72 points = collect(range(grid.limit_lower[dim],stop=grid.limit_upper[dim],length=grid.size[dim])) 81 points = collect(range(grid.limit_lower[dim],stop=grid.limit_upper[dim],length=grid.size[dim]))
73 end 82 end