Fix for linear-normal-inverse-gamma-gaussian when multiplying matrix is not square.
Fix for linear-normal-inverse-gamma-gaussian when multiplying matrix is
not square.
lawmurray/Birch,lawmurray/Birch,lawmurray/Birch