fix chambolle boundary

src/chambolle.jl

@@ -52,7 +52,7 @@ function init(md::ChambolleModel, alg::ChambolleAlgorithm)
     sv = zero(rv)
     p = extend(reshape(reinterpret(SVector{d,Float64}, pv), size(md.g)), StaticKernels.ExtensionNothing())
     r = reshape(reinterpret(Float64, rv), size(md.g))
-    s = extend(reshape(reinterpret(Float64, sv), size(md.g)), StaticKernels.ExtensionConstant(0.))
+    s = extend(reshape(reinterpret(Float64, sv), size(md.g)), StaticKernels.ExtensionReplicate())
 
     z = zero(CartesianIndex{d})
     @inline kf1(pw, gw) = @inbounds divergence(pw) + gw[z]