I'm not completely clear on when the switch over to the array interface in was made, but some years ago there was a change to scipy sparse matrix multiplication such that the * operator performs element-wise multiplication, breaking the code logic here somewhat invisibly. The matmul (@) operator was implemented way back in 0.17.0 so it might be worth switching to that operator to avoid getting caught out, if it's not already affected.