doxygen/html/igi_8cpp-source.html

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
<title>pbrt: igi.cpp Source File</title>
<link href="doxygen.css" rel="stylesheet" type="text/css">
<link href="tabs.css" rel="stylesheet" type="text/css">
</head><body>
<!-- Generated by Doxygen 1.5.6 -->
<div class="navigation" id="top">
  <div class="tabs">
    <ul>
      <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
      <li><a href="annotated.html"><span>Classes</span></a></li>
      <li class="current"><a href="files.html"><span>Files</span></a></li>
    </ul>
  </div>
<h1>igi.cpp</h1><a href="igi_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 
<a name="l00002"></a>00002 <span class="comment">/*</span>
<a name="l00003"></a>00003 <span class="comment">    pbrt source code Copyright(c) 1998-2007 Matt Pharr and Greg Humphreys.</span>
<a name="l00004"></a>00004 <span class="comment"></span>
<a name="l00005"></a>00005 <span class="comment">    This file is part of pbrt.</span>
<a name="l00006"></a>00006 <span class="comment"></span>
<a name="l00007"></a>00007 <span class="comment">    pbrt is free software; you can redistribute it and/or modify</span>
<a name="l00008"></a>00008 <span class="comment">    it under the terms of the GNU General Public License as published by</span>
<a name="l00009"></a>00009 <span class="comment">    the Free Software Foundation; either version 2 of the License, or</span>
<a name="l00010"></a>00010 <span class="comment">    (at your option) any later version.  Note that the text contents of</span>
<a name="l00011"></a>00011 <span class="comment">    the book "Physically Based Rendering" are *not* licensed under the</span>
<a name="l00012"></a>00012 <span class="comment">    GNU GPL.</span>
<a name="l00013"></a>00013 <span class="comment"></span>
<a name="l00014"></a>00014 <span class="comment">    pbrt is distributed in the hope that it will be useful,</span>
<a name="l00015"></a>00015 <span class="comment">    but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
<a name="l00016"></a>00016 <span class="comment">    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the</span>
<a name="l00017"></a>00017 <span class="comment">    GNU General Public License for more details.</span>
<a name="l00018"></a>00018 <span class="comment"></span>
<a name="l00019"></a>00019 <span class="comment">    You should have received a copy of the GNU General Public License</span>
<a name="l00020"></a>00020 <span class="comment">    along with this program.  If not, see &lt;http://www.gnu.org/licenses/&gt;.</span>
<a name="l00021"></a>00021 <span class="comment"></span>
<a name="l00022"></a>00022 <span class="comment"> */</span>
<a name="l00023"></a>00023 
<a name="l00024"></a>00024 <span class="comment">// igi.cpp*</span>
<a name="l00025"></a>00025 <span class="preprocessor">#include "<a class="code" href="pbrt_8h.html">pbrt.h</a>"</span>
<a name="l00026"></a>00026 <span class="preprocessor">#include "<a class="code" href="transport_8h.html">transport.h</a>"</span>
<a name="l00027"></a>00027 <span class="preprocessor">#include "<a class="code" href="scene_8h.html">scene.h</a>"</span>
<a name="l00028"></a>00028 <span class="preprocessor">#include "<a class="code" href="mc_8h.html">mc.h</a>"</span>
<a name="l00029"></a>00029 <span class="preprocessor">#include "<a class="code" href="sampling_8h.html">sampling.h</a>"</span>
<a name="l00030"></a>00030 
<a name="l00031"></a>00031 <span class="comment">// IGI Local Structures</span>
<a name="l00032"></a><a class="code" href="structVirtualLight.html">00032</a> <span class="keyword">struct </span><a class="code" href="structVirtualLight.html">VirtualLight</a> {
<a name="l00033"></a><a class="code" href="structVirtualLight.html#60e8731d5104cbcef3b6c34d6aac229d">00033</a>         <a class="code" href="structVirtualLight.html#60e8731d5104cbcef3b6c34d6aac229d">VirtualLight</a>() { }
<a name="l00034"></a><a class="code" href="structVirtualLight.html#4d39a2d03472d66e1ee5ef56dc97e792">00034</a>         <a class="code" href="structVirtualLight.html#60e8731d5104cbcef3b6c34d6aac229d">VirtualLight</a>(<span class="keyword">const</span> <a class="code" href="classPoint.html">Point</a> &amp;pp, <span class="keyword">const</span> <a class="code" href="classNormal.html">Normal</a> &amp;nn, <span class="keyword">const</span> <a class="code" href="classSpectrum.html">Spectrum</a> &amp;le)
<a name="l00035"></a>00035                 : <a class="code" href="structVirtualLight.html#08758ea29483fb1bbd68f508cba3e941">p</a>(pp), <a class="code" href="structVirtualLight.html#d28596f48afff99d93b9c2f844f6da51">n</a>(nn), <a class="code" href="structVirtualLight.html#ce365dfb2173047408d3252dceb730bd">Le</a>(le) { }
<a name="l00036"></a><a class="code" href="structVirtualLight.html#08758ea29483fb1bbd68f508cba3e941">00036</a>         <a class="code" href="classPoint.html">Point</a> <a class="code" href="structVirtualLight.html#08758ea29483fb1bbd68f508cba3e941">p</a>;
<a name="l00037"></a><a class="code" href="structVirtualLight.html#d28596f48afff99d93b9c2f844f6da51">00037</a>         <a class="code" href="classNormal.html">Normal</a> <a class="code" href="structVirtualLight.html#d28596f48afff99d93b9c2f844f6da51">n</a>;
<a name="l00038"></a><a class="code" href="structVirtualLight.html#ce365dfb2173047408d3252dceb730bd">00038</a>         <a class="code" href="classSpectrum.html">Spectrum</a> <a class="code" href="structVirtualLight.html#ce365dfb2173047408d3252dceb730bd">Le</a>;
<a name="l00039"></a>00039 };
<a name="l00040"></a>00040 
<a name="l00041"></a><a class="code" href="classIGIIntegrator.html">00041</a> <span class="keyword">class </span><a class="code" href="classIGIIntegrator.html">IGIIntegrator</a> : <span class="keyword">public</span> <a class="code" href="classSurfaceIntegrator.html">SurfaceIntegrator</a> {
<a name="l00042"></a>00042 <span class="keyword">public</span>:
<a name="l00043"></a>00043         <span class="comment">// IGIIntegrator Public Methods</span>
<a name="l00044"></a>00044         <a class="code" href="classIGIIntegrator.html#fff5dc7931bb88a71636427a80419571">IGIIntegrator</a>(<span class="keywordtype">int</span> nl, <span class="keywordtype">int</span> ns, <span class="keywordtype">float</span> md, <span class="keywordtype">float</span> rrt, <span class="keywordtype">float</span> is);
<a name="l00045"></a>00045         <a class="code" href="classSpectrum.html">Spectrum</a> <a class="code" href="classIGIIntegrator.html#4eb4c3765e9b190ad61a88c6177f1243">Li</a>(<span class="keyword">const</span> <a class="code" href="classScene.html">Scene</a> *scene, <span class="keyword">const</span> <a class="code" href="classRayDifferential.html">RayDifferential</a> &amp;ray,
<a name="l00046"></a>00046                 <span class="keyword">const</span> <a class="code" href="structSample.html">Sample</a> *sample, <span class="keywordtype">float</span> *alpha) <span class="keyword">const</span>;
<a name="l00047"></a>00047         <span class="keywordtype">void</span> <a class="code" href="classIGIIntegrator.html#bdb96d3633e48f83a3663a29db23e208">RequestSamples</a>(<a class="code" href="structSample.html">Sample</a> *sample, <span class="keyword">const</span> <a class="code" href="classScene.html">Scene</a> *scene);
<a name="l00048"></a>00048         <span class="keywordtype">void</span> <a class="code" href="classIGIIntegrator.html#de14ff50176a55a455f9f2e8ae1f7745">Preprocess</a>(<span class="keyword">const</span> <a class="code" href="classScene.html">Scene</a> *);
<a name="l00049"></a>00049 
<a name="l00050"></a>00050 <span class="keyword">private</span>:
<a name="l00051"></a>00051         <span class="comment">// IGI Private Data</span>
<a name="l00052"></a><a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">00052</a>         <a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> <a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a>, <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>;
<a name="l00053"></a><a class="code" href="classIGIIntegrator.html#246f41da38149de0ad6ca955da445cdb">00053</a>         vector&lt;VirtualLight&gt; *<a class="code" href="classIGIIntegrator.html#246f41da38149de0ad6ca955da445cdb">virtualLights</a>;
<a name="l00054"></a><a class="code" href="classIGIIntegrator.html#177fcdea6b7cf6fd34521a229686af9f">00054</a>         <span class="keyword">mutable</span> <span class="keywordtype">int</span> <a class="code" href="classIGIIntegrator.html#177fcdea6b7cf6fd34521a229686af9f">specularDepth</a>;
<a name="l00055"></a><a class="code" href="classIGIIntegrator.html#aa26af9208abf5180cb38389aeb30228">00055</a>         <span class="keywordtype">int</span> <a class="code" href="classIGIIntegrator.html#aa26af9208abf5180cb38389aeb30228">maxSpecularDepth</a>;
<a name="l00056"></a><a class="code" href="classIGIIntegrator.html#1a530bd523cc3a085556e72552ed2e52">00056</a>         <span class="keywordtype">float</span> <a class="code" href="classIGIIntegrator.html#3cf8c708a42cc273c4b7831e59ce94f2">minDist2</a>, <a class="code" href="classIGIIntegrator.html#1a530bd523cc3a085556e72552ed2e52">rrThreshold</a>, <a class="code" href="classIGIIntegrator.html#64bc2b2651161c53bae80cbcd48b3b7c">indirectScale</a>;
<a name="l00057"></a><a class="code" href="classIGIIntegrator.html#b98e3dd20072c7c484e5e6bef40d8b05">00057</a>         <span class="keywordtype">int</span> <a class="code" href="classIGIIntegrator.html#b98e3dd20072c7c484e5e6bef40d8b05">vlSetOffset</a>;
<a name="l00058"></a>00058 
<a name="l00059"></a><a class="code" href="classIGIIntegrator.html#037e221083d0435bb427618ea0fb8444">00059</a>         <span class="keywordtype">int</span> *<a class="code" href="classIGIIntegrator.html#037e221083d0435bb427618ea0fb8444">lightSampleOffset</a>, <a class="code" href="classIGIIntegrator.html#6061040a08b216cb299a2666e8def6a2">lightNumOffset</a>;
<a name="l00060"></a><a class="code" href="classIGIIntegrator.html#349e0621b3b3635cb218da2db563be9b">00060</a>         <span class="keywordtype">int</span> *<a class="code" href="classIGIIntegrator.html#349e0621b3b3635cb218da2db563be9b">bsdfSampleOffset</a>, *<a class="code" href="classIGIIntegrator.html#7fcb4bd79b08f5de9ee1a0d0b286f8d6">bsdfComponentOffset</a>;
<a name="l00061"></a>00061 };
<a name="l00062"></a>00062 
<a name="l00063"></a>00063 <span class="comment">// IGIIntegrator Implementation</span>
<a name="l00064"></a><a class="code" href="classIGIIntegrator.html#fff5dc7931bb88a71636427a80419571">00064</a> <a class="code" href="classIGIIntegrator.html#fff5dc7931bb88a71636427a80419571">IGIIntegrator::IGIIntegrator</a>(<span class="keywordtype">int</span> nl, <span class="keywordtype">int</span> ns, <span class="keywordtype">float</span> md,
<a name="l00065"></a>00065                 <span class="keywordtype">float</span> rrt, <span class="keywordtype">float</span> is) {
<a name="l00066"></a>00066         <a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a> = <a class="code" href="pbrt_8h.html#8d4117d7c2f2056597da264d3c763ec6">RoundUpPow2</a>((<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a>)nl);
<a name="l00067"></a>00067         <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a> = <a class="code" href="pbrt_8h.html#8d4117d7c2f2056597da264d3c763ec6">RoundUpPow2</a>((<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a>)ns);
<a name="l00068"></a>00068         <a class="code" href="classIGIIntegrator.html#3cf8c708a42cc273c4b7831e59ce94f2">minDist2</a> = md * md;
<a name="l00069"></a>00069         <a class="code" href="classIGIIntegrator.html#1a530bd523cc3a085556e72552ed2e52">rrThreshold</a> = rrt;
<a name="l00070"></a>00070         <a class="code" href="classIGIIntegrator.html#64bc2b2651161c53bae80cbcd48b3b7c">indirectScale</a> = is;
<a name="l00071"></a>00071         <a class="code" href="classIGIIntegrator.html#aa26af9208abf5180cb38389aeb30228">maxSpecularDepth</a> = 5;
<a name="l00072"></a>00072         <a class="code" href="classIGIIntegrator.html#177fcdea6b7cf6fd34521a229686af9f">specularDepth</a> = 0;
<a name="l00073"></a>00073         <a class="code" href="classIGIIntegrator.html#246f41da38149de0ad6ca955da445cdb">virtualLights</a> = <span class="keyword">new</span> vector&lt;VirtualLight&gt;[<a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>];
<a name="l00074"></a>00074 }
<a name="l00075"></a><a class="code" href="classIGIIntegrator.html#bdb96d3633e48f83a3663a29db23e208">00075</a> <span class="keywordtype">void</span> <a class="code" href="classIGIIntegrator.html#bdb96d3633e48f83a3663a29db23e208">IGIIntegrator::RequestSamples</a>(<a class="code" href="structSample.html">Sample</a> *sample,
<a name="l00076"></a>00076                                    <span class="keyword">const</span> <a class="code" href="classScene.html">Scene</a> *scene) {
<a name="l00077"></a>00077         <span class="comment">// Request samples for area light sampling</span>
<a name="l00078"></a>00078         <a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> nLights = scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>.size();
<a name="l00079"></a>00079         <a class="code" href="classIGIIntegrator.html#037e221083d0435bb427618ea0fb8444">lightSampleOffset</a> = <span class="keyword">new</span> <span class="keywordtype">int</span>[nLights];
<a name="l00080"></a>00080         <a class="code" href="classIGIIntegrator.html#349e0621b3b3635cb218da2db563be9b">bsdfSampleOffset</a> = <span class="keyword">new</span> <span class="keywordtype">int</span>[nLights];
<a name="l00081"></a>00081         <a class="code" href="classIGIIntegrator.html#7fcb4bd79b08f5de9ee1a0d0b286f8d6">bsdfComponentOffset</a> = <span class="keyword">new</span> <span class="keywordtype">int</span>[nLights];
<a name="l00082"></a>00082         <span class="keywordflow">for</span> (<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> i = 0; i &lt; nLights; ++i) {
<a name="l00083"></a>00083                 <span class="keyword">const</span> <a class="code" href="classLight.html">Light</a> *light = scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>[i];
<a name="l00084"></a>00084                 <span class="keywordtype">int</span> lightSamples =
<a name="l00085"></a>00085                         scene-&gt;<a class="code" href="classScene.html#69e1a55e5efa7f323049dd06688abb49">sampler</a>-&gt;<a class="code" href="classSampler.html#d2ca20cd81ffb6fc242c3f0cf6c9812f">RoundSize</a>(light-&gt;<a class="code" href="classLight.html#fae2e64a4a14046913c82c175078d8bb">nSamples</a>);
<a name="l00086"></a>00086                 <a class="code" href="classIGIIntegrator.html#037e221083d0435bb427618ea0fb8444">lightSampleOffset</a>[i] = sample-&gt;<a class="code" href="structSample.html#21b7f539ec14f3b5c7bdddf790e1b8a7">Add2D</a>(lightSamples);
<a name="l00087"></a>00087                 <a class="code" href="classIGIIntegrator.html#349e0621b3b3635cb218da2db563be9b">bsdfSampleOffset</a>[i] = sample-&gt;<a class="code" href="structSample.html#21b7f539ec14f3b5c7bdddf790e1b8a7">Add2D</a>(lightSamples);
<a name="l00088"></a>00088                 <a class="code" href="classIGIIntegrator.html#7fcb4bd79b08f5de9ee1a0d0b286f8d6">bsdfComponentOffset</a>[i] = sample-&gt;<a class="code" href="structSample.html#a2b567c5a4160488ba81953b59df59f0">Add1D</a>(lightSamples);
<a name="l00089"></a>00089         }
<a name="l00090"></a>00090         <a class="code" href="classIGIIntegrator.html#6061040a08b216cb299a2666e8def6a2">lightNumOffset</a> = -1;
<a name="l00091"></a>00091         <a class="code" href="classIGIIntegrator.html#b98e3dd20072c7c484e5e6bef40d8b05">vlSetOffset</a> = sample-&gt;<a class="code" href="structSample.html#a2b567c5a4160488ba81953b59df59f0">Add1D</a>(1);
<a name="l00092"></a>00092 }
<a name="l00093"></a><a class="code" href="classIGIIntegrator.html#de14ff50176a55a455f9f2e8ae1f7745">00093</a> <span class="keywordtype">void</span> <a class="code" href="classIGIIntegrator.html#de14ff50176a55a455f9f2e8ae1f7745">IGIIntegrator::Preprocess</a>(<span class="keyword">const</span> <a class="code" href="classScene.html">Scene</a> *scene) {
<a name="l00094"></a>00094         <span class="keywordflow">if</span> (scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>.size() == 0) <span class="keywordflow">return</span>;
<a name="l00095"></a>00095         <span class="comment">// Compute samples for emitted rays from lights</span>
<a name="l00096"></a>00096         <span class="keywordtype">float</span> *lightNum = <span class="keyword">new</span> <span class="keywordtype">float</span>[<a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a> * <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>];
<a name="l00097"></a>00097         <span class="keywordtype">float</span> *lightSamp0 = <span class="keyword">new</span> <span class="keywordtype">float</span>[2 * <a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a> * <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>];
<a name="l00098"></a>00098         <span class="keywordtype">float</span> *lightSamp1 = <span class="keyword">new</span> <span class="keywordtype">float</span>[2 * <a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a> * <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>];
<a name="l00099"></a>00099         <a class="code" href="sampling_8h.html#5f30f3ccf1c6546da244633214a5039c">LDShuffleScrambled1D</a>(<a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a>, <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>, lightNum);
<a name="l00100"></a>00100         <a class="code" href="sampling_8h.html#e40441cee55a8c50ea59ccf08dca8335">LDShuffleScrambled2D</a>(<a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a>, <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>, lightSamp0);
<a name="l00101"></a>00101         <a class="code" href="sampling_8h.html#e40441cee55a8c50ea59ccf08dca8335">LDShuffleScrambled2D</a>(<a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a>, <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>, lightSamp1);
<a name="l00102"></a>00102         <span class="comment">// Precompute information for light sampling densities</span>
<a name="l00103"></a>00103         <span class="keywordtype">int</span> nLights = int(scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>.size());
<a name="l00104"></a>00104         <span class="keywordtype">float</span> *lightPower = (<span class="keywordtype">float</span> *)alloca(nLights * <span class="keyword">sizeof</span>(<span class="keywordtype">float</span>));
<a name="l00105"></a>00105         <span class="keywordtype">float</span> *lightCDF = (<span class="keywordtype">float</span> *)alloca((nLights+1) * <span class="keyword">sizeof</span>(float));
<a name="l00106"></a>00106         <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; nLights; ++i)
<a name="l00107"></a>00107                 lightPower[i] = scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>[i]-&gt;Power(scene).y();
<a name="l00108"></a>00108         <span class="keywordtype">float</span> totalPower;
<a name="l00109"></a>00109         <a class="code" href="mc_8cpp.html#1cb88cade3da8d8ff8474bcfc71cbe81">ComputeStep1dCDF</a>(lightPower, nLights, &amp;totalPower, lightCDF);
<a name="l00110"></a>00110         <span class="keywordflow">for</span> (<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> s = 0; s &lt; <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>; ++s) {
<a name="l00111"></a>00111                 <span class="keywordflow">for</span> (<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> i = 0; i &lt; <a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a>; ++i) {
<a name="l00112"></a>00112                         <span class="comment">// Follow path _i_ from light to create virtual lights</span>
<a name="l00113"></a>00113                         <span class="keywordtype">int</span> sampOffset = s*nLightPaths + i;
<a name="l00114"></a>00114                         <span class="comment">// Choose light source to trace path from</span>
<a name="l00115"></a>00115                         <span class="keywordtype">float</span> lightPdf;
<a name="l00116"></a>00116                         <span class="keywordtype">int</span> lNum = <a class="code" href="pbrt_8h.html#02648670ea7bf29f5782902cd9aa87c1">Floor2Int</a>(<a class="code" href="mc_8cpp.html#22c0752ae44234ec91da37f554d988cb">SampleStep1d</a>(lightPower, lightCDF,
<a name="l00117"></a>00117                                 totalPower, nLights, lightNum[sampOffset], &amp;lightPdf) * nLights);
<a name="l00118"></a>00118                         fprintf(stderr, <span class="stringliteral">"samp %f -&gt; num %d\n"</span>, lightNum[sampOffset], lNum);
<a name="l00119"></a>00119                         <a class="code" href="classLight.html">Light</a> *light = scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>[lNum];
<a name="l00120"></a>00120                         <span class="comment">// Sample ray leaving light source</span>
<a name="l00121"></a>00121                         <a class="code" href="classRayDifferential.html">RayDifferential</a> ray;
<a name="l00122"></a>00122                         <span class="keywordtype">float</span> pdf;
<a name="l00123"></a>00123                         <a class="code" href="classSpectrum.html">Spectrum</a> alpha =
<a name="l00124"></a>00124                                 light-&gt;<a class="code" href="classLight.html#cc1f14dae6074cc9c16b7387020e814c">Sample_L</a>(scene, lightSamp0[2*sampOffset],
<a name="l00125"></a>00125                                                 lightSamp0[2*sampOffset+1],
<a name="l00126"></a>00126                                                 lightSamp1[2*sampOffset],
<a name="l00127"></a>00127                                                 lightSamp1[2*sampOffset+1],
<a name="l00128"></a>00128                                                 &amp;ray, &amp;pdf);
<a name="l00129"></a>00129                         <span class="keywordflow">if</span> (pdf == 0.f || alpha.<a class="code" href="classSpectrum.html#a05b9a07476998192a83fcd015b7619f">Black</a>()) <span class="keywordflow">continue</span>;
<a name="l00130"></a>00130                         alpha /= pdf * lightPdf;
<a name="l00131"></a>00131                         fprintf(stderr, <span class="stringliteral">"initial alpha %f, light # %d\n"</span>, alpha.<a class="code" href="classSpectrum.html#25d4cb8d137b7bb0f85e636ea41461a8">y</a>(), lNum);
<a name="l00132"></a>00132                         <a class="code" href="structIntersection.html">Intersection</a> isect;
<a name="l00133"></a>00133                         <span class="keywordtype">int</span> nIntersections = 0;
<a name="l00134"></a>00134                         <span class="keywordflow">while</span> (scene-&gt;<a class="code" href="classScene.html#36f27fc416cc5b53fe03575e2d15803a">Intersect</a>(ray, &amp;isect) &amp;&amp; !alpha.<a class="code" href="classSpectrum.html#a05b9a07476998192a83fcd015b7619f">Black</a>()) {
<a name="l00135"></a>00135                                 ++nIntersections;
<a name="l00136"></a>00136                                 alpha *= scene-&gt;<a class="code" href="classScene.html#c80ef18b12e54a64a809de677aada2e9">Transmittance</a>(ray);
<a name="l00137"></a>00137                                 <a class="code" href="classVector.html">Vector</a> wo = -ray.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a>;
<a name="l00138"></a>00138                                 <a class="code" href="classBSDF.html">BSDF</a> *bsdf = isect.GetBSDF(ray);
<a name="l00139"></a>00139                                 <span class="comment">// Create virtual light at ray intersection point</span>
<a name="l00140"></a>00140                                 <span class="keyword">static</span> <a class="code" href="classStatsCounter.html">StatsCounter</a> vls(<span class="stringliteral">"IGI Integrator"</span>, <span class="stringliteral">"Virtual Lights Created"</span>); <span class="comment">//NOBOOK</span>
<a name="l00141"></a>00141                                 ++vls; <span class="comment">//NOBOOK</span>
<a name="l00142"></a>00142                                 <a class="code" href="classSpectrum.html">Spectrum</a> Le = alpha * bsdf-&gt;<a class="code" href="classBSDF.html#0aa3132935c2c62d824b3bd45f6939d9">rho</a>(wo) / <a class="code" href="pbrt_8h.html#e71449b1cc6e6250b91f539153a7a0d3">M_PI</a>;
<a name="l00143"></a>00143                                 fprintf(stderr, <span class="stringliteral">"\tmade light with le y %f\n"</span>, Le.<a class="code" href="classSpectrum.html#25d4cb8d137b7bb0f85e636ea41461a8">y</a>());
<a name="l00144"></a>00144                                 <a class="code" href="classIGIIntegrator.html#246f41da38149de0ad6ca955da445cdb">virtualLights</a>[s].push_back(<a class="code" href="structVirtualLight.html">VirtualLight</a>(isect.dg.p, isect.dg.nn, Le));
<a name="l00145"></a>00145                                 <span class="comment">// Sample new ray direction and update weight</span>
<a name="l00146"></a>00146                                 <a class="code" href="classVector.html">Vector</a> wi;
<a name="l00147"></a>00147                                 <span class="keywordtype">float</span> pdf;
<a name="l00148"></a>00148                                 <a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e669730">BxDFType</a> flags;
<a name="l00149"></a>00149                                 <a class="code" href="classSpectrum.html">Spectrum</a> fr = bsdf-&gt;<a class="code" href="classBSDF.html#3498bc3b7484a17abd8429084124a0d4">Sample_f</a>(wo, &amp;wi, <a class="code" href="pbrt_8h.html#4f11994534b888fcf650c04b5d61d9d4">RandomFloat</a>(),
<a name="l00150"></a>00150                                                                  <a class="code" href="pbrt_8h.html#4f11994534b888fcf650c04b5d61d9d4">RandomFloat</a>(), <a class="code" href="pbrt_8h.html#4f11994534b888fcf650c04b5d61d9d4">RandomFloat</a>(),
<a name="l00151"></a>00151                                                                  &amp;pdf, <a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e669730734897d21389ab0737866047d40e7122">BSDF_ALL</a>, &amp;flags);
<a name="l00152"></a>00152                                 <span class="keywordflow">if</span> (fr.<a class="code" href="classSpectrum.html#a05b9a07476998192a83fcd015b7619f">Black</a>() || pdf == 0.f)
<a name="l00153"></a>00153                                         <span class="keywordflow">break</span>;
<a name="l00154"></a>00154                                 <a class="code" href="classSpectrum.html">Spectrum</a> anew = alpha * fr * <a class="code" href="geometry_8h.html#0c07ebf9b1454b90d6646a2f89065bb3">AbsDot</a>(wi, bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#1aa044f6c7b053ca460c10b0cf832c14">nn</a>) / pdf;
<a name="l00155"></a>00155                                 <span class="keywordtype">float</span> r = anew.<a class="code" href="classSpectrum.html#25d4cb8d137b7bb0f85e636ea41461a8">y</a>() / alpha.<a class="code" href="classSpectrum.html#25d4cb8d137b7bb0f85e636ea41461a8">y</a>();
<a name="l00156"></a>00156                                 fprintf(stderr, <span class="stringliteral">"\tr = %f\n"</span>, r);
<a name="l00157"></a>00157                                 <span class="keywordflow">if</span> (<a class="code" href="pbrt_8h.html#4f11994534b888fcf650c04b5d61d9d4">RandomFloat</a>() &gt; r)
<a name="l00158"></a>00158                                         <span class="keywordflow">break</span>;
<a name="l00159"></a>00159                                 alpha = anew / r;
<a name="l00160"></a>00160                                 fprintf(stderr, <span class="stringliteral">"\tnew alpha %f\n"</span>, alpha.<a class="code" href="classSpectrum.html#25d4cb8d137b7bb0f85e636ea41461a8">y</a>());
<a name="l00161"></a>00161                                 ray = <a class="code" href="classRayDifferential.html">RayDifferential</a>(isect.dg.p, wi);
<a name="l00162"></a>00162                         }
<a name="l00163"></a>00163                         <a class="code" href="classBSDF.html#ca0e5890e4ef6d5af745234d9fb4d99e">BSDF::FreeAll</a>();
<a name="l00164"></a>00164                 }
<a name="l00165"></a>00165         }
<a name="l00166"></a>00166         <span class="keyword">delete</span>[] lightNum; <span class="comment">// NOBOOK</span>
<a name="l00167"></a>00167         <span class="keyword">delete</span>[] lightSamp0; <span class="comment">// NOBOOK</span>
<a name="l00168"></a>00168         <span class="keyword">delete</span>[] lightSamp1; <span class="comment">// NOBOOK</span>
<a name="l00169"></a>00169 }
<a name="l00170"></a><a class="code" href="classIGIIntegrator.html#4eb4c3765e9b190ad61a88c6177f1243">00170</a> <a class="code" href="classSpectrum.html">Spectrum</a> <a class="code" href="classIGIIntegrator.html#4eb4c3765e9b190ad61a88c6177f1243">IGIIntegrator::Li</a>(<span class="keyword">const</span> <a class="code" href="classScene.html">Scene</a> *scene,
<a name="l00171"></a>00171                 <span class="keyword">const</span> <a class="code" href="classRayDifferential.html">RayDifferential</a> &amp;ray, <span class="keyword">const</span> <a class="code" href="structSample.html">Sample</a> *sample,
<a name="l00172"></a>00172                    <span class="keywordtype">float</span> *alpha)<span class="keyword"> const </span>{
<a name="l00173"></a>00173         <a class="code" href="classSpectrum.html">Spectrum</a> L(0.);
<a name="l00174"></a>00174         <a class="code" href="structIntersection.html">Intersection</a> isect;
<a name="l00175"></a>00175         <span class="keywordflow">if</span> (scene-&gt;<a class="code" href="classScene.html#36f27fc416cc5b53fe03575e2d15803a">Intersect</a>(ray, &amp;isect)) {
<a name="l00176"></a>00176                 <span class="keywordflow">if</span> (alpha) *alpha = 1.;
<a name="l00177"></a>00177                 <a class="code" href="classVector.html">Vector</a> wo = -ray.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a>;
<a name="l00178"></a>00178                 <span class="comment">// Compute emitted light if ray hit an area light source</span>
<a name="l00179"></a>00179                 L += isect.<a class="code" href="structIntersection.html#a2f623b7b62e37a4e083fb0118d47fa3">Le</a>(wo);
<a name="l00180"></a>00180                 <span class="comment">// Evaluate BSDF at hit point</span>
<a name="l00181"></a>00181                 <a class="code" href="classBSDF.html">BSDF</a> *bsdf = isect.<a class="code" href="structIntersection.html#ff18bb39143b49fa340a3697b9c606b5">GetBSDF</a>(ray);
<a name="l00182"></a>00182                 <span class="keyword">const</span> <a class="code" href="classPoint.html">Point</a> &amp;p = bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#3d5353db1c1a6a93ef6da642d064f46e">p</a>;
<a name="l00183"></a>00183                 <span class="keyword">const</span> <a class="code" href="classNormal.html">Normal</a> &amp;n = bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#1aa044f6c7b053ca460c10b0cf832c14">nn</a>;
<a name="l00184"></a>00184                 L += <a class="code" href="transport_8cpp.html#b9dc243808fac32535dce9d4c4a29699">UniformSampleAllLights</a>(scene, p, n,
<a name="l00185"></a>00185                                             wo, bsdf, sample,
<a name="l00186"></a>00186                                             <a class="code" href="classIGIIntegrator.html#037e221083d0435bb427618ea0fb8444">lightSampleOffset</a>, <a class="code" href="classIGIIntegrator.html#349e0621b3b3635cb218da2db563be9b">bsdfSampleOffset</a>,
<a name="l00187"></a>00187                                             <a class="code" href="classIGIIntegrator.html#7fcb4bd79b08f5de9ee1a0d0b286f8d6">bsdfComponentOffset</a>);
<a name="l00188"></a>00188                 <span class="comment">// Compute indirect illumination with virtual lights</span>
<a name="l00189"></a>00189                 <a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> lSet = min(<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a>(sample-&gt;<a class="code" href="structSample.html#4a1ae5fbae61678e57efaf185a319ceb">oneD</a>[<a class="code" href="classIGIIntegrator.html#b98e3dd20072c7c484e5e6bef40d8b05">vlSetOffset</a>][0] * <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a>),
<a name="l00190"></a>00190                                  nLightSets-1);
<a name="l00191"></a>00191                 <span class="keywordflow">for</span> (<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> i = 0; i &lt; <a class="code" href="classIGIIntegrator.html#246f41da38149de0ad6ca955da445cdb">virtualLights</a>[lSet].size(); ++i) {
<a name="l00192"></a>00192                         <span class="keyword">const</span> <a class="code" href="structVirtualLight.html">VirtualLight</a> &amp;vl = <a class="code" href="classIGIIntegrator.html#246f41da38149de0ad6ca955da445cdb">virtualLights</a>[lSet][i];
<a name="l00193"></a>00193                         <span class="comment">// Add contribution from _VirtualLight_ _vl_</span>
<a name="l00194"></a>00194                         <span class="comment">// Ignore light if it's too close to current point</span>
<a name="l00195"></a>00195                         <span class="keywordtype">float</span> d2 = <a class="code" href="geometry_8h.html#7e7a8e6ca4a481e5f0d1c44312a956d0">DistanceSquared</a>(p, vl.<a class="code" href="structVirtualLight.html#08758ea29483fb1bbd68f508cba3e941">p</a>);
<a name="l00196"></a>00196                         <span class="comment">//if (d2 &lt; .8 * minDist2) continue;</span>
<a name="l00197"></a>00197                         <span class="keywordtype">float</span> distScale = <a class="code" href="pbrt_8h.html#332172a54ac2930e10f7a5c446cb712d">SmoothStep</a>(.8 * <a class="code" href="classIGIIntegrator.html#3cf8c708a42cc273c4b7831e59ce94f2">minDist2</a>, 1.2 * <a class="code" href="classIGIIntegrator.html#3cf8c708a42cc273c4b7831e59ce94f2">minDist2</a>, d2);
<a name="l00198"></a>00198                         <span class="comment">// Compute virtual light's tentative contribution _Llight_</span>
<a name="l00199"></a>00199                         <a class="code" href="classVector.html">Vector</a> wi = <a class="code" href="geometry_8h.html#46448bfa3df1c2de5e3b763f163866a4">Normalize</a>(vl.<a class="code" href="structVirtualLight.html#08758ea29483fb1bbd68f508cba3e941">p</a> - p);
<a name="l00200"></a>00200                         <a class="code" href="classSpectrum.html">Spectrum</a> f = distScale * bsdf-&gt;<a class="code" href="classBSDF.html#137c0efa679347a2787b366effdeeb50">f</a>(wo, wi);
<a name="l00201"></a>00201                         <span class="keywordflow">if</span> (f.<a class="code" href="classSpectrum.html#a05b9a07476998192a83fcd015b7619f">Black</a>()) <span class="keywordflow">continue</span>;
<a name="l00202"></a>00202                         <span class="keywordtype">float</span> G = <a class="code" href="geometry_8h.html#0c07ebf9b1454b90d6646a2f89065bb3">AbsDot</a>(wi, n) * <a class="code" href="geometry_8h.html#0c07ebf9b1454b90d6646a2f89065bb3">AbsDot</a>(wi, vl.<a class="code" href="structVirtualLight.html#d28596f48afff99d93b9c2f844f6da51">n</a>) / d2;
<a name="l00203"></a>00203                         <a class="code" href="classSpectrum.html">Spectrum</a> Llight = <a class="code" href="classIGIIntegrator.html#64bc2b2651161c53bae80cbcd48b3b7c">indirectScale</a> * f * G * vl.<a class="code" href="structVirtualLight.html#ce365dfb2173047408d3252dceb730bd">Le</a> /
<a name="l00204"></a>00204                                 <a class="code" href="classIGIIntegrator.html#246f41da38149de0ad6ca955da445cdb">virtualLights</a>[lSet].size();
<a name="l00205"></a>00205                         Llight *= scene-&gt;<a class="code" href="classScene.html#c80ef18b12e54a64a809de677aada2e9">Transmittance</a>(<a class="code" href="classRay.html">Ray</a>(p, vl.<a class="code" href="structVirtualLight.html#08758ea29483fb1bbd68f508cba3e941">p</a> - p));
<a name="l00206"></a>00206                         <span class="comment">// Possibly skip shadow ray with Russian roulette</span>
<a name="l00207"></a>00207                         <span class="keywordflow">if</span> (Llight.<a class="code" href="classSpectrum.html#25d4cb8d137b7bb0f85e636ea41461a8">y</a>() &lt; <a class="code" href="classIGIIntegrator.html#1a530bd523cc3a085556e72552ed2e52">rrThreshold</a>) {
<a name="l00208"></a>00208                                 <span class="keywordtype">float</span> continueProbability = .1f;
<a name="l00209"></a>00209                                 <span class="keywordflow">if</span> (<a class="code" href="pbrt_8h.html#4f11994534b888fcf650c04b5d61d9d4">RandomFloat</a>() &gt; continueProbability)
<a name="l00210"></a>00210                                         <span class="keywordflow">continue</span>;
<a name="l00211"></a>00211                                 Llight /= continueProbability;
<a name="l00212"></a>00212                         }
<a name="l00213"></a>00213                         <span class="keyword">static</span> <a class="code" href="classStatsCounter.html">StatsCounter</a> vlsr(<span class="stringliteral">"IGI Integrator"</span>, <span class="stringliteral">"Shadow Rays to Virtual Lights"</span>); <span class="comment">//NOBOOK</span>
<a name="l00214"></a>00214                         ++vlsr; <span class="comment">//NOBOOK</span>
<a name="l00215"></a>00215                         <span class="keywordflow">if</span> (!scene-&gt;<a class="code" href="classScene.html#558b5f74407c7479517ef82b0fe3f913">IntersectP</a>(<a class="code" href="classRay.html">Ray</a>(p, vl.<a class="code" href="structVirtualLight.html#08758ea29483fb1bbd68f508cba3e941">p</a> - p, <a class="code" href="pbrt_8h.html#7507e14449e57d80d946ee4509de1b3b">RAY_EPSILON</a>,
<a name="l00216"></a>00216                                         1.f - <a class="code" href="pbrt_8h.html#7507e14449e57d80d946ee4509de1b3b">RAY_EPSILON</a>)))
<a name="l00217"></a>00217                                 L += Llight;
<a name="l00218"></a>00218                 }
<a name="l00219"></a>00219                 <span class="comment">// Trace rays for specular reflection and refraction</span>
<a name="l00220"></a>00220                 <span class="keywordflow">if</span> (<a class="code" href="classIGIIntegrator.html#177fcdea6b7cf6fd34521a229686af9f">specularDepth</a>++ &lt; <a class="code" href="classIGIIntegrator.html#aa26af9208abf5180cb38389aeb30228">maxSpecularDepth</a>) {
<a name="l00221"></a>00221                         <a class="code" href="classVector.html">Vector</a> wi;
<a name="l00222"></a>00222                         <span class="comment">// Trace rays for specular reflection and refraction</span>
<a name="l00223"></a>00223                         <a class="code" href="classSpectrum.html">Spectrum</a> f = bsdf-&gt;<a class="code" href="classBSDF.html#3498bc3b7484a17abd8429084124a0d4">Sample_f</a>(wo, &amp;wi,
<a name="l00224"></a>00224                                                 <a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e669730">BxDFType</a>(<a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e6697304331da0ab7ddceadcfffa81ecf70577e">BSDF_REFLECTION</a> | <a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e669730c94ee1d195ff19883a79d1c08bd773d9">BSDF_SPECULAR</a>));
<a name="l00225"></a>00225                         <span class="keywordflow">if</span> (!f.<a class="code" href="classSpectrum.html#a05b9a07476998192a83fcd015b7619f">Black</a>()) {
<a name="l00226"></a>00226                                 <span class="comment">// Compute ray differential _rd_ for specular reflection</span>
<a name="l00227"></a>00227                                 <a class="code" href="classRayDifferential.html">RayDifferential</a> rd(p, wi);
<a name="l00228"></a>00228                                 rd.<a class="code" href="classRayDifferential.html#87b02390087bc229067e219becef4889">hasDifferentials</a> = <span class="keyword">true</span>;
<a name="l00229"></a>00229                                 rd.<a class="code" href="classRayDifferential.html#15359147d3814d4e592f005e1ae7ef64">rx</a>.<a class="code" href="classRay.html#1014d000992e8ac3b92f150ad2a5cffe">o</a> = p + isect.<a class="code" href="structIntersection.html#27a04b9f68291416c4d864f105a7d991">dg</a>.<a class="code" href="structDifferentialGeometry.html#307387af96b1b7a98358ea448114d697">dpdx</a>;
<a name="l00230"></a>00230                                 rd.<a class="code" href="classRayDifferential.html#a6ab98cac3bffedb13f3c74f0e13f3b4">ry</a>.<a class="code" href="classRay.html#1014d000992e8ac3b92f150ad2a5cffe">o</a> = p + isect.<a class="code" href="structIntersection.html#27a04b9f68291416c4d864f105a7d991">dg</a>.<a class="code" href="structDifferentialGeometry.html#d3e6794a920835d56713b3659437fc00">dpdy</a>;
<a name="l00231"></a>00231                                 <span class="comment">// Compute differential reflected directions</span>
<a name="l00232"></a>00232                                 <a class="code" href="classNormal.html">Normal</a> dndx = bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#f9daf909200f1e446924641fe27029b0">dndu</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#ee263cfa99de799c049e33b552ff0b82">dudx</a> +
<a name="l00233"></a>00233                                         bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#cbcc594e0fff4bc1bf44d82e42cc6390">dndv</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#f0d770bd2370babb99b4eda972a757e2">dvdx</a>;
<a name="l00234"></a>00234                                 <a class="code" href="classNormal.html">Normal</a> dndy = bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#f9daf909200f1e446924641fe27029b0">dndu</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#57f8e45c2aa405802f25f1921ba27edd">dudy</a> +
<a name="l00235"></a>00235                                         bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#cbcc594e0fff4bc1bf44d82e42cc6390">dndv</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#7343b4a399c5001c92fad44696060142">dvdy</a>;
<a name="l00236"></a>00236                                 <a class="code" href="classVector.html">Vector</a> dwodx = -ray.<a class="code" href="classRayDifferential.html#15359147d3814d4e592f005e1ae7ef64">rx</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> - wo, dwody = -ray.<a class="code" href="classRayDifferential.html#a6ab98cac3bffedb13f3c74f0e13f3b4">ry</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> - wo;
<a name="l00237"></a>00237                                 <span class="keywordtype">float</span> dDNdx = <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(dwodx, n) + <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wo, dndx);
<a name="l00238"></a>00238                                 <span class="keywordtype">float</span> dDNdy = <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(dwody, n) + <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wo, dndy);
<a name="l00239"></a>00239                                 rd.<a class="code" href="classRayDifferential.html#15359147d3814d4e592f005e1ae7ef64">rx</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> = wi -
<a name="l00240"></a>00240                                         dwodx + 2 * <a class="code" href="classVector.html">Vector</a>(<a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wo, n) * dndx +
<a name="l00241"></a>00241                                                  dDNdx * n);
<a name="l00242"></a>00242                                 rd.<a class="code" href="classRayDifferential.html#a6ab98cac3bffedb13f3c74f0e13f3b4">ry</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> = wi -
<a name="l00243"></a>00243                                         dwody + 2 * Vector(<a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wo, n) * dndy +
<a name="l00244"></a>00244                                                  dDNdy * n);
<a name="l00245"></a>00245                                 L += scene-&gt;<a class="code" href="classScene.html#79235a24b9d34c5bd08b88f11da4cfb0">Li</a>(rd, sample) * f * <a class="code" href="geometry_8h.html#0c07ebf9b1454b90d6646a2f89065bb3">AbsDot</a>(wi, n);
<a name="l00246"></a>00246                         }
<a name="l00247"></a>00247                         f = bsdf-&gt;<a class="code" href="classBSDF.html#3498bc3b7484a17abd8429084124a0d4">Sample_f</a>(wo, &amp;wi,
<a name="l00248"></a>00248                                            <a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e669730">BxDFType</a>(<a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e669730806674394a395d5dc6383956e4f2733f">BSDF_TRANSMISSION</a> | <a class="code" href="reflection_8h.html#0ab5fa7d70a51348af4188bf1e669730c94ee1d195ff19883a79d1c08bd773d9">BSDF_SPECULAR</a>));
<a name="l00249"></a>00249                         <span class="keywordflow">if</span> (!f.<a class="code" href="classSpectrum.html#a05b9a07476998192a83fcd015b7619f">Black</a>()) {
<a name="l00250"></a>00250                                 <span class="comment">// Compute ray differential _rd_ for specular transmission</span>
<a name="l00251"></a>00251                                 <a class="code" href="classRayDifferential.html">RayDifferential</a> rd(p, wi);
<a name="l00252"></a>00252                                 rd.<a class="code" href="classRayDifferential.html#87b02390087bc229067e219becef4889">hasDifferentials</a> = <span class="keyword">true</span>;
<a name="l00253"></a>00253                                 rd.<a class="code" href="classRayDifferential.html#15359147d3814d4e592f005e1ae7ef64">rx</a>.<a class="code" href="classRay.html#1014d000992e8ac3b92f150ad2a5cffe">o</a> = p + isect.<a class="code" href="structIntersection.html#27a04b9f68291416c4d864f105a7d991">dg</a>.<a class="code" href="structDifferentialGeometry.html#307387af96b1b7a98358ea448114d697">dpdx</a>;
<a name="l00254"></a>00254                                 rd.<a class="code" href="classRayDifferential.html#a6ab98cac3bffedb13f3c74f0e13f3b4">ry</a>.<a class="code" href="classRay.html#1014d000992e8ac3b92f150ad2a5cffe">o</a> = p + isect.<a class="code" href="structIntersection.html#27a04b9f68291416c4d864f105a7d991">dg</a>.<a class="code" href="structDifferentialGeometry.html#d3e6794a920835d56713b3659437fc00">dpdy</a>;
<a name="l00255"></a>00255 
<a name="l00256"></a>00256                                 <span class="keywordtype">float</span> eta = bsdf-&gt;<a class="code" href="classBSDF.html#1ce3e1bc653569c466bba6ca6666f419">eta</a>;
<a name="l00257"></a>00257                                 <a class="code" href="classVector.html">Vector</a> w = -wo;
<a name="l00258"></a>00258                                 <span class="keywordflow">if</span> (<a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wo, n) &lt; 0) eta = 1.f / eta;
<a name="l00259"></a>00259 
<a name="l00260"></a>00260                                 <a class="code" href="classNormal.html">Normal</a> dndx = bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#f9daf909200f1e446924641fe27029b0">dndu</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#ee263cfa99de799c049e33b552ff0b82">dudx</a> + bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#cbcc594e0fff4bc1bf44d82e42cc6390">dndv</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#f0d770bd2370babb99b4eda972a757e2">dvdx</a>;
<a name="l00261"></a>00261                                 <a class="code" href="classNormal.html">Normal</a> dndy = bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#f9daf909200f1e446924641fe27029b0">dndu</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#57f8e45c2aa405802f25f1921ba27edd">dudy</a> + bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#cbcc594e0fff4bc1bf44d82e42cc6390">dndv</a> * bsdf-&gt;<a class="code" href="classBSDF.html#4fbebe2038bb049c50deba997a497f2d">dgShading</a>.<a class="code" href="structDifferentialGeometry.html#7343b4a399c5001c92fad44696060142">dvdy</a>;
<a name="l00262"></a>00262 
<a name="l00263"></a>00263                                 <a class="code" href="classVector.html">Vector</a> dwodx = -ray.<a class="code" href="classRayDifferential.html#15359147d3814d4e592f005e1ae7ef64">rx</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> - wo, dwody = -ray.<a class="code" href="classRayDifferential.html#a6ab98cac3bffedb13f3c74f0e13f3b4">ry</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> - wo;
<a name="l00264"></a>00264                                 <span class="keywordtype">float</span> dDNdx = <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(dwodx, n) + <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wo, dndx);
<a name="l00265"></a>00265                                 <span class="keywordtype">float</span> dDNdy = <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(dwody, n) + <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wo, dndy);
<a name="l00266"></a>00266 
<a name="l00267"></a>00267                                 <span class="keywordtype">float</span> mu = eta * <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(w, n) - <a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wi, n);
<a name="l00268"></a>00268                                 <span class="keywordtype">float</span> dmudx = (eta - (eta*eta*<a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(w,n))/<a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wi, n)) * dDNdx;
<a name="l00269"></a>00269                                 <span class="keywordtype">float</span> dmudy = (eta - (eta*eta*<a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(w,n))/<a class="code" href="geometry_8h.html#ad0d36688d190693f70a5c9ef9342914">Dot</a>(wi, n)) * dDNdy;
<a name="l00270"></a>00270 
<a name="l00271"></a>00271                                 rd.<a class="code" href="classRayDifferential.html#15359147d3814d4e592f005e1ae7ef64">rx</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> = wi + eta * dwodx - <a class="code" href="classVector.html">Vector</a>(mu * dndx + dmudx * n);
<a name="l00272"></a>00272                                 rd.<a class="code" href="classRayDifferential.html#a6ab98cac3bffedb13f3c74f0e13f3b4">ry</a>.<a class="code" href="classRay.html#5bbfc94e178e0a050b51200acc579069">d</a> = wi + eta * dwody - Vector(mu * dndy + dmudy * n);
<a name="l00273"></a>00273                                 L += scene-&gt;<a class="code" href="classScene.html#79235a24b9d34c5bd08b88f11da4cfb0">Li</a>(rd, sample) * f * <a class="code" href="geometry_8h.html#0c07ebf9b1454b90d6646a2f89065bb3">AbsDot</a>(wi, n);
<a name="l00274"></a>00274                         }
<a name="l00275"></a>00275                 }
<a name="l00276"></a>00276                 --<a class="code" href="classIGIIntegrator.html#177fcdea6b7cf6fd34521a229686af9f">specularDepth</a>;
<a name="l00277"></a>00277         }
<a name="l00278"></a>00278         <span class="keywordflow">else</span> {
<a name="l00279"></a>00279                 <span class="comment">// Handle ray with no intersection</span>
<a name="l00280"></a>00280                 <span class="keywordflow">if</span> (alpha) *alpha = 0.;
<a name="l00281"></a>00281                 <span class="keywordflow">for</span> (<a class="code" href="pbrt_8h.html#c319c165d52643e43249fe003e18bdf3">u_int</a> i = 0; i &lt; scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>.size(); ++i)
<a name="l00282"></a>00282                         L += scene-&gt;<a class="code" href="classScene.html#546553c8f864fd2aa709ad95a1d5ecc9">lights</a>[i]-&gt;Le(ray);
<a name="l00283"></a>00283                 <span class="keywordflow">if</span> (alpha &amp;&amp; !L.<a class="code" href="classSpectrum.html#a05b9a07476998192a83fcd015b7619f">Black</a>()) *alpha = 1.;
<a name="l00284"></a>00284                 <span class="keywordflow">return</span> L;
<a name="l00285"></a>00285         }
<a name="l00286"></a>00286         <span class="keywordflow">return</span> L;
<a name="l00287"></a>00287 }
<a name="l00288"></a><a class="code" href="igi_8cpp.html#292dab5329893d80187ce895221b3e6d">00288</a> <span class="keyword">extern</span> <span class="stringliteral">"C"</span> <a class="code" href="pbrt_8h.html#808e08638be3cba36e36759e5b150de0">DLLEXPORT</a> <a class="code" href="classSurfaceIntegrator.html">SurfaceIntegrator</a> *<a class="code" href="bidirectional_8cpp.html#292dab5329893d80187ce895221b3e6d">CreateSurfaceIntegrator</a>(<span class="keyword">const</span> <a class="code" href="classParamSet.html">ParamSet</a> &amp;params)
<a name="l00289"></a>00289 {
<a name="l00290"></a>00290         <span class="keywordtype">int</span> <a class="code" href="classIGIIntegrator.html#389eba70940527b81d7bfbf16cccfa85">nLightPaths</a> = params.<a class="code" href="classParamSet.html#23f24db9c46375e33764910c4fa999d5">FindOneInt</a>(<span class="stringliteral">"nlights"</span>, 64);
<a name="l00291"></a>00291         <span class="keywordtype">int</span> <a class="code" href="classIGIIntegrator.html#ceb019de8d99ce8527ef6370bc085503">nLightSets</a> = params.<a class="code" href="classParamSet.html#23f24db9c46375e33764910c4fa999d5">FindOneInt</a>(<span class="stringliteral">"nsets"</span>, 4);
<a name="l00292"></a>00292         <span class="keywordtype">float</span> minDist = params.<a class="code" href="classParamSet.html#87143fd8e73a75082179a110d42531a4">FindOneFloat</a>(<span class="stringliteral">"mindist"</span>, .1f);
<a name="l00293"></a>00293         <span class="keywordtype">float</span> rrThresh = params.<a class="code" href="classParamSet.html#87143fd8e73a75082179a110d42531a4">FindOneFloat</a>(<span class="stringliteral">"rrthreshold"</span>, .05f);
<a name="l00294"></a>00294         <span class="keywordtype">float</span> is = params.<a class="code" href="classParamSet.html#87143fd8e73a75082179a110d42531a4">FindOneFloat</a>(<span class="stringliteral">"indirectscale"</span>, 1.f);
<a name="l00295"></a>00295         <span class="keywordflow">return</span> <span class="keyword">new</span> <a class="code" href="classIGIIntegrator.html#fff5dc7931bb88a71636427a80419571">IGIIntegrator</a>(nLightPaths, nLightSets, minDist, rrThresh, is);
<a name="l00296"></a>00296 }
</pre></div></div>
<hr size="1"><address style="text-align: right;"><small>Generated on Mon Jul 20 17:31:53 2009 for pbrt by&nbsp;
<a href="http://www.doxygen.org/index.html">
<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.6 </small></address>
</body>
</html>