The mass lumping parts in these two files actually forgot to factorize localF. Therefore, they will cause the effective reaction rate not the same as the real reaction rate!