|
Текст подпрограммы и версий ama4r_c.zip , ama4d_c.zip |
Тексты тестовых примеров tama4r_c.zip , tama4d_c.zip |
Сложение ленточной и ленточной симметричной матриц, заданных в компактной форме.
ama4r_c вычисляет сумму ленточной матрицы А порядка N, заданной в компактной форме, и ленточной симметричной матрицы B того же порядка, заданной также в компактной форме. Результатом этого сложения является ленточная матрица C порядка N с max (NLА, NС) нижними и max (NUА, NС) верхними ко - диагоналями. Матрица C запоминается в компактной форме.
int ama4r_c (real *a, integer *n, integer *nla, integer *nua,
real *b, integer *nc, real *c)
Параметры
| a - | вещественный двумерный массив размера n на (nla + nua + 1), в котором в компактной форме задается ленточная матрица A; |
| n - | заданный порядок матриц A и B (тип: целый); |
|
nla - nua | заданное число нижних и верхних ко - диагоналей матрицы A соответственно (тип: целый); |
| b - | вещественный двумерный массив размера n на (nc + 1), в котором в компактной форме задается ленточная симметричная матрица B; |
| nc - | заданное число ко - диагоналей матрицы B (тип: целый); |
| c - | вещественный двумерный массив размера n на (max (nla, nc) + max (nua, nc) + 1), в котором в компактной форме запоминается ленточная матрица сложения. |
Версии
| ama4d_c - | сложение с повышенной точностью ленточной матрицы и ленточной симметричной матрицы, заданных в компактной форме. |
Вызываемые подпрограммы: нет.
Замечания по использованию
| В подпрограмме ama4d_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[10] /* was [5][2] */ = { 0.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,
1.f,1.f };
/* Local variables */
extern int ama4r_c(float *, int *, int *, int *, float *, int *, float *);
static float c__[20] /* was [5][4] */;
static int n, nc, nla, nua, i__;
n = 5;
nla = 2;
nua = 1;
nc = 1;
ama4r_c(a, &n, &nla, &nua, b, &nc, c__);
for (i__ = 0; i__ <= 18; i__+= 2) {
printf("\n %15.7f %15.7f \n", c__[i__], c__[i__+1]);
}
return 0;
} /* main */
Результаты:
| 0.0 0.0 2.0 2.0 |
| 0.0 2.0 2.0 2.0 |
c__ = | 1.0 2.0 2.0 2.0 |
| 1.0 2.0 2.0 2.0 |
| 1.0 2.0 2.0 0.0 |