MatLDA_Grid function

Fits the JJ-class penalized matrix-normal model for a two-dimensional grid of tuning parameters. Used for tuning parameter selection.