summaryrefslogtreecommitdiffstats
path: root/libdimension/raytrace.c
blob: e91dbb24e4dccec44f1bb6a47e21bec1f671c271 (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
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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
/*************************************************************************
 * Copyright (C) 2008 Tavian Barnes <tavianator@gmail.com>               *
 *                                                                       *
 * This file is part of The Dimension Library.                           *
 *                                                                       *
 * The Dimension Library is free software; you can redistribute it and/  *
 * or modify it under the terms of the GNU Lesser General Public License *
 * as published by the Free Software Foundation; either version 3 of the *
 * License, or (at your option) any later version.                       *
 *                                                                       *
 * The Dimension Library is distributed in the hope that it will be      *
 * useful, but WITHOUT ANY WARRANTY; without even the implied warranty   *
 * of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU  *
 * Lesser General Public License for more details.                       *
 *                                                                       *
 * You should have received a copy of the GNU Lesser General Public      *
 * License along with this program.  If not, see                         *
 * <http://www.gnu.org/licenses/>.                                       *
 *************************************************************************/

#include "dimension.h"
#include <unistd.h> /* For sysconf */

typedef struct {
  dmnsn_scene *scene;
  unsigned int i, n;
} dmnsn_raytrace_thread_payload;

static void *dmnsn_raytrace_scene_thread(void *arg);

void
dmnsn_raytrace_scene(dmnsn_scene *scene)
{
  long n = sysconf(_SC_NPROCESSORS_ONLN);
  unsigned int i;
  pthread_t thread;
  dmnsn_raytrace_thread_payload payload;
  dmnsn_array *threads, *payloads;

  threads  = dmnsn_new_array(sizeof(pthread_t));
  payloads = dmnsn_new_array(sizeof(dmnsn_raytrace_thread_payload));

  if (n <= 0) n = 1;

  payload.scene = scene;
  payload.n = n;
  for (i = 0; i < n; ++i) {
    payload.i = i;
    dmnsn_array_push(payloads, &payload);

    pthread_create(&thread, NULL, &dmnsn_raytrace_scene_thread,
                   dmnsn_array_at(payloads, i));
    dmnsn_array_push(threads, &thread);
  }

  for (i = 0; i < n; ++i) {
    dmnsn_array_get(threads, i, &thread);
    pthread_join(thread, NULL);
  }

  dmnsn_delete_array(payloads);
  dmnsn_delete_array(threads);
}

/* Raytrace a scene */
static void *
dmnsn_raytrace_scene_thread(void *arg)
{
  unsigned int i, j, k, l;
  double t, t_temp;
  dmnsn_object *object;
  dmnsn_line ray, ray_trans;
  dmnsn_raytrace_thread_payload *payload = (dmnsn_raytrace_thread_payload *)arg;
  dmnsn_scene *scene = payload->scene;
  dmnsn_array *intersections;
  dmnsn_color color;
  dmnsn_sRGB sRGB;

  /* Iterate through each pixel */
  for (i = 0; i < scene->canvas->x; ++i) {
    for (j = 0; j < scene->canvas->y; ++j) {
      /* Only do the pixels assigned to this thread */
      if ((j*scene->canvas->x + i)%payload->n == payload->i) {
        /* Set the pixel to the background color */
        color = scene->background;
        t = 0.0;

        /* Get the ray corresponding to the (i,j)th pixel */
        ray = (*scene->camera->ray_fn)(scene->camera, scene->canvas, i, j);

        for (k = 0; k < scene->objects->length; ++k) {
          dmnsn_array_get(scene->objects, k, &object);

          /* Transform the ray according to the object */
          ray_trans = dmnsn_matrix_line_mul(object->trans, ray);

          /* Test for intersections with objects */
          intersections = (*object->intersections_fn)(object, ray_trans);
          for (l = 0; l < intersections->length; ++l) {
            dmnsn_array_get(intersections, l, &t_temp);
            if (t_temp < t || t == 0.0) t = t_temp;
          }
          dmnsn_delete_array(intersections);
        }

        if (t != 0.0) {
          sRGB.R = 1.0 - (t - 2.25)/2.25;
          sRGB.G = sRGB.R;
          sRGB.B = sRGB.R;
          color = dmnsn_color_from_sRGB(sRGB);
        }

        dmnsn_set_pixel(scene->canvas, i, j, color);
      }
    }
  }

  return NULL;
}