joonsrenderer 1.1-java

data/src/main/java/org/sunflow/core/tesselatable/Teapot.java

@@ -0,0 +1,237 @@
+package org.sunflow.core.tesselatable;
+
+public class Teapot extends BezierMesh {
+    // teapot data, from: http://www.cs.ucsb.edu/~cs280/winter2004/hw2/
+
+    private static final float[][] PATCHES = {
+        {-80.00f, 0.00f, 30.00f, -80.00f, -44.80f, 30.00f, -44.80f,
+            -80.00f, 30.00f, 0.00f, -80.00f, 30.00f, -80.00f, 0.00f,
+            12.00f, -80.00f, -44.80f, 12.00f, -44.80f, -80.00f, 12.00f,
+            0.00f, -80.00f, 12.00f, -60.00f, 0.00f, 3.00f, -60.00f,
+            -33.60f, 3.00f, -33.60f, -60.00f, 3.00f, 0.00f, -60.00f,
+            3.00f, -60.00f, 0.00f, 0.00f, -60.00f, -33.60f, 0.00f,
+            -33.60f, -60.00f, 0.00f, 0.00f, -60.00f, 0.00f,},
+        {0.00f, -80.00f, 30.00f, 44.80f, -80.00f, 30.00f, 80.00f, -44.80f,
+            30.00f, 80.00f, 0.00f, 30.00f, 0.00f, -80.00f, 12.00f,
+            44.80f, -80.00f, 12.00f, 80.00f, -44.80f, 12.00f, 80.00f,
+            0.00f, 12.00f, 0.00f, -60.00f, 3.00f, 33.60f, -60.00f,
+            3.00f, 60.00f, -33.60f, 3.00f, 60.00f, 0.00f, 3.00f, 0.00f,
+            -60.00f, 0.00f, 33.60f, -60.00f, 0.00f, 60.00f, -33.60f,
+            0.00f, 60.00f, 0.00f, 0.00f,},
+        {-60.00f, 0.00f, 90.00f, -60.00f, -33.60f, 90.00f, -33.60f,
+            -60.00f, 90.00f, 0.00f, -60.00f, 90.00f, -70.00f, 0.00f,
+            69.00f, -70.00f, -39.20f, 69.00f, -39.20f, -70.00f, 69.00f,
+            0.00f, -70.00f, 69.00f, -80.00f, 0.00f, 48.00f, -80.00f,
+            -44.80f, 48.00f, -44.80f, -80.00f, 48.00f, 0.00f, -80.00f,
+            48.00f, -80.00f, 0.00f, 30.00f, -80.00f, -44.80f, 30.00f,
+            -44.80f, -80.00f, 30.00f, 0.00f, -80.00f, 30.00f,},
+        {0.00f, -60.00f, 90.00f, 33.60f, -60.00f, 90.00f, 60.00f, -33.60f,
+            90.00f, 60.00f, 0.00f, 90.00f, 0.00f, -70.00f, 69.00f,
+            39.20f, -70.00f, 69.00f, 70.00f, -39.20f, 69.00f, 70.00f,
+            0.00f, 69.00f, 0.00f, -80.00f, 48.00f, 44.80f, -80.00f,
+            48.00f, 80.00f, -44.80f, 48.00f, 80.00f, 0.00f, 48.00f,
+            0.00f, -80.00f, 30.00f, 44.80f, -80.00f, 30.00f, 80.00f,
+            -44.80f, 30.00f, 80.00f, 0.00f, 30.00f,},
+        {-56.00f, 0.00f, 90.00f, -56.00f, -31.36f, 90.00f, -31.36f,
+            -56.00f, 90.00f, 0.00f, -56.00f, 90.00f, -53.50f, 0.00f,
+            95.25f, -53.50f, -29.96f, 95.25f, -29.96f, -53.50f, 95.25f,
+            0.00f, -53.50f, 95.25f, -57.50f, 0.00f, 95.25f, -57.50f,
+            -32.20f, 95.25f, -32.20f, -57.50f, 95.25f, 0.00f, -57.50f,
+            95.25f, -60.00f, 0.00f, 90.00f, -60.00f, -33.60f, 90.00f,
+            -33.60f, -60.00f, 90.00f, 0.00f, -60.00f, 90.00f,},
+        {0.00f, -56.00f, 90.00f, 31.36f, -56.00f, 90.00f, 56.00f, -31.36f,
+            90.00f, 56.00f, 0.00f, 90.00f, 0.00f, -53.50f, 95.25f,
+            29.96f, -53.50f, 95.25f, 53.50f, -29.96f, 95.25f, 53.50f,
+            0.00f, 95.25f, 0.00f, -57.50f, 95.25f, 32.20f, -57.50f,
+            95.25f, 57.50f, -32.20f, 95.25f, 57.50f, 0.00f, 95.25f,
+            0.00f, -60.00f, 90.00f, 33.60f, -60.00f, 90.00f, 60.00f,
+            -33.60f, 90.00f, 60.00f, 0.00f, 90.00f,},
+        {80.00f, 0.00f, 30.00f, 80.00f, 44.80f, 30.00f, 44.80f, 80.00f,
+            30.00f, 0.00f, 80.00f, 30.00f, 80.00f, 0.00f, 12.00f,
+            80.00f, 44.80f, 12.00f, 44.80f, 80.00f, 12.00f, 0.00f,
+            80.00f, 12.00f, 60.00f, 0.00f, 3.00f, 60.00f, 33.60f,
+            3.00f, 33.60f, 60.00f, 3.00f, 0.00f, 60.00f, 3.00f, 60.00f,
+            0.00f, 0.00f, 60.00f, 33.60f, 0.00f, 33.60f, 60.00f, 0.00f,
+            0.00f, 60.00f, 0.00f,},
+        {0.00f, 80.00f, 30.00f, -44.80f, 80.00f, 30.00f, -80.00f, 44.80f,
+            30.00f, -80.00f, 0.00f, 30.00f, 0.00f, 80.00f, 12.00f,
+            -44.80f, 80.00f, 12.00f, -80.00f, 44.80f, 12.00f, -80.00f,
+            0.00f, 12.00f, 0.00f, 60.00f, 3.00f, -33.60f, 60.00f,
+            3.00f, -60.00f, 33.60f, 3.00f, -60.00f, 0.00f, 3.00f,
+            0.00f, 60.00f, 0.00f, -33.60f, 60.00f, 0.00f, -60.00f,
+            33.60f, 0.00f, -60.00f, 0.00f, 0.00f,},
+        {60.00f, 0.00f, 90.00f, 60.00f, 33.60f, 90.00f, 33.60f, 60.00f,
+            90.00f, 0.00f, 60.00f, 90.00f, 70.00f, 0.00f, 69.00f,
+            70.00f, 39.20f, 69.00f, 39.20f, 70.00f, 69.00f, 0.00f,
+            70.00f, 69.00f, 80.00f, 0.00f, 48.00f, 80.00f, 44.80f,
+            48.00f, 44.80f, 80.00f, 48.00f, 0.00f, 80.00f, 48.00f,
+            80.00f, 0.00f, 30.00f, 80.00f, 44.80f, 30.00f, 44.80f,
+            80.00f, 30.00f, 0.00f, 80.00f, 30.00f,},
+        {0.00f, 60.00f, 90.00f, -33.60f, 60.00f, 90.00f, -60.00f, 33.60f,
+            90.00f, -60.00f, 0.00f, 90.00f, 0.00f, 70.00f, 69.00f,
+            -39.20f, 70.00f, 69.00f, -70.00f, 39.20f, 69.00f, -70.00f,
+            0.00f, 69.00f, 0.00f, 80.00f, 48.00f, -44.80f, 80.00f,
+            48.00f, -80.00f, 44.80f, 48.00f, -80.00f, 0.00f, 48.00f,
+            0.00f, 80.00f, 30.00f, -44.80f, 80.00f, 30.00f, -80.00f,
+            44.80f, 30.00f, -80.00f, 0.00f, 30.00f,},
+        {56.00f, 0.00f, 90.00f, 56.00f, 31.36f, 90.00f, 31.36f, 56.00f,
+            90.00f, 0.00f, 56.00f, 90.00f, 53.50f, 0.00f, 95.25f,
+            53.50f, 29.96f, 95.25f, 29.96f, 53.50f, 95.25f, 0.00f,
+            53.50f, 95.25f, 57.50f, 0.00f, 95.25f, 57.50f, 32.20f,
+            95.25f, 32.20f, 57.50f, 95.25f, 0.00f, 57.50f, 95.25f,
+            60.00f, 0.00f, 90.00f, 60.00f, 33.60f, 90.00f, 33.60f,
+            60.00f, 90.00f, 0.00f, 60.00f, 90.00f,},
+        {0.00f, 56.00f, 90.00f, -31.36f, 56.00f, 90.00f, -56.00f, 31.36f,
+            90.00f, -56.00f, 0.00f, 90.00f, 0.00f, 53.50f, 95.25f,
+            -29.96f, 53.50f, 95.25f, -53.50f, 29.96f, 95.25f, -53.50f,
+            0.00f, 95.25f, 0.00f, 57.50f, 95.25f, -32.20f, 57.50f,
+            95.25f, -57.50f, 32.20f, 95.25f, -57.50f, 0.00f, 95.25f,
+            0.00f, 60.00f, 90.00f, -33.60f, 60.00f, 90.00f, -60.00f,
+            33.60f, 90.00f, -60.00f, 0.00f, 90.00f,},
+        {-64.00f, 0.00f, 75.00f, -64.00f, 12.00f, 75.00f, -60.00f, 12.00f,
+            84.00f, -60.00f, 0.00f, 84.00f, -92.00f, 0.00f, 75.00f,
+            -92.00f, 12.00f, 75.00f, -100.00f, 12.00f, 84.00f,
+            -100.00f, 0.00f, 84.00f, -108.00f, 0.00f, 75.00f, -108.00f,
+            12.00f, 75.00f, -120.00f, 12.00f, 84.00f, -120.00f, 0.00f,
+            84.00f, -108.00f, 0.00f, 66.00f, -108.00f, 12.00f, 66.00f,
+            -120.00f, 12.00f, 66.00f, -120.00f, 0.00f, 66.00f,},
+        {-60.00f, 0.00f, 84.00f, -60.00f, -12.00f, 84.00f, -64.00f,
+            -12.00f, 75.00f, -64.00f, 0.00f, 75.00f, -100.00f, 0.00f,
+            84.00f, -100.00f, -12.00f, 84.00f, -92.00f, -12.00f,
+            75.00f, -92.00f, 0.00f, 75.00f, -120.00f, 0.00f, 84.00f,
+            -120.00f, -12.00f, 84.00f, -108.00f, -12.00f, 75.00f,
+            -108.00f, 0.00f, 75.00f, -120.00f, 0.00f, 66.00f, -120.00f,
+            -12.00f, 66.00f, -108.00f, -12.00f, 66.00f, -108.00f,
+            0.00f, 66.00f,},
+        {-108.00f, 0.00f, 66.00f, -108.00f, 12.00f, 66.00f, -120.00f,
+            12.00f, 66.00f, -120.00f, 0.00f, 66.00f, -108.00f, 0.00f,
+            57.00f, -108.00f, 12.00f, 57.00f, -120.00f, 12.00f, 48.00f,
+            -120.00f, 0.00f, 48.00f, -100.00f, 0.00f, 39.00f, -100.00f,
+            12.00f, 39.00f, -106.00f, 12.00f, 31.50f, -106.00f, 0.00f,
+            31.50f, -80.00f, 0.00f, 30.00f, -80.00f, 12.00f, 30.00f,
+            -76.00f, 12.00f, 18.00f, -76.00f, 0.00f, 18.00f,},
+        {-120.00f, 0.00f, 66.00f, -120.00f, -12.00f, 66.00f, -108.00f,
+            -12.00f, 66.00f, -108.00f, 0.00f, 66.00f, -120.00f, 0.00f,
+            48.00f, -120.00f, -12.00f, 48.00f, -108.00f, -12.00f,
+            57.00f, -108.00f, 0.00f, 57.00f, -106.00f, 0.00f, 31.50f,
+            -106.00f, -12.00f, 31.50f, -100.00f, -12.00f, 39.00f,
+            -100.00f, 0.00f, 39.00f, -76.00f, 0.00f, 18.00f, -76.00f,
+            -12.00f, 18.00f, -80.00f, -12.00f, 30.00f, -80.00f, 0.00f,
+            30.00f,},
+        {68.00f, 0.00f, 51.00f, 68.00f, 26.40f, 51.00f, 68.00f, 26.40f,
+            18.00f, 68.00f, 0.00f, 18.00f, 104.00f, 0.00f, 51.00f,
+            104.00f, 26.40f, 51.00f, 124.00f, 26.40f, 27.00f, 124.00f,
+            0.00f, 27.00f, 92.00f, 0.00f, 78.00f, 92.00f, 10.00f,
+            78.00f, 96.00f, 10.00f, 75.00f, 96.00f, 0.00f, 75.00f,
+            108.00f, 0.00f, 90.00f, 108.00f, 10.00f, 90.00f, 132.00f,
+            10.00f, 90.00f, 132.00f, 0.00f, 90.00f,},
+        {68.00f, 0.00f, 18.00f, 68.00f, -26.40f, 18.00f, 68.00f, -26.40f,
+            51.00f, 68.00f, 0.00f, 51.00f, 124.00f, 0.00f, 27.00f,
+            124.00f, -26.40f, 27.00f, 104.00f, -26.40f, 51.00f,
+            104.00f, 0.00f, 51.00f, 96.00f, 0.00f, 75.00f, 96.00f,
+            -10.00f, 75.00f, 92.00f, -10.00f, 78.00f, 92.00f, 0.00f,
+            78.00f, 132.00f, 0.00f, 90.00f, 132.00f, -10.00f, 90.00f,
+            108.00f, -10.00f, 90.00f, 108.00f, 0.00f, 90.00f,},
+        {108.00f, 0.00f, 90.00f, 108.00f, 10.00f, 90.00f, 132.00f, 10.00f,
+            90.00f, 132.00f, 0.00f, 90.00f, 112.00f, 0.00f, 93.00f,
+            112.00f, 10.00f, 93.00f, 141.00f, 10.00f, 93.75f, 141.00f,
+            0.00f, 93.75f, 116.00f, 0.00f, 93.00f, 116.00f, 6.00f,
+            93.00f, 138.00f, 6.00f, 94.50f, 138.00f, 0.00f, 94.50f,
+            112.00f, 0.00f, 90.00f, 112.00f, 6.00f, 90.00f, 128.00f,
+            6.00f, 90.00f, 128.00f, 0.00f, 90.00f,},
+        {132.00f, 0.00f, 90.00f, 132.00f, -10.00f, 90.00f, 108.00f,
+            -10.00f, 90.00f, 108.00f, 0.00f, 90.00f, 141.00f, 0.00f,
+            93.75f, 141.00f, -10.00f, 93.75f, 112.00f, -10.00f, 93.00f,
+            112.00f, 0.00f, 93.00f, 138.00f, 0.00f, 94.50f, 138.00f,
+            -6.00f, 94.50f, 116.00f, -6.00f, 93.00f, 116.00f, 0.00f,
+            93.00f, 128.00f, 0.00f, 90.00f, 128.00f, -6.00f, 90.00f,
+            112.00f, -6.00f, 90.00f, 112.00f, 0.00f, 90.00f,},
+        {50.00f, 0.00f, 90.00f, 50.00f, 28.00f, 90.00f, 28.00f, 50.00f,
+            90.00f, 0.00f, 50.00f, 90.00f, 52.00f, 0.00f, 90.00f,
+            52.00f, 29.12f, 90.00f, 29.12f, 52.00f, 90.00f, 0.00f,
+            52.00f, 90.00f, 54.00f, 0.00f, 90.00f, 54.00f, 30.24f,
+            90.00f, 30.24f, 54.00f, 90.00f, 0.00f, 54.00f, 90.00f,
+            56.00f, 0.00f, 90.00f, 56.00f, 31.36f, 90.00f, 31.36f,
+            56.00f, 90.00f, 0.00f, 56.00f, 90.00f,},
+        {0.00f, 50.00f, 90.00f, -28.00f, 50.00f, 90.00f, -50.00f, 28.00f,
+            90.00f, -50.00f, 0.00f, 90.00f, 0.00f, 52.00f, 90.00f,
+            -29.12f, 52.00f, 90.00f, -52.00f, 29.12f, 90.00f, -52.00f,
+            0.00f, 90.00f, 0.00f, 54.00f, 90.00f, -30.24f, 54.00f,
+            90.00f, -54.00f, 30.24f, 90.00f, -54.00f, 0.00f, 90.00f,
+            0.00f, 56.00f, 90.00f, -31.36f, 56.00f, 90.00f, -56.00f,
+            31.36f, 90.00f, -56.00f, 0.00f, 90.00f,},
+        {-50.00f, 0.00f, 90.00f, -50.00f, -28.00f, 90.00f, -28.00f,
+            -50.00f, 90.00f, 0.00f, -50.00f, 90.00f, -52.00f, 0.00f,
+            90.00f, -52.00f, -29.12f, 90.00f, -29.12f, -52.00f, 90.00f,
+            0.00f, -52.00f, 90.00f, -54.00f, 0.00f, 90.00f, -54.00f,
+            -30.24f, 90.00f, -30.24f, -54.00f, 90.00f, 0.00f, -54.00f,
+            90.00f, -56.00f, 0.00f, 90.00f, -56.00f, -31.36f, 90.00f,
+            -31.36f, -56.00f, 90.00f, 0.00f, -56.00f, 90.00f,},
+        {0.00f, -50.00f, 90.00f, 28.00f, -50.00f, 90.00f, 50.00f, -28.00f,
+            90.00f, 50.00f, 0.00f, 90.00f, 0.00f, -52.00f, 90.00f,
+            29.12f, -52.00f, 90.00f, 52.00f, -29.12f, 90.00f, 52.00f,
+            0.00f, 90.00f, 0.00f, -54.00f, 90.00f, 30.24f, -54.00f,
+            90.00f, 54.00f, -30.24f, 90.00f, 54.00f, 0.00f, 90.00f,
+            0.00f, -56.00f, 90.00f, 31.36f, -56.00f, 90.00f, 56.00f,
+            -31.36f, 90.00f, 56.00f, 0.00f, 90.00f,},
+        {8.00f, 0.00f, 102.00f, 8.00f, 4.48f, 102.00f, 4.48f, 8.00f,
+            102.00f, 0.00f, 8.00f, 102.00f, 16.00f, 0.00f, 96.00f,
+            16.00f, 8.96f, 96.00f, 8.96f, 16.00f, 96.00f, 0.00f,
+            16.00f, 96.00f, 52.00f, 0.00f, 96.00f, 52.00f, 29.12f,
+            96.00f, 29.12f, 52.00f, 96.00f, 0.00f, 52.00f, 96.00f,
+            52.00f, 0.00f, 90.00f, 52.00f, 29.12f, 90.00f, 29.12f,
+            52.00f, 90.00f, 0.00f, 52.00f, 90.00f,},
+        {0.00f, 8.00f, 102.00f, -4.48f, 8.00f, 102.00f, -8.00f, 4.48f,
+            102.00f, -8.00f, 0.00f, 102.00f, 0.00f, 16.00f, 96.00f,
+            -8.96f, 16.00f, 96.00f, -16.00f, 8.96f, 96.00f, -16.00f,
+            0.00f, 96.00f, 0.00f, 52.00f, 96.00f, -29.12f, 52.00f,
+            96.00f, -52.00f, 29.12f, 96.00f, -52.00f, 0.00f, 96.00f,
+            0.00f, 52.00f, 90.00f, -29.12f, 52.00f, 90.00f, -52.00f,
+            29.12f, 90.00f, -52.00f, 0.00f, 90.00f,},
+        {-8.00f, 0.00f, 102.00f, -8.00f, -4.48f, 102.00f, -4.48f, -8.00f,
+            102.00f, 0.00f, -8.00f, 102.00f, -16.00f, 0.00f, 96.00f,
+            -16.00f, -8.96f, 96.00f, -8.96f, -16.00f, 96.00f, 0.00f,
+            -16.00f, 96.00f, -52.00f, 0.00f, 96.00f, -52.00f, -29.12f,
+            96.00f, -29.12f, -52.00f, 96.00f, 0.00f, -52.00f, 96.00f,
+            -52.00f, 0.00f, 90.00f, -52.00f, -29.12f, 90.00f, -29.12f,
+            -52.00f, 90.00f, 0.00f, -52.00f, 90.00f,},
+        {0.00f, -8.00f, 102.00f, 4.48f, -8.00f, 102.00f, 8.00f, -4.48f,
+            102.00f, 8.00f, 0.00f, 102.00f, 0.00f, -16.00f, 96.00f,
+            8.96f, -16.00f, 96.00f, 16.00f, -8.96f, 96.00f, 16.00f,
+            0.00f, 96.00f, 0.00f, -52.00f, 96.00f, 29.12f, -52.00f,
+            96.00f, 52.00f, -29.12f, 96.00f, 52.00f, 0.00f, 96.00f,
+            0.00f, -52.00f, 90.00f, 29.12f, -52.00f, 90.00f, 52.00f,
+            -29.12f, 90.00f, 52.00f, 0.00f, 90.00f,},
+        {0.00f, 0.00f, 120.00f, 0.00f, 0.00f, 120.00f, 0.00f, 0.00f,
+            120.00f, 0.00f, 0.00f, 120.00f, 32.00f, 0.00f, 120.00f,
+            32.00f, 18.00f, 120.00f, 18.00f, 32.00f, 120.00f, 0.00f,
+            32.00f, 120.00f, 0.00f, 0.00f, 108.00f, 0.00f, 0.00f,
+            108.00f, 0.00f, 0.00f, 108.00f, 0.00f, 0.00f, 108.00f,
+            8.00f, 0.00f, 102.00f, 8.00f, 4.48f, 102.00f, 4.48f, 8.00f,
+            102.00f, 0.00f, 8.00f, 102.00f,},
+        {0.00f, 0.00f, 120.00f, 0.00f, 0.00f, 120.00f, 0.00f, 0.00f,
+            120.00f, 0.00f, 0.00f, 120.00f, 0.00f, 32.00f, 120.00f,
+            -18.00f, 32.00f, 120.00f, -32.00f, 18.00f, 120.00f,
+            -32.00f, 0.00f, 120.00f, 0.00f, 0.00f, 108.00f, 0.00f,
+            0.00f, 108.00f, 0.00f, 0.00f, 108.00f, 0.00f, 0.00f,
+            108.00f, 0.00f, 8.00f, 102.00f, -4.48f, 8.00f, 102.00f,
+            -8.00f, 4.48f, 102.00f, -8.00f, 0.00f, 102.00f,},
+        {0.00f, 0.00f, 120.00f, 0.00f, 0.00f, 120.00f, 0.00f, 0.00f,
+            120.00f, 0.00f, 0.00f, 120.00f, -32.00f, 0.00f, 120.00f,
+            -32.00f, -18.00f, 120.00f, -18.00f, -32.00f, 120.00f,
+            0.00f, -32.00f, 120.00f, 0.00f, 0.00f, 108.00f, 0.00f,
+            0.00f, 108.00f, 0.00f, 0.00f, 108.00f, 0.00f, 0.00f,
+            108.00f, -8.00f, 0.00f, 102.00f, -8.00f, -4.48f, 102.00f,
+            -4.48f, -8.00f, 102.00f, 0.00f, -8.00f, 102.00f,},
+        {0.00f, 0.00f, 120.00f, 0.00f, 0.00f, 120.00f, 0.00f, 0.00f,
+            120.00f, 0.00f, 0.00f, 120.00f, 0.00f, -32.00f, 120.00f,
+            18.00f, -32.00f, 120.00f, 32.00f, -18.00f, 120.00f, 32.00f,
+            0.00f, 120.00f, 0.00f, 0.00f, 108.00f, 0.00f, 0.00f,
+            108.00f, 0.00f, 0.00f, 108.00f, 0.00f, 0.00f, 108.00f,
+            0.00f, -8.00f, 102.00f, 4.48f, -8.00f, 102.00f, 8.00f,
+            -4.48f, 102.00f, 8.00f, 0.00f, 102.00f,}};
+
+    public Teapot() {
+        super(PATCHES);
+    }
+}

