Project 3: Rishabh Shah

CMakeLists.txt

@@ -72,7 +72,7 @@ if(${CMAKE_SYSTEM_NAME} MATCHES "Darwin")
 endif()
 
 include_directories(.)
-#add_subdirectory(stream_compaction)  # TODO: uncomment if using your stream compaction
+add_subdirectory(stream_compaction)  # TODO: uncomment if using your stream compaction
 add_subdirectory(src)
 
 cuda_add_executable(${CMAKE_PROJECT_NAME}
@@ -82,7 +82,7 @@ cuda_add_executable(${CMAKE_PROJECT_NAME}
 
 target_link_libraries(${CMAKE_PROJECT_NAME}
     src
-    #stream_compaction  # TODO: uncomment if using your stream compaction
+    stream_compaction  # TODO: uncomment if using your stream compaction
     ${CORELIBS}
     )
 

README.md

@@ -3,11 +3,121 @@ CUDA Path Tracer
 
 **University of Pennsylvania, CIS 565: GPU Programming and Architecture, Project 3**
 
-* (TODO) YOUR NAME HERE
-* Tested on: (TODO) Windows 22, i7-2222 @ 2.22GHz 22GB, GTX 222 222MB (Moore 2222 Lab)
+* Rishabh Shah
+* Tested on: Windows 10, i7-6700HQ @ 2.6GHz 16GB, GTX 960M 4096MB (Laptop)
 
-### (TODO: Your README)
+## Overview
+Monte-Carlo Pathtracing is a way of approximating the behavior of light in the physical world. This technique is used in rendering photorealistic images. Because light rays bounce off the surfaces based on the surface properties, trying to calculate this is a computationally intensive process. Movie studios generally use CPU pathtracers due to memory limitations of GPUs, but the algorithm itself is highly parallel.
 
-*DO NOT* leave the README to the last minute! It is a crucial part of the
-project, and we will not be able to grade you without a good README.
+In this project, I have implemented a GPU pathtracer.
 
+## Functionalities Implemented
+
+* Shading:
+    * Lambertian shading for diffuse surfaces
+    * Perfectly specular (Reflective) surfaces
+    * Refractive surfaces using Schlick's approximation. There is a `#define FRESNEL` flag in interaction.cpp, but it needs some fixing.
+
+
+* Stream Compaction using `thrust::partition`
+
+* Sorting paths and intersections based on material
+
+* Caching first bounce, and using it in all subsequent iterations
+
+* Physically-based depth-of-field (by jittering rays within an aperture) (based on PBRT Thin Lens Camera Model)
+
+* Stochastic Sampled Antialiasing.
+
+* Direct Lighting - directly sampling light after the last bounce
+
+* Purely Direct Lighting (based on PBRT) - I implemented this as the first step towards an integrator with Multiple Importance Sampling (MIS). This converges very fast, but does not take into account bouncing of light, so lacks global illumination.
+
+* Russian Roulette path termination
+
+* Full Lighting Integrator - This integrator uses MIS to give weights to Direct Light contribution, and BRDF sampling contribution using a heuristic.
+    * This is currently a work in progress. There is some bug which causes the image to blow-out instead of converging.
+
+
+## Performance Analysis
+
+##### Time comparison of different integration schemes
+
+Direct lighting is only a 1 bounce process. So it is the fastest. Naive integration on the other hand involves light mouncing based on BSDF of the surface, so it is slower. Direct lighting on the last bounce, does not significantly increase the total time cost, but it is difficult to properly weight the contributions of direct and indirect lighting components this way, so not all scenes will look correct.
+
+|   Integrator   |   1000 spp   |   5000 spp |
+|:---------------------|-------:|---------------:|
+|Naive	|61	|306 |
+|Direct	|17	|87	|
+|Naive and Direct on last bounce |70	|347 |
+
+*Plot: performance for a depth of 8*
+
+![](img/integration.png)
+
+
+##### Time comparison of different optimization techniques
+
+Stream compaction should reduce the time taken for rendering, but it does not. Perhaps, the depth here is not enough to make a difference, and in fact the overhead is higher than not sorting.
+
+Caching first bounce effectively reduces 1 bounce from the tracing, so when the number of bounces is low, it definitely makes things faster.
+
+Material sorting was implemented based on the idea that branch divergence would reduce if paths and intersections of same material are contiguous in memory. But it turns out that the sorting overhead is so high that the total computation time is many times higher than not doing it.
+
+Russian roulette path termination is a way of terminating the path without going up to the maximum depth if enough color contribution is already available for the path. This becomes increasingly helpful with larger depths.
+
+|   Optimizattion   |   Depth 8   |   Depth 16 |
+|:---------------------|-------:|---------------:|
+|Without stream compaction	|61	|108 |
+|With stream compaction	|123 |137 |
+|Cashing first bounce	|59	|104 |
+|Material sorting |473 |945 |
+|Russian roulette |60 |104 |
+
+*Plot: performance for 1000 spp*
+
+![](img/optimization.png)
+
+## Sample Renders
+
+#### Naive Integrator
+
+![](img/renders/naive_16_5000_800x800_cornell.2017-10-02_01-19-18z.5000samp.png)
+
+#### Direct Lighting Integrator
+
+![](img/renders/direct_5000_800x800_cornell.2017-10-02_01-32-55z.5000samp.png)
+
+#### Naive Integrator with Direct Lighting at last bounce
+
+![](img/renders/naive_direct_last_bounce_cornell.2017-10-02_01-54-08z.5000samp.png)
+
+#### Anti-Aliasing
+
+This cannot be used with first bounce cache, as we are jittering the rays from the camera before the first bounce for AA.
+
+*Left: with anti-aliasing vs Right: without anti-aliasing*
+
+![](img/renders/nonAA.png)
+![](img/renders/AA.png)
+
+#### Depth of Field
+
+![](img/renders/cornell.2017-10-02_03-23-09z.133samp.png)
+
+#### Full Lighting Integrator with Multiple Importance Sampling
+As this uses both direct and indirect lighting computations with weights, it give global illumination, and converges very quickly. This image below is a nearly converged result in only about 2000 spp, as compared to very noisy image we get by naive integration at 5000 spp. There is some bug due to which, I was losing color. Fixing this issue creates other issues, so I'm still working on it.
+
+*Still Buggy*
+
+![](img/renders/cornell.2017-10-06_22-51-41z.5000samp.png)
+
+#### Bonus
+
+*Specular and Refractive*
+
+![](img/renders/cornell.2017-09-30_02-52-01z.5000samp.png)
+
+*Noise Art*
+
+![](img/renders/cornell.2017-09-25_14-12-54z.153samp.png)

scenes/cornell.txt

@@ -1,12 +1,12 @@
 // Emissive material (light)
 MATERIAL 0
-RGB         1 1 1
+RGB         1 0.5 1
 SPECEX      0
 SPECRGB     0 0 0
 REFL        0
 REFR        0
 REFRIOR     0
-EMITTANCE   5
+EMITTANCE   1
 
 // Diffuse white
 MATERIAL 1
@@ -48,25 +48,68 @@ REFR        0
 REFRIOR     0
 EMITTANCE   0
 
+// Transmissive white - NEW
+// IOR glass = 1.517
+// IOR water = 1.333
+MATERIAL 5
+RGB         .98 .98 .98
+SPECEX      0
+SPECRGB     .98 .98 .98
+REFL        0
+REFR        1
+REFRIOR     1.517
+EMITTANCE   0
+
+// Transmissive yellow
+MATERIAL 6
+RGB         .98 .98 .2
+SPECEX      0
+SPECRGB     .98 .98 .2
+REFL        0
+REFR        1
+REFRIOR     1.517
+EMITTANCE   0
+
+// Diffuse yellow
+MATERIAL 7
+RGB         .98 .98 .2
+SPECEX      0
+SPECRGB     0 0 0
+REFL        0
+REFR        0
+REFRIOR     0
+EMITTANCE   0
+
+MATERIAL 8
+RGB         0.5 1 1
+SPECEX      0
+SPECRGB     0 0 0
+REFL        0
+REFR        0
+REFRIOR     0
+EMITTANCE   5
+
 // Camera
 CAMERA
 RES         800 800
 FOVY        45
 ITERATIONS  5000
-DEPTH       8
+DEPTH       16
 FILE        cornell
+FOCALLENGTH 5.0
+LENSRAD     0.2
 EYE         0.0 5 10.5
 LOOKAT      0 5 0
 UP          0 1 0
 
 
 // Ceiling light
 OBJECT 0
-cube
+sphere
 material 0
-TRANS       0 10 0
+TRANS       2 7 -2
 ROTAT       0 0 0
-SCALE       3 .3 3
+SCALE       2 2 2
 
 // Floor
 OBJECT 1
@@ -112,6 +155,30 @@ SCALE       .01 10 10
 OBJECT 6
 sphere
 material 4
-TRANS       -1 4 -1
+TRANS       -3 2 0
+ROTAT       0 0 0
+SCALE       2.5 2.5 2.5
+
+// Sphere
+OBJECT 7
+sphere
+material 5
+TRANS       3 2 0
+ROTAT       0 0 0
+SCALE       2.5 2.5 2.5
+
+// Sphere
+OBJECT 8
+sphere
+material 7
+TRANS       0 2 0
+ROTAT       0 0 0
+SCALE       2.5 2.5 2.5
+
+// light yellow
+OBJECT 9
+sphere
+material 8
+TRANS       -2 7 -2
 ROTAT       0 0 0
-SCALE       3 3 3
+SCALE       2 2 2