--- a/DiffOps/test/runtests.jl	Wed Jun 26 17:53:15 2019 +0200
+++ b/DiffOps/test/runtests.jl	Wed Jun 26 17:54:32 2019 +0200
@@ -1,4 +1,37 @@
+using Test
 using DiffOps
-using Test
+using Grids
+using SbpOperators
+using RegionIndices
+using LazyTensors
 
 @test_broken false
+
+@testset "BoundaryValue" begin
+    op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
+    g = EquidistantGrid((3,3), (0.0, 0.0), (1.0,1.0))
+
+    e_w = BoundaryValue(op, g, CartesianBoundary{1,Lower}())
+    e_e = BoundaryValue(op, g, CartesianBoundary{1,Upper}())
+    e_s = BoundaryValue(op, g, CartesianBoundary{2,Lower}())
+    e_n = BoundaryValue(op, g, CartesianBoundary{2,Upper}())
+
+    v = [
+        1 2 3;
+        4 5 6;
+        7 8 9.0;
+        10 11 12;
+    ]
+
+    @test e_w  isa TensorMapping{T,2,1} where T
+    @test e_w' isa TensorMapping{T,1,2} where T
+
+    @test domain_size(e_w, (3,2)) == (2,)
+    @test domain_size(e_e, (3,2)) == (2,)
+    @test domain_size(e_s, (3,2)) == (3,)
+    @test domain_size(e_n, (3,2)) == (3,)
+
+
+    @test collect(e_w'*v) == [1,4,7.0]
+
+end