data/src/main/java/org/sunflow/image/Bitmap.java

@@ -0,0 +1,15 @@
+package org.sunflow.image;
+
+public abstract class Bitmap {
+
+    protected static final float INV255 = 1.0f / 255;
+    protected static final float INV65535 = 1.0f / 65535;
+
+    public abstract int getWidth();
+
+    public abstract int getHeight();
+
+    public abstract Color readColor(int x, int y);
+
+    public abstract float readAlpha(int x, int y);
+}

data/src/main/java/org/sunflow/image/BitmapReader.java

@@ -0,0 +1,39 @@
+package org.sunflow.image;
+
+import java.io.IOException;
+
+/**
+ * This is a very simple interface, designed to handle loading of bitmap data.
+ */
+public interface BitmapReader {
+
+    /**
+     * Load the specified filename. This method should throw exception if it
+     * encounters any errors. If the file is valid but its contents are not
+     * (invalid header for example), a {@link BitmapFormatException} may be
+     * thrown. It is an error for this method to return
+     * <code>null</code>.
+     *
+     * @param filename image filename to load
+     * @param isLinear if this is <code>true</code>, the bitmap is assumed to be
+     * already in linear space. This can be usefull when reading greyscale
+     * images for bump mapping for example. HDR formats can ignore this flag
+     * since they usually always store data in linear form.
+     * @return a new {@link Bitmap} object
+     * @throws java.io.IOException
+     * @throws org.sunflow.image.BitmapReader.BitmapFormatException
+     */
+    public Bitmap load(String filename, boolean isLinear) throws IOException, BitmapFormatException;
+
+    /**
+     * This exception can be used internally by bitmap readers to signal they
+     * have encountered a valid file but which contains invalid content.
+     */
+    @SuppressWarnings("serial")
+    public static final class BitmapFormatException extends Exception {
+
+        public BitmapFormatException(String message) {
+            super(message);
+        }
+    }
+}

