Hi Johannes,
Thanks for pointing this out. This is actually Tao's code, so it'd be better to talk to him before changing it.
Michel, would you like to ask Tao about this?
Best, Houtan
On 28 Jun 2013, at 5:44 PM, Johannes Pfeifer wrote:
Dear Houtan,
in csminwel.m of \contrib\ms-sbvar\TZcode\MatlabFiles in line 183 after
gcliff=((f2-f1)/((norm(x2-x1))^2))*(x2-x1);
thee following line is missing:
if(size(x0,2)>1), gcliff=gcliff'; end
Otherwise the code crashes due to non-conformable matrix multiplication. This line is actually contained in our csminwel1.
Best,
Johannes
-- Johannes Pfeifer Haußerstr. 29 72076 Tübingen Tel.: +49-(0)7071-6396184 Mobil.: +49-(0)170-6936820 Germany
Thank you Houtan and Johannes.
Yes, let's add that line.
Thanks for pointing that out.
Tao
--------- On Jun 28, 2013, at 12:34 PM, Houtan Bastani wrote:
Hi Johannes,
Thanks for pointing this out. This is actually Tao's code, so it'd be better to talk to him before changing it.
Michel, would you like to ask Tao about this?
Best, Houtan
On 28 Jun 2013, at 5:44 PM, Johannes Pfeifer wrote:
Dear Houtan,
in csminwel.m of \contrib\ms-sbvar\TZcode\MatlabFiles in line 183 after
gcliff=((f2-f1)/((norm(x2-x1))^2))*(x2-x1);
thee following line is missing:
if(size(x0,2)>1), gcliff=gcliff'; end
Otherwise the code crashes due to non-conformable matrix multiplication. This line is actually contained in our csminwel1.
Best,
Johannes
-- Johannes Pfeifer Haußerstr. 29 72076 Tübingen Tel.: +49-(0)7071-6396184 Mobil.: +49-(0)170-6936820 Germany
_______________________________________________ Dev mailing list Dev@dynare.org https://www.dynare.org/cgi-bin/mailman/listinfo/dev