|
Текст подпрограммы и версий ( Фортран ) ammdr.zip , ammdd.zip , ammdc.zip |
Тексты тестовых примеров ( Фортран ) tammdr.zip , tammdd.zip , tammdc.zip |
|
Текст подпрограммы и версий ( Си ) ammdr_c.zip , ammdd_c.zip , ammdc_c.zip |
Тексты тестовых примеров ( Си ) tammdr_c.zip , tammdd_c.zip , tammdc_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) ammdr_p.zip , ammde_p.zip , ammdc_p.zip |
Тексты тестовых примеров ( Паскаль ) tammdr_p.zip , tammde_p.zip , tammdc_p.zip |
Умножение ленточных матриц, заданных в компактной форме.
АММDR вычисляет произведение ленточной матрицы А порядка N с NLА нижними и NUА верхними ко - диагоналями на ленточную матрицу B порядка N с NLВ нижними и NUВ верхними ко - диагоналями. Ленточные матрицы А и B задаются в компактной форме. Результатом этого произведения является ленточная матрица C порядка N с NLС нижними и NUС верхними ко - диагоналями, где NLС = min (N - 1, NLА + NLВ) и NUС = min (N - 1, NUА + NUВ). Ленточная матрица C запоминается также в компактной форме представления.
SUBROUTINE AMMDR (A, N, NLA, NUA, B, NLB, NUB, C)
Параметры
| A - | вещественный двумерный массив размера N* (NLА + NUА + 1), в котором задается ленточная матрица А в компактной форме; |
| N - | заданный порядок матриц А и B (тип: целый); |
|
NLA - NUA | заданное число нижних и верхних ко - диагоналей матрицы А соответственно (тип: целый); |
| B - | вещественный двумерный массив размера N* (NLВ + NUВ + 1), в котором задается ленточная матрица в компактной форме; |
|
NLB - NUB | заданное число нижних и верхних ко - диагоналей матрицы B соответственно (тип: целый); |
| C - | вещественный двумерный массив размера N* (NLС + NUС + 1), в котором запоминается ленточная матрица произведения в компактной форме; здесь NLС = min (N - 1, NLА+NLВ) и NUС = min (N - 1, NUА+NUВ). |
Версии
| AMMDD - | умножение с повышенной точностью ленточных матриц, заданных в компактной форме. |
| AMMDC - | умножение ленточных комплексных матриц, заданных в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
|
В подпрограмме АММDD параметры А, B и C должны иметь тип DОUВLЕ РRЕСISIОN. В подпрограмме АММDС параметры А, B и C должны иметь тип СОМРLЕХ. |
DIMENSION A(5, 4), B(5, 4), C(5, 7)
DATA A /2*0.0, 3*1.0, 0.0, 13*1.0, 0.0/
DATA B /0.0, 13*1.0, 0.0, 3*1.0, 2*0.0/
N = 5
NLA = 2
NUA = 1
NLB = 1
NUB = 2
CALL AMMDR (A, N, NLA, NUA, B, NLB, NUB, C)
Результаты:
| 0. 0. 0. 2. 2. 2. 1. |
| 0. 0. 2. 3. 3. 2. 1. |
C = | 0. 2. 3. 4. 3. 2. 0. |
| 1. 2. 3. 4. 3. 0. 0. |
| 1. 2. 3. 3. 0. 0. 0. |