Normalization of bndtobnd and minor stuff

src/Solvers/Propagators.jl

@@ -150,7 +150,7 @@ end
 """
     function bndtobnd(bndpro, U, dpar, dws, lp)
 
-Returns the boundary to boundary propagator of the Schrodinger functional given that bndpro is the propagator from t = 0 to the bulk, given by the function bndpropagator.
+Returns the boundary to boundary propagator of the Schrodinger functional given that bndpro is the propagator from t = 0 to the bulk, given by the function bndpropagator!.
 
 """
 function bndtobnd(bndpro::AbstractArray, U::AbstractArray, dpar::DiracParam, dws::DiracWorkspace, lp::SpaceParm{4,6,1,D}) where {D}
@@ -160,7 +160,7 @@ function bndtobnd(bndpro::AbstractArray, U::AbstractArray, dpar::DiracParam, dws
         r=Int64(CUDA.blockIdx().x)
 
         if point_time((b, r), lp) == lp.iL[end]
-            psi[b,r] = gdagpmul(Pgamma{4,1},U[b,4,r],bndpro[b,r])
+            psi[b,r] = gdagpmul(Pgamma{4,1},U[b,4,r],bndpro[b,r])/2
         else
             psi[b,r] = 0.0*bndpro[b,r]
         end