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