I set out on this project with the goal of writing a software path tracer that could render Himalayan Salt Lamps using subsurface scattering, but got sidetracked making performance improvements. Vectorization with numpy led to a 300x speed improvement, and I learned a lot about caches and data locality along the way.