AWbd

description:

[Q_hat,bs,A_new]=AWbd(A) is the Murota-Kanno-Kojima-Kojima-Maehara numerical algorithm for the Artin-Wedderburn block-diagonal decomposition of matrix *-algebras.

arguments:

A is a matrix where each row represents a symmetric square matrix A_i

output:

Q_hat is a matrix so all the matrices Q_hat'*A_i*Q_hat are simultaneously block diagonalized

bs is a row of block sizes

A_new is a matrix where each row represents simultaneously block diagonalized (A_new_i=Q_hat'*A_i*Q_hat)

possible usage:

AWbd(A)