This function computes the expression y = a*x1*x2+y. The column
dimension of x1 must match the row dimension of x2.
The values stored in the matrices are assumed to be double-precision
floating-point values. This code is suitable for
dense matrices. No optimizations are made for sparse matrices.
@param x1 Pointer to r1 by c1 input matrix.
@param a Scalar value.
@param r1 Number of rows in x1.
@param c1 Number of columns in x1. Also number of rows in x2.
@param x2 Pointer to c1 by c2 input matrix.
@param c2 Number of columns in x2.
@param y Pointer to r1 by c2 output matrix.
- Algorithm:
- DSPF_dp_mat_mul_gemm.c is the natural C equivalent of the optimized intrinsic C code withoutrestrictions. Note that the intrinsic C code is optimized and restrictions may apply.
- Assumptions:
- The arrays x1, x2 and y are stored in distinct arrays. That is, in-place processing is not allowed.
r1, c2 are assumed to be even. c1 is assumed to be an integer multiple of 4.
- Implementation Notes:
- Interruptibility : The code is interruptible.
Endian support : supports little endian mode.