Scythe-1.0.3
template<matrix_order RO, matrix_style RS, typename T , matrix_order PO, matrix_style PS>
Matrix<T, RO, RS> scythe::crossprod ( const Matrix< T, PO, PS > &  A)

Fast caclulation of $A'A$.

This function calculates $A'A$ efficiently, traversing the matrices in storage order where possible, and avoiding the use of the temporary matrix objects.

Scythe will use LAPACK/BLAS routines to compute the cross-product of column-major matrices of double-precision floating point numbers if LAPACK/BLAS is available and you compile your program with the SCYTHE_LAPACK flag enabled.

Parameters:
AThe Matrix to return the cross product of.
See also:
diag (const Matrix<T, PO, PS>& M)

References scythe::Matrix_base< ORDER, STYLE >::cols(), and scythe::Matrix_base< ORDER, STYLE >::rows().