Mercurial > repos > public > sbplib_julia

changeset 820:b4acd25943f4 refactor/sbp_operators_method_signatures

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove some more types and type parameters
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 13 Jan 2022 08:27:36 +0100
parents f378cc64bafd
children 9b2910a6411d
files src/SbpOperators/boundaryops/boundary_operator.jl src/SbpOperators/volumeops/derivatives/second_derivative.jl src/SbpOperators/volumeops/laplace/laplace.jl
diffstat 3 files changed, 7 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/boundaryops/boundary_operator.jl	Thu Jan 13 07:42:49 2022 +0100
+++ b/src/SbpOperators/boundaryops/boundary_operator.jl	Thu Jan 13 08:27:36 2022 +0100
@@ -9,7 +9,7 @@
 of `IdentityMappings` in orthogonal coordinate directions, e.g for `Dim=3`,
 the boundary restriction operator in the y-direction direction is `Ix⊗op⊗Iz`.
 """
-function boundary_operator(grid::EquidistantGrid{Dim,T}, closure_stencil::Stencil{T}, boundary::CartesianBoundary) where {Dim,T}
+function boundary_operator(grid::EquidistantGrid, closure_stencil, boundary::CartesianBoundary)
     #TODO:Check that dim(boundary) <= Dim?
 
     # Create 1D boundary operator
@@ -18,8 +18,8 @@
     op = BoundaryOperator(restrict(grid, d), closure_stencil, r)
 
     # Create 1D IdentityMappings for each coordinate direction
-    one_d_grids = restrict.(Ref(grid), Tuple(1:Dim))
-    Is = IdentityMapping{T}.(size.(one_d_grids))
+    one_d_grids = restrict.(Ref(grid), Tuple(1:dimension(grid)))
+    Is = IdentityMapping{eltype(grid)}.(size.(one_d_grids))
 
     # Formulate the correct outer product sequence of the identity mappings and
     # the boundary operator
--- a/src/SbpOperators/volumeops/derivatives/second_derivative.jl	Thu Jan 13 07:42:49 2022 +0100
+++ b/src/SbpOperators/volumeops/derivatives/second_derivative.jl	Thu Jan 13 08:27:36 2022 +0100
@@ -1,6 +1,5 @@
 """
-    second_derivative(grid::EquidistantGrid{Dim}, inner_stencil, closure_stencils, direction)
-    second_derivative(grid::EquidistantGrid{1}, inner_stencil, closure_stencils)
+    second_derivative(grid::EquidistantGrid, inner_stencil, closure_stencils, direction)
 
 Creates the second-derivative operator `D2` as a `TensorMapping`
 
@@ -12,7 +11,7 @@
 one-dimensional operator with the `IdentityMapping`s in orthogonal coordinate dirrections.
 Also see the documentation of `SbpOperators.volume_operator(...)` for more details.
 """
-function second_derivative(grid::EquidistantGrid{Dim}, inner_stencil, closure_stencils, direction) where Dim
+function second_derivative(grid::EquidistantGrid, inner_stencil, closure_stencils, direction)
     h_inv = inverse_spacing(grid)[direction]
     return SbpOperators.volume_operator(grid, scale(inner_stencil,h_inv^2), scale.(closure_stencils,h_inv^2), even, direction)
 end
--- a/src/SbpOperators/volumeops/laplace/laplace.jl	Thu Jan 13 07:42:49 2022 +0100
+++ b/src/SbpOperators/volumeops/laplace/laplace.jl	Thu Jan 13 08:27:36 2022 +0100
@@ -11,9 +11,9 @@
 multi-dimensional `grid`, `Δ` is the sum of multi-dimensional `second_derivative`s
 where the sum is carried out lazily.
 """
-function laplace(grid::EquidistantGrid{Dim}, inner_stencil, closure_stencils) where Dim
+function laplace(grid::EquidistantGrid, inner_stencil, closure_stencils)
     Δ = second_derivative(grid, inner_stencil, closure_stencils, 1)
-    for d = 2:Dim
+    for d = 2:dimension(grid)
         Δ += second_derivative(grid, inner_stencil, closure_stencils, d)
     end
     return Δ