when the phonon expansion order is quite large, a segfault occurred due to the sice of an array. That array size was increased. In addition, additional output is written to the standard output stream to inform the user of calculation progress.
when the beta grid contains very small values, the trans subroutine calculates a delta value of roughly that size. For a value of for example 1e-13, an inifite loop can occur because the updated beta value (beta+delta) remains the same when beta is large enough. When that happens, we now increase the delta value to compensate.
This PR fixes a few issues in LEAPR:
test 80 was added to verify this last issue.
PR #314 must be merged prior to this one.