Текст подпрограммы и версий
amz6r_p.zip , amz6e_p.zip
Тексты тестовых примеров
tamz6r_p.zip , tamz6e_p.zip

Подпрограмма:  AMZ6R (модуль AMZ6R_p)

Назначение

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

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

Симметричная ленточная матрица  А порядка  N с NС ко - диагоналями, заданная в компактной форме, преобразуется в компактную форму представления симметричной матрицы.

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

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

Параметры

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

Версии

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

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

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

 

Массив  А и вектор  B должны занимать разное место в памяти.

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

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

Unit tamz6r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMZ6R_p;

function tamz6r: String; 

implementation

function tamz6r: String;
var
N,NC,_i :Integer;
B :Array [0..9] of Real;
const
A :Array [0..7] of Real = ( 0.0,2.1,3.2,4.3,1.1,2.2,3.3,4.4 );
begin
Result := '';
N := 4;
NC := 1;
AMZ6R(A,N,NC,B);
Result := Result + #$0D#$0A;
for _i:=0 to 9 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('tamz6r',Result);  { вывод результатов в файл tamz6r.res }
exit;
end;

end.


Результаты:    В = ( 1.1, 2.1, 2.2, 0., 3.2, 3.3, 0., 0., 4.3, 4.4 )