Commit71fef4b

committed

Check eigenvalues' real part using NumPy

Faster than converting to a pure Python list.

1 parent9a33509 commit71fef4bCopy full SHA for 71fef4b

File tree

-3

lines changed

-3

lines changed

Lines changed: 2 additions & 2 deletions

Original file line number	Diff line number	Diff line change
`@@ -149,7 +149,7 @@ def modred(sys, ELIM, method='matchdc'):`
`149`	`149`
`150`	`150`
`151`	`151`	`#Check system is stable`
`152`		`-ifany(e.real>=0.0foreinnp.linalg.eigvals(sys.A)):`
	`152`	`+ifnp.any(np.linalg.eigvals(sys.A).real>=0.0):`
`153`	`153`	`raiseValueError("Oops, the system is unstable!")`
`154`	`154`
`155`	`155`	`ELIM=np.sort(ELIM)`
`@@ -251,7 +251,7 @@ def balred(sys, orders, method='truncate'):`
`251`	`251`	`dico='C'`
`252`	`252`
`253`	`253`	`#Check system is stable`
`254`		`-ifany(e.real>=0.0foreinnp.linalg.eigvals(sys.A)):`
	`254`	`+ifnp.any(np.linalg.eigvals(sys.A).real>=0.0):`
`255`	`255`	`raiseValueError("Oops, the system is unstable!")`
`256`	`256`
`257`	`257`	`ifmethod=='matchdc':`

Lines changed: 1 addition & 1 deletion

Original file line number	Diff line number	Diff line change
`@@ -355,7 +355,7 @@ def gram(sys,type):`
`355`	`355`
`356`	`356`	`#TODO: Check system is stable, perhaps a utility in ctrlutil.py`
`357`	`357`	`# or a method of the StateSpace class?`
`358`		`-ifany(e.real>=0.0foreinnp.linalg.eigvals(sys.A)):`
	`358`	`+ifnp.any(np.linalg.eigvals(sys.A).real>=0.0):`
`359`	`359`	`raiseValueError("Oops, the system is unstable!")`
`360`	`360`
`361`	`361`	`iftype=='c':`

Comments

(0)