Текст подпрограммы и версий ( Фортран )
amm8r.zip , amm8d.zip
Тексты тестовых примеров ( Фортран )
tamm8r.zip , tamm8d.zip
Текст подпрограммы и версий ( Си )
amm8r_c.zip , amm8d_c.zip
Тексты тестовых примеров ( Си )
tamm8r_c.zip , tamm8d_c.zip
Текст подпрограммы и версий ( Паскаль )
amm8r_p.zip , amm8e_p.zip
Тексты тестовых примеров ( Паскаль )
tamm8r_p.zip , tamm8e_p.zip

Подпрограмма:  AMM8R

Назначение

Умножение симметричной матрицы, заданной в компактной форме, на ленточную матрицу, заданную в компактной форме.

Математическое описание

АММ8R вычисляет произведение симметричной матрицы А порядка N на ленточную матрицу B порядка N с NLВ нижними и NUВ верхними ко - диагоналями. Матрицы А и B задаются в компактной форме представления. Результатом этого произведения является квадратная матрица C порядка N.

Использование

    SUBROUTINE  AMM8R (A, N, B, NLB, NUB, C) 

Параметры

A - вещественный вектор длины N * (N + 1)/2, в котором задается симметричная матрица А в компактной форме;
N - заданный порядок матриц А и B (тип: целый);
B - вещественный двумерный массив размера N * (NLВ + NUВ + 1), в котором задается ленточная матрица в компактной форме;
NL, NUB - заданное число нижних и верхних ко - диагоналей матрицы B соответственно (тип: целый);
C - вещественный двумерный массив размера N*N, в котором запоминается матрица произведения.

Версии

AMM8D - умножение с повышенной точностью симметричной матрицы, заданной в компактной форме, на ленточную матрицу, заданную в компактной форме.

Вызываемые подпрограммы: нет

Замечания по использованию

  В подпрограмме АММ8D параметры А, B и C должны иметь тип DОUВLЕ РRЕСISIОN.

Пример использования

          DIMENSION  A(15), B(5, 4), C(5, 5)
          DATA  A /15*1.0/
          DATA  B /2*0.0, 3*1.0, 0.0, 13*1.0, 0.0/
          N = 5
          NLB = 2
          NUB = 1
          CALL  AMM8R (A, N, B, NLB, NUB, C)

Результаты:

                     |  3.  4.  4.  3.  2. |
                     |  3.  4.  4.  3.  2. |
           C  =   |  3.  4.  4.  3.  2. |
                     |  3.  4.  4.  3.  2. |
                     |  3.  4.  4.  3.  2. |