|
Текст подпрограммы и версий ama6r_c.zip , ama6d_c.zip |
Тексты тестовых примеров tama6r_c.zip , tama6d_c.zip |
Сложение симметричной матрицы и ленточной симметричной, заданных в компактной форме.
ama6r_c вычисляет сумму симметричной матрицы А порядка N, заданной в компактной форме, и ленточной симметричной матрицы B того же порядка, заданной также в компактной форме. Результатом этого сложения является результирующая симметричная матрица C порядка N. Матрица C запоминается в компактной форме.
int ama6r_c (real *a, integer *n, real *b, integer *nc,
real *c)
Параметры
| a - | вещественный вектор длины n на (n + 1)/2, в котором задается в компактной форме симметричная матрица A; |
| n - | заданный порядок матриц A и B (тип: целый); |
| b - | вещественный двумерный массив размера n на (nc + 1), в котором в компактной форме задается ленточная симметричная матрица B; |
| nc - | заданное число ко - диагоналей матрицы B (тип: целый); |
| c - | вещественный вектор длины n на (n + 1)/2, в котором в компактной форме запоминается результирующая ленточная матрица сложения. |
Версии
| ama6d_c - | сложение с повышенной точностью симметричной матрицы и ленточной симметричной, заданных в компактной форме. |
Вызываемые подпрограммы: нет.
Замечания по использованию
| В подпрограмме ama6d_c параметры a, b и c должны иметь тип double. |
int main(void)
{
/* Initialized data */
static float a[10] = { 1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f };
static float b[8] /* was [4][2] */ = { 0.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f };
/* Local variables */
extern int ama6r_c(float *, int *, float *, int *, float *);
static float c__[10];
static int n, nc, i__;
n = 4;
nc = 1;
ama6r_c(a, &n, b, &nc, c__);
for (i__ = 1; i__ <= 10; ++i__) {
printf("\n %15.7f \n", c__[i__-1]);
}
return 0;
} /* main */
Результаты:
| 2.0 |
| 2.0 |
| 2.0 |
| 1.0 |
c__ = | 2.0 |
| 2.0 |
| 1.0 |
| 1.0 |
| 2.0 |
| 2.0 |