|
Текст подпрограммы и версий ama5r_c.zip , ama5d_c.zip |
Тексты тестовых примеров tama5r_c.zip , tama5d_c.zip |
Сложение симметричной матрицы и ленточной, заданных в компактной форме.
ama5r_c вычисляет сумму симметричной матрицы А порядка N, заданной в компактной форме, и ленточной матрицы B того же порядка, заданной также в компактной форме. Результатом этого сложения является квадратная матрица C порядка N.
int ama5r_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), в котором в компактной форме задается ленточная матрица B; |
|
nlb - nub | заданное число нижних и верхних ко - диагоналей матрицы B соответственно (тип: целый); |
| c - | вещественный двумерный массив размера n на n, в котором запоминается результирующая матрица сложения. |
Версии
| ama5d_c - | сложение с повышенной точностью симметричной матрицы и ленточной, заданных в компактной форме. |
Вызываемые подпрограммы: нет.
Замечания по использованию
| В подпрограмме ama5d_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 ama5r_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;
ama5r_c(a, &n, b, &nlb, &nub, c__);
for (i__ = 0; i__ <= 22; i__+= 2) {
printf("\n %15.7f %15.7f \n", c__[i__], c__[i__+1]);
}
printf("\n %15.7f \n", c__[24]);
return 0;
} /* main */
Результаты:
| 2.0 2.0 1.0 1.0 1.0 |
| 2.0 2.0 2.0 1.0 1.0 |
c__ = | 2.0 2.0 2.0 2.0 1.0 |
| 1.0 2.0 2.0 2.0 2.0 |
| 1.0 1.0 2.0 2.0 2.0 |