AZEN-SGG/2nd_Sem_Bogachev
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Точно работает

Browse source
This commit is contained in:
AZEN-SGG 2025-04-08 13:52:44 +03:00
parent 48ab5717a0
commit f22449e1e1
1 changed files with 19 additions and 18 deletions
Download patch file Download diff file Expand all files Collapse all files

37
2025.04.04/dist/Krivoruchenko_SK/solve.c vendored
View file

@ -32,8 +32,8 @@ int t14_solve(int n, double * restrict A, double * restrict X, int * restrict c)
} }
} }
printf("\n------- K = %d -------\n", k); // printf("\n------- K = %d -------\n", k);
printf("Maximum = %lf i = %d j = %d\n", maximum, max_i, max_j); // printf("Maximum = %lf i = %d j = %d\n", maximum, max_i, max_j);
// Если максимальный по модулю элемент равен нулю, значит матрица вырождена // Если максимальный по модулю элемент равен нулю, значит матрица вырождена
if (fabs(maximum) < eps) if (fabs(maximum) < eps)
@ -81,19 +81,19 @@ int t14_solve(int n, double * restrict A, double * restrict X, int * restrict c)
} }
} }
printf("BEFORE GAUSS\n"); // printf("BEFORE GAUSS\n");
printf("Original matrix:\n"); // printf("Original matrix:\n");
print_matrix(A, n, n); // print_matrix(A, n, n);
printf("Inverse matrix:\n"); // printf("Inverse matrix:\n");
print_matrix(X, n, n); // print_matrix(X, n, n);
gauss_inverse(n, k, A, X); gauss_inverse(n, k, A, X);
printf("AFTER GAUSS\n"); // printf("AFTER GAUSS\n");
printf("Original matrix:\n"); // printf("Original matrix:\n");
print_matrix(A, n, n); // print_matrix(A, n, n);
printf("Inverse matrix:\n"); // printf("Inverse matrix:\n");
print_matrix(X, n, n); // print_matrix(X, n, n);
} }
gauss_back_substitution(n, A, X); gauss_back_substitution(n, A, X);
@ -136,6 +136,7 @@ int t14_solve(int n, double * restrict A, double * restrict X, int * restrict c)
return 0; return 0;
} }
// Прямой ход Го ----- йда // Прямой ход Го ----- йда
void gauss_inverse(const int n, const int k, double * restrict A, double * restrict X) void gauss_inverse(const int n, const int k, double * restrict A, double * restrict X)
{ {
@ -147,7 +148,7 @@ void gauss_inverse(const int n, const int k, double * restrict A, double * restr
for (int ij = kk+1; ij < kn+n; ij++) for (int ij = kk+1; ij < kn+n; ij++)
A[ij] *= inv_akk; A[ij] *= inv_akk;
for (int ij = kn; ij <= kk; ij++) for (int ij = kn; ij < kn + n; ij++)
X[ij] *= inv_akk; X[ij] *= inv_akk;
for (int i = k+1; i < n; ++i) for (int i = k+1; i < n; ++i)
@ -182,11 +183,11 @@ void gauss_back_substitution(const int n, double * restrict A, double * restrict
X[in + j] -= X[kn + j] * aik; X[in + j] -= X[kn + j] * aik;
} }
printf("\n------- K = %d -------\n", k); // printf("\n------- K = %d -------\n", k);
printf("Original matrix:\n"); // printf("Original matrix:\n");
print_matrix(A, n, n); // print_matrix(A, n, n);
printf("Inverse matrix:\n"); // printf("Inverse matrix:\n");
print_matrix(X, n, n); // print_matrix(X, n, n);
} }
} }