From 75331d7a4f9dc4ab0b4862598bd4d71140a5a5ce Mon Sep 17 00:00:00 2001
From: Ferhat Mihoubi <ferhat.mihoubi@univ-evry.fr>
Date: Fri, 9 Mar 2012 11:45:10 +0100
Subject: [PATCH] The global oo_.steady_state should also be updated in the
homotopy loop
---
matlab/homotopy1.m | 2 +-
matlab/homotopy2.m | 2 +-
matlab/homotopy3.m | 2 +-
3 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/matlab/homotopy1.m b/matlab/homotopy1.m
index 8b8305d82..f0913b86b 100644
--- a/matlab/homotopy1.m
+++ b/matlab/homotopy1.m
@@ -77,5 +77,5 @@ for i=1:step_nbr+1
oo_.exo_steady_state(values(ix,2)) = points(ix,i);
oo_.exo_det_steady_state(values(ixd,2)) = points(ixd,i);
- steady_(M_,options_,oo_);
+ oo_.steady_state = steady_(M_,options_,oo_);
end
diff --git a/matlab/homotopy2.m b/matlab/homotopy2.m
index 3c972a0ee..726a8aac0 100644
--- a/matlab/homotopy2.m
+++ b/matlab/homotopy2.m
@@ -100,6 +100,6 @@ for i = 1:nv
disp([ 'HOMOTOPY mode 2: lauching solver with ' deblank(varname) ' = ' num2str(v) ' ...'])
- steady_(M_,options_,oo_);
+ oo_.steady_state = steady_(M_,options_,oo_);
end
end
diff --git a/matlab/homotopy3.m b/matlab/homotopy3.m
index 9075a4a10..5359680bf 100644
--- a/matlab/homotopy3.m
+++ b/matlab/homotopy3.m
@@ -91,7 +91,7 @@ while iter < step_nbr
old_ss = oo_.steady_state;
try
- steady_(M_,options_,oo_);
+ oo_.steady_state = steady_(M_,options_,oo_);
if length([kplus; kminus]) == nv
return
--
GitLab