|
Текст подпрограммы и версий ammbr_c.zip , ammbd_c.zip |
Тексты тестовых примеров tammbr_c.zip , tammbd_c.zip |
Умножение симметричной ленточной матрицы, заданной в компактной форме, на симметричную матрицу, заданную в компактной форме.
ammbr_c вычисляет произведение симметричной ленточной матрицы А порядка N с NС ко - диагоналями на симметричную матрицу B порядка N. Матрицы А и B задаются в компактной форме представления. Результатом этого произведения является квадратная матрица C порядка N.
int ammbr_c (real *a, integer *n, integer *nc, real *b,
real *c)
Параметры
| a - | вещественный двумерный массив размера n* (nc + 1), в котором задается симметричная ленточная матрица A в компактной форме; |
| n - | заданный порядок матриц A и B (тип: целый); |
| nc - | заданное число ко - диагоналей матрицы A (тип: целый); |
| b - | вещественный вектор длины n* (n + 1)/2, в котором задается симметричная матрица B в компактной форме; |
| c - | вещественный двумерный массив размера n*n, в котором запоминается матрица произведения. |
Версии
| ammbd_c - | умножение с повышенной точностью симметричной ленточной матрицы, заданной в компактной форме, на симметричную матрицу, заданную в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
| В подпрограмме ammbd_c параметры a, b и c должны иметь тип double. |
int main(void)
{
/* Initialized data */
static float a[10] /* was [5][2] */ = { 0.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,
1.f,1.f };
static float b[15] = { 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__[25] /* was [5][5] */;
static int n;
extern int ammbr_c(float *, int *, int *, float *, float *);
static int nc, i__;
n = 5;
nc = 1;
ammbr_c(a, &n, &nc, b, c__);
for (i__ = 0; i__ <= 21; i__+=3) {
printf("\n %20.12e %20.12e %20.12e \n",
c__[i__], c__[i__+1], c__[i__+2]);
}
printf("\n %20.12e \n", c__[24]);
return 0;
} /* main */
Результаты:
| 2 2 2 2 2 |
| 3 3 3 3 3 |
c__ = | 3 3 3 3 3 |
| 3 3 3 3 3 |
| 2 2 2 2 2 |