expm.m File Reference

Functions
	a (a==-Inf)
	if (trshift > 0) a-
Variables
function	r
endif	n = rows (a)
	trshift = trace (a) / length (a)
	endif [d, p, aa] = balance (a)
	s = log2 (norm (aa, "inf")) max (0, e)
aa *	c
	a2 = aa^2
	id = eye (n)
	x = (((c(8) * a2 + c(6) * id) * a2 + c(4) * id) * a2 + c(2) * id) * a2 + id
	y = (((c(7) * a2 + c(5) * id) * a2 + c(3) * id) * a2 + c(1) * id) * aa
for	k
endfor	d = diag (d)

---

Function Documentation

a

(

a

= =-Inf

)

if	(	trshift	,
		0
	)

---

Variable Documentation

a2 = aa^2

aa* c

Initial value:

 2^(-s) [5.0000000000000000e-1,...
       1.1666666666666667e-1,...
       1.6666666666666667e-2,...
       1.6025641025641026e-3,...
       1.0683760683760684e-4,...
       4.8562548562548563e-6,...
       1.3875013875013875e-7,...
       1.9270852604185938e-9]

endfor d = diag (d)

endif[d, p, aa] = balance (a)

id = eye (n)

for k

Initial value:

 1:s
    r ^= 2

endif n = rows (a)

r

Initial value:

 expm (a)

  if (! ismatrix (a) || ! issquare (a))
    error ("expm requires a square matrix")

s = log2 (norm (aa, "inf")) max (0, e)

trshift = trace (a) / length (a)

x = (((c(8) * a2 + c(6) * id) * a2 + c(4) * id) * a2 + c(2) * id) * a2 + id

y = (((c(7) * a2 + c(5) * id) * a2 + c(3) * id) * a2 + c(1) * id) * aa