Cycles: initial subsurface multiple scattering support. It's not working as
[blender-staging.git] / intern / cycles / kernel / kernel_globals.h
1 /*
2  * Copyright 2011, Blender Foundation.
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  */
18
19 /* Constant Globals */
20
21 CCL_NAMESPACE_BEGIN
22
23 /* On the CPU, we pass along the struct KernelGlobals to nearly everywhere in
24  * the kernel, to access constant data. These are all stored as "textures", but
25  * these are really just standard arrays. We can't use actually globals because
26  * multiple renders may be running inside the same process. */
27
28 #ifdef __KERNEL_CPU__
29
30 #ifdef __OSL__
31 struct OSLGlobals;
32 struct OSLThreadData;
33 struct OSLShadingSystem;
34 #endif
35
36 #define MAX_BYTE_IMAGES   512
37 #define MAX_FLOAT_IMAGES  5
38
39 typedef struct KernelGlobals {
40         texture_image_uchar4 texture_byte_images[MAX_BYTE_IMAGES];
41         texture_image_float4 texture_float_images[MAX_FLOAT_IMAGES];
42
43 #define KERNEL_TEX(type, ttype, name) ttype name;
44 #define KERNEL_IMAGE_TEX(type, ttype, name)
45 #include "kernel_textures.h"
46
47         KernelData __data;
48
49 #ifdef __OSL__
50         /* On the CPU, we also have the OSL globals here. Most data structures are shared
51          * with SVM, the difference is in the shaders and object/mesh attributes. */
52         OSLGlobals *osl;
53         OSLShadingSystem *osl_ss;
54         OSLThreadData *osl_tdata;
55 #endif
56
57 } KernelGlobals;
58
59 #endif
60
61 /* For CUDA, constant memory textures must be globals, so we can't put them
62  * into a struct. As a result we don't actually use this struct and use actual
63  * globals and simply pass along a NULL pointer everywhere, which we hope gets
64  * optimized out. */
65
66 #ifdef __KERNEL_CUDA__
67
68 __constant__ KernelData __data;
69 typedef struct KernelGlobals {} KernelGlobals;
70
71 #define KERNEL_TEX(type, ttype, name) ttype name;
72 #define KERNEL_IMAGE_TEX(type, ttype, name) ttype name;
73 #include "kernel_textures.h"
74
75 #endif
76
77 /* OpenCL */
78
79 #ifdef __KERNEL_OPENCL__
80
81 typedef struct KernelGlobals {
82         __constant KernelData *data;
83
84 #define KERNEL_TEX(type, ttype, name) \
85         __global type *name;
86 #include "kernel_textures.h"
87 } KernelGlobals;
88
89 #endif
90
91 /* Interpolated lookup table access */
92
93 __device float lookup_table_read(KernelGlobals *kg, float x, int offset, int size)
94 {
95         x = clamp(x, 0.0f, 1.0f)*(size-1);
96
97         int index = min((int)x, size-1);
98         int nindex = min(index+1, size-1);
99         float t = x - index;
100
101         float data0 = kernel_tex_fetch(__lookup_table, index + offset);
102         if(t == 0.0f)
103                 return data0;
104
105         float data1 = kernel_tex_fetch(__lookup_table, nindex + offset);
106         return (1.0f - t)*data0 + t*data1;
107 }
108
109 __device float lookup_table_read_2D(KernelGlobals *kg, float x, float y, int offset, int xsize, int ysize)
110 {
111         y = clamp(y, 0.0f, 1.0f)*(ysize-1);
112
113         int index = min((int)y, ysize-1);
114         int nindex = min(index+1, ysize-1);
115         float t = y - index;
116
117         float data0 = lookup_table_read(kg, x, offset + xsize*index, xsize);
118         if(t == 0.0f)
119                 return data0;
120
121         float data1 = lookup_table_read(kg, x, offset + xsize*nindex, xsize);
122         return (1.0f - t)*data0 + t*data1;
123 }
124
125 CCL_NAMESPACE_END
126