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
>