Текст подпрограммы и версий
amm2r_p.zip , amm2e_p.zip 		Тексты тестовых примеров
tamm2r_p.zip , tamm2e_p.zip

Подпрограмма:  AMM2R (модуль AMM2R_p)

Назначение

Умножение прямоугольной матрицы на симметричную матрицу, заданную в компактной форме.

Математическое описание

AMM2R вычисляет произведение прямоугольной матрицы А размера N на М на симметричную матрицу B порядка М, заданную в компактной форме. Результатом этого произведения является матрица C размера N на М.

Использование 

procedure AMM2R(var A :Array of Real; N :Integer; M :Integer;
                var B :Array of Real; var C :Array of Real); 

   Параметры
A - вещественный двумерный массив размера N*М, в котором задается матрица A;
N - заданное число строк матрицы А (тип: целый);
M - заданное число столбцов матрицы А (тип: целый);
B - вещественный вектор длины М * (M + 1) /2, в котором задается симметричная матрица B в компактной форме;
C - вещественный двумерный массив размера N*М, в котором запоминается матрица произведения. 
   Версии
AMM2E - умножение с расширенной (Extended) точностью прямоугольной матрицы на симметричную, заданную в компактной форме. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
  В подпрограмме АММ2E параметры А, B и C должны иметь тип Extended.

Пример использования 

Unit tamm2r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMM2R_p;

function tamm2r: String; 

implementation

function tamm2r: String;
var
N,M,_i :Integer;
C :Array [0..11] of Real;
const
A :Array [0..11] 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 
);
B :Array [0..9] of Real = ( 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 := 4;
AMM2R(A,N,M,B,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('tamm2r',Result);  { вывод результатов в файл tamm2r.res }
exit;
end;

end.


Результаты:

                    |  4.   4.   4.   4.  |
          S  =    |  4.   4.   4.   4.  |
                    |  4.   4.   4.   4.  |