Текст подпрограммы и версий
ammer_p.zip , ammee_p.zip , ammec_p.zip 		Тексты тестовых примеров
tammer_p.zip , tammee_p.zip , tammec_p.zip

Подпрограмма:  AMMER (модуль AMMER_p)

Назначение

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

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

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

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

procedure AMMER(var A :Array of Real; N :Integer; NLA :Integer;
                NUA :Integer; var B :Array of Real; M :Integer;
                var C :Array of Real); 

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

В подпрограмме АММЕE параметры А, И и С должны иметь тип Extended.

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

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

Unit tammer_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMMER_p;

function tammer: String; 

implementation

function tammer: String;
var
N,NLA,NUA,M,_i :Integer;
C :Array [0..14] 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..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;
NLA := 2;
NUA := 1;
M := 3;
AMMER(A,N,NLA,NUA,B,M,C);
Result := Result + #$0D#$0A;
for _i:=0 to 14 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('tammer',Result);  { вывод результатов в файл tammer.res }
exit;
end;

end.


Результаты:

                |  2.  2.  2. |
                |  3.  3.  3. |
     C  =    |  4.  4.  4. |
                |  4.  4.  4. |
                |  3.  3.  3. |