Efficient numerics for thermally-assisted trap-limited conduction in chalcogenides