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