1. Hd_i == Hd_j so only compute one of them; 2. add permut array to work_set
add new Q and Q_bar computing routine (v2) and update coordinate descent routine accordingly
accelerate updateLBFGS routine by using blas lv2
fix some memory leaks
print out to file LHAC_log
track elapsed time instead of cputime
ability to shrink the gamma in lbfgs
add sufficient decrease
print out to an additional log file (more brief)
maintenance
track # chol factorizations
added function to read input from text files
bug fixes
maintenance
prepare to add functionality to time gradient/function evaluations and subproblem solving