@@ -1042,13 +1042,9 @@ function solve_time1_bellman(model::Model{TR}, μgrid::AbstractArray) where {TR
10421042 cf[s_, sprime] = y -> splc(y)
10431043 nf[s_, sprime] = y -> spln(y)
10441044 xprimef[s_, sprime] = y -> splx(y)
1045- # cf[s_, sprime] = LinInterp(x[end:-1:1], c[:, sprime][end:-1:1])
1046- # nf[s_, sprime] = LinInterp(x[end:-1:1], n[:, sprime][end:-1:1])
1047- # xprimef[s_, sprime] = LinInterp(x[end:-1:1], xprimes[:, sprime][end:-1:1])
10481045 end
10491046 splV = CubicSpline(V[end:-1:1], x[end:-1:1])
10501047 Vf[s_] = y -> splV(y)
1051- # Vf[s_] = LinInterp(x[end:-1:1], V[end:-1:1])
10521048 end
10531049
10541050 policies = [cf, nf, xprimef]
@@ -1096,7 +1092,6 @@ function fit_policy_function(T::BellmanEquation_Recursive,
10961092 end
10971093 splV = CubicSpline(PFvec[1,:], xgrid)
10981094 Vf[s_] = y -> splV(y)
1099- # Vf[s_] = LinInterp(xgrid, PFvec[1, :])
11001095 for sprime=1:S
11011096 splc = CubicSpline(PFvec[1 + sprime, :], xgrid)
11021097 spln = CubicSpline(PFvec[1 + S + sprime, :], xgrid)
0 commit comments