Re: Speedup/Refaktoring loop? Extended...

[Dmitri A. Sergatskov]

On Wed, Aug 27, 2008 at 9:10 AM, Andreas Romeyke <address@hidden> wrote:
> Hello,
>
>>     a=transpose(U) * X;
>>     b=transpose(U) * U * V;
>>     V=V .* a ./ b;
>>     a=X * transpose(V);
>>     b=U * V * transpose(V);
>>     U=U .* a ./ b;
>
> BTW, could anyone describe why Matlab 7.3 runs 4 times faster than
> Octave 3.02 in this piece of code? AFAIK both programs uses the
> Altlas/Lapack libraries. And the code is fully vectorized here.
>
> Could you explain the difference?

Have you tried to ask Mathworks people? After all they do have access to both
Matlab and Octave source code (and we do not).

Also, have you tried to replace transpose(X) with X.'?

>
> Bye Andreas
>

Sincerely,

Dmitri.
--