2
* Adapted from Open Shading Language with this license:
4
* Copyright (c) 2009-2010 Sony Pictures Imageworks Inc., et al.
7
* Modifications Copyright 2011, Blender Foundation.
9
* Redistribution and use in source and binary forms, with or without
10
* modification, are permitted provided that the following conditions are
12
* * Redistributions of source code must retain the above copyright
13
* notice, this list of conditions and the following disclaimer.
14
* * Redistributions in binary form must reproduce the above copyright
15
* notice, this list of conditions and the following disclaimer in the
16
* documentation and/or other materials provided with the distribution.
17
* * Neither the name of Sony Pictures Imageworks nor the names of its
18
* contributors may be used to endorse or promote products derived from
19
* this software without specific prior written permission.
20
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
21
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
22
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
23
* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
24
* OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
25
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
26
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
27
* DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
28
* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
29
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
30
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
33
#include <OpenImageIO/fmath.h>
35
#include <OSL/genclosure.h>
37
#include "osl_closures.h"
43
class BSSRDFCubicClosure : public BSSRDFClosure {
50
static inline T pow3 (const T &x) { return x * x * x; }
53
static inline T pow5 (const T &x) { T x2 = x * x; return x2 * x2 * x; }
55
BSSRDFCubicClosure() { }
59
// pre-compute some terms
61
for (int i = 0; i < 3; i++) {
62
m_scale[i] = m_radius[i] > 0 ? 4 / pow5 (m_radius[i]) : 0;
63
m_max_radius = std::max (m_max_radius, m_radius[i]);
67
bool mergeable (const ClosurePrimitive *other) const {
68
const BSSRDFCubicClosure *comp = (const BSSRDFCubicClosure *)other;
69
return m_radius == comp->m_radius && BSSRDFClosure::mergeable(other);
72
size_t memsize () const { return sizeof(*this); }
74
const char *name () const { return "bssrdf_cubic"; }
76
void print_on (std::ostream &out) const
78
out << name() << " ((" << m_radius[0] << ", " << m_radius[1] << ", " << m_radius[2] << "), ("
79
<< m_scale[0] << ", " << m_scale[1] << ", " << m_scale[2] << "))";
82
Color3 eval (float r) const
84
return Color3 ((r < m_radius.x) ? pow3 (m_radius.x - r) * m_scale.x : 0,
85
(r < m_radius.y) ? pow3 (m_radius.y - r) * m_scale.y : 0,
86
(r < m_radius.z) ? pow3 (m_radius.z - r) * m_scale.z : 0);
89
float max_radius() const
97
ClosureParam closure_bssrdf_cubic_params[] = {
98
CLOSURE_COLOR_PARAM (BSSRDFCubicClosure, m_radius),
99
CLOSURE_STRING_KEYPARAM ("label"),
100
CLOSURE_FINISH_PARAM(BSSRDFCubicClosure) };
102
CLOSURE_PREPARE(closure_bssrdf_cubic_prepare, BSSRDFCubicClosure)