* Used BFGS to optimize
* Changed default sub-gaussian parameter to alpha=4
Fixed bug in post-optimization phase
Slight bugfix in source generation
Changed defaults
Now using row-major layout
Set build type to release
Optimized std::pow by casting a to integer