convert 3D matrix into diagonal block matrix

function M = mat3D2blk(A, d)

%   input:
%       A is a 3D d*d*n matrix 
%   output:
%       convert [A1;A2;...;An] into blkdiag([A1;A2;...;An])
%                   A1                       A1
%                   A2                          A2
%                   .                               .
%                   An                                 An
%   created by seamanj on 23/01/2017 in NCCA@UK


M = reshape(permute(A,[1 3 2]), [],d);
M = mat2cell(M,d*ones((size(M,1)/d),1),d);
M = blkdiag(M{:});