Efficient inference in matrix-variate Gaussian models with iid observation noise