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