First of all congrats for the great work :)
I write this issue just to let you know that the code breaks in all the calls to projection_2d , binding_2d and unbinding_2d as the dim argument is missing (probably this is due to the recent update to the hrr package.
I suppose that adding dim=(-2,-1) will solve the problem.
Is that solution correct?