|
|
Blitz Support : |
From: Matthias Schillinger (matthias_at_[hidden])
Date: 2003-08-29 10:21:30
Hi Fernando,
declare j as thirdIndex and it should work, I think.
See also the user guide on page 57.
Regards,
Matthias
> // 2d version
> void mat_ten_inner0(Array<double,2>& M,
> Array<double,2>& T, Array<double,2>& U ) {
> firstIndex i0;
> secondIndex j;
>
> secondIndex i1;
> // We want for(i0,i1) U(i0,i1) = sum(...,j): j is the sum index
> U = sum(M(i0,j)*T(i1,j),j);
> }