pub type PSolveFn = Option<unsafe extern "C" fn(P_data: *mut c_void, r: N_Vector, z: N_Vector, tol: realtype, lr: c_int) -> c_int>;