/**
* <p>
* Performs the following operation:<br>
* <br>
* c = c + α * a * b<br>
* c<sub>ij</sub> = c<sub>ij</sub> + α * ∑<sub>k=1:n</sub> { a<sub>ik</sub> * b<sub>kj</sub>}
* </p>
*
* @param realAlpha real component of scaling factor.
* @param imgAlpha 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 multAdd(double realAlpha, double imgAlpha, ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
{
if (b.numCols >= EjmlParameters.CMULT_COLUMN_SWITCH)
{
MatrixMatrixMult_ZDRM.multAdd_reorder(realAlpha, imgAlpha, a, b, c);
}
else
{
MatrixMatrixMult_ZDRM.multAdd_small(realAlpha, imgAlpha, 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 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 * b<br>
* c<sub>ij</sub> = c<sub>ij</sub> + ∑<sub>k=1:n</sub> { a<sub>ik</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 multAdd(ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
{
if (b.numCols >= EjmlParameters.MULT_COLUMN_SWITCH)
{
MatrixMatrixMult_ZDRM.multAdd_reorder(a, b, c);
}
else
{
MatrixMatrixMult_ZDRM.multAdd_small(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);
}
}
/**
* <p>Performs the following operation:<br>
* <br>
* c = a<sup>H</sup> * b <br>
* <br>
* c<sub>ij</sub> = ∑<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 multTransA(ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
{
if (a.numCols >= EjmlParameters.CMULT_COLUMN_SWITCH ||
b.numCols >= EjmlParameters.CMULT_COLUMN_SWITCH)
{
MatrixMatrixMult_ZDRM.multTransA_reorder(a, b, c);
}
else
{
MatrixMatrixMult_ZDRM.multTransA_small(a, b, c);
}
}
/**
* <p>
* Performs the following operation:<br>
* <br>
* c = c + α * a<sup>H</sup> * b<br>
* c<sub>ij</sub> =c<sub>ij</sub> + α * ∑<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);
}
}
/**
* <p>
* Performs the following operation:<br>
* <br>
* c = c + α * a * b<sup>H</sup><br>
* c<sub>ij</sub> = c<sub>ij</sub> + α * ∑<sub>k=1:n</sub> { a<sub>ik</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 multAddTransB(double realAlpha, double imagAlpha, ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
{
// TODO add a matrix vectory multiply here
MatrixMatrixMult_ZDRM.multAddTransB(realAlpha, imagAlpha, a, b, c);
}
/**
* <p>
* Performs the following operation:<br>
* <br>
* c = c + a * b<sup>H</sup> <br>
* c<sub>ij</sub> = c<sub>ij</sub> + ∑<sub>k=1:n</sub> { a<sub>ik</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 multAddTransB(ZMatrixRMaj a, ZMatrixRMaj b, ZMatrixRMaj c)
{
MatrixMatrixMult_ZDRM.multAddTransB(a, b, c);
}
