Trying also truncating state vector for speed.
Trying to use symmetry properties of the exp matrices
Fixed error
Fixed error
Lapac was not faster
Added comment
Merge branch 'single_precision'