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