Multiplying a 2d matrix with each slice of 3d matrix
38 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Søren sønderby
el 1 de Ag. de 2014
Comentada: Ray M
el 13 de Oct. de 2021
What is the fastest way multiply a 2d matrix with each slice of a 3d matrix?
x = rand(1000,200);
F = rand(100,1000,10);
b = zeros(100,200,10);
for i = 1:10
b(:,:,i) = F(:,:,i)*x;
end
2 comentarios
Shravankumar P
el 2 de Ag. de 2014
I would like to have the transposed F. i.e., F=F'; is it possible.
Respuesta aceptada
Andrei Bobrov
el 2 de Ag. de 2014
f1 = size(F);
x1 = size(x);
b = reshape(sum(bsxfun(@times,reshape(F,[f1(1),1,f1(2:3)]),...
reshape(x',1,x1(2),[])),3),f1(1),x1(2),[]);
2 comentarios
Ray M
el 11 de Oct. de 2021
Great Solution, the solution you have provided is equivalent to:
for i = 1:whatever
b(:,:,i) = F(:,:,i) * x;
end
How would you produce the trasnpose multiplication effect (assuming correct sizes)?
for i = 1:whatever
b(:,:,i) = x * F(:,:,i);
end
Más respuestas (4)
Steven Lord
el 11 de Oct. de 2021
In release R2020b we introduced the pagemtimes function for this purpose.
3 comentarios
Ray M
el 11 de Oct. de 2021
Editada: Ray M
el 11 de Oct. de 2021
I am not sure if this meant to be a reply to my question under Andrei Bobrov solution, but it does not appear so. I am simply interested in performing a single matrix multiplication by slices of a 3D matrix with the use of bsxfun as follows:
mat_mxn * 3Dmat_nxkxl
and
3Dmat_mxnxk * mat_nxh
Thanks
Steven Lord
el 12 de Oct. de 2021
@Image Analyst Your bsxfun call is calling times not mtimes. times can work with implicit expansion.
A = int16(magic(4));
B = repmat(A, 1, 1, 3);
C = A.*B
mtimes isn't fully defined for integer arrays. If mtimes doesn't work pagemtimes probably shouldn't.
D = A*A % errors
Edric Ellis
el 4 de Ag. de 2014
1 comentario
Nils Melchert
el 19 de Mayo de 2020
Editada: Nils Melchert
el 19 de Mayo de 2020
Is there the possibility to get a minimal example for exactly this use case? I am struggling with the same thing.
Image Analyst
el 2 de Ag. de 2014
You can do it easily if the number of rows and columns in your 3D and 2D match, which they don't in your example:
rows = 1000;
columns = 200;
slices = 10;
x = rand(rows, columns);
F = rand(rows, columns, slices);
b = zeros(rows, columns, slices);
for slice = 1 : slices
b(:,:, slice) = F(:,:, slice) .* x; % Use dot star, not just star.
end
If the number of rows and columns are different you need to make some decisions about exactly where you want to multiply, if one is smaller than the other, or one extends out past the other.
1 comentario
Image Analyst
el 2 de Ag. de 2014
Alternate way. Not sure which is faster:
% Mask the 3D image called "image3D" with 2D image called "mask".
masked3DImage = bsxfun(@times, image3D, cast(mask, class(image3D)));
Ver también
Categorías
Más información sobre GPU Computing en Help Center y File Exchange.
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!