data/src/main/java/org/sunflow/image/BitmapWriter.java

@@ -0,0 +1,79 @@
+package org.sunflow.image;
+
+import java.io.IOException;
+
+/**
+ * This interface is used to represents an image output format. The methods are
+ * tile oriented so that tiled image formats may be optimally supported. Note
+ * that if the header is declared with a 0 tile size, the image will not be
+ * written with identical sized tiles. The image should either be buffered so it
+ * can be written all at once on close, or an eror should be thrown. The bitmap
+ * writer should be designed so that it is thread safe. Specifically, this means
+ * that the tile writing method can be called by several threads.
+ */
+public interface BitmapWriter {
+
+    /**
+     * This method will be called before writing begins. It is used to set
+     * common attributes to file writers. Currently supported keywords include:
+     * <ul>
+     * <li>"compression"</li>
+     * <li>"channeltype": "byte", "short", "half", "float"</li>
+     * </ul>
+     * Note that this method should not fail if its input is not supported or
+     * invalid. It should gracefully ignore the error and keep its default
+     * state.
+     *
+     * @param option
+     * @param value
+     */
+    public abstract void configure(String option, String value);
+
+    /**
+     * Open a handle to the specified file for writing. If the writer buffers
+     * the image and writes it on close, then the filename should be stored.
+     *
+     * @param filename filename to write the bitmap to
+     * @throws IOException thrown if an I/O error occurs
+     */
+    public abstract void openFile(String filename) throws IOException;
+
+    /**
+     * Write the bitmap header. This may be defered if the image is buffered for
+     * writing all at once on close. Note that if tile size is positive, data
+     * sent to this class is guarenteed to arrive in tiles of that size (except
+     * at borders). Otherwise, it should be assumed that the data is random, and
+     * that it may overlap. The writer should then either throw an error or
+     * start buffering data manually.
+     *
+     * @param width image width
+     * @param height image height
+     * @param tileSize tile size or 0 if the image will not be sent in tiled
+     * form
+     * @throws IOException thrown if an I/O error occurs
+     * @throws UnsupportedOperationException thrown if this writer does not
+     * support writing the image with the supplied tile size
+     */
+    public abstract void writeHeader(int width, int height, int tileSize) throws IOException, UnsupportedOperationException;
+
+    /**
+     * Write a tile of data. Note that this method may be called by more than
+     * one thread, so it should be made thread-safe if possible.
+     *
+     * @param x tile x coordinate
+     * @param y tile y coordinate
+     * @param w tile width
+     * @param h tile height
+     * @param color color data
+     * @param alpha alpha data
+     * @throws IOException thrown if an I/O error occurs
+     */
+    public abstract void writeTile(int x, int y, int w, int h, Color[] color, float[] alpha) throws IOException;
+
+    /**
+     * Close the file, this completes the bitmap writing process.
+     *
+     * @throws IOException thrown if an I/O error occurs
+     */
+    public abstract void closeFile() throws IOException;
+}

