diff options
| author | Gloria Mundi <gloria@gloria-mundi.eu> | 2024-11-16 01:24:14 +0100 |
|---|---|---|
| committer | Gloria Mundi <gloria@gloria-mundi.eu> | 2024-11-16 01:24:14 +0100 |
| commit | 98567ec798aa8ca2cfbcb85c774dd470f30e30d4 (patch) | |
| tree | 5113d5cc24d1ad5f93810b6442ce584a36950dc8 /math/gauss.cpp | |
| parent | ad3856a6b766087df0036de0b556f4700a6498c9 (diff) | |
| parent | 8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff) | |
mzuenni tests
Diffstat (limited to 'math/gauss.cpp')
| -rw-r--r-- | math/gauss.cpp | 36 |
1 files changed, 0 insertions, 36 deletions
diff --git a/math/gauss.cpp b/math/gauss.cpp deleted file mode 100644 index 3e3b7aa..0000000 --- a/math/gauss.cpp +++ /dev/null @@ -1,36 +0,0 @@ -void normalLine(int line) { - double factor = mat[line][line]; - for (double& x : mat[line]) x /= factor; -} - -void takeAll(int n, int line) { - for (int i = 0; i < n; i++) { - if (i == line) continue; - double diff = mat[i][line]; - for (int j = 0; j <= n; j++) { - mat[i][j] -= diff * mat[line][j]; -}}} - -int gauss(int n) { - vector<bool> done(n, false); - for (int i = 0; i < n; i++) { - int swappee = i; // Sucht Pivotzeile für bessere Stabilität. - for (int j = 0; j < n; j++) { - if (done[j]) continue; - if (abs(mat[j][i]) > abs(mat[i][i])) swappee = j; - } - swap(mat[i], mat[swappee]); - if (abs(mat[i][i]) > EPS) { - normalLine(i); - takeAll(n, i); - done[i] = true; - }} - // Ab jetzt nur checks bzgl. Eindeutigkeit/Existenz der Lösung. - for (int i = 0; i < n; i++) { - bool allZero = true; - for (int j = i; j < n; j++) allZero &= abs(mat[i][j]) <= EPS; - if (allZero && abs(mat[i][n]) > EPS) return INCONSISTENT; - if (allZero && abs(mat[i][n]) <= EPS) return MULTIPLE; - } - return UNIQUE; -} |
