Pack all the vectors together
Replaced ViennaCL by Eigen (small problem size...)
* Used BFGS to optimize
* Changed default sub-gaussian parameter to alpha=4