Текст подпрограммы и версий
amz2r_p.zip , amz2e_p.zip , amz2c_p.zip
Тексты тестовых примеров
tamz2r_p.zip , tamz2e_p.zip , tamz2c_p.zip

Подпрограмма:  AMZ2R (модуль AMZ2R_p)

Назначение

Преобразование полной ленточной матрицы в компактную форму представления.

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

AMZ2R преобразует полную ленточную матрицу  А порядка  N с NUС верхними и NLС нижними ко - диагоналями в компактную форму представления. При этом следует иметь в виду, что экономия памяти происходит только тогда, когда (NUС + 1 + NLС) меньше  N.

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

procedure AMZ2R(var A :Array of Real; N :Integer; NUC :Integer;
                NLC :Integer; var B :Array of Real); 

Параметры

A - вещественный двумерный массив размера N на N, в котором задается ленточная матрица в полной форме;
N - заданный порядок матрицы (тип: целый);
NUC - заданное число верхних ко - диагоналей (тип:целый);
NLC - число нижних ко - диагоналей (тип: целый);
B - вещественный двумерный массив размера N на (NUС + 1 + NLС), в котором запоминается исходная ленточная матрица в компактной форме представления.

Версии

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

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

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

 

Массив  B может занимать то же место в памяти, что и массив  А (это означает, что первый элемент массива  B может совпадать с первым элементом массива  А).

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

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

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

Unit tamz2r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMZ2R_p;

function tamz2r: String; 

implementation

function tamz2r: String;
var
N,NLC,NUC,_i :Integer;
B :Array [0..19] of Real;
const
A :Array [0..24] of Real = ( 1.1,2.1,3.1,0.0,0.0,1.2,2.2,3.2,4.2,0.0,0.0,2.3,
3.3,4.3,5.3,0.0,0.0,3.4,4.4,5.4,0.0,0.0,0.0,4.5,
5.5 );
begin
Result := '';
N := 5;
NLC := 2;
NUC := 1;
AMZ2R (A,N,NUC,NLC,B);
Result := Result + #$0D#$0A;
for _i:=0 to 19 do
 begin
  Result := Result + Format('%20.16f ',[B[_i]]);
  if ( ((_i+1) mod 2)=0 )
   then Result := Result + #$0D#$0A;
 end;
Result := Result + #$0D#$0A;
UtRes('tamz2r',Result);  { вывод результатов в файл tamz2r.res }
exit;
end;

end.


Результаты:

       B  =  ( 0., 0., 3.1, 4.2, 5.3, 0., 2.1, 3.2, 4.3, 5.4, 1.1, 2.2, 3.3, 4.4, 
                 5.5, 1.2, 2.3, 3.4, 4.5, 0. )