crossprod(x) <- t(x) %*% y