Текст подпрограммы и версий
ammdr_p.zip , ammde_p.zip , 		Тексты тестовых примеров
tammdr_p.zip , tammde_p.zip ,

Подпрограмма:  AMMDR (модуль AMMDR_p)

Назначение

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

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

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

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

procedure AMMDR(var A :Array of Real; N :Integer; NLA :Integer;
                NUA :Integer; var B :Array of Real; NLB :Integer;
                NUB :Integer; var C :Array of Real); 

   Параметры
A - вещественный двумерный массив размера N* (NLА + NUА + 1), в котором задается ленточная матрица А в компактной форме;
N - заданный порядок матриц А и B (тип: целый);
   NLA -
   NUA   		заданное число нижних и верхних ко - диагоналей матрицы А соответственно (тип: целый);
B - вещественный двумерный массив размера N* (NLВ + NUВ + 1), в котором задается ленточная матрица в компактной форме;
   NLB -
   NUB   		заданное число нижних и верхних ко - диагоналей матрицы B соответственно (тип: целый);
C - вещественный двумерный массив размера N* (NLС + NUС + 1), в котором запоминается ленточная матрица произведения в компактной форме; здесь NLС = min (N - 1, NLА+NLВ) и NUС = min (N - 1, NUА+NUВ). 
   Версии
AMMDE - умножение с расширенной (Extended) точностью ленточных матриц, заданных в компактной форме.
AMMDC - умножение ленточных комплексных матриц, заданных в компактной форме. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
 

В подпрограмме АММDE параметры А, B и C должны иметь тип Extended.

В подпрограмме АММDС параметры А, B и C должны иметь тип Complex.

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

Unit tammdr_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMMDR_p;

function tammdr: String; 

implementation

function tammdr: String;
var
N,NLA,NUA,NLB,NUB,_i :Integer;
C :Array [0..34] of Real;
const
A :Array [0..19] of Real = ( 0.0,0.0,1.0,1.0,1.0,0.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,0.0 );
B :Array [0..19] of Real = ( 0.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,0.0,1.0,1.0,1.0,0.0,0.0 );
begin
Result := '';
N := 5;
NLA := 2;
NUA := 1;
NLB := 1;
NUB := 2;
AMMDR(A,N,NLA,NUA,B,NLB,NUB,C);
Result := Result + #$0D#$0A;
for _i:=0 to 34 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('tammdr',Result);  { вывод результатов в файл tammdr.res }
exit;
end;

end.


Результаты:

               |  0.  0.  0.  2.  2.  2.  1. |
               |  0.  0.  2.  3.  3.  2.  1. |
    C  =    |  0.  2.  3.  4.  3.  2.  0. |
               |  1.  2.  3.  4.  3.  0.  0. |
               |  1.  2.  3.  3.  0.  0.  0. |