I have a dxd matrix that I would like to multiply with a dx1xn matrix. The idea is to get back another dx1xn matrix, so I would essentially like to multiply each dxd matrix with each dx1 matrix n times along that third axis. I have no idea how to go about this though. I have tried np.apply_along_axis and np.apply_over_axes, but the first only selects one axis rather than the dx1 matrices, and the second doesn’t do what I want either.
<code>k = 2
s = 36
d = k + s - 1
n = 1119
T = np.zeros(shape=(d, d)
alpha_filt = np.zeros(shape=(d, 1, n))
# Doesn't work obviously
np.matmul(T, alpha_filt)
</code>
<code>k = 2
s = 36
d = k + s - 1
n = 1119
T = np.zeros(shape=(d, d)
alpha_filt = np.zeros(shape=(d, 1, n))
# Doesn't work obviously
np.matmul(T, alpha_filt)
</code>
k = 2
s = 36
d = k + s - 1
n = 1119
T = np.zeros(shape=(d, d)
alpha_filt = np.zeros(shape=(d, 1, n))
# Doesn't work obviously
np.matmul(T, alpha_filt)