each value of A multiplied by B

Call signature
('ij,kl->ijkl', A, B)
NumPy equivalent
A[:, :, None, None] * B