vec3.hpp

Go to the documentation of this file.

#pragma once

 
#include <cmath>
 
#include "orrery/core/types.hpp"
 
namespace orrery::core {

struct Vec3 {
    Real x{};
    Real y{};
    Real z{};
 
    constexpr Vec3& operator+=(Vec3 other) noexcept {
        x += other.x;
        y += other.y;
        z += other.z;
        return *this;
    }
 
    constexpr Vec3& operator-=(Vec3 other) noexcept {
        x -= other.x;
        y -= other.y;
        z -= other.z;
        return *this;
    }
 
    constexpr Vec3& operator*=(Real scale) noexcept {
        x *= scale;
        y *= scale;
        z *= scale;
        return *this;
    }
 
    constexpr Vec3& operator/=(Real divisor) noexcept {
        x /= divisor;
        y /= divisor;
        z /= divisor;
        return *this;
    }

    [[nodiscard]] constexpr bool operator==(const Vec3& other) const = default;
};
 
[[nodiscard]] constexpr Vec3 operator+(Vec3 a, Vec3 b) noexcept {
    return {a.x + b.x, a.y + b.y, a.z + b.z};
}
 
[[nodiscard]] constexpr Vec3 operator-(Vec3 a, Vec3 b) noexcept {
    return {a.x - b.x, a.y - b.y, a.z - b.z};
}
 
[[nodiscard]] constexpr Vec3 operator-(Vec3 v) noexcept {
    return {-v.x, -v.y, -v.z};
}
 
[[nodiscard]] constexpr Vec3 operator*(Vec3 v, Real scale) noexcept {
    return {v.x * scale, v.y * scale, v.z * scale};
}
 
[[nodiscard]] constexpr Vec3 operator*(Real scale, Vec3 v) noexcept {
    return v * scale;
}
 
[[nodiscard]] constexpr Vec3 operator/(Vec3 v, Real divisor) noexcept {
    return {v.x / divisor, v.y / divisor, v.z / divisor};
}
 
[[nodiscard]] constexpr Real dot(Vec3 a, Vec3 b) noexcept {
    return (a.x * b.x) + (a.y * b.y) + (a.z * b.z);
}
 
[[nodiscard]] constexpr Vec3 cross(Vec3 a, Vec3 b) noexcept {
    return {(a.y * b.z) - (a.z * b.y), (a.z * b.x) - (a.x * b.z), (a.x * b.y) - (a.y * b.x)};
}

[[nodiscard]] constexpr Real squared_norm(Vec3 v) noexcept {
    return dot(v, v);
}

[[nodiscard]] inline Real norm(Vec3 v) noexcept {
    return std::sqrt(squared_norm(v));
}
 
} // namespace orrery::core