Re: Matrix multiplication using subarrays?


Richard Maine wrote:
I think you are misdiagnosing the problem. Since you haven't given
enough data for us to do the diagnosis, there isn't much help I can
offer.

The problem turns out to be that I had the dimensions of A reversed in
the actual pgm. Hence, my diagnosis of the proximate cause was correct
(I had the advantage of a numerical example), but I hadn't given you
the data you needed to make this diagnosis. However, you provided a
key bit of help, namely, that the submatrix should behave as expected.

As an aside, since you are presumably using f90/f95 anyway (given the
array slice notation and the MATMUL intrinsic), I often find it more
convenient to use an allocatable array for cases like this.

This is the kind of wisdom I was looking for.

.