Computing A(:, m)*B(m, :) without a for loop

1 visualización (últimos 30 días)
Nicolas Perez
Nicolas Perez el 12 de Nov. de 2018
Comentada: Nicolas Perez el 12 de Nov. de 2018
Provided I have two matrices A and B of dimensions NxM and MxN I want to compute a matrix C of dimensions NxNxM such that each NxN slice is the result of the prduct A(:, m)*B(m, :). I know this can be obtained in the following way:
for m=1:M
C(:, :, m) = A(:, m)*B(m, :);
end
Is there a way of obtaining the same result without using a for loop (with the aim of making the computation faster)?

Respuesta aceptada

Matt J
Matt J el 12 de Nov. de 2018
Ar=reshape(A,N,1,[]);
Br=reshape(B.',1,N,[]);
C=Ar.*Br;

Más respuestas (1)

madhan ravi
madhan ravi el 12 de Nov. de 2018
C(:, :, 1:M) = A(:, 1:M)*B(1:M, :);
  1 comentario
Nicolas Perez
Nicolas Perez el 12 de Nov. de 2018
The output of A(:, 1:M)*B(1:M, :); is the sum of C(:, :, 1:M) along the third dimension (an NxN matrix). Thus, you get the "Assignment has fewer non-singleton rhs dimensions than non-singleton subscripts" error.

Iniciar sesión para comentar.

Categorías

Más información sobre Chemistry en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by