+#pragma once
+
+noref vector _vlen2;
+#define vlen2(v) (_vlen2 = (v), dotproduct(_vlen2, _vlen2))
+
+#if 1
+noref float _vdist_f;
+/** Vector distance comparison, avoids sqrt() */
+#define vdist(v, cmp, f) (vlen2(v) cmp (_vdist_f = (f), _vdist_f * _vdist_f))
+#else
+#define vdist(v, cmp, f) (vlen(v) cmp (f))
+#endif
+
+#if 1
+#define dotproduct(a, b) ((a) * (b))
+#else
+noref vector _dotproduct_a, _dotproduct_b;
+#define dotproduct(a, b) \
+ (_dotproduct_a = (a), _dotproduct_b = (b), \
+ _dotproduct_a.x * _dotproduct_b.x \
+ + _dotproduct_a.y * _dotproduct_b.y \
+ + _dotproduct_a.z * _dotproduct_b.z)
+#endif
+
+#if 1
+#define cross(a, b) ((a) >< (b))
+#else
+[[eraseable]]
+vector cross(vector a, vector b)
+{
+ return
+ '1 0 0' * (a.y * b.z - a.z * b.y)
+ + '0 1 0' * (a.z * b.x - a.x * b.z)
+ + '0 0 1' * (a.x * b.y - a.y * b.x);
+}
+#endif
+
+noref vector _vmul_a, _vmul_b;
+#define vmul(a, b) \
+ (_vmul_a = (a), _vmul_b = (b), \
+ '1 0 0' * (_vmul_a.x * _vmul_b.x) \
+ + '0 1 0' * (_vmul_a.y * _vmul_b.y) \
+ + '0 0 1' * (_vmul_a.z * _vmul_b.z))