Infeasibility not detected

[danielpmorton]

Hey Preston,

These matrices below are flagged as infeasible by OSQP but not by JaxOSQP, which returns a solution that does not satisfy the constraints

(Interestingly, this actually works well for my problem, but I figured it was weird behavior)

    P = np.array([[2.0, 0.0, 0.0], [0.0, 2.0, 0.0], [0.0, 0.0, 2.0]])
    q = np.array([0.00414928, 0.03791134, 0.01410465])
    A = np.array(
        [
            [-0.76501865, 2.82754636, -2.62377954],
            [1.0, -0.0, -0.0],
            [-0.0, 1.0, -0.0],
            [-0.0, -0.0, 1.0],
            [-1.0, -0.0, -0.0],
            [-0.0, -1.0, -0.0],
            [-0.0, -0.0, -1.0],
        ]
    )
    l = np.array(
        [
            -np.inf,
            -np.inf,
            -np.inf,
            -np.inf,
            -np.inf,
            -np.inf,
            -np.inf,
        ]
    )
    u = np.array(
        [
            -86.72295701,
            0.29792536,
            0.28104433,
            0.29294767,
            0.30207464,
            0.31895567,
            0.30705233,
        ]
    )