Текст подпрограммы и версий 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. )