Index: tests/vector/.expect/vec2.txt
===================================================================
--- tests/vector/.expect/vec2.txt	(revision 11881953ca54519ac6075af439004447ec006cda)
+++ 	(revision )
@@ -1,31 +1,0 @@
-ctor(x,y):<1,2>
-copy ctor:<1,2>
-assignment:<3,4.2>
-move assignment:<1,2>
-zero-init:<0,0>
-zero-assign:<0,0>
-fill-ctor:<1.23,1.23>
-?-?:<0.02,0.43>
-?-=?:<3.23,-6.57>
--?:<-3.23,6.57>
-?+?:<2.3,2.45>
-?+=?:<2.3,2.45>
-v*s:<4.5,8.25>
-s*v:<4.5,8.25>
-?*=?:<4.5,8.25>
-?/?:<0.666667,-0.0333333>
-?/=?:<0.666667,-0.0333333>
-dot_1:0.
-dot_2:8.6
-length:5.
-length_squared:25.
-distance:6.5146
-normalize:<0.882353,-0.470588>
-normalize_2:<1,0>
-project:<5,0>
-project_2:<5.5,5.5>
-reflect:<-5,6>
-refract:<5,6>
-refract:<3.7594,4.5592>
-faceforward_nochange:<4,5.5>
-faceforward_flip:<-4,-5.5>
Index: tests/vector/.expect/vec3_float.txt
===================================================================
--- tests/vector/.expect/vec3_float.txt	(revision 1dd929e5b6f43607cda0fc28bc6feab1dc9ec530)
+++ tests/vector/.expect/vec3_float.txt	(revision 1dd929e5b6f43607cda0fc28bc6feab1dc9ec530)
@@ -0,0 +1,30 @@
+ctor(x,y):<1.,2.,3.>
+copy ctor:<1.,2.,3.>
+assignment:<3.,4.2,-2.>
+move assignment:<1.,2.,3.>
+zero-init:<0.,0.,0.>
+zero-assign:<0.,0.,0.>
+fill-ctor:<1.23,1.23,1.23>
+?-?:<0.02,0.43,-0.999998>
+?-=?:<0.,0.,0.>
+-?:<-0.,-0.,-0.>
+?+?:<2.3,2.45,-9.2>
+?+=?:<2.3,2.45,-9.2>
+v*s:<4.5,8.25,300.9>
+s*v:<4.5,8.25,300.9>
+?*=?:<4.5,8.25,300.9>
+?/?:<0.666667,-0.0333333,15.>
+?/=?:<0.666667,-0.0333333,15.>
+dot_1:0.
+dot_2:37.8
+length:3.74166
+length_squared:14.
+distance:7.48599
+normalize:<0.872964,-0.465581,0.145494>
+project:<0.,0.,0.>
+project_2:<5.88034,0.588034,2.35214>
+reflect:<5.28917,6.,-1.57084>
+refract:<0.613403,0.736083,0.286214>
+refract:<0.461205,0.82471,-0.327329>
+faceforward_nochange:<4.,5.5,2.>
+faceforward_flip:<-4.,-5.5,-2.>
Index: tests/vector/glm_equivalents/vec2_float.cc
===================================================================
--- tests/vector/glm_equivalents/vec2_float.cc	(revision 1dd929e5b6f43607cda0fc28bc6feab1dc9ec530)
+++ tests/vector/glm_equivalents/vec2_float.cc	(revision 1dd929e5b6f43607cda0fc28bc6feab1dc9ec530)
@@ -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/glm_equivalents/vec3_float.cc
===================================================================
--- tests/vector/glm_equivalents/vec3_float.cc	(revision 1dd929e5b6f43607cda0fc28bc6feab1dc9ec530)
+++ tests/vector/glm_equivalents/vec3_float.cc	(revision 1dd929e5b6f43607cda0fc28bc6feab1dc9ec530)
@@ -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/vec2.cfa
===================================================================
--- tests/vector/vec2.cfa	(revision 11881953ca54519ac6075af439004447ec006cda)
+++ 	(revision )
@@ -1,93 +1,0 @@
-#include "../../libcfa/src/vector.hfa"
-#include <fstream.hfa>
-
-int main(void) {
-    vec2 v1 = {1.,2.};
-    sout | "ctor(x,y):" | v1;
-
-    vec2 v2 = v1;
-    sout | "copy ctor:" | v2;
-
-    v2 = (vec2){3., 4.2};
-    sout | "assignment:" | v2;
-
-    v2 = v1;
-    sout | "move assignment:" | v2;
-
-    vec2 v3 = 0;
-    sout | "zero-init:" | v3;
-
-    v1 = 0;
-    sout | "zero-assign:" | v1;
-
-    vec2 v4 = {1.23};
-    sout | "fill-ctor:" | v4;
-
-    v1 = (vec2){1.23, 3.43};
-    sout | "?-?:" | (v1 - (vec2){1.21,3});
-
-    v1 -= (vec2){-2, 10};
-    sout | "?-=?:" | v1;
-
-    v1 = -v1;
-    sout | "-?:" | v1;
-
-    v1 = (vec2){1.5, 2.75};
-    sout | "?+?:" | (v1 + (vec2){0.8, -0.3});
-
-    v1 += (vec2){0.8, -0.3};
-    sout | "?+=?:" | v1;
-
-    v1 = (vec2){1.5, 2.75};
-    sout | "v*s:" | v1 * 3;
-
-    sout | "s*v:" | 3 * v1;
-
-    v1 *= 3;
-    sout | "?*=?:" | v1;
-
-    v1 = (vec2){2, -0.1};
-    sout | "?/?:" | (v1 / 3);
-
-    v1 /= 3;
-    sout | "?/=?:" | v1;
-
-    v1 = (vec2){2, 3};
-    v2 = (vec2){-3, 2};
-    sout | "dot_1:" | dot(v1, v2);
-
-    v2 = (vec2){1.3, 2};
-    sout | "dot_2:" | dot(v1, v2);
-
-    v1 = (vec2){4, 3};
-    sout | "length:" | length(v1);
-
-    sout | "length_squared:" | length_squared(v1);
-
-    v2 = (vec2){6, -3.2};
-    sout | "distance:" | distance(v1, v2);
-
-    sout | "normalize:" | normalize(v2);
-
-    v1 = (vec2){1,0};
-    sout | "normalize_2:" | normalize(v1);
-
-    sout | "project:" | project((vec2){5,6}, (vec2){1, 0});
-    sout | "project_2:" | project((vec2){5,6}, (vec2){1, 1});
-
-    v1 = (vec2){5,6};
-    v2 = (vec2){1,0};
-    sout | "reflect:" | reflect(v1,v2);
-
-    v2 = (vec2){0,-1};
-    sout | "refract:" | refract(v1,v2,1);
-    sout | "refract:" | refract(v1,v2,1./1.33);
-
-    vec2 geometric_normal = {5,6};
-    vec2 perturbed_normal = {4,5.5};
-    vec2 eyeline = {-1,0.002f};
-    sout | "faceforward_nochange:" | faceforward(perturbed_normal, eyeline, geometric_normal);
-
-    eyeline = (vec2){1,0.002f};
-    sout | "faceforward_flip:" | faceforward(perturbed_normal, eyeline, geometric_normal);
-}
Index: tests/vector/vec2_float.cfa
===================================================================
--- tests/vector/vec2_float.cfa	(revision 11881953ca54519ac6075af439004447ec006cda)
+++ tests/vector/vec2_float.cfa	(revision 1dd929e5b6f43607cda0fc28bc6feab1dc9ec530)
@@ -81,6 +81,6 @@
 
     v2 = (vec2(float)){0,-1};
-    sout | "refract:" | refract(v1,v2,1);
-    sout | "refract:" | refract(v1,v2,1.f/1.33f);
+    sout | "refract:" | refract(v1,normalize(v2),1);
+    sout | "refract:" | refract(v1,normalize(v2),1.f/1.33f);
 
     vec2(float) geometric_normal = {5,6};