Текст подпрограммы и версий amm4r_p.zip , amm4e_p.zip , amm4c_p.zip |
Тексты тестовых примеров tamm4r_p.zip , tamm4e_p.zip , tamm4c_p.zip |
Умножение прямоугольной матрицы на ленточную матрицу, заданную в компактной форме.
AMM4R вычисляет произведение прямоугольной матрицы А размера N*М на ленточную матрицу B порядка М с МLВ нижними и МUВ верхними ко - диагоналями. Матрица B задается в компактной форме представления. Результатом этого произведения является прямоугольная матрица C размера N*М.
procedure AMM4R(var A :Array of Real; N :Integer; M :Integer; var B :Array of Real; MLB :Integer; MUB :Integer; var C :Array of Real);
Параметры
A - | вещественный двумерный массив размера N*М, в котором задается матрица A; |
N - | заданное число строк матрицы А (тип: целый); |
M - | заданное число столбцов матрицы А (тип: целый); |
B - | вещественный двумерный массив размера М * (МLВ + МUВ + 1), в котором задается ленточная матрица B в компактной форме; |
MLB - | заданное число нижних ко - диагоналей матрицы B (тип: целый); |
MUB - | заданное число верхних ко - диагоналей матрицы B (тип: целый); |
C - | вещественный двумерный массив размера N*M, в котором запоминается матрица произведения. |
Версии
AMM4E - | умножение с расширенной (Extended) точностью прямоугольной матрицы на ленточную матрицу, заданную в компактной форме. |
AMM4C - | умножение прямоугольной комплексной матрицы на ленточную комплексную матрицу, заданную в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
В подпрограмме АММ4E параметры А, B и C должны иметь тип Extended. В подпрограмме AMM4C параметры А, B и С должны иметь тип Complex. |
Unit tamm4r_p; interface uses SysUtils, Math, { Delphi } LStruct, Lfunc, UtRes_p, AMM4R_p; function tamm4r: String; implementation function tamm4r: String; var N,M,MLB,MUB,_i :Integer; C :Array [0..14] of Real; const A :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 ); B :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 ); begin Result := ''; N := 3; M := 5; MUB := 1; MLB := 2; AMM4R(A,N,M,B,MLB,MUB,C); Result := Result + #$0D#$0A; for _i:=0 to 14 do begin Result := Result + Format('%20.16f ',[C[_i]]); if ( ((_i+1) mod 3)=0 ) then Result := Result + #$0D#$0A; end; Result := Result + #$0D#$0A; UtRes('tamm4r',Result); { вывод результатов в файл tamm4r.res } exit; end; end. Результаты: | 3. 4. 4. 3. 2. | C = | 3. 4. 4. 3. 2. | | 3. 4. 4. 3. 2. |