improved math

1 files changed, 8 insertions, 10 deletions

diff --git a/math/gauss.cpp b/math/gauss.cpp
index 0c7ab03..3e3b7aa 100644
--- a/math/gauss.cpp
+++ b/math/gauss.cpp
@@ -1,8 +1,7 @@
-void normalLine(int n, int line) {
+void normalLine(int line) {
 	double factor = mat[line][line];
-	for (int i = 0; i <= n; i++) {
-		mat[line][i] /= factor;
-}}
+	for (double& x : mat[line]) x /= factor;
+}
 
 void takeAll(int n, int line) {
 	for (int i = 0; i < n; i++) {
@@ -21,18 +20,17 @@ int gauss(int n) {
 			if (abs(mat[j][i]) > abs(mat[i][i])) swappee = j;
 		}
 		swap(mat[i], mat[swappee]);
-		if (abs(mat[i][i]) > EPSILON) {
-			normalLine(n, i);
+		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++)
-			if (abs(mat[i][j]) > EPSILON) allZero = false;
-		if (allZero && abs(mat[i][n]) > EPSILON) return INCONSISTENT;
-		if (allZero && abs(mat[i][n]) < EPSILON) return MULTIPLE;
+		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;
 }