Implement symmetric_product() to produce a symmetric matrix: C = alpha * X @ X.T + beta * C

Implement symmetric_product() to produce a symmetric matrix: C = alpha * X @ X.T + beta * C #60696