/* Copyright 1994-2004 The MathWorks, Inc.
 *
 * File: rt_matmultandinccc_sgl.c     $Revision: 1.3.4.1 $
 *
 * Abstract:
 *      Real-Time Workshop support routine for matrix multiplication
 *      and increment for two complex single precision float operands
 *
 */

#include "rtlibsrc.h"

/*
 * Function: rt_MatMultAndIncCC_Sgl
 * Abstract:
 *      2-input matrix multiply and increment function
 *      Input 1: Complex, single-precision
 *      Input 2: Complex, single-precision
 */
#ifdef CREAL_T
void rt_MatMultAndIncCC_Sgl(creal32_T       *y,
                            const creal32_T *A,
                            const creal32_T *B,
                            const int       dims[3])
{
  int k;
  for(k=dims[2]; k-- > 0; ) {
    const creal32_T *A1 = A;
    int i;
    for(i=dims[0]; i-- > 0; ) {
      const creal32_T *A2 = A1++;
      const creal32_T *B1 = B;
      creal32_T acc;
      int j;
      acc.re = (real32_T)0.0;
      acc.im = (real32_T)0.0;
      for(j=dims[1]; j-- > 0; ) {

        creal32_T A2_val = *A2;
        creal32_T B1_val = *B1++;
        acc.re += CMULT_RE(A2_val, B1_val);
        acc.im += CMULT_IM(A2_val, B1_val);
        A2 += dims[0];
      }
      y->re += acc.re;
      y->im += acc.im;
      y++;
    }
    B += dims[1];
  }
}
#endif
/* [EOF] rt_matmultandinccc_sgl.c */