stupid bug

constants/juobs_const.jl

@@ -101,36 +101,11 @@ rm(beta::Float64)         = RM[b_values .== beta][1]
 ba_bp(beta::Float64)      = BA_BP[b_values .== beta][1]
 Z(beta::Float64)          = ZZ[b_values .== beta][1]
 
-
-@doc raw""" 
-  ca(beta:Float64; verbose::Bool = true, acceptance:Float64=10^-2)
-
-it compute the ca coefficient of the axial vector current improvement
-As of now, check whether g_0^2 < 1.8+acceptance and if true give ca according to arXiv:1502.04999
-otherwise it gives 0
-The default value of acceptance has no real justification at the moment
-If verbose = true, it prints a warning if the g_0^2 is > 1.8  
-
-"""
-function ca(beta::Float64; verbose::Bool = true, acceptance::Float64 = 10^-2)
+function ca(beta::Float64)
   g02 = 6/beta; 
-  if (1.8-g02)>-acceptance # for g_0^2 about 1.8 and below it work. I decided on a 10^-2 acceptance range
   p0 = 9.2056
   p1 = -13.9847
   e= exp(p0+p1/g02)
   return -0.006033*g02*(1+e)
-  end
-  if verbose
-    @warn raw"""
-    no valid expression for g_0^2>1.8
-    your g_0^2 = $(g02)
-    return value set to 0. 
-    If you want to use this formula, then modify the acceptance value by calling
-    ```@example
-    ca(beta, acceptance=<your value>) 
-    ``` 
-    """
-  end
-  return 0.0
 end
 

src/juobs_obs.jl

@@ -327,7 +327,7 @@ function mpcac_ren(a0p::Corr, pp::Corr, plat::Vector{Int64}; ca::Float64=0.0,ba_
         error("mu or kappa values does not match")
     end
 end
-
+ =#
 ## Decay constants
 @doc raw"""
     dec_const(a0p::Vector{uwreal}, pp::Vector{uwreal}, plat::Vector{Int64}, m::uwreal, y0::Int64; ca::Float64=0.0, pl::Bool=true, data::Bool=false, wpm::Union{Dict{Int64,Vector{Float64}},Dict{String,Vector{Float64}}, Nothing}=nothing)