/**
  * <p>
  * Performs the following operation:<br>
  * <br>
  * c = c + a<sup>H</sup> * b<br>
  * c<sub>ij</sub> = c<sub>ij</sub> + &sum;<sub>k=1:n</sub> { a<sub>ki</sub> * b<sub>kj</sub>}
  * </p>
  *
  * @param a The left matrix in the multiplication operation. Not modified.
  * @param b The right matrix in the multiplication operation. Not modified.
  * @param c Where the results of the operation are stored. Modified.
  */
 public static void multAddTransA(ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
 {
     if (a.numCols >= EjmlParameters.CMULT_COLUMN_SWITCH ||
         b.numCols >= EjmlParameters.CMULT_COLUMN_SWITCH)
     {
         MatrixMatrixMult_ZDRM.multAddTransA_reorder(a, b, c);
     }
     else
     {
         MatrixMatrixMult_ZDRM.multAddTransA_small(a, b, c);
     }
 }
 /**
  * <p>
  * Performs the following operation:<br>
  * <br>
  * c = c + &alpha; * a<sup>H</sup> * b<br>
  * c<sub>ij</sub> =c<sub>ij</sub> +  &alpha; * &sum;<sub>k=1:n</sub> { a<sub>ki</sub> * b<sub>kj</sub>}
  * </p>
  *
  * @param realAlpha Real component of scaling factor.
  * @param imagAlpha Imaginary component of scaling factor.
  * @param a The left matrix in the multiplication operation. Not modified.
  * @param b The right matrix in the multiplication operation. Not modified.
  * @param c Where the results of the operation are stored. Modified.
  */
 public static void multAddTransA(double realAlpha, double imagAlpha, ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
 {
     // TODO add a matrix vectory multiply here
     if (a.numCols >= EjmlParameters.CMULT_COLUMN_SWITCH ||
         b.numCols >= EjmlParameters.CMULT_COLUMN_SWITCH)
     {
         MatrixMatrixMult_ZDRM.multAddTransA_reorder(realAlpha, imagAlpha, a, b, c);
     }
     else
     {
         MatrixMatrixMult_ZDRM.multAddTransA_small(realAlpha, imagAlpha, a, b, c);
     }
 }