Текст подпрограммы и версий ammfr_p.zip , ammfe_p.zip |
Тексты тестовых примеров tammfr_p.zip , tammfe_p.zip |
Умножение ленточной матрицы, заданной в компактной форме, на симметричную матрицу, заданную в компактной форме.
АММFR вычисляет произведение ленточной матрицы А порядка N с NLА нижними и NUА верхними ко - диагоналями на симметричную матрицу B порядка N. Матрицы А и B задаются в компактной форме представления. Результатом этого произведения является квадратная матрица С.
procedure AMMFR(var A :Array of Real; N :Integer; NLA :Integer; NUA :Integer; var B :Array of Real; var C :Array of Real);
Параметры
A - | вещественный двумерный массив размера N * (NLА+NUА+1), в котором задается ленточная матрица в компактной форме; |
N - | заданный порядок матриц А и B (тип: целый); |
NLA - NUA | заданное число нижних и верхних ко - диагоналей матрицы А соответственно (тип: целый); |
B - | вещественный вектор длины N * (N+1)/2, в котором задается симметричная матрица B в компактной форме; |
C - | вещественный двумерный массив размера N*N, в котором запоминается матрица произведения. |
Версии
AMMFE - | умножение с расширенной (Extended) точностью ленточной матрицы, заданной в компактной форме, на симметричную матрицу, заданную в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
В подпрограмме АММFE параметры А, B и C должны иметь тип Extended. |
Unit tammfr_p; interface uses SysUtils, Math, { Delphi } LStruct, Lfunc, UtRes_p, AMMFR_p; function tammfr: String; implementation function tammfr: String; var N,NLA,NUA,_i :Integer; C :Array [0..24] of Real; const A :Array [0..19] of Real = ( 0.0,0.0,1.0,1.0,1.0,0.0,1.0,1.0,1.0,1.0,1.0,1.0, 1.0,1.0,1.0,1.0,1.0,1.0,1.0,0.0 ); B :Array [0..14] of Real = ( 1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0, 1.0,1.0,1.0 ); begin Result := ''; N := 5; NLA := 2; NUA := 1; AMMFR(A,N,NLA,NUA,B,C); Result := Result + #$0D#$0A; for _i:=0 to 24 do begin Result := Result + Format('%20.16f ',[C[_i]]); if ( ((_i+1) mod 1)=0 ) then Result := Result + #$0D#$0A; end; Result := Result + #$0D#$0A; UtRes('tammfr',Result); { вывод результатов в файл tammfr.res } exit; end; end. Результаты: | 2. 2. 2. 2. 2. | | 3. 3. 3. 3. 3. | C = | 4. 4. 4. 4. 4. | | 4. 4. 4. 4. 4. | | 3. 3. 3. 3. 3. |