|
Текст подпрограммы и версий amm3r_c.zip , amm3d_c.zip |
Тексты тестовых примеров tamm3r_c.zip , tamm3d_c.zip |
Умножение прямоугольной матрицы на симметричную ленточную матрицу, заданную в компактной форме.
amm3r_c вычисляет произведение прямоугольной матрицы А размера N на М на симметричную ленточную матрицу B порядка М с NС ко - диагоналями, заданную в компактной форме. Результатом этого произведения является прямоугольная матрица C размера N на М.
int amm3r_c (real *a, integer *n, integer *m, real *b,
integer *nc, real *c)
Параметры
| a - | вещественный двумерный массив размера n*m, в котором задается матрица A; |
| n - | заданное число строк матрицы A (тип: целый); |
| m - | заданное число столбцов матрицы A (тип: целый); |
| b - | вещественный двумерный массив размера m * (nc + 1), в котором задается симметричная ленточная матрица B в компактной форме; |
| nc - | заданное число ко - диагоналей матрицы B (тип: целый); |
| c - | вещественный двумерный массив размера n*m, в котором запоминается матрица произведения. |
Версии
| amm3d_c - | умножение с повышенной точностью прямоугольной матрицы на симметричную ленточную матрицу, заданную в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
| В подпрограмме amm3d_c параметры a, b и c должны иметь тип double. |
int main(void)
{
/* Initialized data */
static float a[12] /* was [3][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 };
static float b[8] /* was [4][2] */ = { 0.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f };
/* Local variables */
extern int amm3r_c(float *, int *, int *, float *, int *, float *);
static float c__[12] /* was [3][4] */;
static int m, n, nc, i__;
n = 3;
m = 4;
nc = 1;
amm3r_c(a, &n, &m, b, &nc, c__);
for (i__ = 0; i__ <= 9; i__+= 3) {
printf("\n %20.12e %20.12e %20.12e \n",
c__[i__], c__[i__+1], c__[i__+2]);
}
return 0;
} /* main */
Результаты:
| 2. 3. 3. 2. |
c__ = | 2. 3. 3. 2. |
| 2. 3. 3. 2. |