/**
  * <p>
  * Performs the following operation:<br>
  * <br>
  * c = a<sup>T</sup> * b<sup>T</sup><br>
  * 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 multTransAB(ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
 {
     if (a.numCols >= EjmlParameters.CMULT_TRANAB_COLUMN_SWITCH)
     {
         MatrixMatrixMult_ZDRM.multTransAB_aux(a, b, c, null);
     }
     else
     {
         MatrixMatrixMult_ZDRM.multTransAB(a, b, c);
     }
 }
 /**
  * <p>
  * Performs the following operation:<br>
  * <br>
  * c = &alpha; * a<sup>H</sup> * b<sup>H</sup><br>
  * 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 multTransAB(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.multTransAB_aux(realAlpha, imagAlpha, a, b, c, null);
     }
     else
     {
         MatrixMatrixMult_ZDRM.multTransAB(realAlpha, imagAlpha, a, b, c);
     }
 }