| 1 |
#include "gbasic.h"
|
| 2 |
|
| 3 |
|
| 4 |
int gb_epsorder(g3Float a,g3Float b)
|
| 5 |
{
|
| 6 |
a = a - b;
|
| 7 |
return gb_signum(a);
|
| 8 |
}
|
| 9 |
|
| 10 |
int gb_epseq(g3Float a,g3Float b)
|
| 11 |
{
|
| 12 |
return (fabs(a - b) < GB_EPSILON);
|
| 13 |
}
|
| 14 |
|
| 15 |
int gb_inrange(g3Float a,g3Float rb,g3Float re)
|
| 16 |
{
|
| 17 |
if (((a >= rb) && (a <= re)) || ((a <= rb) && (a >= re)))
|
| 18 |
return 1;
|
| 19 |
return 0;
|
| 20 |
}
|
| 21 |
|
| 22 |
int gb_signum(g3Float a)
|
| 23 |
{
|
| 24 |
return ((gb_epseq(a,0.0)) ? 0 : ((a < 0.0) ? -1 : 1));
|
| 25 |
}
|
| 26 |
|
| 27 |
|
| 28 |
g3Float gb_max(g3Float a,g3Float b)
|
| 29 |
{
|
| 30 |
return (a < b) ? b : a;
|
| 31 |
}
|
| 32 |
|
| 33 |
g3Float gb_min(g3Float a,g3Float b)
|
| 34 |
{
|
| 35 |
return (a < b) ? a : b;
|
| 36 |
}
|
| 37 |
|
| 38 |
int gb_getroots(g3Float* r1,g3Float* r2,g3Float a,g3Float b,g3Float c)
|
| 39 |
{
|
| 40 |
g3Float d;
|
| 41 |
if (gb_epseq(a,0.0)) {
|
| 42 |
if (gb_epseq(b,0.0))
|
| 43 |
return 0;
|
| 44 |
*r1 = *r2 = -c/b;
|
| 45 |
return 1;
|
| 46 |
}
|
| 47 |
d = b*b - 4*a*c;
|
| 48 |
if (d < -GB_EPSILON)
|
| 49 |
return 0;
|
| 50 |
if (gb_epseq(d,0.0)) {
|
| 51 |
*r1 = *r2 = 0.5*b/a;
|
| 52 |
return 1;
|
| 53 |
}
|
| 54 |
d = sqrt(d);
|
| 55 |
*r1 = 0.5*(b + d)/a;
|
| 56 |
*r2 = 0.5*(b - d)/a;
|
| 57 |
return 2;
|
| 58 |
}
|
