First version with accumulating stats
Use cholfact() rather than factorize(Hermitian()), as indicated in https://groups.google.com/d/msg/julia-users/QBXmLGQp3Mw/IPXaxh-6Fk0J