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 single-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_sp_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. In-place processing is not allowed.
r1, c2 are assumed to be an integer multiple of 4. c1 is assumed to be even.
- Implementation Notes:
- Interruptibility : The code is interruptible.
Endian support : supports Little endian mode.