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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
#include <u.h>
#include <libc.h>
#include <geometry.h>
uvlong nanosec(void);
double min(double, double);
double dppd(Point2, Point2);
double dppda(Point2, Point2);
double dppd3(Point3, Point3);
double dppd3a(Point3, Point3);
Point2 Pt2b(double, double, double);
Point3 xvec3(Point3, Point3);
double hsubpd(double, double);
double fma(double, double, double);
double
fmin(double a, double b)
{
return a<b? a: b;
}
void
main(int argc, char *argv[])
{
uvlong t0, t1;
double a, b, r;
Point2 p0, p1;
Point3 p0t, p1t, pr;
GEOMfmtinstall();
ARGBEGIN{default:sysfatal("shit");}ARGEND
if(argc != 2)
sysfatal("shit");
a = strtod(argv[0], nil);
b = strtod(argv[1], nil);
t0 = nanosec();
r = fmin(a, b);
t1 = nanosec();
print("fmin(%g, %g) = %g\ttook %lludns\n", a, b, r, t1-t0);
t0 = nanosec();
r = min(a, b);
t1 = nanosec();
print("min(%g, %g) = %g\ttook %lludns\n", a, b, r, t1-t0);
print("\n");
p0 = Pt2b(a, 1, 1);
p1 = Pt2b(b, 3, 1);
t0 = nanosec();
r = dppd(p0, p1);
t1 = nanosec();
print("dppd(%v, %v) = %g\ttook %lludns\n", p0, p1, r, t1-t0);
t0 = nanosec();
r = dotvec2(p0, p1);
t1 = nanosec();
print("dotvec2(%v, %v) = %g\ttook %lludns\n", p0, p1, r, t1-t0);
t0 = nanosec();
r = dppda(p0, p1);
t1 = nanosec();
print("dppda(%v, %v) = %g\ttook %lludns\n", p0, p1, r, t1-t0);
print("\n");
p0t = Pt3(a, 1, 9, 1);
p1t = Pt3(b, 3, 4, 1);
t0 = nanosec();
r = dppd3(p0t, p1t);
t1 = nanosec();
print("dppd3(%V, %V) = %g\ttook %lludns\n", p0t, p1t, r, t1-t0);
t0 = nanosec();
r = dotvec3(p0t, p1t);
t1 = nanosec();
print("dotvec3(%V, %V) = %g\ttook %lludns\n", p0t, p1t, r, t1-t0);
t0 = nanosec();
r = dppd3a(p0t, p1t);
t1 = nanosec();
print("dppd3a(%V, %V) = %g\ttook %lludns\n", p0t, p1t, r, t1-t0);
print("\n");
t0 = nanosec();
r = hsubpd(a, b);
t1 = nanosec();
print("hsubpd(%g, %g) = %g\ttook %lludns\n", a, b, r, t1-t0);
print("\n");
p0t = Pt3(a, 1, 9, 1);
p1t = Pt3(b, 3, 4, 1);
t0 = nanosec();
pr = xvec3(p0t, p1t);
t1 = nanosec();
print("xvec3(%V, %V) = %V\ttook %lludns\n", p0t, p1t, pr, t1-t0);
t0 = nanosec();
pr = crossvec3(p0t, p1t);
t1 = nanosec();
print("crossvec3(%V, %V) = %V\ttook %lludns\n", p0t, p1t, pr, t1-t0);
print("\n");
t0 = nanosec();
r = fma(a, b, 21);
t1 = nanosec();
print("fma(%g, %g, 21) = %g\ttook %lludns\n", a, b, r, t1-t0);
exits(nil);
}
|
