|
Текст подпрограммы и версий amazr_c.zip , amazd_c.zip |
Тексты тестовых примеров tamazr_c.zip , tamazd_c.zip |
Сложение квадратной матрицы с симметричной матрицей, заданной в компактной форме.
amazr_c вычисляет сумму квадратной матрицы А порядка N и симметричной матрицы B того же порядка, заданной в компактной форме. Результатом этого сложения является квадратная матрица C порядка N.
int amazr_c (real *a, integer *n, real *b, real *c)
Параметры
| a - | вещественный двумерный массив размера n на n, в котором задается матрица A; |
| n - | заданный порядок матриц A и B (тип: целый); |
| b - | вещественный вектор длины n на (n + 1)/2, в котором задается в компактной форме симметричная матрица B; |
| c - | вещественный двумерный массив размера n на n, в котором запоминается результирующая матрица сложения. |
Версии
| amazd_c - | сложение с повышенной точностью квадратной матрицы с симметричной матрицей, заданной в компактной форме. |
Вызываемые подпрограммы: нет.
Замечания по использованию
| В подпрограмме amazd_c параметры a, b, и c должны иметь тип double. |
int main(void)
{
/* Initialized data */
static float a[16] /* was [4][4] */ = { 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,1.f };
static float b[10] = { 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__[16] /* was [4][4] */;
static int n, i__;
extern int amazr_c(float *, int *, float *, float *);
n = 4;
amazr_c(a, &n, b, c__);
for (i__ = 0; i__ <= 14; i__+= 2) {
printf("\n %24.12f %24.12f \n", c__[i__], c__[i__+1]);
}
return 0;
} /* main */
Результаты:
| 2.0 2.0 2.0 2.0 |
| 2.0 2.0 2.0 2.0 |
c__ = | 2.0 2.0 2.0 2.0 |
| 2.0 2.0 2.0 2.0 |