Mercurial > repos > public > sbplib_julia

comparison src/Grids/equidistant_grid.jl @ 1696:29b96fc75bee feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge feature/grids/manifolds
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 28 Aug 2024 10:50:15 +0200
parents a4c52ae93b11
children 03894fd7b132
comparison
equal deleted inserted replaced
1693:c7eee3952dcd 1696:29b96fc75bee
115 gs = map(equidistant_grid, limit_lower, limit_upper, dims) 115 gs = map(equidistant_grid, limit_lower, limit_upper, dims)
116 return TensorGrid(gs...) 116 return TensorGrid(gs...)
117 end 117 end
118 118
119 """ 119 """
120 equidistant_grid(limit_lower::T, limit_upper::T, size::Int) 120 equidistant_grid(limit_lower::Number, limit_upper::Number, size::Int)
121 121
122 Constructs a 1D equidistant grid. 122 Constructs a 1D equidistant grid.
123 """ 123 """
124 function equidistant_grid(limit_lower::T, limit_upper::T, size::Int) where T 124 function equidistant_grid(limit_lower::Number, limit_upper::Number, size::Int)
125 if any(size .<= 0) 125 if size <= 0
126 throw(DomainError("size must be postive")) 126 throw(DomainError("size must be postive"))
127 end 127 end
128 128
129 if any(limit_upper.-limit_lower .<= 0) 129 if limit_upper-limit_lower <= 0
130 throw(DomainError("side length must be postive")) 130 throw(DomainError("side length must be postive"))
131 end 131 end
132
132 return EquidistantGrid(range(limit_lower, limit_upper, length=size)) # TBD: Should it use LinRange instead? 133 return EquidistantGrid(range(limit_lower, limit_upper, length=size)) # TBD: Should it use LinRange instead?
133 end 134 end
134 135
135 136
136 equidistant_grid(hb::HyperBox, dims::Vararg{Int}) = equidistant_grid(hb.a, hb.b, dims...) 137 equidistant_grid(hb::HyperBox, dims::Vararg{Int}) = equidistant_grid(hb.a, hb.b, dims...)