Skip to content
Snippets Groups Projects
Commit d30d0418 authored by MichelJuillard's avatar MichelJuillard
Browse files

modify FastLapackInterface API

parent ad1cfae1
Branches master
No related tags found
No related merge requests found
Pipeline #6634 failed
using LinearAlgebra
using FastLapackInterface
using FastLapackInterface.LinSolveAlgo
using FastLapackInterface.QrAlgo
using PolynomialMatrixEquations
#using SolveEyePlusMinusAkronB: EyePlusAtKronBWs, generalized_sylvester_solver!
......@@ -145,7 +143,7 @@ struct BackwardLinearRationalExpectationsWs
colsUE = backward_nbr .+ findall(in(backward_indices), forward_indices)
linsolve_static_ws = LinSolveWs(static_nbr)
AGplusB = Matrix{Float64}(undef, endogenous_nbr, endogenous_nbr)
AGplusB_linsolve_ws = LinSolveAlgo.LinSolveWs(endogenous_nbr)
AGplusB_linsolve_ws = LinSolveWs(endogenous_nbr)
# if m.serially_correlated_exogenous
# eye_plus_at_kron_b_ws = EyePlusAtKronBWs(ma, mb, mc, 1)
# else
......
......@@ -324,7 +324,7 @@ function make_lu_AGplusB!(AGplusB::AbstractMatrix{Float64},
mul!(ws.temp7, A, ws.temp3)
vAG = view(AGplusB, :, ws.backward_indices)
vAG .+= ws.temp7
LinSolveAlgo.lu!(AGplusB, ws.AGplusB_linsolve_ws)
FastLapackInterface.lu!(AGplusB, ws.AGplusB_linsolve_ws)
end
function solve_for_derivatives_with_respect_to_shocks!(results::LinearRationalExpectationsResults,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment