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

Do smth

Browse source
This commit is contained in:
AZEN-SGG 2025-05-09 12:21:22 +03:00
parent 031a771f7e
commit c10e3d3bb1
4 changed files with 68 additions and 13 deletions
Download patch file Download diff file Expand all files Collapse all files

4
2025.05.02/03Ex/solve.c
View file

@ -46,9 +46,7 @@ int t3_solve (
memcpy(&bits, &y, sizeof(bits));
sgn_c = (bits >> 63) & 1;
if (is_equal(a, c) || is_equal(c, b))
it = m+1;
else if (is_eps(y, eps))
if (is_eps(y, eps))
break;
else if (sgn_c == sgn_a)
{

50
2025.05.02/04Ex/comp.h Normal file
View file

@ -0,0 +1,50 @@
#ifndef COMP_H
#define COMP_H
#include <float.h>
#include <math.h>
static inline double * fpmax (double *pa, double *pb, double fa, double fb, double *max_f_p)
{
if ((fa - fb) > DBL_EPSILON)
{
*max_f_p = fa;
return pa;
} else
{
*max_f_p = fb;
return pb;
}
}
static inline double * fp_abs_max (double *pa, double *pb, double *fa, double *fb, double **max_f_p)
{
if ((fabs(*fa) - fabs(*fb)) > DBL_EPSILON)
{
*max_f_p = fa;
return pa;
} else
{
*max_f_p = fb;
return pb;
}
}
static inline int is_equal (const double a, const double b)
{
double diff = a - b;
double max_val = (a > b) ? a : b;
return ((diff < 0) ? -diff : diff) < (DBL_EPSILON * max_val);
}
static inline int is_null (const double a)
{
return ((a < 0) ? -a : a) < DBL_EPSILON;
}
static inline int is_eps (const double a, const double eps)
{
return (((a < 0) ? -a : a) - eps) < DBL_EPSILON;
}
#endif

16
2025.05.02/04Ex/solve.c
View file

@ -1,4 +1,5 @@
#include "solve.h"
#include "comp.h"
#include <math.h>
#include <float.h>
@ -11,15 +12,15 @@ int t4_solve (
int it = 0;
double c = DBL_MAX, y, y_a = f(a), y_b = f(b);
if (fabs(y_a) - eps < DBL_EPSILON)
if (is_eps(y_a, eps))
{
*x = a;
return 1;
} if (fabs(y_b) - eps < DBL_EPSILON)
} if (is_eps(y_b, eps))
{
*x = b;
return 1;
} if (fabs(fabs(y_b) - fabs(y_a)) < DBL_EPSILON)
} if (is_equal(fabs(y_b), fabs(y_a)))
{
*x = a;
return -1;
@ -30,17 +31,20 @@ int t4_solve (
c = b - ((b - a) / (y_b - y_a)) * y_b;
y = f(c);
if (fabs(y) - eps < DBL_EPSILON)
if (is_eps(y, eps))
break;
else if (fabs(y_a) - fabs(y_b) > DBL_EPSILON)
{
a = c;
y_a = y;
} else if (fabs(y_b) - fabs(y_a) > DBL_EPSILON)
} else
// if (fabs(y_b) - fabs(y_a) > DBL_EPSILON)
{
b = c;
y_b = y;
} else
}
if (is_equal(y_a, y_b))
it = m+1;
}

11
2025.05.02/tests/main_test.sh
View file

@ -29,18 +29,21 @@ rm -f $outlog $errlog
echo "Тест запущен..."
i=2
for (( k = 3 ; k < 7; k++ )); do
echo "------- K = $k -------"
for (( a = -100 ; a < -40 ; a++ )); do
for (( b = -9 ; b < 10 ; b++ )); do
echo "./a0$prog.out "$(echo "$a / 10" | bc -l)" "$(echo "$b / 10" | bc -l)" 1e-16 $iter $k"
./a0$prog.out "$(echo "$a / 10" | bc -l)" "$(echo "$b / 10" | bc -l)" 1e-16 $iter $k
echo "$i $(./a0$prog.out "$(echo "$a / 10" | bc -l)" "$(echo "$b / 10" | bc -l)" 1e-16 $iter $k)"
((i+=2))
done
done
for (( a = -9 ; a < 10 ; a++ )); do
for (( b = 11 ; b < 100 ; b++ )); do
echo "./a0$prog.out "$(echo "$a / 10" | bc -l)" "$(echo "$b / 10" | bc -l)" 1e-16 $iter $k"
./a0$prog.out "$(echo "$a / 10" | bc -l)" "$(echo "$b / 10" | bc -l)" 1e-16 $iter $k
echo "$i ./a0$prog.out "$(echo "$a / 10" | bc -l)" "$(echo "$b / 10" | bc -l)" 1e-16 $iter $k"
echo "$i $(./a0$prog.out "$(echo "$a / 10" | bc -l)" "$(echo "$b / 10" | bc -l)" 1e-16 $iter $k)"
((i+=2))
done
done
done >$outlog 2>$errlog