Convolution of NumPy arrays of arbitrary dimension for Cauchy product of multivariate power series

Question

I\'m trying to implement the idea I have suggested here, for Cauchy product of multivariate finite power series (i.e. polynomials) represented as NumPy ndarrays. numpy