//Autor: Artur Czekalski (Sator)  www.epokaY.net/artur  artur@epokaY.net
//Mnożenie macierzy o wyrazach typu single
program Project1;
{$APPTYPE CONSOLE}
uses SysUtils, windows, Math;
//---------------------------------------------------------------------------
const M :Integer = 1000; //liczba wierszy w A
const N :Integer = 700; //liczba kolumn w A i wierszy w B
const K :Integer = 500; //liczba kolumn w B
var
 A: array [0..1000*700-1] of single; //A[i][j]=A[i*N+j]
 B: array [0..700*500-1] of single;
 C: array [0..1000*500-1] of single;
 i, j: Integer;
 w: single;
 czas: Cardinal;
//---------------------------------------------------------------------------
Procedure MnozMacierze(M: Integer; N: Integer; K: Integer;
 var A : array of single; var B : array of single; var C : array of single);
//A[M][N]; B[N][K]; C[M][K];  A[i][j] = A[i*(l.kolumn)+j]; i-wiersz; j-kolumna
//C[i][j] = Suma po 0<x<N : A[i][x]*B[x][j]
var
 S: single;
 i, j, x, w, z: Integer;
begin
 For i:=0 To M-1 Do
  begin z := i*K;
   For j:=0 To K-1 Do
    begin S := 0.0;  w := i*N;
     For x:=0 To N-1 Do
      S := S + A[w+x] * B[x*K+j];
     C[z+j] := S;
    end;
  end;
end;
//---------------------------------------------------------------------------
begin
 w := 0.0;
 //---A[]=
 For i:=0 To M-1 Do //wiersze
  begin
   for j:=0 To N-1 Do //kolukmy
    A[i*N+j] := w;  w := w + 0.01;
  end;
 //---B[]=
 For i:=0 To N-1 Do //wiersze
  begin
   For j:=0 To K-1 Do //kolukmy
    B[i*K+j] := w;  w := w - 0.01;
  end;
 //---
 czas := GetTickCount;
 MnozMacierze(M, N, K, A, B, C);
 czas := GetTickCount - czas;
 //---Spr. C[]
 w := 0.0;
 For i:=0 To M*K-1 Do
  w := w + C[i];
 writeln('Suma wyrazow C[]=', w:0:10);
 //---
 writeln;  writeln('Delphi 7: Czas=', czas);
 readln;
end.