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