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)