/**
* <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> + ∑<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 + α * a<sup>H</sup> * b<sup>H</sup><br>
* c<sub>ij</sub> = c<sub>ij</sub> + α * ∑<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);
}
}