     const float3 H           = optix::normalize( w_out + w_in );
     const float  inv_alpha_u = 1.0f / alpha_u;
     const float  inv_alpha_v = 1.0f / alpha_v;
-    const float  spec_coeff = cos_in * legion::ONE_DIV_PI *
-                              0.25f * inv_alpha_u * inv_alpha_v / sqrtf( cos_in*cos_out ) *
-                              wardExpTerm( inv_alpha_u, inv_alpha_v, H, onb );
+    const float  spec_coeff = 
+        cos_in * legion::ONE_DIV_PI *
+        0.25f * inv_alpha_u * inv_alpha_v / sqrtf( cos_in*cos_out ) *
+        wardExpTerm( inv_alpha_u, inv_alpha_v, H, onb );
 
     return spec_coeff * spec_reflectance;
 }

        float x_coeff   = cos_phi_h / alpha_u;
        float y_coeff   = sin_phi_h / alpha_v;

        float theta_h   = atanf( sqrtf( -logf( xi2 ) / ( x_coeff*x_coeff + y_coeff*y_coeff ) ) );
        float sin_theta_h = sinf( theta_h );
        float cos_theta_h = 1.0f - sin_theta_h*sin_theta_h;
        if (cos_theta_h > 0.0f)

        else if (xi1 < 0.5f)
        {
            xi1 = 4.0f * (xi1 - 0.25f);
            phi_h = legion::PI - atanf(alpha_v * tanf(0.5f * legion::PI * xi1) /
                    alpha_u);
        }
        else if (xi1 < 0.75f)
        {
            xi1 = 4.0f * (xi1 - 0.5f);
            phi_h = legion::PI + atanf(alpha_v * tanf(0.5f * legion::PI * xi1) /
                    alpha_u);
        }
        else
        {
            xi1 = 4.0f * (xi1 - 0.75f);
            phi_h = legion::PI*2.0f - 
                    atanf(alpha_v * tanf( 0.5f * legion::PI * xi1) / alpha_u);
        }

        float cos_phi_h = cosf( phi_h );
