Cycles: initial subsurface multiple scattering support. It's not working as

well as I would like, but it works, just add a subsurface scattering node and
you can use it like any other BSDF.

It is using fully raytraced sampling compatible with progressive rendering
and other more advanced rendering algorithms we might used in the future, and
it uses no extra memory so it's suitable for complex scenes.

Disadvantage is that it can be quite noisy and slow. Two limitations that will
be solved are that it does not work with bump mapping yet, and that the falloff
function used is a simple cubic function, it's not using the real BSSRDF
falloff function yet.

The node has a color input, along with a scattering radius for each RGB color
channel along with an overall scale factor for the radii.

There is also no GPU support yet, will test if I can get that working later.

Node Documentation:
http://wiki.blender.org/index.php/Doc:2.6/Manual/Render/Cycles/Nodes/Shaders#BSSRDF

Implementation notes:
http://wiki.blender.org/index.php/Dev:2.6/Source/Render/Cycles/Subsurface_Scattering

1 files changed, 53 insertions, 0 deletions

diff --git a/intern/cycles/kernel/svm/svm_closure.h b/intern/cycles/kernel/svm/svm_closure.h
index b5bd2b42cb4..a37646beb2e 100644
--- a/intern/cycles/kernel/svm/svm_closure.h
+++ b/intern/cycles/kernel/svm/svm_closure.h
@@ -306,6 +306,59 @@ __device void svm_node_closure_bsdf(KernelGlobals *kg, ShaderData *sd, float *st
 			}
 			break;
 		}
+#ifdef __SUBSURFACE__
+		case CLOSURE_BSSRDF_ID: {
+			ShaderClosure *sc = &sd->closure[sd->num_closure];
+			float3 weight = sc->weight * mix_weight;
+			float sample_weight = fabsf(average(weight));
+
+			if(sample_weight > 1e-5f && sd->num_closure+2 < MAX_CLOSURE) {
+				/* radius * scale */
+				float3 radius = stack_load_float3(stack, data_node.w)*param1;
+				/* index of refraction */
+				float eta = fmaxf(param2, 1.0f + 1e-5f);
+
+				/* create one closure per color channel */
+				if(fabsf(weight.x) > 0.0f) {
+					sc->weight = make_float3(weight.x, 0.0f, 0.0f);
+					sc->sample_weight = sample_weight;
+					sc->data0 = radius.x;
+					sc->data1 = eta;
+					sc->N = N;
+					sd->flag |= bssrdf_setup(sc);
+
+					sd->num_closure++;
+					sc++;
+				}
+
+				if(fabsf(weight.y) > 0.0f) {
+					sc->weight = make_float3(0.0f, weight.y, 0.0f);
+					sc->sample_weight = sample_weight;
+					sc->data0 = radius.y;
+					sc->data1 = eta;
+					sc->N = N;
+					sd->flag |= bssrdf_setup(sc);
+
+					sd->num_closure++;
+					sc++;
+				}
+
+				if(fabsf(weight.z) > 0.0f) {
+					sc->weight = make_float3(0.0f, 0.0f, weight.z);
+					sc->sample_weight = sample_weight;
+					sc->data0 = radius.z;
+					sc->data1 = eta;
+					sc->N = N;
+					sd->flag |= bssrdf_setup(sc);
+
+					sd->num_closure++;
+					sc++;
+				}
+			}
+
+			break;
+		}
+#endif
 		default:
 			break;
 	}