Re: searching for the elegant code (matrix indexing and more)

[Juan Pablo Carbajal]

On Wed, Jan 22, 2014 at 5:11 PM, oxy <address@hidden> wrote:
> ... you probably realized a small mistake. Parameter E is also
> an index to the matrix sig:
>
> dim=11;
>  for A=(1:dim)
>      for B=(1:dim)
>          for C=(1:dim)
>              for D=(1:dim)
>                  for E=(1:dim)
>                      sig(A, B, C, D, E) = A*B^2 - C*D^3 +10*E;
>                  endfor
>              endfor
>          endfor
>      endfor
>  endfor
> _______________________________________________
> Help-octave mailing list
> address@hidden
> https://mailman.cae.wisc.edu/listinfo/help-octave

You need to implement tensor products. I never tried with rank-5, you
can for sure reduce the thing t 3 loops (so work with rank-2) but I
can imagine a clever use of kron, permute and reshape can do the
trick. However, it is problably more efficient if you just write the
loops in C and make and oct file :D