Efficient Numerics for Thermally-Assisted trap-limited Conduction in Chalcogenides