Index: tests/vector_math/glm_equivalents/vec2_float.cc
===================================================================
--- tests/vector_math/glm_equivalents/vec2_float.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
+++ tests/vector_math/glm_equivalents/vec2_float.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
@@ -0,0 +1,28 @@
+#include <glm/glm.hpp>
+#include <iostream>
+
+using namespace glm;
+
+std::ostream& operator<<(std::ostream& os, glm::tvec2<float>& v) {
+    os << "<" << v.x << "," << v.y << ">";
+    return os;
+}
+std::ostream& operator<<(std::ostream& os, glm::tvec2<float>&& v) {
+    os << "<" << v.x << "," << v.y << ">";
+    return os;
+}
+
+float length_squared(glm::tvec2<float> v) {
+    return glm::length(v) * glm::length(v);
+}
+
+tvec2<float> project(glm::tvec2<float> u, glm::tvec2<float> v) {
+    return normalize(v) * dot(u, normalize(v));
+}
+
+int main() {
+    tvec2<float> v1 = {5.f,6.f};
+    tvec2<float> v2 = {0.f,-1.f};
+    std::cout << "refract:" << refract(v1,normalize(v2),1.f) << std::endl;
+    std::cout << "refract:" << refract(v1,normalize(v2),1.f/1.33f) << std::endl;
+}
Index: tests/vector_math/glm_equivalents/vec2_int.cc
===================================================================
--- tests/vector_math/glm_equivalents/vec2_int.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
+++ tests/vector_math/glm_equivalents/vec2_int.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
@@ -0,0 +1,120 @@
+#include <glm/glm.hpp>
+#include <iostream>
+
+using namespace glm;
+
+std::ostream& operator<<(std::ostream& os, glm::tvec2<int>& v) {
+    os << "<" << v.x << "," << v.y << ">";
+    return os;
+}
+std::ostream& operator<<(std::ostream& os, glm::tvec2<int>&& v) {
+    os << "<" << v.x << "," << v.y << ">";
+    return os;
+}
+
+int length_squared(glm::tvec2<float> v) {
+    return dot(v, v);
+}
+
+// tvec2<int> project(glm::tvec2<int> u, glm::tvec2<int> v) {
+//     return normalize(v) * dot(u, normalize(v));
+// }
+
+int main(void) {
+    tvec2<int> v1 = {1,2};
+    std::cout << "ctor(x,y):" << v1 << std::endl;
+
+    tvec2<int> v2 = v1;
+    std::cout << "copy ctor:" << v2 << std::endl;
+
+    v2 = (tvec2<int>){3, 4};
+    std::cout << "assignment:" << v2 << std::endl;
+
+    v2 = v1;
+    std::cout << "move assignment:" << v2 << std::endl;
+
+
+                // NOT SUPPORTED BY GLM
+    // tvec2<int> v3 = 0;
+    // std::cout << "zero-init:" << v3 << std::endl;
+
+    // v1 = 0;
+    // std::cout << "zero-assign:" << v1 << std::endl;
+
+    // tvec2<int> v4 = {1.23f};
+    // std::cout << "fill-ctor:" << v4 << std::endl;
+
+
+    v1 = (tvec2<int>){1, 3};
+    std::cout << "?-?:" << (v1 - (tvec2<int>){1003, -14}) << std::endl;
+
+    v1 -= (tvec2<int>){1003, -14};
+    std::cout << "?-=?:" << v1 << std::endl;
+
+    v1 = -v1;
+    std::cout << "-?:" << v1 << std::endl;
+
+    v1 = (tvec2<int>){1, 3};
+    std::cout << "?+?:" << (v1 + (tvec2<int>){1003, -14}) << std::endl;
+
+    v1 += (tvec2<int>){1003, -14};
+    std::cout << "?+=?:" << v1 << std::endl;
+
+    v1 = (tvec2<int>){15, 275};
+    std::cout << "v*s:" << v1 * 3 << std::endl;
+
+    std::cout << "s*v:" << 3 * v1 << std::endl;
+
+    v1 *= 3;
+    std::cout << "?*=?:" << v1 << std::endl;
+
+    v1 = (tvec2<int>){21, -10};
+    std::cout << "?/?:" << (v1 / 3) << std::endl;
+
+    v1 /= 3;
+    std::cout << "?/=?:" << v1 << std::endl;
+
+    //      FORCE GLM TO FLOAT, BUT WORKS
+    tvec2<float> v1_f = (tvec2<float>){2.f, 3};
+    tvec2<float> v2_f = (tvec2<float>){-3.f, 2.f};
+    std::cout << "dot_1:" << dot(v1_f, v2_f) << std::endl;
+
+    v2_f = (tvec2<float>){13, 2};
+    std::cout << "dot_2:" << dot(v1_f, v2_f) << std::endl;
+
+
+    //      NOT IN GLM BUT
+    //      USES DOT, BUT THAT'S FLOATING-POINT ONLY!
+    std::cout << "length_squared:" << length_squared(v1_f) << std::endl;
+
+    //      FLOATING-POINT ONLY GLM
+    // v1 = (tvec2<int>){100, -101};
+    // v2 = (tvec2<int>){6, 3};
+    // std::cout << "distance:" << distance(v1, v2) << std::endl;
+    //
+    // std::cout << "normalize:" << normalize(v2) << std::endl;
+    //
+    // v1 = (tvec2<int>){1,0};
+    // std::cout << "normalize_2:" << normalize(v1) << std::endl;
+    //
+    // std::cout << "project:" << project((tvec2<int>){5,6}, (tvec2<int>){1, 0}) << std::endl;
+    // std::cout << "project_2:" << project((tvec2<int>){5,6}, (tvec2<int>){1, 1}) << std::endl;
+    //
+    // v1 = (tvec2<int>){5,6};
+    // v2 = (tvec2<int>){1,0};
+    // std::cout << "reflect:" << reflect(v1,v2) << std::endl;
+    //
+    // v2 = (tvec2<int>){0,-1};
+    // std::cout << "refract:" << refract(v1,normalize(v2),1) << std::endl;
+    // std::cout << "refract:" << refract(v1,normalize(v2),1.f/1.33f) << std::endl;
+    //
+    //
+    //      FLOATING-POINT ONLY IN GLM BUT JUST USES DOT SO...
+    tvec2<float> geometric_normal = {5,6};
+    tvec2<float> perturbed_normal = {4,5};
+    tvec2<float> eyeline = {-1,0};
+    std::cout << "faceforward_nochange:" << faceforward(perturbed_normal, eyeline, geometric_normal) << std::endl;
+
+    eyeline = (tvec2<int>){1,0};
+    std::cout << "faceforward_flip:" << faceforward(perturbed_normal, eyeline, geometric_normal) << std::endl;
+}
Index: tests/vector_math/glm_equivalents/vec3_float.cc
===================================================================
--- tests/vector_math/glm_equivalents/vec3_float.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
+++ tests/vector_math/glm_equivalents/vec3_float.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
@@ -0,0 +1,117 @@
+#include <glm/glm.hpp>
+#include <iostream>
+#include <ostream>
+
+using namespace glm;
+
+std::ostream& operator<<(std::ostream& os, glm::tvec3<float>& v) {
+    os << "<" << v.x << "," << v.y << "," << v.z << ">";
+    return os;
+}
+std::ostream& operator<<(std::ostream& os, glm::tvec3<float>&& v) {
+    os << "<" << v.x << "," << v.y << "," << v.z << ">";
+    return os;
+}
+
+float length_squared(glm::tvec3<float> v) {
+    return glm::length(v) * glm::length(v);
+}
+
+tvec3<float> project(glm::tvec3<float> u, glm::tvec3<float> v) {
+    return normalize(v) * dot(u, normalize(v));
+}
+
+int main(void) {
+    tvec3<float> v1 = {1.f,2.f,3.f};
+    std::cout << "ctor(x,y):" << v1 << std::endl;
+
+    tvec3<float> v2 = v1;
+    std::cout << "copy ctor:" << v2 << std::endl;
+
+    v2 = (tvec3<float>){3.f, 4.2f, -2.f};
+    std::cout << "assignment:" << v2 << std::endl;
+
+    v2 = v1;
+    std::cout << "move assignment:" << v2 << std::endl;
+
+    // tvec3<float> v3 = 0;
+    // std::cout << "zero-init:" << v3 << std::endl;
+
+    // v1 = 0;
+    // std::cout << "zero-assign:" << v1 << std::endl;
+    //
+    // tvec3<float> v4 = {1.23f};
+    // std::cout << "fill-ctor:" << v4 << std::endl;
+    //
+    v1 = (tvec3<float>){1.23f, 3.43f, 0.000002f};
+    std::cout << "?-?:" << (v1 - (tvec3<float>){1.21f,3,1}) << std::endl;
+
+    v1 -= (tvec3<float>){1.23f, 3.43f, 0.000002f};
+    std::cout << "?-=?:" << v1 << std::endl;
+
+    v1 = -v1;
+    std::cout << "-?:" << v1 << std::endl;
+
+
+
+
+
+    v1 = (tvec3<float>){1.5f, 2.75f, -14.2f};
+    std::cout << "?+?:" << (v1 + (tvec3<float>){0.8f, -0.3f, 5}) << std::endl;
+
+    v1 += (tvec3<float>){0.8f, -0.3f, 5};
+    std::cout << "?+=?:" << v1 << std::endl;
+
+    v1 = (tvec3<float>){1.5f, 2.75f, 100.3f};
+    std::cout << "v*s:" << v1 * 3.f << std::endl;
+
+    std::cout << "s*v:" << 3.f * v1 << std::endl;
+
+    v1 *= 3;
+    std::cout << "?*=?:" << v1 << std::endl;
+
+    v1 = (tvec3<float>){2, -0.1f, 45};
+    std::cout << "?/?:" << (v1 / 3.f) << std::endl;
+
+    v1 /= 3.f;
+    std::cout << "?/=?:" << v1 << std::endl;
+
+
+    v1 = (tvec3<float>){4,2,3};
+    v2 = (tvec3<float>){0,-3,2};
+    std::cout << "dot_1:" << dot(v1, v2) << std::endl;
+
+    v2 = (tvec3<float>){1.3f, -2, 12.2};
+    std::cout << "dot_2:" << dot(v1, v2) << std::endl;
+
+
+    v1 = (tvec3<float>){1,2,3};
+    std::cout << "length:" << length(v1) << std::endl;
+
+    std::cout << "length_squared:" << length_squared(v1) << std::endl;
+
+    v2 = (tvec3<float>){6, -3.2f, 1};
+    std::cout << "distance:" << distance(v1, v2) << std::endl;
+
+    std::cout << "normalize:" << normalize(v2) << std::endl;
+
+    std::cout << "project:" << project((tvec3<float>){5,6,0}, (tvec3<float>){0, 0, 1}) << std::endl;
+    std::cout << "project_2:" << project((tvec3<float>){5,6,3.2f}, (tvec3<float>){10, 1, 4}) << std::endl;
+
+    v1 = (tvec3<float>){5,6,2.333f};
+    v2 = (tvec3<float>){1,0,-13.5f};
+    std::cout << "reflect:" << reflect(v1,normalize(v2)) << std::endl;
+
+    v2 = (tvec3<float>){0,-1,2};
+    std::cout << "refract:" << refract(normalize(v1),normalize(v2),1.f) << std::endl;
+    std::cout << "refract:" << refract(normalize(v1),normalize(v2),1.f/1.33f) << std::endl;
+
+    tvec3<float> geometric_normal = {5,6,1};
+    tvec3<float> perturbed_normal = {4,5.5f,2};
+    tvec3<float> eyeline = {-1,0.002f,-1.0345f};
+    std::cout << "faceforward_nochange:" << faceforward(perturbed_normal, eyeline, geometric_normal) << std::endl;
+
+    eyeline = (tvec3<float>){1,0.002f,-1.0345f};
+    std::cout << "faceforward_flip:" << faceforward(perturbed_normal, eyeline, geometric_normal) << std::endl;
+}
+
Index: tests/vector_math/glm_equivalents/vec4_float.cc
===================================================================
--- tests/vector_math/glm_equivalents/vec4_float.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
+++ tests/vector_math/glm_equivalents/vec4_float.cc	(revision 9b026f1dbc3026997339356f6745d77521eb933c)
@@ -0,0 +1,114 @@
+#include <glm/glm.hpp>
+#include <iostream>
+#include <ostream>
+
+using namespace glm;
+
+std::ostream& operator<<(std::ostream& os, glm::tvec4<float>& v) {
+    os << "<" << v.x << "," << v.y << "," << v.z << "," << v.w << ">";
+    return os;
+}
+std::ostream& operator<<(std::ostream& os, glm::tvec4<float>&& v) {
+    os << "<" << v.x << "," << v.y << "," << v.z << "," << v.w << ">";
+    return os;
+}
+
+float length_squared(glm::tvec4<float> v) {
+    return glm::length(v) * glm::length(v);
+}
+
+tvec4<float> project(glm::tvec4<float> u, glm::tvec4<float> v) {
+    return normalize(v) * dot(u, normalize(v));
+}
+
+int main(void) {
+    tvec4<float> v1 = {1.f,2.f,3.f,4.f};
+    std::cout << "ctor(x,y):" << v1 << std::endl;
+
+    tvec4<float> v2 = v1;
+    std::cout << "copy ctor:" << v2 << std::endl;
+
+    v2 = (tvec4<float>){3.f,4.2f,-2.f,-100.f};
+    std::cout << "assignment:" << v2 << std::endl;
+
+    v2 = v1;
+    std::cout << "move assignment:" << v2 << std::endl;
+
+    // tvec4<float> v3 = 0;
+    // std::cout << "zero-init:" << v3 << std::endl;
+
+    // v1 = 0;
+    // std::cout << "zero-assign:" << v1 << std::endl;
+    //
+    // tvec4<float> v4 = {1.23f};
+    // std::cout << "fill-ctor:" << v4 << std::endl;
+    //
+    v1 = (tvec4<float>){1.23f, 3.43f, 0.000002f, -0.000002f};
+    std::cout << "?-?:" << (v1 - (tvec4<float>){1.21f,3,1,-0.000001f}) << std::endl;
+
+    v1 -= (tvec4<float>){1.21f,3,1,-0.000001f};
+    std::cout << "?-=?:" << v1 << std::endl;
+
+    v1 = -v1;
+    std::cout << "-?:" << v1 << std::endl;
+
+
+    v1 = (tvec4<float>){1.5f, 2.75f, -14.2f, -13.5f};
+    std::cout << "?+?:" << (v1 + (tvec4<float>){0.8f, -0.3f, 5, 1}) << std::endl;
+
+    v1 += (tvec4<float>){0.8f, -0.3f, 5, 1};
+    std::cout << "?+=?:" << v1 << std::endl;
+
+    v1 = (tvec4<float>){1.5f, 2.75f, 100.3f, -2.f};
+    std::cout << "v*s:" << v1 * 3.f << std::endl;
+
+    std::cout << "s*v:" << 3.f * v1 << std::endl;
+
+    v1 *= 3;
+    std::cout << "?*=?:" << v1 << std::endl;
+
+    v1 = (tvec4<float>){2, -0.1f, 45, -15};
+    std::cout << "?/?:" << (v1 / 3.f) << std::endl;
+
+    v1 /= 3.f;
+    std::cout << "?/=?:" << v1 << std::endl;
+
+
+    v1 = (tvec4<float>){4,3,2,1};
+    v2 = (tvec4<float>){0,4,-3,-6};
+    std::cout << "dot_1:" << dot(v1, v2) << std::endl;
+
+    v2 = (tvec4<float>){1.3f, -2, 12.2, 24.1345};
+    std::cout << "dot_2:" << dot(v1, v2) << std::endl;
+
+
+    v1 = (tvec4<float>){1,2,3,4};
+    std::cout << "length:" << length(v1) << std::endl;
+
+    std::cout << "length_squared:" << length_squared(v1) << std::endl;
+
+    v2 = (tvec4<float>){6, -3.2f, 1, 12};
+    std::cout << "distance:" << distance(v1, v2) << std::endl;
+
+    std::cout << "normalize:" << normalize(v2) << std::endl;
+
+    std::cout << "project:" << project((tvec4<float>){5,6,0,-1}, (tvec4<float>){0, 0, 0, 1}) << std::endl;
+    std::cout << "project_2:" << project((tvec4<float>){5,6,3.2f,-16.75f}, (tvec4<float>){10, 1, 4, 2}) << std::endl;
+
+    v1 = (tvec4<float>){5,6,2.333f,1};
+    v2 = (tvec4<float>){1,0,-13.5f,2};
+    std::cout << "reflect:" << reflect(v1,normalize(v2)) << std::endl;
+
+    v2 = (tvec4<float>){0,-1,2,0.5f};
+    std::cout << "refract:" << refract(normalize(v1),normalize(v2),1.f) << std::endl;
+    std::cout << "refract:" << refract(normalize(v1),normalize(v2),1.f/1.33f) << std::endl;
+
+    tvec4<float> geometric_normal = {5,6,1,2};
+    tvec4<float> perturbed_normal = {4,5.5f,2,2};
+    tvec4<float> eyeline = {-1,0.002f,-1.0345f,-2};
+    std::cout << "faceforward_nochange:" << faceforward(perturbed_normal, eyeline, geometric_normal) << std::endl;
+
+    eyeline = (tvec4<float>){1,0.002f,-1.0345f,-2};
+    std::cout << "faceforward_flip:" << faceforward(perturbed_normal, eyeline, geometric_normal) << std::endl;
+}
+