|
Текст подпрограммы и версий amz4r_p.zip , amz4e_p.zip |
Тексты тестовых примеров tamz4r_p.zip , tamz4e_p.zip |
Преобразование полной симметричной ленточной матрицы в компактную форму представления.
AMZ4R преобразует полную симметричную ленточную матрицу А порядка N с NС ко - диагоналями в компактную форму представления. Тем самым, минимизируется расход памяти, требуемой на хранение симметричной ленточной матрицы.
procedure AMZ4R(var A :Array of Real; N :Integer; NC :Integer;
var B :Array of Real);
Параметры
| A - | вещественный двумерный массив размера N на N, в котором задается симметричная ленточная матрица в полной форме; |
| N - | заданный порядок матрицы (тип: целый); |
| NC - | заданное число ко - диагоналей (тип: целый); |
| B - | вещественный двумерный массив размера N на (NС + 1), в котором запоминается исходная симметричная ленточная матрица в компактной форме представления. |
Версии
| AMZ4E - | преобразование полной симметричной ленточной матрицы, заданной с расширенной (Extended) точностью, в компактную форму представления. |
Вызываемые подпрограммы: нет.
Замечания по использованию
|
Массив B может занимать то же место в памяти, что и массив А (это означает, что первый элемент массива B может совпадать с первым элементом массива А). В подпрограмме АМZ4E параметры А и B должны иметь тип Extended. |
Unit tamz4r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMZ4R_p;
function tamz4r: String;
implementation
function tamz4r: String;
var
NC,N,_i :Integer;
B :Array [0..7] of Real;
const
A :Array [0..15] of Real = ( 1.1,2.1,0.0,0.0,1.2,2.2,3.2,0.0,0.0,2.3,3.3,4.3,
0.0,0.0,3.4,4.4 );
begin
Result := '';
NC := 1;
N := 4;
AMZ4R(A,N,NC,B);
Result := Result + #$0D#$0A;
for _i:=0 to 7 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('tamz4r',Result); { вывод результатов в файл tamz4r.res }
exit;
end;
end.
Результаты: В = ( 0., 2.1, 3.2, 4.3, 1.1, 2.2, 3.3, 4.4 )