BUGS: Many fixes, notably a "safe" log function to handle zeros, and a correction to the entrop calculation for the matrix-variate posteriors