performance enhancements to decay by reducing function calls and use of std::vector::at in LMatrix
merging robert's changes from use of at() to operator[] in the LMatrix class