Division by zero in hybrid.reference.HybridizedPopulationAnnealing about dwave-hybrid

[randomir]

MCVE boils down to:

import dimod, hybrid

bqm = dimod.BQM({0: -1, 1: 0.01}, {}, 0, 'BINARY')
workflow = hybrid.PopulationAnnealing(num_reads=1, num_iter=2, num_sweeps=1)
workflow.run(hybrid.State.from_problem(bqm)).result()

It fails with:

/home/radomir/work/dwave-hybrid/hybrid/reference/pa.py:65: RuntimeWarning: overflow encountered in exp
  w = np.exp(-beta * ss.record.energy)
/home/radomir/work/dwave-hybrid/hybrid/reference/pa.py:66: RuntimeWarning: invalid value encountered in true_divide
  p = w / sum(w)
Traceback (most recent call last):
...
    idx = self.random.choice(len(ss), len(ss), p=p)
  File "mtrand.pyx", line 928, in numpy.random.mtrand.RandomState.choice
ValueError: probabilities contain NaN