Текст подпрограммы и версий
amm4r_p.zip , amm4e_p.zip , amm4c_p.zip
Тексты тестовых примеров
tamm4r_p.zip , tamm4e_p.zip , tamm4c_p.zip

Подпрограмма:  AMM4R (модуль AMM4R_p)

Назначение

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

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

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

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

procedure AMM4R(var A :Array of Real; N :Integer; M :Integer;
                var B :Array of Real; MLB :Integer; MUB :Integer;
                var C :Array of Real); 

Параметры

A - вещественный двумерный массив размера N*М, в котором задается матрица A;
N - заданное число строк матрицы А (тип: целый);
M - заданное число столбцов матрицы А (тип: целый);
B - вещественный двумерный массив размера М * (МLВ + МUВ + 1), в котором задается ленточная матрица B в компактной форме;
MLB - заданное число нижних ко - диагоналей матрицы B (тип: целый);
MUB - заданное число верхних ко - диагоналей матрицы B (тип: целый);
C - вещественный двумерный массив размера N*M, в котором запоминается матрица произведения.

Версии

AMM4E - умножение с расширенной (Extended) точностью прямоугольной матрицы на ленточную матрицу, заданную в компактной форме.
AMM4C - умножение прямоугольной комплексной матрицы на ленточную комплексную матрицу, заданную в компактной форме.

Вызываемые подпрограммы: нет

Замечания по использованию

 

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

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

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

Unit tamm4r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMM4R_p;

function tamm4r: String; 

implementation

function tamm4r: String;
var
N,M,MLB,MUB,_i :Integer;
C :Array [0..14] of Real;
const
A :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 );
B :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 );
begin
Result := '';
N := 3;
M := 5;
MUB := 1;
MLB := 2;
AMM4R(A,N,M,B,MLB,MUB,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('tamm4r',Result);  { вывод результатов в файл tamm4r.res }
exit;
end;

end.


Результаты:

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