Added semi-auto verbose option 'trace=TRUE/FALSE' based on SQUAREM to print current progress, which is helpful for monitoring large data. Fixed a bug that might cause NaN due to negative value in log() due to rounding