summaryrefslogtreecommitdiff
path: root/math/gauss.cpp
diff options
context:
space:
mode:
authorGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
committerGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
commit98567ec798aa8ca2cfbcb85c774dd470f30e30d4 (patch)
tree5113d5cc24d1ad5f93810b6442ce584a36950dc8 /math/gauss.cpp
parentad3856a6b766087df0036de0b556f4700a6498c9 (diff)
parent8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff)
mzuenni tests
Diffstat (limited to 'math/gauss.cpp')
-rw-r--r--math/gauss.cpp36
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;
-}