|
Текст подпрограммы и версий amm2r_c.zip , amm2d_c.zip |
Тексты тестовых примеров tamm2r_c.zip , tamm2d_c.zip |
Умножение прямоугольной матрицы на симметричную матрицу, заданную в компактной форме.
amm2r_c вычисляет произведение прямоугольной матрицы А размера N на М на симметричную матрицу B порядка М, заданную в компактной форме. Результатом этого произведения является матрица C размера N на М.
int amm2r_c (real *a, integer *n, integer *m, real *b,
real *c)
Параметры
| a - | вещественный двумерный массив размера n*m, в котором задается матрица A; |
| n - | заданное число строк матрицы A (тип: целый); |
| m - | заданное число столбцов матрицы A (тип: целый); |
| b - | вещественный вектор длины m * (m + 1) /2, в котором задается симметричная матрица B в компактной форме; |
| c - | вещественный двумерный массив размера n*m, в котором запоминается матрица произведения. |
Версии
| amm2d_c - | умножение с повышенной точностью прямоугольной матрицы на симметричную, заданную в компактной форме. |
Вызываемые подпрограммы: нет
Замечания по использованию
| В подпрограмме amm2d_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[10] = { 1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f,1.f };
/* Local variables */
extern int amm2r_c(float *, int *, int *, float *, float *);
static float c__[12] /* was [3][4] */;
static int m, n, i__;
n = 3;
m = 4;
amm2r_c(a, &n, &m, b, 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 */
Результаты:
| 4. 4. 4. 4. |
c__ = | 4. 4. 4. 4. |
| 4. 4. 4. 4. |