|
Текст подпрограммы и версий ammfr_c.zip , ammfd_c.zip |
Тексты тестовых примеров tammfr_c.zip , tammfd_c.zip |
Умножение ленточной матрицы, заданной в компактной форме, на симметричную матрицу, заданную в компактной форме.
ammfr_c вычисляет произведение ленточной матрицы А порядка N с NLA нижними и NUA верхними ко - диагоналями на симметричную матрицу B порядка N. Матрицы А и B задаются в компактной форме представления. Результатом этого произведения является квадратная матрица С.
int ammfr_c (real *a, integer *n, integer *nla, integer *nua,
real *b, real *c)
Параметры
| a - | вещественный двумерный массив размера n * (nla+nua+1), в котором задается ленточная матрица в компактной форме; |
| n - | заданный порядок матриц A и B (тип: целый); |
|
nla - nua | заданное число нижних и верхних ко - диагоналей матрицы A соответственно (тип: целый); |
| b - | вещественный вектор длины n * (n+1)/2, в котором задается симметричная матрица B в компактной форме; |
| c - | вещественный двумерный массив размера n*n, в котором запоминается матрица произведения. |
Версии
| ammfd_c - | умножение с повышенной точностью ленточной матрицы, заданной в компактной форме, на симметричную матрицу, заданную в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
| В подпрограмме ammfd_c параметры a, b и c должны иметь тип double. |
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] = { 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 ammfr_c(float *, int *, int *, int *, float *, float *);
static int nla, nua, i__;
n = 5;
nla = 2;
nua = 1;
ammfr_c(a, &n, &nla, &nua, b, c__);
for (i__ = 1; i__ <= 25; ++i__) {
printf("\n %16.7e \n", c__[i__-1]);
}
return 0;
} /* main */
Результаты:
| 2. 2. 2. 2. 2. |
| 3. 3. 3. 3. 3. |
c__ = | 4. 4. 4. 4. 4. |
| 4. 4. 4. 4. 4. |
| 3. 3. 3. 3. 3. |