36 lines
608 B
C
36 lines
608 B
C
#include "solve.h"
|
|
#include "math.h"
|
|
|
|
#define eps 1e-9
|
|
|
|
int compare(const double a, const double b)
|
|
{
|
|
if (a - b > eps) return 1;
|
|
else if (b - a > eps) return -1;
|
|
return 0;
|
|
}
|
|
|
|
double t5_solve(const double *A, const double *b, int n, int m)
|
|
{
|
|
double temp, summa, maximum = 0;
|
|
int i, j, k;
|
|
|
|
if (n <= 0 || m <= 0) return 0;
|
|
|
|
for (i = 0; i < n; i++)
|
|
{
|
|
summa = 0;
|
|
for (k = 0; k < n; k++) {
|
|
temp = 0;
|
|
for (j = 0; j < m; j++)
|
|
temp += A[i * m + j] * b[j * n + k];
|
|
|
|
if (k == i) temp--;
|
|
summa += fabs(temp);
|
|
}
|
|
if (compare(summa, maximum) == 1) maximum = summa;
|
|
}
|
|
|
|
return maximum;
|
|
}
|
|
|