In place matrix multiplication?

["chris.danx" <spamoff.danx@ntlworld.com>]

Hi,

Does anyone know how to do this?  At first I thought it'd be something 
along the lines of...

    procedure Multiply (A : in out Matrix;
                        B : in     Matrix) is
       T : Scalar;
       Z : Scalar := Create(0);
    begin
       for rCount in Matrix_Size'range loop
          for cCount in Matrix_Size'range loop
             T := Z;
             for iCount in Matrix'range loop
                T := T + A.Mx(rCount, iCount) * B.Mx(iCount, cCount);
             end loop;
             A.Mx(rCount, iCount) := T;
          end loop;
       end loop;
    end Multiply;

But this is obviously wrong.  Is there an algorithm to do this?  I 
googled for it, but haven't found anything yet.  I gather it's possible 
because some libraries do it.  They just don't say how they do so.



Cheers,
Chris