--- a/test/SbpOperators/stencil_test.jl	Wed Jan 19 11:08:43 2022 +0100
+++ b/test/SbpOperators/stencil_test.jl	Wed Jan 19 21:49:34 2022 +0100
@@ -1,6 +1,8 @@
 using Test
 using Sbplib.SbpOperators
 import Sbplib.SbpOperators.Stencil
+import Sbplib.SbpOperators.NestedStencil
+import Sbplib.SbpOperators.scale
 
 @testset "Stencil" begin
     s = Stencil((-2,2), (1.,2.,2.,3.,4.))
@@ -29,3 +31,27 @@
         @test convert(Stencil{Rational}, Stencil(1.,2.,3.,4.,5.; center=2)) == Stencil(1//1,2//1,3//1,4//1,5//1; center=2)
     end
 end
+
+@testset "NestedStencil" begin
+
+    @testset "Constructors" begin
+        s1 = Stencil(-1, 1, 0; center = 1)
+        s2 = Stencil(-1, 0, 1; center = 2)
+        s3 = Stencil( 0,-1, 1; center = 3)
+
+        ns = NestedStencil(CenteredStencil(s1,s2,s3))
+        @test ns isa NestedStencil{Int,3}
+
+        @test CenteredNestedStencil(s1,s2,s3) == ns
+
+        @test NestedStencil(s1,s2,s3, center = 2) == ns
+        @test NestedStencil(s1,s2,s3, center = 1) == NestedStencil(Stencil(s1,s2,s3, center=1))
+
+        @test NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=2) == ns
+
+
+        @testset "Error handling" begin
+
+        end
+    end
+end