I have programmed up the LSCF, MLFD, and p-LSCF methods and am having trouble with the final piece p-MLFD. The dimensions of the matrices in the early part of the paper don't seem to be right and I think the definition of the X and Y matrices are swapped. The definitions of the R, S, T, V, & W matrices are mostly correct. I could be messing up the un-vec operator of delta vec(D) or delta vec(No). I believe the un-vec of a column array should be a 2 row or column matrix with one row/column being every other value in the vec array.
Has anyone tried to code this up?