Mercurial > repos > public > sbplib_julia

changeset 795:4ff2dad2dd7f operator_storage_array_of_table

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in constructor for ConstantInteriorScalingOperator when constructing on a grid
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 25 Jul 2021 14:38:04 +0200
parents e0c5eeb3cbbb
children f682e4fe3cef
files src/SbpOperators/volumeops/constant_interior_scaling_operator.jl test/SbpOperators/volumeops/constant_interior_scaling_operator_test.jl
diffstat 2 files changed, 7 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/constant_interior_scaling_operator.jl	Fri Jul 23 10:50:55 2021 +0200
+++ b/src/SbpOperators/volumeops/constant_interior_scaling_operator.jl	Sun Jul 25 14:38:04 2021 +0200
@@ -19,7 +19,7 @@
 end
 
 function ConstantInteriorScalingOperator(grid::EquidistantGrid{1}, interior_weight, closure_weights)
-    return ConstantInteriorScalingOperator(interior_weight, Tuple(closure_weights), size(grid))
+    return ConstantInteriorScalingOperator(interior_weight, Tuple(closure_weights), size(grid)[1])
 end
 
 closure_size(::ConstantInteriorScalingOperator{T,N}) where {T,N} = N
--- a/test/SbpOperators/volumeops/constant_interior_scaling_operator_test.jl	Fri Jul 23 10:50:55 2021 +0200
+++ b/test/SbpOperators/volumeops/constant_interior_scaling_operator_test.jl	Sun Jul 25 14:38:04 2021 +0200
@@ -3,6 +3,7 @@
 using Sbplib.LazyTensors
 using Sbplib.SbpOperators
 import Sbplib.SbpOperators: ConstantInteriorScalingOperator
+using Sbplib.Grids
 
 @testset "ConstantInteriorScalingOperator" begin
     @test ConstantInteriorScalingOperator(1, (2,3), 10) isa ConstantInteriorScalingOperator{Int,2}
@@ -27,4 +28,9 @@
     @test domain_size(a) == (7,)
 
     @test_throws DomainError ConstantInteriorScalingOperator(4,(2,3), 3)
+
+    @testset "Grid constructor" begin
+        g = EquidistantGrid(11, 0., 2.)
+        @test ConstantInteriorScalingOperator(g, 3., (.1,.2)) isa ConstantInteriorScalingOperator{Float64}
+    end
 end