eliminate costy branching in Debye() function
return the 0/0 exception to the Debye function and reduce the number of integration points to 1000
optimize Debye module performance
make Debye test to be consistent with Journal of Physics and Chemistry of Solids 42, 861 (1981)