Текст подпрограммы и версий
amm9r_p.zip , amm9e_p.zip 		Тексты тестовых примеров
tamm9r_p.zip , tamm9e_p.zip

Подпрограмма:  AMM9R (модуль AMM9R_p)

Назначение

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

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

AMM9R вычисляет произведение симметричной ленточной матрицы А порядка N с NСА ко - диагоналями на симметричную ленточную матрицу B порядка N с NСВ ко - диагоналями. Матрицы А и B задаются в компактной форме представления. Результатом этого произведения является ленточная матрица C порядка N с М верхними и М нижними ко - диагоналями, где М = min (N - 1, NСА + NСВ). Матрица C запоминается в компактной форме.

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

procedure AMM9R(var A :Array of Real; N :Integer; NCA :Integer;
                var B :Array of Real; NCB :Integer;
                var C :Array of Real); 

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

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

Unit TAMM9R_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMM9R_p;

function TAMM9R: String; 

implementation

function TAMM9R: String;
var
N,NCA,NCB,_i :Integer;
C :Array [0..24] of Real;
const
A :Array [0..9] of Real = ( 0.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 = ( 0.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0 );
begin
Result := '';
N := 5;
NСА := 1;
NСВ := 1;
AMM9R(A,N,NCA,B,NCB,C);
Result := Result + #$0D#$0A;
for _i:=0 to 24 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('TAMM9R',Result);  { вывод результатов в файл TAMM9R.res }
exit;
end;

end.

Результаты:

               |  0  0  2  2  1  |
               |  0  2  3  2  1  |
     C  =   |  1  2  3  2  1  |
               |  1  2  3  2  0  |
               |  1  2  2  0  0  |