|
Текст подпрограммы и версий amm1r_p.zip , amm1e_p.zip , |
Тексты тестовых примеров tamm1r_p.zip , tamm1e_p.zip , |
Умножение прямоугольных матриц.
AMM1R вычисляет произведение прямоугольной матрицы А размера N на М на прямоугольную матрицу B размера М на L. Результатом этого произведения является прямоугольная матрица C размера N на L.
procedure AMM1R(var A :Array of Real; N :Integer; M :Integer;
var B :Array of Real; L :Integer;
var C :Array of Real);
Параметры
| A - | вещественный двумерный массив размера N*М, в котором задается матрица A; |
| N - | заданное число строк матрицы А (тип: целый); |
| M - | заданное число столбцов матрицы А (тип: целый); |
| B - | вещественный двумерный массив размера М*L, в котором задается матрица B; |
| L - | заданное число столбцов матрицы B (тип: целый); |
| C - | вещественный двумерный массив размера N*L, в котором запоминается матрица произведения. |
Версии
| AMM1E - | умножение прямоугольных матриц с повыщенной точностью. |
| AMM1C - | умножение прямоугольных комплексных матриц. |
Вызываемые подпрограммы: нет
Замечания по использованию
|
В подпрограмме АММ1E параметры А, B и C должны иметь тип Extended. В подпрограмме AMM1C параметры А, B и C должны иметь тип Complex. |
Unit tamm1r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMM1R_p;
function tamm1r: String;
implementation
function tamm1r: String;
var
N,M,L,_i :Integer;
C :Array [0..11] of Real;
const
A :Array [0..14] of Real = ( 1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,
1.0,1.0,1.0 );
B :Array [0..19] of Real = ( 1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,
1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0 );
begin
Result := '';
N := 3;
M := 5;
L := 4;
AMM1R(A,N,M,B,L,C);
Result := Result + #$0D#$0A;
for _i:=0 to 11 do
begin
Result := Result + Format('%20.16f ',[C[_i]]);
if ( ((_i+1) mod 3)=0 )
then Result := Result + #$0D#$0A;
end;
Result := Result + #$0D#$0A;
UtRes('tamm1r',Result); { вывод результатов в файл tamm1r.res }
exit;
end;
end.
Результаты:
| 5. 5. 5. 5. |
C = | 5. 5. 5. 5. |
| 5. 5. 5. 5. |