|
12 | 12 | import pathlib |
13 | 13 |
|
14 | 14 | import numpy as np |
| 15 | +import unittest |
15 | 16 |
|
16 | 17 | python_path = pathlib.Path(".").absolute().parent.parent / "python" |
17 | 18 | os.sys.path.insert(1, str(python_path)) |
|
23 | 24 |
|
24 | 25 | print(" TEST OSQP ".center(LINE_WIDTH, "-")) |
25 | 26 |
|
26 | | -problem, xs_init, us_init = create_clqr_problem() |
27 | 27 |
|
28 | | -ddp1 = CSQP(problem, "CustomOSQP") |
29 | | -ddp2 = CSQP(problem, "OSQP") |
30 | | - |
31 | | -ddp1.with_callbacks = True |
32 | | -ddp2.with_callbacks = True |
33 | | - |
34 | | -max_qp_iters = 10000 |
35 | | -ddp1.max_qp_iters = max_qp_iters |
36 | | -ddp2.max_qp_iters = max_qp_iters |
37 | | - |
38 | | -eps_abs = 1e-8 |
39 | | -eps_rel = 0.0 |
40 | | -ddp1.eps_abs = eps_abs |
41 | | -ddp2.eps_abs = eps_abs |
42 | | - |
43 | | -converged = ddp1.solve(xs_init, us_init, 1) |
44 | | -converged = ddp2.solve(xs_init, us_init, 1) |
45 | | - |
46 | | - |
47 | | -assert ddp1.qp_iters == ddp2.qp_iters |
48 | | - |
49 | | -set_tol = 1e-8 |
50 | | -assert np.linalg.norm(np.array(ddp1.xs) - np.array(ddp2.xs)) < set_tol, "Test failed" |
51 | | -assert np.linalg.norm(np.array(ddp1.us) - np.array(ddp2.us)) < set_tol, "Test failed" |
52 | | -assert np.linalg.norm(np.array(ddp1.lag_mul) - np.array(ddp2.lag_mul)) < set_tol, ( |
53 | | - "Test failed" |
54 | | -) |
55 | | -for t in range(len(ddp1.y)): |
56 | | - assert np.linalg.norm(ddp1.y[t] - ddp2.y[t]) < set_tol, "Test failed" |
| 28 | +class TestCLQROSQP(unittest.TestCase): |
| 29 | + def setUp(self): |
| 30 | + self.problem, self.xs_init, self.us_init = create_clqr_problem() |
| 31 | + self.ddp1 = CSQP(self.problem, "CustomOSQP") |
| 32 | + self.ddp2 = CSQP(self.problem, "OSQP") |
| 33 | + self.ddp1.with_callbacks = True |
| 34 | + self.ddp2.with_callbacks = True |
| 35 | + max_qp_iters = 10000 |
| 36 | + self.ddp1.max_qp_iters = max_qp_iters |
| 37 | + self.ddp2.max_qp_iters = max_qp_iters |
| 38 | + eps_abs = 1e-8 |
| 39 | + self.ddp1.eps_abs = eps_abs |
| 40 | + self.ddp2.eps_abs = eps_abs |
| 41 | + |
| 42 | + # @unittest.skip("Skipping this test as it seems that the \"OSQP\" solver does not converge as of 12/09/2025.") |
| 43 | + def test_osqp_match(self): |
| 44 | + self.ddp1.solve(self.xs_init, self.us_init, 1) |
| 45 | + self.ddp2.solve(self.xs_init, self.us_init, 1) |
| 46 | + set_tol = 1e-8 |
| 47 | + self.assertEqual(self.ddp1.qp_iters, self.ddp2.qp_iters) |
| 48 | + self.assertTrue( |
| 49 | + np.linalg.norm(np.array(self.ddp1.xs) - |
| 50 | + np.array(self.ddp2.xs)) < set_tol, |
| 51 | + "Test failed: xs mismatch", |
| 52 | + ) |
| 53 | + self.assertTrue( |
| 54 | + np.linalg.norm(np.array(self.ddp1.us) - |
| 55 | + np.array(self.ddp2.us)) < set_tol, |
| 56 | + "Test failed: us mismatch", |
| 57 | + ) |
| 58 | + self.assertTrue( |
| 59 | + np.linalg.norm(np.array(self.ddp1.lag_mul) - |
| 60 | + np.array(self.ddp2.lag_mul)) |
| 61 | + < set_tol, |
| 62 | + "Test failed: lag_mul mismatch", |
| 63 | + ) |
| 64 | + for t in range(len(self.ddp1.y)): |
| 65 | + self.assertTrue( |
| 66 | + np.linalg.norm(self.ddp1.y[t] - self.ddp2.y[t]) < set_tol, |
| 67 | + f"Test failed: y mismatch at t={t}", |
| 68 | + ) |
| 69 | + |
| 70 | + |
| 71 | +if __name__ == "__main__": |
| 72 | + unittest.main() |
0 commit comments