diff --git a/src/DynareModel.jl b/src/DynareModel.jl index f9d84fd8383fad7bf6b693a02a7bc84b79ede838..651a6bfb50aea3960a60f37f50a311884ed140d8 100644 --- a/src/DynareModel.jl +++ b/src/DynareModel.jl @@ -55,6 +55,7 @@ struct AuxVars orig_lead_lag::Int eq_nbr::Int orig_expr::String + unary_op::String end struct PredVars diff --git a/src/DynareSolvers.jl b/src/DynareSolvers.jl index 038f3e6b36985b54571281325546688d476be1a4..1f83cc538eee96b0b2a7194f0cc1aa9ccd10e9dc 100644 --- a/src/DynareSolvers.jl +++ b/src/DynareSolvers.jl @@ -189,6 +189,7 @@ function trustregion(f!::Function, j!::Function, x0::Vector{Float64}, factor::Fl # Initialize scale parameter. scale, scale0 = one(Float64), one(Float64) # Newton iterations + δ = 0.0 while iter<=maxiter && info==0 # Compute columns norm for the Jacobian matrix. @inbounds for i=1:n diff --git a/src/SteadyState.jl b/src/SteadyState.jl index 6de7f8eb812401d9d540e809647bb70305b8993f..898bdf5241ccfadcda666d26e555591d02659dd0 100644 --- a/src/SteadyState.jl +++ b/src/SteadyState.jl @@ -19,6 +19,7 @@ module SteadyState # along with Dynare. If not, see <http://www.gnu.org/licenses/>. ## +using LinearAlgebra using DynareSolvers using Printf @@ -120,7 +121,7 @@ function display_steady_state(model::Model, oo::Output) print(aa, " ", model.endo[i].name, "\t=\t", @sprintf("%.6f", tmp)) println(aa) end - print(STDOUT, aa) + print(aa) end function display_steady_state(model::Model, oo::Output, ys::Vector{Float64}) @@ -135,7 +136,7 @@ function display_steady_state(model::Model, oo::Output, ys::Vector{Float64}) print(aa, " ", model.endo[i].name, "\t=\t", @sprintf("%.6f", tmp)) println(aa) end - print(STDOUT, aa) + print(aa) end function issteadystate(model::Model, oo::Output, ys::Vector{Float64}) diff --git a/src/solvers/CyclicReduction.jl b/src/solvers/CyclicReduction.jl index e8da52bb48e8984c65a2802bdc9414e73cf4f523..3a48d8eafd07fba9dfbb8c6faae7e1d8b9791dbb 100644 --- a/src/solvers/CyclicReduction.jl +++ b/src/solvers/CyclicReduction.jl @@ -1,5 +1,7 @@ module CyclicReduction +using LinearAlgebra + import ...DynLinAlg.LinSolveAlgo: LinSolveWS, linsolve_core! import Base.LinAlg.BLAS: scal!, gemm!