Multidimensional matrix multiplication elementwise

For sake of simplicity, let’s just say, I have a 10 by 31 by 31 matrix. I want the 10 th pageof the 31 by 31 matrix to be multiplied with a 5139 by 31 matrix giving me a 5139 by 31 matrix as result.i want hadamard or element by element product between the square and rectangular matrix.

Please suggest.

[James Tursa]

How can you have an element-by-element product between the square and rectangular matrix? Can you give a specific example, say with a 2x3x3 matrix and a 5x3 matrix, showing the full inputs and desired output?

Side Note: In general, it is best to have your 2D pages as the first two dimensions, not the last two. That way the 2D pages are contiguous in memory.

[Raisul Islam]

So the objective is to take 31 variables and apply armairf. For 10 step ahead forecast you get 10 by 31 by 31 as the output. Now the original 31 variables have 5139 data each. So the residual matrix is 5139 by 31. I want to multiply the irf and residual. If I get 5139 rows for 31 variables I can plot the indices.

[Raisul Islam]

It’s important to achieve the 5139 rows for all variables. No matter how I change the dimension. Now I know 31 by 31 contains about 961 data or so. But by having that multiplied with 10 gets 9610 data. I need 5139 by 31. Which is much higher. Maybe I can try reshaping with repeat. I don’t know.

[James Tursa]

For this small example:

x = reshape(1:18,2,3,3);
y = reshape(1:15,5,3);

what would be the desired result, ? (i.e., show us what the elements of the result would be for this specific case)

result = ...?