Текст подпрограммы и версий
ama2r_c.zip , ama2d_c.zip
Тексты тестовых примеров
tama2r_c.zip , tama2d_c.zip

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

Назначение

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

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

ama2r_c вычисляет сумму квадратной матрицы  А порядка  N и ленточной симметричной матрицы  B того же порядка, заданной в компактной форме. Результатом этого сложения является квадратная матрица  C порядка  N.

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

    int ama2r_c (real *a, integer *n, real *b, integer *nc,
        real *c)

Параметры

a - вещественный двумерный массив размера n на n, в котором задается матрица  A;
n - заданный порядок матриц  A и  B (тип: целый);
b - вещественный двумерный массив размера n на (nc + 1), в котором в компактной форме задается ленточная симметричная матрица  B;
nc - заданное число ко - диагоналей матрицы  B (тип: целый);
c - вещественный двумерный массив размера n на n, в котором запоминается результирующая матрица сложения.

Версии

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

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

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

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

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

int main(void)
{
    /* Initialized data */
    static float a[16] /* was [4][4] */ = { 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,1.f };
    static float b[8] /* was [4][2] */ = { 0.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f };

    /* Local variables */
    extern int ama2r_c(float *, int *, float *, int *, float *);
    static float c__[16] /* was [4][4] */;
    static int n, nc, i__;

    n = 4;
    nc = 1;
    ama2r_c(a, &n, b, &nc, c__);

    for (i__ = 0; i__ <= 14; i__+= 2) {
         printf("\n  %16.7f %16.7f \n", c__[i__], c__[i__+1]);
    }
    return 0;
} /* main */


Результаты:

                      | 2.0   2.0   1.0   1.0 |
                      | 2.0   2.0   2.0   1.0 |
         c__  =   | 1.0   2.0   2.0   2.0 |
                      | 1.0   1.0   2.0   2.0 |