Текст подпрограммы и версий
ammer_c.zip , ammed_c.zip , ammec_c.zip 		Тексты тестовых примеров
tammer_c.zip , tammed_c.zip , tammec_c.zip

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

Назначение

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

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

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

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

    int ammer_c (real *a, integer *n, integer *nla, integer *nua,
            real *b, integer *m, real *c)

   Параметры
a - вещественный двумерный массив размера n* (nla+nua+1), в котором задается ленточная матрица A в компактной форме;
n - заданный порядок матрицы A (тип: целый);
   nla -
   nua   		заданное число нижних и верхних ко - диагоналей матрицы A соответственно (тип: целый);
b - вещественный двумерный массив размера n*m, в котором задается матрица B;
m - заданное число столбцов матрицы B (тип: целый);
c - вещественный двумерный массив размера n*m, в котором запоминается матрица произведения. 
   Версии
ammed_c - умножение с повышенной точностью ленточной матрицы, заданной в компактной форме, на прямоугольную матрицу.
ammec_c - умножение ленточной комплексной матрицы, заданной в компактной форме, на прямоугольную комплексную матрицу. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
 

В подпрограмме ammed_c параметры a, b и c должны иметь тип double.

В подпрограмме ammec_c параметры a, b и c должны иметь тип complex.

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

int main(void)
{
    /* Initialized data */
    static float a[20] /* was [5][4] */ = { 0.f,0.f,1.f,1.f,1.f,0.f,1.f,1.f,
                            1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,0.f };
    static float b[15] /* was [5][3] */ = { 1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,
                                            1.f,1.f,1.f,1.f,1.f,1.f,1.f };
    /* Local variables */
    static float c__[15] /* was [5][3] */;
    static int m, n;
    extern int ammer_c(float *, int *, int *, int *, float *, int *, float *);
    static int nla, nua, i__;

    n = 5;
    nla = 2;
    nua = 1;
    m = 3;
    ammer_c(a, &n, &nla, &nua, b, &m, c__);

    for (i__ = 0; i__ <= 10; i__+=5) {
         printf("\n %15.7e %15.7e %15.7e %15.7e %15.7e \n",
                c__[i__], c__[i__+1], c__[i__+2], c__[i__+3], c__[i__+4]);
    }
    return 0;
} /* main */

Результаты:

                   |  2.  2.  2. |
                   |  3.  3.  3. |
     c__  =    |  4.  4.  4. |
                   |  4.  4.  4. |
                   |  3.  3.  3. |