Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

su2 matrix exponentiation

Hi, does anyone know if sage math is able to successfully calculate exp(G) where G belongs to su(2) ? The sagemath commands are given below :

var('phi theta', domain='real')

u = vector([cos(phi)*sin(theta),sin(phi)*sin(theta),cos(theta)])
u.norm().simplify_trig()

sigma_x = matrix([[0,1],[1,0]])
sigma_y = matrix([[0,-i],[i,0]])
sigma_z = matrix([[1,0],[0,-1]])
sigma_x, sigma_y, sigma_z

var('t',domain='real')

G = -i*t*(sigma_x*u[0]+sigma_y*u[1]+sigma_z*u[2])
G.trace(), (G.det()/t**2).simplify_full(), G.is_hermitian()

 exp(G)

returns the following error code :

TypeError: ECL says: Error executing code in Maxima: Unable to find the spectral representation

Thank you for helping

Epi

su2 matrix exponentiation

Hi, does anyone know if sage math is able to successfully calculate exp(G) where G belongs to su(2) ? The sagemath commands are given below :

var('phi theta', domain='real')

u = vector([cos(phi)*sin(theta),sin(phi)*sin(theta),cos(theta)])
u.norm().simplify_trig()

sigma_x = matrix([[0,1],[1,0]])
sigma_y = matrix([[0,-i],[i,0]])
sigma_z = matrix([[1,0],[0,-1]])
sigma_x, sigma_y, sigma_z

var('t',domain='real')

G = -i*t*(sigma_x*u[0]+sigma_y*u[1]+sigma_z*u[2])
G.trace(), (G.det()/t**2).simplify_full(), G.is_hermitian()

 exp(G)

exp(G) returns the following error code :

TypeError: ECL says: Error executing code in Maxima: Unable to find the spectral representation

Thank you for helping

Epi

su2 matrix exponentiation

Hi, Hi,

does anyone know if sage math is able to successfully calculate exp(G) where G belongs to su(2) ? The sagemath commands are given below :

var('phi theta', domain='real')

u = vector([cos(phi)*sin(theta),sin(phi)*sin(theta),cos(theta)])
u.norm().simplify_trig()

sigma_x = matrix([[0,1],[1,0]])
sigma_y = matrix([[0,-i],[i,0]])
sigma_z = matrix([[1,0],[0,-1]])
sigma_x, sigma_y, sigma_z

var('t',domain='real')

G = -i*t*(sigma_x*u[0]+sigma_y*u[1]+sigma_z*u[2])
G.trace(), (G.det()/t**2).simplify_full(), G.is_hermitian()

 exp(G)

exp(G) returns the following error code :

TypeError: ECL says: Error executing code in Maxima: Unable to find the spectral representation

Thank you for helpinghelping.

Epi