1 files changed, 106 insertions, 0 deletions
diff --git a/src/core/vector.cpp b/src/core/vector.cpp
new file mode 100644
index 0000000..51d8e2e
--- /dev/null
+++ b/src/core/vector.cpp
@@ -0,0 +1,106 @@
+#include "vector.hpp"
+
+#include <math.h>
+#include <stdexcept>
+
+Vec3d::Vec3d() {
+    set(0, 0, 0);
+}
+
+Vec3d::Vec3d(double x, double y, double z) {
+    set(x, y, z);
+}
+
+void Vec3d::set(double x, double y, double z) {
+    m_x = x;
+    m_y = y;
+    m_z = z;
+}
+
+void Vec3d::normalize() {
+    auto len = length();
+    if (len == 0) {
+        throw std::runtime_error("Normalizing zero vector");
+    }
+
+    m_x /= len;
+    m_y /= len;
+    m_z /= len;
+}
+
+double Vec3d::length() const {
+    return sqrt(m_x * m_x + m_y * m_y + m_z * m_z);
+}
+
+double Vec3d::dot(const Vec3d &vec) const {
+    return m_x * vec.m_x + m_y * vec.m_y + m_z * vec.m_z;
+}
+
+Vec3d Vec3d::cross(const Vec3d &vec) const {
+    return Vec3d(
+        m_y * vec.m_z - m_z * vec.m_y,
+        m_z * vec.m_x - m_x * vec.m_z,
+        m_x * vec.m_y - m_y * vec.m_x
+        );
+}
+
+Vec3d Vec3d::operator+(const Vec3d &vec) const {
+    return Vec3d(
+            m_x + vec.m_x, 
+            m_y + vec.m_y,
+            m_z + vec.m_z
+            );
+}
+
+Vec3d& Vec3d::operator+=(const Vec3d &vec) {
+    m_x += vec.m_x;
+    m_y += vec.m_y;
+    m_z += vec.m_z;
+    return *this;
+}
+
+Vec3d Vec3d::operator-(const Vec3d &vec) const {
+    return Vec3d(
+            m_x - vec.m_x, 
+            m_y - vec.m_y,
+            m_z - vec.m_z
+            );
+}
+
+Vec3d Vec3d::operator-() const {
+    return Vec3d(
+            -m_x,
+            -m_y,
+            -m_z
+            );
+}
+
+
+Vec3d Vec3d::operator*(double op) const {
+    return Vec3d(
+            m_x * op,
+            m_y * op,
+            m_z * op
+            );
+}
+
+Vec3d Vec3d::operator*(const Vec3d &vec) const {
+    return Vec3d(
+            m_x * vec.m_x,
+            m_y * vec.m_y,
+            m_z * vec.m_z
+            );
+}
+
+Vec3d Vec3d::operator/(double op) const {
+    return Vec3d(
+            m_x / op,
+            m_y / op,
+            m_z / op
+            );
+}
+
+std::ostream& operator<<(std::ostream &out, const Vec3d &v){
+    out << "[ " << v.m_x << ", " << v.m_y << ", " << v.m_z << " ]";
+    return out;
+}