-
Notifications
You must be signed in to change notification settings - Fork 0
/
qurt.c
64 lines (64 loc) · 1.07 KB
/
qurt.c
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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
#include "common.h"
int qurt(double a[], double x1[], double x2[])
{
double d, w1, w2;
if(a[0] == 0.0)
return (999);
d = a[1] * a[1] - 4 * a[0] * a[2];
w1 = 2.0 * a[0];
w2 = sqrt(fabs(d));
if(d > 0.0)
{
x1[0] = (-a[1] + w2) / w1;
x1[1] = 0.0;
x2[0] = (-a[1] - w2) / w1;
x2[1] = 0.0;
return (0);
}
else if(d == 0.0)
{
x1[0] = -a[1] / w1;
x1[1] = 0.0;
x2[0] = x1[0];
x2[1] = 0.0;
return (0);
}
else
{
w2 /= w1;
x1[0] = -a[1] / w1;
x1[1] = w2;
x2[0] = x1[0];
x2[1] = -w2;
return (0);
}
}
int main()
{
double a[3], x1[2], x2[2];
int result;
a[0] = 1.75;
a[1] = -3.2;
a[2] = 2.45;
qurt(a, x1, x2);
result = *(int *) &x1[0];
a[0] = 1.5;
a[1] = -2.5;
a[2] = 1.5;
qurt(a, x1, x2);
result += *(int *) &x1[1];
a[0] = 1.8;
a[1] = -4.275;
a[2] = 8.31;
qurt(a, x1, x2);
result -= *(int *) &x1[1];
if(result != -1776094907)
{
puts("qurt: fail\n");
}
else
{
puts("qurt: success\n");
}
return 0;
}