From d8e9745762b56217b83d87f550d3993c47fe1c75 Mon Sep 17 00:00:00 2001 From: RJDennis Date: Sat, 29 Apr 2023 19:55:19 +0100 Subject: [PATCH] Update solution_functions.jl --- src/solution_functions.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/solution_functions.jl b/src/solution_functions.jl index c939d93..60f6b7d 100644 --- a/src/solution_functions.jl +++ b/src/solution_functions.jl @@ -685,7 +685,7 @@ function solve_fourth_order_det(model::REModel,scheme::PerturbationScheme) first_derivs[i,:] .= first_d(point) #second_d(x) = ForwardDiff.hessian(model_equations[i],x,ForwardDiff.HessianConfig(model_equations[i],x,ForwardDiff.Chunk{1}()))[:,1:2*nv] - second_d(x) = ForwardDiff.jacobian(first_d,x,ForwardDiff.JacobianConfig(first_d,x,ForwardDiff.Chunk{1}()))[1:2*n,1:2*n] + second_d(x) = ForwardDiff.jacobian(first_d,x,ForwardDiff.JacobianConfig(first_d,x,ForwardDiff.Chunk{1}()))[1:2*nv,1:2*nv] second_derivs[i,:] .= vec(second_d(point)) third_d(x) = ForwardDiff.jacobian(second_d,x,ForwardDiff.JacobianConfig(second_d,x,ForwardDiff.Chunk{1}()))[:,1:2*nv]