|
Текст подпрограммы и версий ammdr_c.zip , ammdd_c.zip , ammdc_c.zip |
Тексты тестовых примеров tammdr_c.zip , tammdd_c.zip , tammdc_c.zip |
Умножение ленточных матриц, заданных в компактной форме.
ammdr_c вычисляет произведение ленточной матрицы А порядка N с NLA нижними и NUA верхними ко - диагоналями на ленточную матрицу B порядка N с NLB нижними и NUB верхними ко - диагоналями. Ленточные матрицы А и B задаются в компактной форме. Результатом этого произведения является ленточная матрица C порядка N с NLC нижними и NUC верхними ко - диагоналями, где NLС = min (N - 1, NLА + NLВ) и NUС = min (N - 1, NUА + NUВ). Ленточная матрица C запоминается также в компактной форме представления.
int ammdr_c (real *a, integer *n, integer *nla, integer *nua,
real *b, integer *nlb, integer *nub, real *c)
Параметры
| a - | вещественный двумерный массив размера n* (nla + nua + 1), в котором задается ленточная матрица A в компактной форме; |
| n - | заданный порядок матриц A и B (тип: целый); |
|
nla - nua | заданное число нижних и верхних ко - диагоналей матрицы A соответственно (тип: целый); |
| b - | вещественный двумерный массив размера n* (nlb + nub + 1), в котором задается ленточная матрица в компактной форме; |
|
nlb - nub | заданное число нижних и верхних ко - диагоналей матрицы B соответственно (тип: целый); |
| c - | вещественный двумерный массив размера n* (nlc + nuc + 1), в котором запоминается ленточная матрица произведения в компактной форме; здесь nlc = min (n - 1, nla+nlb) и nuc = min (n - 1, nua+nub). |
Версии
| ammdd_c - | умножение с повышенной точностью ленточных матриц, заданных в компактной форме. |
| ammdc_c - | умножение ленточных комплексных матриц, заданных в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
|
В подпрограмме ammdd_c параметры a, b и c должны иметь тип double. В подпрограмме ammdc_c параметры a, b и c должны иметь тип complex. |
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[20] /* was [5][4] */ = { 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,1.f,1.f,1.f,0.f,0.f };
/* Local variables */
static float c__[35] /* was [5][7] */;
static int n;
extern int ammdr_c(float *, int *, int *, int *, float *, int *,
int *, float *);
static int nla, nlb, nua, nub, i__;
n = 5;
nla = 2;
nua = 1;
nlb = 1;
nub = 2;
ammdr_c(a, &n, &nla, &nua, b, &nlb, &nub, c__);
for (i__ = 0; i__ <= 30; i__+=5) {
printf("\n %15.7e %15.7e %15.7e %15.7e %15.7e \n",
c__[i__], c__[i__+1], c__[i__+2], c__[i__+3], c__[i__+4]);
}
return 0;
} /* main */
Результаты:
| 0. 0. 0. 2. 2. 2. 1. |
| 0. 0. 2. 3. 3. 2. 1. |
c__ = | 0. 2. 3. 4. 3. 2. 0. |
| 1. 2. 3. 4. 3. 0. 0. |
| 1. 2. 3. 3. 0. 0. 0. |