Added proper symmetry check in CustomKernel and updated QuadraticTimeMMD

- added method `bool SGMatrix<T>::is_symmetric()` and unit-tests
- added an extra parameter in CCustomKernel, m_is_symmetric, to preserve the lhs_equals_rhs flag of the original kernel
- updated QuadraticTimeMMD with newly added CKernel::sum methods.