--- a/SbpOperators/test/runtests.jl	Thu Jun 27 16:05:26 2019 +0200
+++ b/SbpOperators/test/runtests.jl	Thu Jun 27 16:10:25 2019 +0200
@@ -1,4 +1,23 @@
 using SbpOperators
 using Test
 
-@test_broken false
+@testset "apply_quadrature" begin
+    op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
+    h = 0.5
+
+    @test apply_quadrature(op, h, 1.0, 10, 100) == h
+
+    N = 10
+    qc = op.quadratureClosure
+    q = h.*(qc..., ones(N-2*closuresize(op))..., reverse(qc)...)
+    @assert length(q) == N
+
+    for i ∈ 1:N
+        @test apply_quadrature(op, h, 1.0, i, N) == q[i]
+    end
+
+    v = [2.,3.,2.,4.,5.,4.,3.,4.,5.,4.5]
+    for i ∈ 1:N
+        @test apply_quadrature(op, h, v[i], i, N) == q[i]*v[i]
+    end
+end