Dear Maple Support Team and Community,
I would like to report what seems to be a bug in Maple's MatrixExponential
function in the LinearAlgebra
package.
Consider the following 8×8 skew-Hermitian matrix B:
B := Matrix(8, 8, [[0, I, 0, -I, 0, I, 0, -I],
[-I, 0, I, 0, -I, 0, I, 0],
[0, -I, 0, I, 0, -I, 0, I],
[I, 0, -I, 0, I, 0, -I, 0],
[0, I, 0, -I, 0, I, 0, -I],
[-I, 0, I, 0, -I, 0, I, 0],
[0, -I, 0, I, 0, -I, 0, I],
[I, 0, -I, 0, I, 0, -I, 0]]);
Now compute the matrix exponential exp(I*t*B) using MatrixExponential
H := LinearAlgebra:-MatrixExponential(I * t * B);
However, this result is inconsistent with the matrix exponential computed using MATLAB and SageMath, both of which give the same result, different from Maple's.
I would appreciate it if someone from the development team could look into this.