data/src/main/java/org/sunflow/image/BlackbodySpectrum.java

@@ -0,0 +1,16 @@
+package org.sunflow.image;
+
+public class BlackbodySpectrum extends SpectralCurve {
+
+    private final float temp;
+
+    public BlackbodySpectrum(float temp) {
+        this.temp = temp;
+    }
+
+    @Override
+    public float sample(float lambda) {
+        double wavelength = lambda * 1e-9;
+        return (float) ((3.74183e-16 * Math.pow(wavelength, -5.0)) / (Math.exp(1.4388e-2 / (wavelength * temp)) - 1.0));
+    }
+}

data/src/main/java/org/sunflow/image/ChromaticitySpectrum.java

@@ -0,0 +1,55 @@
+package org.sunflow.image;
+
+/**
+ * This spectral curve represents a given (x,y) chromaticity pair as explained
+ * in the sun/sky paper (section A.5)
+ */
+public final class ChromaticitySpectrum extends SpectralCurve {
+
+    private static final float[] S0_AMPLITUDES = {0.04f, 6.0f, 29.6f, 55.3f,
+        57.3f, 61.8f, 61.5f, 68.8f, 63.4f, 65.8f, 94.8f, 104.8f, 105.9f,
+        96.8f, 113.9f, 125.6f, 125.5f, 121.3f, 121.3f, 113.5f, 113.1f,
+        110.8f, 106.5f, 108.8f, 105.3f, 104.4f, 100.0f, 96.0f, 95.1f,
+        89.1f, 90.5f, 90.3f, 88.4f, 84.0f, 85.1f, 81.9f, 82.6f, 84.9f,
+        81.3f, 71.9f, 74.3f, 76.4f, 63.3f, 71.7f, 77.0f, 65.2f, 47.7f,
+        68.6f, 65.0f, 66.0f, 61.0f, 53.3f, 58.9f, 61.9f};
+    private static final float[] S1_AMPLITUDES = {0.02f, 4.5f, 22.4f, 42.0f,
+        40.6f, 41.6f, 38.0f, 42.4f, 38.5f, 35.0f, 43.4f, 46.3f, 43.9f,
+        37.1f, 36.7f, 35.9f, 32.6f, 27.9f, 24.3f, 20.1f, 16.2f, 13.2f,
+        8.6f, 6.1f, 4.2f, 1.9f, 0.0f, -1.6f, -3.5f, -3.5f, -5.8f, -7.2f,
+        -8.6f, -9.5f, -10.9f, -10.7f, -12.0f, -14.0f, -13.6f, -12.0f,
+        -13.3f, -12.9f, -10.6f, -11.6f, -12.2f, -10.2f, -7.8f, -11.2f,
+        -10.4f, -10.6f, -9.7f, -8.3f, -9.3f, -9.8f};
+    private static final float[] S2_AMPLITUDES = {0.0f, 2.0f, 4.0f, 8.5f, 7.8f,
+        6.7f, 5.3f, 6.1f, 3.0f, 1.2f, -1.1f, -0.5f, -0.7f, -1.2f, -2.6f,
+        -2.9f, -2.8f, -2.6f, -2.6f, -1.8f, -1.5f, -1.3f, -1.2f, -1.0f,
+        -0.5f, -0.3f, 0.0f, 0.2f, 0.5f, 2.1f, 3.2f, 4.1f, 4.7f, 5.1f, 6.7f,
+        7.3f, 8.6f, 9.8f, 10.2f, 8.3f, 9.6f, 8.5f, 7.0f, 7.6f, 8.0f, 6.7f,
+        5.2f, 7.4f, 6.8f, 7.0f, 6.4f, 5.5f, 6.1f, 6.5f};
+    private static final RegularSpectralCurve KS0_SPECTRUM = new RegularSpectralCurve(S0_AMPLITUDES, 300, 830);
+    private static final RegularSpectralCurve KS1_SPECTRUM = new RegularSpectralCurve(S1_AMPLITUDES, 300, 830);
+    private static final RegularSpectralCurve KS2_SPECTRUM = new RegularSpectralCurve(S2_AMPLITUDES, 300, 830);
+    private static final XYZColor S0_XYZ = KS0_SPECTRUM.toXYZ();
+    private static final XYZColor S1_XYZ = KS1_SPECTRUM.toXYZ();
+    private static final XYZColor S2_XYZ = KS2_SPECTRUM.toXYZ();
+    private final float M1, M2;
+
+    public ChromaticitySpectrum(float x, float y) {
+        M1 = (-1.3515f - 1.7703f * x + 5.9114f * y) / (0.0241f + 0.2562f * x - 0.7341f * y);
+        M2 = (0.03f - 31.4424f * x + 30.0717f * y) / (0.0241f + 0.2562f * x - 0.7341f * y);
+    }
+
+    @Override
+    public float sample(float lambda) {
+        return KS0_SPECTRUM.sample(lambda) + M1 * KS1_SPECTRUM.sample(lambda) + M2 * KS2_SPECTRUM.sample(lambda);
+    }
+
+    public static final XYZColor get(float x, float y) {
+        float M1 = (-1.3515f - 1.7703f * x + 5.9114f * y) / (0.0241f + 0.2562f * x - 0.7341f * y);
+        float M2 = (0.03f - 31.4424f * x + 30.0717f * y) / (0.0241f + 0.2562f * x - 0.7341f * y);
+        float X = S0_XYZ.getX() + M1 * S1_XYZ.getX() + M2 * S2_XYZ.getX();
+        float Y = S0_XYZ.getY() + M1 * S1_XYZ.getY() + M2 * S2_XYZ.getY();
+        float Z = S0_XYZ.getZ() + M1 * S1_XYZ.getZ() + M2 * S2_XYZ.getZ();
+        return new XYZColor(X, Y, Z);
+    }
+}

