test_basic()1 static char* test_basic()
2 {
3 //A matrix data
4 QDLDL_int Ap[] = {0, 1, 2, 4, 5, 6, 8, 10, 12, 14, 17};
5 QDLDL_int Ai[] = {0, 1, 1, 2, 3, 4, 1, 5, 0, 6, 3, 7, 6, 8, 1, 2, 9};
6 QDLDL_float Ax[] = {1.0, 0.460641, -0.121189, 0.417928, 0.177828,
7 0.1, -0.0290058, -1.0, 0.350321, -0.441092, -0.0845395,
8 -0.316228, 0.178663, -0.299077, 0.182452, -1.56506, -0.1};
9 QDLDL_int An = 10;
10
11 // RHS and solution to Ax = b
12 QDLDL_float b[] = {1,2,3,4,5,6,7,8,9,10};
13 QDLDL_float xsol[] = {10.2171, 3.9416, -5.69096, 9.28661, 50.0, -6.11433,
14 -26.3104, -27.7809, -45.8099, -3.74178};
15
16 //x replaces b during solve
17 int status = ldl_factor_solve(An,Ap,Ai,Ax,b);
18
19 mu_assert("Factorisation failed", status >= 0);
20 mu_assert("Solve accuracy failed", vec_diff_norm(b,xsol,An) < QDLDL_TESTS_TOL);
21
22 return 0;
23 }
24