performance enhancements to decay by reducing function calls and use of std::vector::at in LMatrix
hdf5back speedups by using bracket operator instead of vector::at and by hoisting some operations out of an inner loop