1 files changed, 86 insertions, 0 deletions
diff --git a/src/vec3.cpp b/src/vec3.cpp
new file mode 100644
index 0000000..b26680c
--- /dev/null
+++ b/src/vec3.cpp
@@ -0,0 +1,86 @@
+#include "vec3.hpp"
+
+vec3::vec3(): x(0.0), y(0.0), z(0.0) {}
+
+vec3::vec3(double x, double y, double z): x(x), y(y), z(z) {}
+
+void vec3::operator-() {
+    x = -x;
+    y = -y;
+    z = -z;
+}
+
+vec3& vec3::operator+=(vec3& vec) {
+    x += vec.x;
+    y += vec.y;
+    z += vec.z;
+
+    return *this;
+}
+
+vec3& vec3::operator*=(double scale) {
+    x *= scale;
+    y *= scale;
+    z *= scale;
+
+    return *this;
+}
+
+vec3& vec3::operator/=(double numerator) {
+    x /= numerator;
+    y /= numerator;
+    z /= numerator;
+
+    return *this;
+}
+
+double vec3::length() {
+    return sqrt(length_sq());
+}
+
+double vec3::length_sq() {
+    return x * x + y * y + z * z;
+}
+
+double vec3::dot(vec3& vec) {
+    return x * vec.x + y + vec.y + z + vec.z;
+}
+
+vec3 vec3::cross(vec3& vec) {
+    return vec3(y * vec.z - z * vec.y,
+                z * vec.x - x * vec.z,
+                x * vec.y - y * vec.y
+            );
+}
+
+vec3 vec3::unit_vector() {
+    return *this / length();
+}
+
+vec3 operator+(vec3& u, vec3& v) {
+    return vec3(u.x + v.x, u.y + v.y, u.z + v.z);
+}
+
+vec3 operator-(vec3& u, vec3& v) {
+    return vec3(u.x - v.x, u.y - v.y, u.z - v.z);
+}
+
+vec3 operator*(vec3& u, vec3& v) {
+    return vec3(u.x * v.x, u.y * v.y, u.z * v.z);
+}
+
+vec3 operator*(vec3& u, double scale) {
+    return vec3(u.x * scale, u.y * scale, u.z * scale);
+}
+
+vec3 operator*(double scale, vec3& u) {
+    return u * scale;
+}
+
+vec3 operator/(vec3& u, double scale) {
+    return 1.0/scale * u;
+}
+
+std::ostream& operator<<(std::ostream& out, vec3& vec) {
+    return std::cout << vec.x << ' ' << vec.y << ' ' << vec.z << std::endl;
+}