device 0.0.4

A generative engine
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

use super::*;

pub(crate) const CENTERED_RGB: Mat4f = matrix!(
  2.0, 0.0, 0.0, -1.0;
  0.0, 2.0, 0.0, -1.0;
  0.0, 0.0, 2.0, -1.0;
  0.0, 0.0, 0.0,  1.0;
);

pub(crate) const CENTERED_RGB_INVERSE: Mat4f = matrix!(
  0.5, 0.0, 0.0, 0.5;
  0.0, 0.5, 0.0, 0.5;
  0.0, 0.0, 0.5, 0.5;
  0.0, 0.0, 0.0, 1.0;
);

pub(crate) const YCGCO: Mat3f = matrix!(
   0.25, 0.5,  0.25;
   0.50, 0.0, -0.50;
  -0.25, 0.5, -0.25;
);

pub(crate) const YCGCO_INVERSE: Mat3f = matrix!(
  1.0,  1.0, -1.0;
  1.0,  0.0,  1.0;
  1.0, -1.0, -1.0;
);

pub(crate) const YIQ: Mat3f = matrix!(
  0.2990,  0.5870,  0.1140;
  0.5959, -0.2746, -0.3213;
  0.2115, -0.5227,  0.3112;
);

pub(crate) const YIQ_INVERSE: Mat3f = matrix!(
  1.0,  0.956,  0.619;
  1.0, -0.272, -0.647;
  1.0, -1.106,  1.703;
);

pub(crate) fn invert() -> Mat4f {
  Mat3f::from_diagonal(&vector!(-1.0, -1.0, -1.0))
    .to_homogeneous()
    .append_translation(&vector!(1.0, 1.0, 1.0))
}

pub(crate) fn rotate_hue_blaster(r: f32) -> Mat4f {
  Transformation3 {
    space: Space::Blaster,
    rotation: UnitQuaternion::from_axis_angle(&Vec3f::x_axis(), r),
    ..default()
  }
  .response(1.0)
}

pub(crate) fn rotate_hue_yiq(r: f32) -> Mat4f {
  Transformation3 {
    space: Space::Yiq,
    rotation: UnitQuaternion::from_axis_angle(&Vec3f::x_axis(), r),
    ..default()
  }
  .response(1.0)
}

pub(crate) fn saturate(s: f32) -> Mat4f {
  Transformation3 {
    space: Space::Ycgco,
    scaling: Vec3f::new(1.0, s, s),
    ..default()
  }
  .response(1.0)
}

#[cfg(test)]
mod tests {
  use super::*;

  #[test]
  fn centered_rgb() {
    assert_eq!(CENTERED_RGB * CENTERED_RGB_INVERSE, Mat4f::identity());
  }

  #[test]
  fn ycgco() {
    assert_eq!(YCGCO * YCGCO_INVERSE, Mat3f::identity());
  }

  #[test]
  fn yiq() {
    assert!((YIQ * YIQ_INVERSE).is_identity(0.001));
  }
}