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 of the original kernel
- updated QuadraticTimeMMD with newly added CKernel::sum methods.