From 7aa93c068492ff8b9b7e7769635234a55e1158dd Mon Sep 17 00:00:00 2001
From: adjemian <adjemian@ac1d8469-bf42-47a9-8791-bf33cf982152>
Date: Mon, 19 Dec 2005 15:18:19 +0000
Subject: [PATCH] Change the call to resol.m (without systematically using
ys_).
git-svn-id: https://www.dynare.org/svn/dynare/dynare_v3@563 ac1d8469-bf42-47a9-8791-bf33cf982152
---
matlab/EvalWelfare.m | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/matlab/EvalWelfare.m b/matlab/EvalWelfare.m
index 58ff19fe2c..fabb40d9a4 100644
--- a/matlab/EvalWelfare.m
+++ b/matlab/EvalWelfare.m
@@ -47,6 +47,7 @@ deep = get_posterior_parameters('posterior_mean');
% [2.3] Compute the posterior distribution of Welfare:
hfid = waitbar(0,'Posterior welfare distribution...');
compt = 0;
+ys = ys_;
for i=1:B
linea = 1+floor(rand*TotalNumberOfDraws);
tmp = find(METRO(:,3)<linea);
@@ -62,9 +63,10 @@ for i=1:B
DEEP = x2(linee,:);
deep(subindx) = DEEP(subindx);
set_parameters(deep);
- [dr,info] = resol(ys_,0);
+ [dr,info] = resol(ys,0);
if ~info(1)
- WelfDistribution(i) = dr.ys(kk(windx))+.5*dr.ghs2(windx);
+ ys = dr.ys;
+ WelfDistribution(i) = ys(kk(windx))+.5*dr.ghs2(windx);
else
WelfDistribution(i) = Inf;
compt = compt+1;
--
GitLab