The function mc_sandwich is just an auxiliar
function to compute product matrix in the sandwich form
bord1 * middle * bord2. An special case appears when
computing the derivative of the covariance matrix with respect to
the power parameter. Always the bord1 and bord2 should be
diagonal matrix. If it is not true, this product is too slow.