intern/cycles/kernel/osl/bssrdf.cpp

   1 /*
   2  * Adapted from Open Shading Language with this license:
   3  *
   4  * Copyright (c) 2009-2010 Sony Pictures Imageworks Inc., et al.
   5  * All Rights Reserved.
   6  *
   7  * Modifications Copyright 2011, Blender Foundation.
   8  *
   9  * Redistribution and use in source and binary forms, with or without
  10  * modification, are permitted provided that the following conditions are
  11  * met:
  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.
  31  */
  32
  33 #include <OpenImageIO/fmath.h>
  34
  35 #include <OSL/genclosure.h>
  36
  37 #include "osl_closures.h"
  38
  39 CCL_NAMESPACE_BEGIN
  40
  41 using namespace OSL;
  42
  43 class BSSRDFCubicClosure : public BSSRDFClosure {
  44 public:
  45         Color3 m_radius;
  46         Color3 m_scale;
  47         float m_max_radius;
  48
  49         template <typename T>
  50         static inline T pow3(const T &x) { return x * x * x; }
  51
  52         template <typename T>
  53         static inline T pow5(const T &x) { T x2 = x * x; return x2 * x2 * x; }
  54
  55         BSSRDFCubicClosure() {}
  56
  57         void setup()
  58         {
  59                 // pre-compute some terms
  60                 m_max_radius = 0;
  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]);
  64                 }
  65         }
  66
  67         bool mergeable(const ClosurePrimitive *other) const {
  68                 const BSSRDFCubicClosure *comp = (const BSSRDFCubicClosure *)other;
  69                 return m_radius == comp->m_radius && BSSRDFClosure::mergeable(other);
  70         }
  71
  72         size_t memsize() const { return sizeof(*this); }
  73
  74         const char *name() const { return "bssrdf_cubic"; }
  75
  76         void print_on(std::ostream &out) const
  77         {
  78                 out << name() << " ((" << m_radius[0] << ", " << m_radius[1] << ", " << m_radius[2] << "), ("
  79                     << m_scale[0] << ", " << m_scale[1] << ", " << m_scale[2] << "))";
  80         }
  81
  82         Color3 eval(float r) const
  83         {
  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);
  87         }
  88
  89         float max_radius() const
  90         {
  91                 return m_max_radius;
  92         }
  93 };
  94
  95
  96
  97 ClosureParam closure_bssrdf_cubic_params[] = {
  98         CLOSURE_COLOR_PARAM(BSSRDFCubicClosure, m_radius),
  99         CLOSURE_STRING_KEYPARAM("label"),
 100         CLOSURE_FINISH_PARAM(BSSRDFCubicClosure)
 101 };
 102
 103 CLOSURE_PREPARE(closure_bssrdf_cubic_prepare, BSSRDFCubicClosure)
 104
 105 CCL_NAMESPACE_END
 106