illuminator 0.1.3

A simple yet powerful rendering tool
Documentation 
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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96

use std::fmt::Debug;

use crate::{
    core::{spherical::xyz2spherical, tensor::Vec3f, vec::Vector},
    raycast::{bounds::Bounds3f, primitive::Primitive, *},
};

#[derive(Clone)]
pub struct Sphere {
    pub cnt: Vec3f,
    pub r: f32,
}

impl Sphere {
    pub fn new(cnt: Vec3f, r: f32) -> Sphere {
        Sphere { cnt, r }
    }

    /// y-up, z-forward
    /// return local spherical coordinates of sphere
    pub fn worldpos2sphere(&self, pos: Vec3f) -> Vec3f {
        let v = pos - self.cnt;
        xyz2spherical(v)
    }

    pub fn intersect(&self, ray_src: Vec3f, ray_dir: Vec3f) -> Option<f32> {
        // Solve t^2*d.d + 2*t*(o-p).d + (o-p).(o-p)-R^2 = 0
        let op = ray_src - self.cnt;
        let a = ray_dir.dot(ray_dir);
        let b = op.dot(ray_dir);
        let c = op.dot(op) - self.r * self.r;
        let det = b * b - c * a;
        if det < 0. {
            None
        } else {
            let det = det.sqrt();
            if -b - det >= 0. {
                return Some((-b - det) / a);
            } else if -b + det >= 0. {
                return Some((-b + det) / a);
            }
            None
        }
    }
}

impl Raycast for Sphere {
    fn raycast(&self, ray: &Ray) -> Option<Hit> {
        if let Some(t) = self.intersect(ray.org, ray.dir) {
            return Some(Hit { t });
        }
        None
    }
}

impl Debug for Sphere {
    fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
        write!(f, "{} {}", self.cnt, self.r)
    }
}

impl Primitive for Sphere {
    fn bounds(&self) -> bounds::Bounds3f {
        let r = Vec3f::vec([self.r; 3]);
        let min = self.cnt - r;
        let max = self.cnt + r;
        Bounds3f { min, max }
    }
}

#[test]
fn test_sphere() {
    let s = Sphere::new(Vec3f::vec([0.; 3]), 1.);

    for i in 0..11 {
        let y = i as f32 / 10.;
        let x = (1. - y * y).sqrt();
        let z = 0.;

        let org = Vec3f::vec([x, y, z]) * 2.;
        let dir = org * -1.;
        let ray = Ray::new(org, dir);

        let hit = s.raycast(&ray).unwrap();
        assert_eq!(hit.position(&ray)[0], x);
        assert_eq!(hit.position(&ray)[1], y);
    }

    let b = s.bounds();
    assert_eq!(b.min[0], -1.);
    assert_eq!(b.min[1], -1.);
    assert_eq!(b.min[2], -1.);
    assert_eq!(b.max[0], 1.);
    assert_eq!(b.max[1], 1.);
    assert_eq!(b.max[2], 1.);
}