blob: 719f573e0792f4ee72e544659fc1a2d40b7e9b53 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
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 < ssize(mat[i]); 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 j = i; // Sucht Pivotzeile für bessere Stabilität.
for (int k = 0; k < n; k++) {
if (!done[k] && abs(mat[k][i]) > abs(mat[i][i])) j = k;
}
swap(mat[i], mat[j]);
if (abs(mat[i][i]) > EPS) {
normalLine(i);
takeAll(n, i);
done[i] = true;
}}
for (int i = 0; i < n; i++) { // gauss fertig, prüfe Lösung
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;
}
|
