Las exponenciales matriciales de las matrices oblicuas-hermitianas son baratas de calcular:
Suponga que es su matriz oblicua - ermitaña , luego i A es hermitiana y, a través de zheevd y sus amigos, puede obtener la descomposiciónUNAi A
i A = UΛ UH,
donde es la matriz unitaria de vector propio y Λ es real y diagonal. Entonces, trivialmente,UΛ
A = U( - i Λ ) UH.
Una vez que tienes y Λ , es fácil calcularUΛ
Exp( A ) = exp( U( - i Λ ) UH) = UExp( - i Λ ) UH
exponiendo primero los valores propios, estableciendo mediante zcopy , realizando B : = B exp ( - i Λ ) ejecutando zscal en cada columna con un valor propio exponencial, y finalmente estableciendo su resultado enB : = UB : = B exp( - i Λ )
Exp( A ) : = B UH
a través de zgemm .