Текст подпрограммы и версий ( Фортран )
ammfr.zip , ammfd.zip 		Тексты тестовых примеров ( Фортран )
tammfr.zip , tammfd.zip
Текст подпрограммы и версий ( Си )
ammfr_c.zip , ammfd_c.zip 		Тексты тестовых примеров ( Си )
tammfr_c.zip , tammfd_c.zip
Текст подпрограммы и версий ( Паскаль )
ammfr_p.zip , ammfe_p.zip 		Тексты тестовых примеров ( Паскаль )
tammfr_p.zip , tammfe_p.zip

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

Назначение

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

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

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

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

    SUBROUTINE  AMMFR (A, N, NLA, NUA, B, C) 
   Параметры
A - вещественный двумерный массив размера N * (NLА+NUА+1), в котором задается ленточная матрица в компактной форме;
N - заданный порядок матриц А и B (тип: целый);
   NLA -
   NUA   		заданное число нижних и верхних ко - диагоналей матрицы А соответственно (тип: целый);
B - вещественный вектор длины N * (N+1)/2, в котором задается симметричная матрица B в компактной форме;
C - вещественный двумерный массив размера N*N, в котором запоминается матрица произведения. 
   Версии
AMMFD - умножение с повышенной точностью ленточной матрицы, заданной в компактной форме, на симметричную матрицу, заданную в компактной форме. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
  В подпрограмме АММFD параметры А, B и C должны иметь тип DОUВLЕ РRЕСISIОN.

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

          DIMENSION  A(5, 4), B(15), C(5, 5)
          DATA  A /2*0.0, 3*1.0, 0.0, 13*1.0, 0.0/
          DATA  B /15*1.0/
          N = 5
          NLA = 2
          NUA = 1
          CALL  AMMFR (A, N, NLA, NUA, B, C)

Результаты:

               |  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. |