summaryrefslogtreecommitdiff
path: root/src/math.zig
blob: a2cbcd29da61cdca7ed36b5a05f46c261b0aaa63 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
const std = @import("std");
pub const tan = std.math.tan;
pub const cos = std.math.cos;
pub const sin = std.math.sin;
pub const rad = std.math.degreesToRadians;

pub const Matrix = struct {
    rows: [4]@Vector(4, f32),

    pub fn lookAt(eye: @Vector(3, f32), target: @Vector(3, f32), arbitrary_up: @Vector(3, f32)) Matrix {
        const forward = normalize(eye - target);
        const right = normalize(cross(arbitrary_up, forward));
        const up = cross(forward, right);

        const view = [_]@Vector(4, f32){
            @Vector(4, f32){ right[0], right[1], right[2], 0.0 },
            @Vector(4, f32){ up[0], up[1], up[2], 0.0 },
            @Vector(4, f32){ forward[0], forward[1], forward[2], 0.0 },
            @Vector(4, f32){ 0.0, 0.0, 1.0, eye[2] },
        };

        return Matrix{
            .rows = view,
        };
    }

    pub fn perspective(fov: f32, aspect: f32, near: f32, far: f32) Matrix {
        const projection = [_]@Vector(4, f32){
            @Vector(4, f32){ 1.0 / (aspect * tan(fov / 2.0)), 0.0, 0.0, 0.0 },
            @Vector(4, f32){ 0.0, 1.0 / tan(fov / 2.0), 0.0, 0.0 },
            @Vector(4, f32){ 0.0, 0.0, -((far + near) / (far - near)), -((2 * far * near) / (far - near)) },
            @Vector(4, f32){ 0.0, 0.0, -1.0, 1.0 },
        };

        return Matrix{
            .rows = projection,
        };
    }

    pub fn identity() Matrix {
        const view = [_]@Vector(4, f32){
            @Vector(4, f32){ 1.0, 0.0, 0.0, 0.0 },
            @Vector(4, f32){ 0.0, 1.0, 0.0, 0.0 },
            @Vector(4, f32){ 0.0, 0.0, 1.0, 0.0 },
            @Vector(4, f32){ 0.0, 0.0, 0.0, 1.0 },
        };

        return Matrix{
            .rows = view,
        };
    }
};

pub fn dot(a: @Vector(3, f32), b: @Vector(3, f32)) f32 {
    return @reduce(.Add, a * b);
}

pub fn cross(a: @Vector(3, f32), b: @Vector(3, f32)) @Vector(3, f32) {
    return @Vector(3, f32){ a[1] * b[2] - a[2] * b[1], a[2] * b[0] - a[0] * b[2], a[0] * b[1] - a[1] * b[0] };
}

pub fn normalize(a: @Vector(3, f32)) @Vector(3, f32) {
    return a / @as(@Vector(3, f32), @splat(@sqrt(dot(a, a))));
}