Commitbf6c446

committed

Added singularity check to the residualization method

* Changed the unit test in matlab_test.py so that no error is thrown

1 parent1fe639c commitbf6c446Copy full SHA for bf6c446

File tree

-3

lines changed

-3

lines changed

Lines changed: 2 additions & 2 deletions

Original file line number	Diff line number	Diff line change
`@@ -181,8 +181,8 @@ def modred(sys, ELIM, method='matchdc'):`
`181`	`181`	`# if matchdc, residualize`
`182`	`182`
`183`	`183`	`# Check if the matrix A22 is invertible`
`184`		`-#if np.linalg.matrix_rank(A22) != len(ELIM):`
`185`		`-# raise ValueError("Matrix A22 is singular to working precision.")`
	`184`	`+ifnp.linalg.matrix_rank(A22)!=len(ELIM):`
	`185`	`+raiseValueError("Matrix A22 is singular to working precision.")`
`186`	`186`
`187`	`187`	`# Now precompute A22\A21 and A22\B2 (A22I = inv(A22))`
`188`	`188`	`# We can solve two linear systems in one pass, since the`

Lines changed: 1 addition & 1 deletion

Original file line number	Diff line number	Diff line change
`@@ -363,7 +363,7 @@ def testBalred(self):`
`363`	`363`
`364`	`364`	`deftestModred(self):`
`365`	`365`	`modred(self.siso_ss1, [1])`
`366`		`-modred(self.siso_ss2*self.siso_ss3, [1,2])`
	`366`	`+modred(self.siso_ss2*self.siso_ss3, [0,1])`
`367`	`367`	`modred(self.siso_ss3, [1],'matchdc')`
`368`	`368`	`modred(self.siso_ss3, [1],'truncate')`
`369`	`369`

Comments

(0)