Текст подпрограммы и версий
ammbr_p.zip , ammbe_p.zip 		Тексты тестовых примеров
tammbr_p.zip , tammbe_p.zip

Подпрограмма:  AMMBR (модуль AMMBR_p)

Назначение

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

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

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

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

procedure AMMBR(var A :Array of Real; N :Integer; NC :Integer;
                var B :Array of Real; var C :Array of Real); 

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

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

Unit tammbr_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMMBR_p;

function tammbr: String; 

implementation

function tammbr: String;
var
N,NC,_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..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 );
begin
Result := '';
N := 5;
NC := 1;
AMMBR(A,N,NC,B,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('tammbr',Result);  { вывод результатов в файл tammbr.res }
exit;
end;

end.


Результаты:

               |  2  2  2  2  2  |
               |  3  3  3  3  3  |
    C  =    |  3  3  3  3  3  |
               |  3  3  3  3  3  |
               |  2  2  2  2  2  |