Skip to content

Conversation

@chris-maes
Copy link
Contributor

When a basis had no singletons we were incorrectly setting the rank to m, even if the right-looking LU reported that the basis was rank deficient. This caused the column pointers and row indices in the L factor to be invalid. (Since the L factor itself was invalid).

We now correctly handle rank deficiency in the basis when there are no singletons.

@chris-maes chris-maes requested a review from a team as a code owner October 8, 2025 20:41
@chris-maes chris-maes requested review from hlinsen and rg20 October 8, 2025 20:41
@copy-pr-bot
Copy link

copy-pr-bot bot commented Oct 8, 2025

This pull request requires additional validation before any workflows can run on NVIDIA's runners.

Pull request vetters can view their responsibilities here.

Contributors can view more details about this message here.

@chris-maes chris-maes added bug Something isn't working doc Improvements or additions to documentation non-breaking Introduces a non-breaking change dual_simplex and removed doc Improvements or additions to documentation labels Oct 8, 2025
@chris-maes chris-maes added this to the 25.10 milestone Oct 8, 2025
@chris-maes
Copy link
Contributor Author

/ok to test 1792314

for (i_t p = col_start; p < col_end; ++p) {
const i_t i = this->i[p];
if (row_marker[i] == j) { printf("CSC error: repeated row index %d in column %d\n", i, j); }
if (i < 0 || i >= this->m) { printf("CSC error: row index %d not in range [0, %d]\n", i, this->m-1); }
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We should probably use cuopt_assert for these checks

@chris-maes
Copy link
Contributor Author

/ok to test 40200ea

@chris-maes
Copy link
Contributor Author

/merge

@rapids-bot rapids-bot bot merged commit 905b4e4 into NVIDIA:branch-25.10 Oct 9, 2025
88 checks passed
@chris-maes chris-maes deleted the factorize_basis_fix branch October 9, 2025 01:23
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

bug Something isn't working dual_simplex non-breaking Introduces a non-breaking change

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants