10 static float clamp(float val, float min, float max);
12 // Simple wrappers around library routines
13 static float abs(float f);
14 static float sqrt(float f);
15 static float ceil(float f);
16 static float sin(float f);
17 static float cos(float f);
18 static float tan(float f);
19 static float atan2(float y, float x);
21 // Takes two args and runs afoul of the Koenig rules.
22 static float pow(double base, double exp);
24 // double variants of the above
25 static double abs(double f);
26 static double sqrt(double f);
27 static double ceil(double f);
28 static double sin(double f);
29 static double cos(double f);
30 static double tan(double f);
31 static double atan2(double y, double x);
33 // Some 3D vector stuff. In all cases, it is permissible for the
34 // "out" vector to be the same as one of the inputs.
35 static void set3(float* v, float* out);
36 static float dot3(float* a, float* b);
37 static void cross3(float* a, float* b, float* out);
38 static void mul3(float scalar, float* v, float* out);
39 static void add3(float* a, float* b, float* out);
40 static void sub3(float* a, float* b, float* out);
41 static float mag3(float* v);
42 static void unit3(float* v, float* out);
44 // Matrix array convention: 0 1 2
48 // Multiply two matrices
49 static void mmul33(float* a, float* b, float* out);
52 static void vmul33(float* m, float* v, float* out);
54 // Multiply the vector by the matrix transpose. Or pre-multiply the
55 // matrix by v as a row vector. Same thing.
56 static void tmul33(float* m, float* v, float* out);
59 static void invert33(float* m, float* out);
61 // Transpose matrix (for an orthonormal orientation matrix, this
62 // is the same as the inverse).
63 static void trans33(float* m, float* out);
65 // Generates an orthonormal basis:
66 // xOut becomes the unit vector in the direction of x
67 // yOut is perpendicular to xOut in the x/y plane
68 // zOut becomes the unit vector: (xOut cross yOut)
69 static void ortho33(float* x, float* y,
70 float* xOut, float* yOut, float* zOut);