From 42081aec886db22fc499a53e3909c8fed35e5390 Mon Sep 17 00:00:00 2001
From: Johannes Pfeifer <jpfeifer@gmx.de>
Date: Wed, 8 May 2013 18:18:34 +0200
Subject: [PATCH] Filter out stale loglinear option

With a previously set loglinear option, both k_order_pert and
stochastic_solvers at order>1 transformed the first order solution, but
kept the higher order terms untransformed. Now, both functions output an
error.
---
 matlab/k_order_pert.m       | 4 ++++
 matlab/stochastic_solvers.m | 3 +++
 2 files changed, 7 insertions(+)

diff --git a/matlab/k_order_pert.m b/matlab/k_order_pert.m
index 1622fc2f3..c6a9d38b0 100644
--- a/matlab/k_order_pert.m
+++ b/matlab/k_order_pert.m
@@ -27,6 +27,10 @@ endo_nbr = M.endo_nbr;
 exo_nbr = M.exo_nbr;
 nspred = M.nspred;
 
+if order>1 && options.loglinear 
+   error('The loglinear-option currently only works at order 1') 
+end
+
 switch(order)
   case 1
     [err, g_1] = k_order_perturbation(dr,M,options);
diff --git a/matlab/stochastic_solvers.m b/matlab/stochastic_solvers.m
index 3c6efeedb..c420bc99c 100644
--- a/matlab/stochastic_solvers.m
+++ b/matlab/stochastic_solvers.m
@@ -271,5 +271,8 @@ if options_.loglinear == 1
     dr.ghx = repmat(1./dr.ys(k1),1,size(dr.ghx,2)).*dr.ghx.* ...
              repmat(dr.ys(k1(klag(:,1)))',size(dr.ghx,1),1);
     dr.ghu = repmat(1./dr.ys(k1),1,size(dr.ghu,2)).*dr.ghu;
+    if options_.order>1
+       error('Loglinear options currently only works at order 1') 
+    end
 end