diff --git a/src/lib/Parametric_Contractor.jl b/src/lib/Parametric_Contractor.jl index 0703a4b..3ea24aa 100644 --- a/src/lib/Parametric_Contractor.jl +++ b/src/lib/Parametric_Contractor.jl @@ -416,7 +416,6 @@ function PIn_NewtonGS(X0::Vector{Interval{T}},P::Vector{Interval{T}}, hj!::Function,h!::Function, opt::Vector{Any},Eflag::Bool, Iflag::Bool,eDflag::Bool) where {T<:AbstractFloat} - println(" ---------- begin in place newton ------------- ") # unpacks option file kmax::Int64 = opt[1] etol::Float64 = opt[2] @@ -599,7 +598,6 @@ function PIn_NewtonGS(X0::Vector{Interval{T}},P::Vector{Interval{T}}, Eflag = true end Xtemp = copy(X) - println(" ---------- end in place newton ------------- ") return X,Xtemp,Eflag,Iflag,eDflag,inclusionLow,inclusionHigh end diff --git a/src/lib/Parametric_Utility.jl b/src/lib/Parametric_Utility.jl index 8087cdc..b87aa3e 100644 --- a/src/lib/Parametric_Utility.jl +++ b/src/lib/Parametric_Utility.jl @@ -178,43 +178,29 @@ function extProcess(N::Interval{T},X::Interval{T},Mii::Interval{T}, Ntemp::Interval{T} = copy(N) M::Interval{T} = (B+S1+S2)+Interval(-rtol,rtol) if (M.lo<=0 && M.hi>=0) - println("branch1") return 0, Interval(-Inf,Inf), Ntemp end if (v == 1) - println("Mii: $Mii") k,IML::Interval{T},IMR::Interval{T} = extDivide(Mii) - println("k: $k") if (k == 1) - println("branch2") return 0, (mid(X)-M*IML), Ntemp elseif (k == 2) - println("branch3") return 0, (mid(X)-M*IMR), Ntemp elseif (k == 3) NR = mid(X)-M*IMR NL = mid(X)-M*IML - print("NR: $NR") - print("NL: $NL") if (~isdisjoint(NL,X) && isdisjoint(NR,X)) - println("branch4") return 0, NL, Ntemp elseif (~isdisjoint(NR,X) && isdisjoint(NL,X)) - println("branch5") return 0, NR, Ntemp elseif (~isdisjoint(NL,X) && ~isdisjoint(NR,X)) - println("branch6") N = NL Ntemp = NR - print("NR: $NR") - print("NL: $NL") return 1, NL, NR else - println("branch7") return -1, N, Ntemp end end end - println("branch8") return 0, N, Ntemp end