Remove DEBUG constant (#389)

Author: devmotion

src/AdvancedHMC.jl

@@ -1,7 +1,5 @@
 module AdvancedHMC
 
-const DEBUG = convert(Bool, parse(Int, get(ENV, "DEBUG_AHMC", "0")))
-
 using Statistics: mean, var, middle
 using LinearAlgebra:
     Symmetric, UpperTriangular, mul!, ldiv!, dot, I, diag, cholesky, UniformScaling

src/adaptation/Adaptation.jl

@@ -5,7 +5,7 @@ using LinearAlgebra: LinearAlgebra
 using Statistics: Statistics
 using SimpleUnPack: @unpack, @pack!
 
-using ..AdvancedHMC: DEBUG, AbstractScalarOrVec
+using ..AdvancedHMC: AbstractScalarOrVec
 using DocStringExtensions
 
 """

src/adaptation/stepsize.jl

@@ -102,7 +102,7 @@ function adapt_stepsize!(
     da::NesterovDualAveraging{T},
     α::AbstractScalarOrVec{<:T},
 ) where {T<:AbstractFloat}
-    DEBUG && @debug "Adapting step size..." α
+    @debug "Adapting step size..." α
 
     # Clip average MH acceptance probability
     if α isa AbstractVector
@@ -124,7 +124,7 @@ function adapt_stepsize!(
     x_bar = (one(T) - η_x) * x_bar .+ η_x * x
 
     ϵ = exp.(x)
-    DEBUG && @debug "Adapting step size..." "new ϵ = $ϵ" "old ϵ = $(da.state.ϵ)"
+    @debug "Adapting step size..." new_ϵ = ϵ old_ϵ = da.state.ϵ
 
     # TODO: we might want to remove this when all other numerical issues are correctly handelled
     if any(isnan.(ϵ)) || any(isinf.(ϵ))

src/trajectory.jl

@@ -791,7 +791,7 @@ function find_good_stepsize(
         ϵ′ = direction == 1 ? d * ϵ : 1 / d * ϵ
         z′, H′ = A(h, z, ϵ)
         ΔH = H - H′
-        DEBUG && @debug "Crossing step" direction H′ ϵ "α = $(min(1, exp(ΔH)))"
+        @debug "Crossing step" direction H′ ϵ α = min(1, exp(ΔH))
         if (direction == 1) && !(ΔH > log(a_cross))
             break
         elseif (direction == -1) && !(ΔH < log(a_cross))
@@ -815,7 +815,7 @@ function find_good_stepsize(
         ϵ_mid = middle(ϵ, ϵ′)
         z′, H′ = A(h, z, ϵ_mid)
         ΔH = H - H′
-        DEBUG && @debug "Bisection step" H′ ϵ_mid "α = $(min(1, exp(ΔH)))"
+        @debug "Bisection step" H′ ϵ_mid α = min(1, exp(ΔH))
         if (exp(ΔH) > a_max)
             ϵ = ϵ_mid
         elseif (exp(ΔH) < a_min)