version 1.1.0
* Allows passing in of a custom matrix multiplication function and a custom square root of the diagonal of a matrix inverse thereby giving the opportunity to use CUDA / GPUs
* Includes a fast algorithm for inference of only one coefficient using conjugate gradient descent

version 1.0.1
* Removed OMP dependency to please CRAN

version 1.0.0
* Initial Release