Add blank lines for readability

Author: hollasch

src/InOneWeekend/sphere.h

@@ -40,6 +40,7 @@ bool sphere::hit(const ray& r, double t_min, double t_max, hit_record& rec) cons
 
     if (discriminant > 0) {
         auto root = sqrt(discriminant);
+
         auto temp = (-half_b - root)/a;
         if (temp < t_max && temp > t_min) {
             rec.t = temp;
@@ -49,6 +50,7 @@ bool sphere::hit(const ray& r, double t_min, double t_max, hit_record& rec) cons
             rec.mat_ptr = mat_ptr;
             return true;
         }
+
         temp = (-half_b + root) / a;
         if (temp < t_max && temp > t_min) {
             rec.t = temp;

src/TheNextWeek/aarect.h

@@ -79,53 +79,62 @@ bool xy_rect::hit(const ray& r, double t0, double t1, hit_record& rec) const {
     auto t = (k-r.origin().z()) / r.direction().z();
     if (t < t0 || t > t1)
         return false;
+
     auto x = r.origin().x() + t*r.direction().x();
     auto y = r.origin().y() + t*r.direction().y();
     if (x < x0 || x > x1 || y < y0 || y > y1)
         return false;
+
     rec.u = (x-x0)/(x1-x0);
     rec.v = (y-y0)/(y1-y0);
     rec.t = t;
     vec3 outward_normal = vec3(0, 0, 1);
     rec.set_face_normal(r, outward_normal);
     rec.mat_ptr = mp;
     rec.p = r.at(t);
+
     return true;
 }
 
 bool xz_rect::hit(const ray& r, double t0, double t1, hit_record& rec) const {
     auto t = (k-r.origin().y()) / r.direction().y();
     if (t < t0 || t > t1)
         return false;
+
     auto x = r.origin().x() + t*r.direction().x();
     auto z = r.origin().z() + t*r.direction().z();
     if (x < x0 || x > x1 || z < z0 || z > z1)
         return false;
+
     rec.u = (x-x0)/(x1-x0);
     rec.v = (z-z0)/(z1-z0);
     rec.t = t;
     vec3 outward_normal = vec3(0, 1, 0);
     rec.set_face_normal(r, outward_normal);
     rec.mat_ptr = mp;
     rec.p = r.at(t);
+
     return true;
 }
 
 bool yz_rect::hit(const ray& r, double t0, double t1, hit_record& rec) const {
     auto t = (k-r.origin().x()) / r.direction().x();
     if (t < t0 || t > t1)
         return false;
+
     auto y = r.origin().y() + t*r.direction().y();
     auto z = r.origin().z() + t*r.direction().z();
     if (y < y0 || y > y1 || z < z0 || z > z1)
         return false;
+
     rec.u = (y-y0)/(y1-y0);
     rec.v = (z-z0)/(z1-z0);
     rec.t = t;
     vec3 outward_normal = vec3(1, 0, 0);
     rec.set_face_normal(r, outward_normal);
     rec.mat_ptr = mp;
     rec.p = r.at(t);
+
     return true;
 }
 

src/TheNextWeek/constant_medium.h

@@ -55,6 +55,7 @@ bool constant_medium::hit(const ray& r, double t_min, double t_max, hit_record&
 
             if (rec1.t >= rec2.t)
                 return false;
+
             if (rec1.t < 0)
                 rec1.t = 0;
 

src/TheNextWeek/moving_sphere.h

@@ -59,8 +59,11 @@ bool moving_sphere::hit(const ray& r, double t_min, double t_max, hit_record& re
     auto a = dot(r.direction(), r.direction());
     auto b = dot(oc, r.direction());
     auto c = dot(oc, oc) - radius*radius;
+
     auto discriminant = b*b - a*c;
+
     if (discriminant > 0) {
+
         auto temp = (-b - sqrt(discriminant))/a;
         if (temp < t_max && temp > t_min) {
             rec.t = temp;
@@ -70,6 +73,7 @@ bool moving_sphere::hit(const ray& r, double t_min, double t_max, hit_record& re
             rec.mat_ptr = mat_ptr;
             return true;
         }
+
         temp = (-b + sqrt(discriminant))/a;
         if (temp < t_max && temp > t_min) {
             rec.t = temp;
@@ -80,6 +84,7 @@ bool moving_sphere::hit(const ray& r, double t_min, double t_max, hit_record& re
             return true;
         }
     }
+
     return false;
 }
 

src/TheRestOfYourLife/cos_cubed.cc

@@ -18,15 +18,18 @@
 
 int main() {
     int N = 1000000;
+
     auto sum = 0.0;
     for (int i = 0; i < N; i++) {
         auto r1 = random_double();
         auto r2 = random_double();
         auto x = cos(2*pi*r1)*2*sqrt(r2*(1-r2));
         auto y = sin(2*pi*r1)*2*sqrt(r2*(1-r2));
         auto z = 1 - r2;
+
         sum += z*z*z / (1.0/(2.0*pi));
     }
+
     std::cout << std::fixed << std::setprecision(12);
     std::cout << "PI/2 = " << pi/2 << '\n';
     std::cout << "Estimate = " << sum/N << '\n';

src/TheRestOfYourLife/cos_density.cc

@@ -23,17 +23,20 @@ inline vec3 random_cosine_direction() {
     auto phi = 2*pi*r1;
     auto x = cos(phi)*sqrt(r2);
     auto y = sin(phi)*sqrt(r2);
+
     return vec3(x, y, z);
 }
 
 
 int main() {
     int N = 1000000;
+
     auto sum = 0.0;
     for (int i = 0; i < N; i++) {
         vec3 v = random_cosine_direction();
         sum += v.z()*v.z()*v.z() / (v.z()/(pi));
     }
+
     std::cout << std::fixed << std::setprecision(12);
     std::cout << "PI/2 = " << pi/2 << '\n';
     std::cout << "Estimate = " << sum/N << '\n';

src/TheRestOfYourLife/integrate_x_sq.cc

@@ -25,11 +25,13 @@ int main() {
     int inside_circle = 0;
     int inside_circle_stratified = 0;
     int N = 1;
+
     auto sum = 0.0;
     for (int i = 0; i < N; i++) {
-            auto x = pow(random_double(0,8), 1./3.);
-            sum += x*x / pdf(x);
+        auto x = pow(random_double(0,8), 1./3.);
+        sum += x*x / pdf(x);
     }
+
     std::cout << std::fixed << std::setprecision(12);
     std::cout << "I = " << sum/N << '\n';
 }

src/TheRestOfYourLife/pi.cc

@@ -21,14 +21,18 @@ int main() {
     int inside_circle = 0;
     int inside_circle_stratified = 0;
     int sqrt_N = 10000;
+
     for (int i = 0; i < sqrt_N; i++) {
         for (int j = 0; j < sqrt_N; j++) {
             auto x = random_double(-1,1);
             auto y = random_double(-1,1);
+
             if (x*x + y*y < 1)
                 inside_circle++;
+
             x = 2*((i + random_double()) / sqrt_N) - 1;
             y = 2*((j + random_double()) / sqrt_N) - 1;
+
             if (x*x + y*y < 1)
                 inside_circle_stratified++;
         }

src/TheRestOfYourLife/sphere_plot.cc

@@ -22,6 +22,7 @@ int main() {
         auto x = cos(2*pi*r1)*2*sqrt(r2*(1-r2));
         auto y = sin(2*pi*r1)*2*sqrt(r2*(1-r2));
         auto z = 1 - 2*r2;
+
         std::cout << x << " " << y << " " << z << '\n';
     }
 }