|
Текст подпрограммы и версий amz2r_p.zip , amz2e_p.zip , amz2c_p.zip |
Тексты тестовых примеров tamz2r_p.zip , tamz2e_p.zip , tamz2c_p.zip |
Преобразование полной ленточной матрицы в компактную форму представления.
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. )