Per the title, I am trying to replicate the result from Scipy when multiplying a csr_matrix with an nD array, n>2, specifically where one dimension being multiplied is intentionally a singleton. In ...