|
Текст подпрограммы и версий amm9r_c.zip , amm9d_c.zip |
Тексты тестовых примеров tamm9r_c.zip , tamm9d_c.zip |
Умножение симметричных ленточных матриц, заданных в компактной форме.
amm9r_c вычисляет произведение симметричной ленточной матрицы А порядка N с NCA ко - диагоналями на симметричную ленточную матрицу B порядка N с NCB ко - диагоналями. Матрицы А и B задаются в компактной форме представления. Результатом этого произведения является ленточная матрица C порядка N с М верхними и М нижними ко - диагоналями, где М = min (N - 1, NСА + NСВ). Матрица C запоминается в компактной форме.
int amm9r_c (real *a, integer *n, integer *nca, real *b,
integer *ncb, real *c)
Параметры
| a - | вещественный двумерный массив размера n * (nca + 1), в котором задается симметричная ленточная матрица A в компактной форме; |
| n - | заданный порядок матриц A и B (тип: целый); |
| nca - | заданное число ко - диагоналей матрицы A (тип: целый); |
| b - | вещественный двумерный массив размера n * (ncb + 1), в котором задается симметричная ленточная матрица B в компактной форме; |
| ncb - | заданное число ко - диагоналей матрицы B (тип: целый); |
| c - | вещественный двумерный массив размера n * (min (2n - 1, 2 (nca + ncb) + 1)), в котором запоминается ленточная матрица произведения в компактной форме. |
Версии
| amm9d_c - | умножение с повышенной точностью симметричных ленточных матриц, заданных в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
| В подпрограмме amm9d_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[10] /* was [5][2] */ = { 0.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,
1.f,1.f };
/* Local variables */
extern int amm9r_c(float *, int *, int *, float *, int *, float *);
static float c__[25] /* was [5][5] */;
static int n, nca, ncb, i__;
n = 5;
nca = 1;
ncb = 1;
amm9r_c(a, &n, &nca, b, &ncb, 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 */
Результаты:
| 0 0 2 2 1 |
| 0 2 3 2 1 |
c__ = | 1 2 3 2 1 |
| 1 2 3 2 0 |
| 1 2 2 0 0 |