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
|
#include "config.h"
#include "point.h"
#include "vector.h"
flt dotproduct(struct vector * a, struct vector * b)
{
return a->dx * b->dx + a->dy * b->dy + a->dz * b->dz;
}
void between(struct point * p, struct point * q,
/*out*/ struct vector * v)
{
v->dx = q->x - p->x;
v->dy = q->y - p->y;
v->dz = q->z - p->z;
}
void opposite(struct vector * v,
/*out*/ struct vector * w)
{
w->dx = - v->dx;
w->dy = - v->dy;
w->dz = - v->dz;
}
void point_along(struct point * p, struct vector * v,
flt ac,
/*out*/ struct point * q)
{
q->x = p->x + v->dx * ac;
q->y = p->y + v->dy * ac;
q->z = p->z + v->dz * ac;
}
void product(struct vector * a, struct vector * b,
/*out*/ struct vector * v)
{
v->dx = a->dy * b->dz - a->dz * b->dy;
v->dy = a->dz * b->dx - a->dx * b->dz;
v->dz = a->dx * b->dy - a->dy * b->dx;
}
flt vlength2(struct vector * a)
{
return a->dx * a->dx + a->dy * a->dy + a->dz * a->dz;
}
flt vlength(struct vector * a)
{
return sqrt(vlength2(a));
}
void vscale(struct vector * a, flt s,
/*out*/ struct vector * v)
{
v->dx = a->dx * s;
v->dy = a->dy * s;
v->dz = a->dz * s;
}
void vnormalize(struct vector * a, /*out*/ struct vector * v)
{
vscale(a, 1 / vlength(a), v);
}
void vsub(struct vector * a, struct vector * b,
/*out*/ struct vector * v)
{
v->dx = a->dx - b->dx;
v->dy = a->dy - b->dy;
v->dz = a->dz - b->dz;
}
|