test_rank_deficient()1 static char* test_rank_deficient()
2 {
3   //A matrix data
4   QDLDL_int Ap[]  = {0, 1, 3};
5   QDLDL_int Ai[]  = {0, 0, 1};
6   QDLDL_float Ax[] = {1.0, 1.0, 1.0};
7   QDLDL_int An = 2;
8 
9   // RHS for Ax = b (should fail to solve)
10   QDLDL_float b[]    = {1,1};
11 
12   //x replaces b during solve
13   int status = ldl_factor_solve(An,Ap,Ai,Ax,b);
14 
15   mu_assert("Rank deficiency not detected", status < 0);
16 
17   return 0;
18 }
19