/**
  * <p>
  * Performs the following operation:<br>
  * <br>
  * c = c + a<sup>H</sup> * b<sup>H</sup><br>
  * c<sub>ij</sub> = c<sub>ij</sub> + &sum;<sub>k=1:n</sub> { a<sub>ki</sub> * b<sub>jk</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 multAddTransAB(ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
 {
     if (a.numCols >= EjmlParameters.CMULT_TRANAB_COLUMN_SWITCH)
     {
         MatrixMatrixMult_ZDRM.multAddTransAB_aux(a, b, c, null);
     }
     else
     {
         MatrixMatrixMult_ZDRM.multAddTransAB(a, b, c);
     }
 }
 /**
  * <p>
  * Performs the following operation:<br>
  * <br>
  * c = c + &alpha; * a<sup>H</sup> * b<sup>H</sup><br>
  * c<sub>ij</sub> = c<sub>ij</sub> + &alpha; * &sum;<sub>k=1:n</sub> { a<sub>ki</sub> * b<sub>jk</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 multAddTransAB(double realAlpha, double imagAlpha, ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
 {
     // TODO add a matrix vectory multiply here
     if (a.numCols >= EjmlParameters.CMULT_TRANAB_COLUMN_SWITCH)
     {
         MatrixMatrixMult_ZDRM.multAddTransAB_aux(realAlpha, imagAlpha, a, b, c, null);
     }
     else
     {
         MatrixMatrixMult_ZDRM.multAddTransAB(realAlpha, imagAlpha, a, b, c);
     }
 }