>>167
// return normalize(vec3( d(pos+dx)-d(pos-dx), d(pos+dy)-d(pos-dy), d(pos+dz)-d(pos-dz) ));
が分かりやすい。dx,dy,dzはepsilon。
面がこっちを向いているほどd(pos+dx)-d(pos-dx)が大きくなるよね？
d()はMandelBoxとHitするpositionまでの距離を求めるだけ。
iterationがあるから最後にpowするとか、式は結構変わってるけど。