def testCyclic(self):
a = numpy.eye(2)
b = -1.0 * numpy.eye(2)
r = numpy.eye(2)
q = numpy.matrix([[1.0, 0.0], [0.0, 0.5]])
ds = DareSolver(a, b, q, r)
ds.use_cyclic = True
x = ds.solve()
self.assertAlmostEqual(1.6180, x[0, 0], 3)
self.assertAlmostEqual(1.0, x[1, 1], 3)
self.assertAlmostEqual(0.0, x[0, 1], 3)
self.assertAlmostEqual(0.0, x[1, 0], 3)
def testCyclic(self): a = numpy.eye(2) b = -1.0*numpy.eye(2) r = numpy.eye(2) q = numpy.matrix([[1.0,0.0],[0.0,0.5]]) ds = DareSolver(a,b,q,r) ds.use_cyclic = True x = ds.solve() self.assertAlmostEqual(1.6180,x[0,0],3) self.assertAlmostEqual(1.0,x[1,1],3) self.assertAlmostEqual(0.0,x[0,1],3) self.assertAlmostEqual(0.0,x[1,0],3)
