Hi there, code author here. I've been lurking here for years, but couldn't resist an invitation like this.
I've always enjoyed Fabien's analyses so it's great fun to see what he makes of my own code. Anyway, a few clarifications from what I remember of this. (I wrote this '09, and I don't have my old notes in front of me right now):
* The `n` value returned by the trace function, `T()`, is the surface normal. It returns this whether it hit the plane or a sphere.
* The `r` vector in `S()` is the reflection vector off whatever was hit.
* The mystery `c` point in the `main()` function is the offset from the eye point (ignoring the lens perturbation `t`) to the corner of the focal plane. Note below that we're nominally tracing rays for pixel (x,y) in the direction of ax+by+c. At the image midpoint, this should be `g`.
* "although I suspect this to be a side effect: This is not how soft-shadows are done." True. There are true soft shadows in this, but this isn't where they're computed. That's the job of the randomization where `l` is computed in `S()`.
Anyway, please feel free to ask any questions about it.
I've always enjoyed Fabien's analyses so it's great fun to see what he makes of my own code. Anyway, a few clarifications from what I remember of this. (I wrote this '09, and I don't have my old notes in front of me right now):
* The `n` value returned by the trace function, `T()`, is the surface normal. It returns this whether it hit the plane or a sphere.
* The `r` vector in `S()` is the reflection vector off whatever was hit.
* The mystery `c` point in the `main()` function is the offset from the eye point (ignoring the lens perturbation `t`) to the corner of the focal plane. Note below that we're nominally tracing rays for pixel (x,y) in the direction of ax+by+c. At the image midpoint, this should be `g`.
* "although I suspect this to be a side effect: This is not how soft-shadows are done." True. There are true soft shadows in this, but this isn't where they're computed. That's the job of the randomization where `l` is computed in `S()`.
Anyway, please feel free to ask any questions about it.