ProjetMAM/Pivot
Archived
4
1
Fork 0
Code Issues Pull Requests Actions Releases 2 Activity

fix triangular system

Browse Source
This commit is contained in:
Simon Pribylski
2024-03-05 22:23:29 +01:00
parent 911f016bb7
commit 432fa99f71
1 changed files with 12 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/Solver.cpp
View File

@@ -36,7 +36,18 @@ VectAffine Solver::TriangularSystem() const {
Vect noyau = solver.Kernel(); Vect noyau = solver.Kernel();
Matrix origin = mat.SubMatrix(0, mat.GetColumnCount() - 1, mat.GetRawCount(), 1); Matrix origin = mat.SubMatrix(0, mat.GetColumnCount() - 1, mat.GetRawCount(), 1);
return {noyau, origin}; // on rajoute des 0 si il faut
Matrix fullOrigin {mat.GetColumnCount() - 1, 1};
for (int i = 0; i < mat.GetRawCount(); i++) {
fullOrigin.at(i, 0) = origin.at(i, 0);
}
for (int i = mat.GetRawCount(); i < mat.GetColumnCount() - 1; i++) {
fullOrigin.at(i, 0) = 0;
}
return {noyau, fullOrigin};
} }
std::size_t Solver::Rank() const { std::size_t Solver::Rank() const {