diff --git a/src/tests/vec4_tests.cpp b/src/tests/vec4_tests.cpp
new file mode 100644
index 0000000..5186e30
--- /dev/null
+++ b/src/tests/vec4_tests.cpp
@@ -0,0 +1,178 @@
+#include "tenvex.h"
+
+#include "gtest/gtest.h"
+
+namespace tnvx
+{
+
+TEST(vec4, add)
+{
+	vec4 a = { 1.0f, 2.0f, 3.0f, 4.0f };
+	vec4 b = { 4.0f, 3.0f, 2.0f, 1.0f };
+	vec4 check{ 5.0f, 5.0f, 5.0f, 5.0f };
+	vec4 result = a + b;
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, sub)
+{
+	vec4 a = { 5.0f, 5.0f, 5.0f, 5.0f };
+	vec4 b = { 4.0f, 3.0f, 2.0f, 1.0f };
+	vec4 check = { 1.0f, 2.0f, 3.0f, 4.0f };
+	vec4 result = a - b;
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, mult_a)
+{
+	vec4 a = { 1.0f, 2.0f, 0.0f, 0.0f };
+	float b = 3.0f;
+	vec4 check = { 3.0f, 6.0f, 0.0f, 0.0f };
+	vec4 result = a * b;
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, mult_b)
+{
+	float a = 2.0f;
+	vec4 b = { 1.0f, 0.0f, 2.0f, 0.0f };
+	vec4 check = { 2.0f, 0.0f, 4.0f, 0.0f };
+	vec4 result = a * b;
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, dot)
+{
+	vec4 a = { 1.0f, 0.0f, 3.0f, 0.0f };
+	vec4 b = { 1.0f, 2.0f, 3.0f, 1.0f };
+	float check = 10.0f;
+	float result = dot(a, b);
+
+	EXPECT_FLOAT_EQ(check, result);
+}
+
+TEST(vec4, cross)
+{
+	vec4 a = { 1.0f, 0.0f, 0.0f, 0.0f };
+	vec4 b = { 0.0f, 1.0f, 0.0f, 0.0f };
+	vec4 check = { 0.0f, 0.0f, 1.0f, 0.0f };
+	vec4 result = cross(a, b);
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, norm_a)
+{
+	vec4 a = { 1.0f, 0.0f, 0.0f, 0.0f };
+	vec4 check = { 1.0f, 0.0f, 0.0f, 0.0f };
+	vec4 result = norm(a);
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, norm_b)
+{
+	vec4 a = { 1.0f, 0.0f, 1.0f, 0.0f };
+	vec4 check = { 0.707106f, 0.0f, 0.707106f, 0.0f };
+	vec4 result = norm(a);
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, magnitude_a)
+{
+	vec4 a = { 1.0f, 0.0f, 0.0f, 0.0f };
+	float check = 1.0f;
+	float result = magnitude(a);
+
+	EXPECT_FLOAT_EQ(check, result);
+}
+
+TEST(vec4, magnitude_b)
+{
+	vec4 a = { 1.0f, 0.0f, 1.0f, 0.0f };
+	float check = 1.4142135f;
+	float result = magnitude(a);
+
+	EXPECT_FLOAT_EQ(check, result);
+}
+
+TEST(vec4, expression_a)
+{
+	vec4 a = { 1.0f, 1.0f, 1.0f, 0.0f };
+	vec4 b = { 0.0f, 1.0f, 0.0f, 0.0f };
+	vec4 c = { 1.15f, 0.0f, 11.0f, 0.0f };
+
+	vec4 check = { 3.45f, 0.0f, 33.0f, 0.0f };
+
+	vec4 result = dot(a, b) * 3.0f * c;
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, expression_b)
+{
+	vec4 a = { 1.0f, 1.0f, 1.0f, 0.0f };
+	vec4 b = { 0.0f, 1.0f, 0.0f, 0.0f };
+	vec4 c = { 1.15f, 0.0f, 11.0f, 0.0f };
+
+	vec4 check = { 33.0f, 0.0f, -3.45f, 0.0f };
+
+	vec4 result = dot(a, b) * 3.0f * cross(b, c);
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, expression_c)
+{
+	vec4 a = { 1.0f, 2.0f, 3.0f, 0.0f };
+	vec4 b = { 4.0f, 5.0f, 6.0f, 0.0f };
+	vec4 c = { 1.0f, 0.0f, 1.0f, 0.0f };
+
+	vec4 check = { 7.243f, 5.657f, 10.071f, 0.0f };
+
+	vec4 result = norm(a + b * 2.0f) * dot(b, c) + c * 3.0f;
+
+	EXPECT_TRUE(approx_eq(check, result));
+}
+
+TEST(vec4, cross_anticommutative)
+{
+	vec4 a = { 1.0f, 0.0f, 0.0f, 0.0f };
+	vec4 b = { 0.0f, 1.0f, 0.0f, 0.0f };
+	EXPECT_TRUE(approx_eq(cross(a, b), cross(b, a) * -1.0f));
+}
+
+TEST(vec4, dot_commutative)
+{
+	vec4 a = { 1.0f, 2.0f, 3.0f, 0.0f };
+	vec4 b = { 4.0f, 5.0f, 6.0f, 0.0f };
+	EXPECT_FLOAT_EQ(dot(a, b), dot(b, a));
+}
+
+TEST(vec4, norm_magnitude)
+{
+	vec4 a = { 3.0f, 1.0f, 4.0f, 0.0f };
+	EXPECT_NEAR(magnitude(norm(a)), 1.0f, 1e-5f);
+}
+
+TEST(vec4, dot_perpendicular)
+{
+	vec4 a = { 1.0f, 0.0f, 0.0f, 0.0f };
+	vec4 b = { 0.0f, 1.0f, 0.0f, 0.0f };
+	EXPECT_FLOAT_EQ(dot(a, b), 0.0f);
+}
+
+TEST(vec4, cross_parallel)
+{
+	vec4 a = { 1.0f, 0.0f, 0.0f, 0.0f };
+	vec4 b = { 2.0f, 0.0f, 0.0f, 0.0f };
+	vec4 check = { 0.0f, 0.0f, 0.0f, 0.0f };
+	EXPECT_TRUE(approx_eq(cross(a, b), check));
+}
+
+}