powered by
Function to calculate dvech(X t(X))
vdtcrossprod(X, a, d, p)
vector
matrix
indices that belong to the diagonal of the matrix
dimension of the matrix
dimension of the vectorized matrix