Given a square matrix Z, return a symmetric matrix S = Z + Zᵀ - diag(Z) (so S[i,j] == S[j,i]).