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