data/src/main/java/org/sunflow/image/Color.java

@@ -0,0 +1,374 @@
+package org.sunflow.image;
+
+import org.sunflow.math.MathUtils;
+
+public final class Color {
+
+    private float r, g, b;
+    public static final RGBSpace NATIVE_SPACE = RGBSpace.SRGB;
+    public static final Color BLACK = new Color(0, 0, 0);
+    public static final Color WHITE = new Color(1, 1, 1);
+    public static final Color RED = new Color(1, 0, 0);
+    public static final Color GREEN = new Color(0, 1, 0);
+    public static final Color BLUE = new Color(0, 0, 1);
+    public static final Color YELLOW = new Color(1, 1, 0);
+    public static final Color CYAN = new Color(0, 1, 1);
+    public static final Color MAGENTA = new Color(1, 0, 1);
+    public static final Color GRAY = new Color(0.5f, 0.5f, 0.5f);
+
+    public static Color black() {
+        return new Color();
+    }
+
+    public static Color white() {
+        return new Color(1, 1, 1);
+    }
+    private static final float[] EXPONENT = new float[256];
+
+    static {
+        EXPONENT[0] = 0;
+        for (int i = 1; i < 256; i++) {
+            float f = 1.0f;
+            int e = i - (128 + 8);
+            if (e > 0) {
+                for (int j = 0; j < e; j++) {
+                    f *= 2.0f;
+                }
+            } else {
+                for (int j = 0; j < -e; j++) {
+                    f *= 0.5f;
+                }
+            }
+            EXPONENT[i] = f;
+        }
+    }
+
+    public Color() {
+    }
+
+    public Color(float gray) {
+        r = g = b = gray;
+    }
+
+    public Color(float r, float g, float b) {
+        this.r = r;
+        this.g = g;
+        this.b = b;
+    }
+
+    public Color toNonLinear() {
+        r = NATIVE_SPACE.gammaCorrect(r);
+        g = NATIVE_SPACE.gammaCorrect(g);
+        b = NATIVE_SPACE.gammaCorrect(b);
+        return this;
+    }
+
+    public Color toLinear() {
+        r = NATIVE_SPACE.ungammaCorrect(r);
+        g = NATIVE_SPACE.ungammaCorrect(g);
+        b = NATIVE_SPACE.ungammaCorrect(b);
+        return this;
+    }
+
+    public Color(Color c) {
+        r = c.r;
+        g = c.g;
+        b = c.b;
+    }
+
+    public Color(int rgb) {
+        r = ((rgb >> 16) & 0xFF) / 255.0f;
+        g = ((rgb >> 8) & 0xFF) / 255.0f;
+        b = (rgb & 0xFF) / 255.0f;
+    }
+
+    public Color copy() {
+        return new Color(this);
+    }
+
+    public final Color set(float r, float g, float b) {
+        this.r = r;
+        this.g = g;
+        this.b = b;
+        return this;
+    }
+
+    public final Color set(Color c) {
+        r = c.r;
+        g = c.g;
+        b = c.b;
+        return this;
+    }
+
+    public final Color setRGB(int rgb) {
+        r = ((rgb >> 16) & 0xFF) / 255.0f;
+        g = ((rgb >> 8) & 0xFF) / 255.0f;
+        b = (rgb & 0xFF) / 255.0f;
+        return this;
+    }
+
+    public final Color setRGBE(int rgbe) {
+        float f = EXPONENT[rgbe & 0xFF];
+        r = f * ((rgbe >>> 24) + 0.5f);
+        g = f * (((rgbe >> 16) & 0xFF) + 0.5f);
+        b = f * (((rgbe >> 8) & 0xFF) + 0.5f);
+        return this;
+    }
+
+    public final boolean isBlack() {
+        return r <= 0 && g <= 0 && b <= 0;
+    }
+
+    public final float getLuminance() {
+        return (0.2989f * r) + (0.5866f * g) + (0.1145f * b);
+    }
+
+    public final float getMin() {
+        return MathUtils.min(r, g, b);
+    }
+
+    public final float getMax() {
+        return MathUtils.max(r, g, b);
+    }
+
+    public final float getAverage() {
+        return (r + g + b) / 3.0f;
+    }
+
+    public final float[] getRGB() {
+        return new float[]{r, g, b};
+    }
+
+    public final int toRGB() {
+        int ir = (int) (r * 255 + 0.5);
+        int ig = (int) (g * 255 + 0.5);
+        int ib = (int) (b * 255 + 0.5);
+        ir = MathUtils.clamp(ir, 0, 255);
+        ig = MathUtils.clamp(ig, 0, 255);
+        ib = MathUtils.clamp(ib, 0, 255);
+        return (ir << 16) | (ig << 8) | ib;
+    }
+
+    public final int toRGBA(float a) {
+        int ir = (int) (r * 255 + 0.5);
+        int ig = (int) (g * 255 + 0.5);
+        int ib = (int) (b * 255 + 0.5);
+        int ia = (int) (a * 255 + 0.5);
+        ir = MathUtils.clamp(ir, 0, 255);
+        ig = MathUtils.clamp(ig, 0, 255);
+        ib = MathUtils.clamp(ib, 0, 255);
+        ia = MathUtils.clamp(ia, 0, 255);
+        return (ia << 24) | (ir << 16) | (ig << 8) | ib;
+    }
+
+    public final int toRGBE() {
+        // encode the color into 32bits while preserving HDR using Ward's RGBE
+        // technique
+        float v = MathUtils.max(r, g, b);
+        if (v < 1e-32f) {
+            return 0;
+        }
+
+        // get mantissa and exponent
+        float m = v;
+        int e = 0;
+        if (v > 1.0f) {
+            while (m > 1.0f) {
+                m *= 0.5f;
+                e++;
+            }
+        } else if (v <= 0.5f) {
+            while (m <= 0.5f) {
+                m *= 2.0f;
+                e--;
+            }
+        }
+        v = (m * 255.0f) / v;
+        int c = (e + 128);
+        c |= ((int) (r * v) << 24);
+        c |= ((int) (g * v) << 16);
+        c |= ((int) (b * v) << 8);
+        return c;
+    }
+
+    public final Color constrainRGB() {
+        // clamp the RGB value to a representable value
+        float w = -MathUtils.min(0, r, g, b);
+        if (w > 0) {
+            r += w;
+            g += w;
+            b += w;
+        }
+        return this;
+    }
+
+    public final boolean isNan() {
+        return Float.isNaN(r) || Float.isNaN(g) || Float.isNaN(b);
+    }
+
+    public final boolean isInf() {
+        return Float.isInfinite(r) || Float.isInfinite(g) || Float.isInfinite(b);
+    }
+
+    public final Color add(Color c) {
+        r += c.r;
+        g += c.g;
+        b += c.b;
+        return this;
+    }
+
+    public static final Color add(Color c1, Color c2) {
+        return Color.add(c1, c2, new Color());
+    }
+
+    public static final Color add(Color c1, Color c2, Color dest) {
+        dest.r = c1.r + c2.r;
+        dest.g = c1.g + c2.g;
+        dest.b = c1.b + c2.b;
+        return dest;
+    }
+
+    public final Color madd(float s, Color c) {
+        r += (s * c.r);
+        g += (s * c.g);
+        b += (s * c.b);
+        return this;
+    }
+
+    public final Color madd(Color s, Color c) {
+        r += s.r * c.r;
+        g += s.g * c.g;
+        b += s.b * c.b;
+        return this;
+    }
+
+    public final Color sub(Color c) {
+        r -= c.r;
+        g -= c.g;
+        b -= c.b;
+        return this;
+    }
+
+    public static final Color sub(Color c1, Color c2) {
+        return Color.sub(c1, c2, new Color());
+    }
+
+    public static final Color sub(Color c1, Color c2, Color dest) {
+        dest.r = c1.r - c2.r;
+        dest.g = c1.g - c2.g;
+        dest.b = c1.b - c2.b;
+        return dest;
+    }
+
+    public final Color mul(Color c) {
+        r *= c.r;
+        g *= c.g;
+        b *= c.b;
+        return this;
+    }
+
+    public static final Color mul(Color c1, Color c2) {
+        return Color.mul(c1, c2, new Color());
+    }
+
+    public static final Color mul(Color c1, Color c2, Color dest) {
+        dest.r = c1.r * c2.r;
+        dest.g = c1.g * c2.g;
+        dest.b = c1.b * c2.b;
+        return dest;
+    }
+
+    public final Color mul(float s) {
+        r *= s;
+        g *= s;
+        b *= s;
+        return this;
+    }
+
+    public static final Color mul(float s, Color c) {
+        return Color.mul(s, c, new Color());
+    }
+
+    public static final Color mul(float s, Color c, Color dest) {
+        dest.r = s * c.r;
+        dest.g = s * c.g;
+        dest.b = s * c.b;
+        return dest;
+    }
+
+    public final Color div(Color c) {
+        r /= c.r;
+        g /= c.g;
+        b /= c.b;
+        return this;
+    }
+
+    public static final Color div(Color c1, Color c2) {
+        return Color.div(c1, c2, new Color());
+    }
+
+    public static final Color div(Color c1, Color c2, Color dest) {
+        dest.r = c1.r / c2.r;
+        dest.g = c1.g / c2.g;
+        dest.b = c1.b / c2.b;
+        return dest;
+    }
+
+    public final Color exp() {
+        r = (float) Math.exp(r);
+        g = (float) Math.exp(g);
+        b = (float) Math.exp(b);
+        return this;
+    }
+
+    public final Color opposite() {
+        r = 1 - r;
+        g = 1 - g;
+        b = 1 - b;
+        return this;
+    }
+
+    public final Color clamp(float min, float max) {
+        r = MathUtils.clamp(r, min, max);
+        g = MathUtils.clamp(g, min, max);
+        b = MathUtils.clamp(b, min, max);
+        return this;
+    }
+
+    public static final Color blend(Color c1, Color c2, float b) {
+        return blend(c1, c2, b, new Color());
+    }
+
+    public static final Color blend(Color c1, Color c2, float b, Color dest) {
+        dest.r = (1.0f - b) * c1.r + b * c2.r;
+        dest.g = (1.0f - b) * c1.g + b * c2.g;
+        dest.b = (1.0f - b) * c1.b + b * c2.b;
+        return dest;
+    }
+
+    public static final Color blend(Color c1, Color c2, Color b) {
+        return blend(c1, c2, b, new Color());
+    }
+
+    public static final Color blend(Color c1, Color c2, Color b, Color dest) {
+        dest.r = (1.0f - b.r) * c1.r + b.r * c2.r;
+        dest.g = (1.0f - b.g) * c1.g + b.g * c2.g;
+        dest.b = (1.0f - b.b) * c1.b + b.b * c2.b;
+        return dest;
+    }
+
+    public static final boolean hasContrast(Color c1, Color c2, float thresh) {
+        if (Math.abs(c1.r - c2.r) / (c1.r + c2.r) > thresh) {
+            return true;
+        }
+        if (Math.abs(c1.g - c2.g) / (c1.g + c2.g) > thresh) {
+            return true;
+        }
+        return (Math.abs(c1.b - c2.b) / (c1.b + c2.b) > thresh);
+    }
+
+    @Override
+    public String toString() {
+        return String.format("(%.3f, %.3f, %.3f)", r, g, b);
+    }
+}