# | Line 102 | Line 102 | double omega; /* light source size */ | |
---|---|---|
102 | /* + source if flat */ | |
103 | if (np->specfl & SP_FLAT) | |
104 | dtmp += omega/(4.0*PI); | |
105 | < | /* delta */ |
105 | > | /* half vector */ |
106 | vtmp[0] = ldir[0] - np->rp->rdir[0]; | |
107 | vtmp[1] = ldir[1] - np->rp->rdir[1]; | |
108 | vtmp[2] = ldir[2] - np->rp->rdir[2]; | |
109 | d2 = DOT(vtmp, np->pnorm); | |
110 | < | d2 = 2.0 - 2.0*d2/sqrt(DOT(vtmp,vtmp)); |
110 | > | d2 *= d2; |
111 | > | d2 = (DOT(vtmp,vtmp) - d2) / d2; |
112 | /* gaussian */ | |
113 | dtmp = exp(-d2/dtmp)/(4.*PI*dtmp); | |
114 | /* worth using? */ |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |