mingzailao
/
mitsuba
mirror of https://github.com/Quartenia/mitsuba.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

added the empty space bonus

metadata
Wenzel Jakob 2010-10-10 13:23:04 +02:00
parent 1058cfe7a7
commit 87b82cc5e1
1 changed files with 14 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
include/mitsuba/render/gkdtree.h
View File

@ -275,10 +275,10 @@ public:
GenericKDTree() : m_root(NULL) { GenericKDTree() : m_root(NULL) {
m_traversalCost = 15; m_traversalCost = 15;
m_intersectionCost = 20; m_intersectionCost = 20;
m_emptySpaceBonus = 0.8f; m_emptySpaceBonus = 0.9f;
m_clip = true; m_clip = true;
m_stopPrims = 2; m_stopPrims = 1;
m_maxBadRefines = 3; m_maxBadRefines = 2;
m_exactDepth = 1; m_exactDepth = 1;
m_maxDepth = 100; m_maxDepth = 100;
} }
@ -1037,6 +1037,7 @@ protected:
/* Calculate a score using the surface area heuristic */ /* Calculate a score using the surface area heuristic */
if (EXPECT_TAKEN(pos >= nodeAABB.min[axis] && pos <= nodeAABB.max[axis])) { if (EXPECT_TAKEN(pos >= nodeAABB.min[axis] && pos <= nodeAABB.max[axis])) {
size_type nL = numLeft[axis], nR = numRight[axis];
Float tmp = nodeAABB.max[axis]; Float tmp = nodeAABB.max[axis];
aabb.max[axis] = pos; aabb.max[axis] = pos;
Float pLeft = invSA * aabb.getSurfaceArea(); Float pLeft = invSA * aabb.getSurfaceArea();
@ -1046,22 +1047,26 @@ protected:
Float pRight = invSA * aabb.getSurfaceArea(); Float pRight = invSA * aabb.getSurfaceArea();
aabb.min[axis] = tmp; aabb.min[axis] = tmp;
Float sahCostPlanarLeft = m_traversalCost + m_intersectionCost Float sahCostPlanarLeft = m_traversalCost + m_intersectionCost
* (pLeft * (numLeft[axis] + numPlanar) + pRight * numRight[axis]); * (pLeft * (nL + numPlanar) + pRight * nR);
Float sahCostPlanarRight = m_traversalCost + m_intersectionCost Float sahCostPlanarRight = m_traversalCost + m_intersectionCost
* (pLeft * numLeft[axis] + pRight * (numRight[axis] + numPlanar)); * (pLeft * nL + pRight * (nR + numPlanar));
if (nL + numPlanar == 0 || nR == 0)
sahCostPlanarLeft *= m_emptySpaceBonus;
if (nL == 0 || nR + numPlanar == 0)
sahCostPlanarRight *= m_emptySpaceBonus;
if (sahCostPlanarLeft < bestSplit.sahCost || sahCostPlanarRight < bestSplit.sahCost) { if (sahCostPlanarLeft < bestSplit.sahCost || sahCostPlanarRight < bestSplit.sahCost) {
bestSplit.pos = pos; bestSplit.pos = pos;
bestSplit.axis = axis; bestSplit.axis = axis;
if (sahCostPlanarLeft < sahCostPlanarRight) { if (sahCostPlanarLeft < sahCostPlanarRight) {
bestSplit.sahCost = sahCostPlanarLeft; bestSplit.sahCost = sahCostPlanarLeft;
bestSplit.numLeft = numLeft[axis] + numPlanar; bestSplit.numLeft = nL + numPlanar;
bestSplit.numRight = numRight[axis]; bestSplit.numRight = nR;
bestSplit.planarLeft = true; bestSplit.planarLeft = true;
} else { } else {
bestSplit.sahCost = sahCostPlanarRight; bestSplit.sahCost = sahCostPlanarRight;
bestSplit.numLeft = numLeft[axis]; bestSplit.numLeft = nL;
bestSplit.numRight = numRight[axis] + numPlanar; bestSplit.numRight = nR + numPlanar;
bestSplit.planarLeft = false; bestSplit.planarLeft = false;
} }
} }