GPU support for breakup

[abulenok]

[github-actions]

Stale issue message

[slayoo]

Here's a recent error caught on CI:

IndexError: index 2 is out of bounds for axis 0 with size 2


i = 0, j = array([1]), k = array([2]), idx = array([0, 1])
is_first_in_pair = array([False, False]), prob_like = array([1.], dtype=float32)

IIUC, happens at the end of:

E               def Commons_pair_indices(
E                     i,
E                     j,
E                     k,
E                     idx,
E                     is_first_in_pair,
E                     prob_like
E                 ):
E                     if (prob_like[i] == 0):
E                         return True;
E                     
E                     offset = 1 - is_first_in_pair[2 * i];
E                     j[0] = idx[2 * i + offset];
E                     k[0] = idx[2 * i + 1 + offset];

[github-actions]

Stale issue message

[github-actions]

Stale issue message