|
Текст подпрограммы и версий ammcr_c.zip , ammcd_c.zip |
Тексты тестовых примеров tammcr_c.zip , tammcd_c.zip |
Умножение симметричной ленточной матрицы, заданной в компактной форме, на ленточную матрицу, заданную в компактной форме.
ammcr_c вычисляет произведение симметричной ленточной матрицы А порядка N с NCA ко - диагоналями на ленточную матрицу B порядка N с NLB нижними и NUB верхними ко - диагоналями. Матрицы А и B задаются в компактной форме представления. Результатом этого произведения является ленточная матрица C порядка N с NLC нижними и NUC верхними ко - диагоналями, где NLС = min (N - 1, NСА + NLВ) и NUС = min (N - 1, NСА + NUВ). Матрица C запоминается в компактной форме.
int ammcr_c (real *a, integer *n, integer *nca, real *b,
integer *nlb, integer *nub, real *c)
Параметры
| a - | вещественный двумерный массив размера n* (nca + 1), в котором задается симметричная ленточная матрица A в компактной форме; |
| n - | заданный порядок матриц A и B (тип: целый); |
| nca - | заданное число ко - диагоналей матрицы A (тип: целый); |
| b - | вещественный двумерный массив размера n* (nlb + nub + 1), в котором задается ленточная матрица B в компактной форме; |
|
nlb - nub | заданное число нижних и верхних ко - диагоналей матрицы B соответственно (тип: целый); |
| c - | вещественный двумерный массив размера n* (nlc + nuc + 1), в котором запоминается ленточная матрица произведения в компактной форме; здесь nlc = min (n - 1, nca + nlb) и nuc = min (n - 1, nca + nub). |
Версии
| ammcd_c - | умножение с повышенной точностью симметричной ленточной матрицы, заданной в компактной форме, на ленточную матрицу, заданную в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
| В подпрограмме ammcd_c параметры a, b и c должны иметь тип double. |
int main(void)
{
/* Initialized data */
static float a[10] /* was [5][2] */ = { 0.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 */
static float c__[30] /* was [5][6] */;
static int n;
extern int ammcr_c(float *, int *, int *, float *, int *, int *, float *);
static int nca, nlb, nub, i__;
n = 5;
nca = 1;
nlb = 2;
nub = 1;
ammcr_c(a, &n, &nca, b, &nlb, &nub, c__);
for (i__ = 0; i__ <= 27; i__+=3) {
printf("\n %20.12e %20.12e %20.12e \n",
c__[i__], c__[i__+1], c__[i__+2]);
}
return 0;
} /* main */
Результаты:
| 0. 0. 0. 2. 2. 1. |
| 0. 0. 3. 3. 2. 1. |
c__ = | 0. 2. 3. 3. 2. 1. |
| 1. 2. 3. 3. 2. 0. |
| 1. 2. 2. 2. 0. 0. |