launch

@@ -1,13 +1,12 @@
 #!/usr/bin/env bash
 
 executable="$1"
-options="$2"
 prime="prime-run"
 
 if command -v $prime &> /dev/null
 then
-    echo "launching $executable $options with nvidia prime offloading"
-    $prime $executable $options
+    echo "launching $executable with nvidia prime offloading"
+    $prime $executable
 else
-    $executable $options
+    $executable
 fi

makefile

@@ -46,11 +46,7 @@ clean:
 	-rm -r $(BIN_DIR)
 	-rm */*.o
 
-# TODO: extract run-* out to like, python scripts? arguments can be implemented with argparse to
-# manage simple interface with C binary
 run: $(TARGET)
-	$(LAUNCH) $(TARGET) 1.0
-run-converge: $(TARGET) 
-	$(LAUNCH) $(TARGET) 0.0
+	$(LAUNCH) $(TARGET)
 
 .PHONY: run clean

shader/include/constants.glsl

@@ -1,7 +1,7 @@
-const float INF = 30.0;
+const float INF = 50.0;
 const float PI = 3.14159;
 const float E = 2.71828;
-const int BOUNCES = 5;
+const int BOUNCES = 4;
 
 // materials
 const int MAT_SKY = -1;

shader/include/depth.glsl

@@ -1,16 +0,0 @@
-float getLinearDepth(float distance)
-{
-    return distance/INF;
-}
-
-float getLogarithmicDepth(float distance)
-{
-    // n roughly correlates to steepness of log curve
-    // TODO: what does this mean in mathematical terms??
-    float n = 4;
-    float f = INF;
-    float z = distance;
-    
-    // logarithmic depth
-    return max(0,log(z*pow(E,n)/f)/n);
-}

shader/include/intersect.glsl

@@ -33,7 +33,7 @@ void intersectPlane(Ray ray, inout RayHit bestHit, vec3 p, vec3 normal)
             bestHit.distance = t;
             bestHit.position = ray.origin + t*ray.direction;
             bestHit.normal = normal;
-            bestHit.albedo = vec3(1.0,.4,.4);
+            bestHit.albedo = vec3(1.0,1.0,1.0);
             bestHit.material = MAT_LAMBERT;
         }
     }

shader/include/ray.glsl

@@ -3,7 +3,6 @@ struct Ray
     vec3 origin;
     vec3 direction;
     vec3 energy;
-    float distance;
 };
 Ray createRay(vec3 origin, vec3 direction)
 {
@@ -12,7 +11,6 @@ Ray createRay(vec3 origin, vec3 direction)
     ray.origin = origin;
     ray.direction = direction;
     ray.energy = vec3(1.0,1.0,1.0);
-    ray.distance = 0;
 
     return ray;
 }

shader/include/scene.glsl

@@ -1,17 +1,19 @@
-RayHit trace(inout Ray ray)
+RayHit trace(Ray ray)
 {
     RayHit hit = createRayHit();
 
-    // floor
+    Sphere s;
+    s.cr = vec4(0.0,0.0,0.0,5.0);
+    s.albedo = vec3(1.0,0.0,0.0);
+    s.material = MAT_CHROME;
+    intersectSphere(ray, hit, s);
+
     intersectPlane(ray, hit, vec3(0.0,-1.5,0.0),vec3(0.0,1.0,0.0));
 
-    // spheres
     for (int i = 0; i < _activeSpheres; i++)
     {
         intersectSphere(ray, hit, _spheres[i]);
     }
 
-    ray.distance += hit.distance * float(hit.distance < INF);
-
     return hit;
 }

shader/root/rt.glsl

@@ -24,7 +24,6 @@ uniform vec3 _skyColor = vec3(0.68,0.85,0.9);
 // scene.glsl includes scene trace function
 #include scene.glsl     
 #include lighting.glsl
-#include depth.glsl
 
 vec3 shade(inout Ray ray, RayHit hit)
 {
@@ -78,27 +77,16 @@ void main()
     int bounces = (1-sky) * BOUNCES; 
     pixel.xyz = mix(pixel.xyz, _skyColor, sky);
 
-    // sample
-    int samples = 2;
-    for (int i = 0; i < samples; i++)
-    {
-        float sampleDepth = 0;
-
-        // trace the ray's path around the scene
+    // trace the rays path around the scene
     for (int j = 0; j < bounces; j++)
     {
         RayHit hit = trace(ray);
 
-            depth = getLogarithmicDepth(ray.distance);
-
         pixel.xyz += ray.energy * shade(ray, hit);
 
         if (length(ray.energy) < 0.001) break;
     }
 
-        depth += sampleDepth / float(samples);
-    }
-
     pixel.xyz = mix(pixel.xyz, vec3(1.0), depth);
 
     // output to a specific pixel in the image

shader/root/rtpre.glsl

@@ -18,7 +18,6 @@ layout (rgba32f, binding = 3) writeonly uniform image2D g1_output;
 #include camera.glsl
 #include image.glsl
 #include scene.glsl
-#include depth.glsl
 
 void main()
 {
@@ -35,7 +34,18 @@ void main()
     Ray ray = createCameraRay(uv);
     RayHit hit = trace(ray);
 
-    float depth = getLogarithmicDepth(hit.distance);
+    // n roughly correlates to steepness of log curve
+    // TODO: what does this mean in mathematical terms??
+    float n = 4;
+    float f = INF;
+    float z = hit.distance;
+    
+    float depth;
+    // linear depth
+    //depth = z/f;
+
+    // logarithmic depth
+    depth = max(0,log(z*pow(E,n)/f)/n);
 
     // pack normal into texture
     pixel.xyz = hit.normal*0.5+0.5;

src/cam.c

@@ -2,8 +2,10 @@
 
 const float FOVY = 90.0;
 
-void updateCameraUniforms(GLuint shaderProgram, float aspect, float t)
+void updateCameraUniforms(GLuint shaderProgram, float aspect)
 {
+    float t = now();
+
     // wobble up dir
     vec3 cdir, cright, cup;
     vec3 up = {0.1*sin(t),1.0,0.05*cos(0.5*t)};

src/cam.h

@@ -6,4 +6,4 @@
 
 #include "clock.h"
 
-void updateCameraUniforms(GLuint shaderProgram, float aspect, float t);
+void updateCameraUniforms(GLuint shaderProgram, float aspect);

src/clock.c

@@ -4,8 +4,3 @@ float now()
 { 
     return (float)SDL_GetTicks() / 1000.0; 
 }
-
-float nowScaled(struct Epoch t)
-{ 
-    return t.scale * now();
-}

src/clock.h

@@ -2,12 +2,5 @@
 
 #include <SDL2/SDL.h>
 
-struct Epoch
-{
-    float scale;
-};
-
-// real seconds since program start
+// seconds since program start
 float now();
-// simulation seconds since program start
-float nowScaled(struct Epoch e);

src/gfx.c

@@ -174,6 +174,16 @@ void setVertexAttributes()
     // color
     glVertexAttribPointer(1, 3, GL_FLOAT, GL_FALSE, stride, (void*)(3*sizeof(float)));
     glEnableVertexAttribArray(1);
+
+
+
+
+
+
+
+
+
+
     // uv 
     glVertexAttribPointer(2, 2, GL_FLOAT, GL_FALSE, stride, (void*)(6*sizeof(float)));
     glEnableVertexAttribArray(2);

src/main.c

@@ -10,13 +10,12 @@
 const int WIDTH = 420;
 const int HEIGHT = 420;
 
-void updateUniforms(GLuint shaderProgram, float t);
+void updateUniforms(GLuint shaderProgram);
 
 SDL_Window *window;
 
 struct Shaders shaders;
 struct Textures textures;
-struct Epoch epoch;
 
 void initialise()
 {
@@ -33,21 +32,8 @@ void initialise()
     createTextures(WIDTH, HEIGHT, shaders, &textures);
 }
 
-void parseArgs(int argc, char* argv[], struct Epoch* e)
+int main()
 {
-    // check we have the right number of arguments
-    if (argc != 2)
-    {
-        fprintf(stderr, "usage: oglc TIMESPEED\n");
-    }
-
-    sscanf(argv[1], "%f", &(e->scale));
-}
-
-int main(int argc, char* argv[])
-{
-    parseArgs(argc, argv, &epoch);
-
     initialise();
 
     float start = now();
@@ -55,7 +41,6 @@ int main(int argc, char* argv[])
     for (frames = 0; !checkQuit(); frames++)
     {
         GLuint shader;
-        float t = nowScaled(epoch);
 
         // prepass
         // TODO: write output to different texture than main output
@@ -63,17 +48,18 @@ int main(int argc, char* argv[])
         glActiveTexture(GL_TEXTURE0);
         glBindTexture(GL_TEXTURE_2D, textures.g0);
         glUseProgram(shader);
-        updateUniforms(shader, t);
+        updateUniforms(shader);
         glDispatchCompute((GLuint)WIDTH, (GLuint)HEIGHT, 1);
 
         // make sure we're finished writing to the texture before trying to read it
         glMemoryBarrier(GL_ALL_BARRIER_BITS);
 
+        // dispatch compute shaders
         shader = shaders.lighting;
         glActiveTexture(GL_TEXTURE0);
         glBindTexture(GL_TEXTURE_2D, textures.target);
         glUseProgram(shader);
-        updateUniforms(shader, t);
+        updateUniforms(shader);
         //int loc = glGetUniformLocation(shader, "_g0");
         //glUniform1i(loc, textures.g0);
         glDispatchCompute((GLuint)WIDTH, (GLuint)HEIGHT, 1);
@@ -81,18 +67,6 @@ int main(int argc, char* argv[])
         // make sure we're finished writing to the texture before trying to read it
         glMemoryBarrier(GL_SHADER_IMAGE_ACCESS_BARRIER_BIT);
 
-        // TODO: frame blending
-        //
-        // before we get to our normal rendering pass we need to blend our samples
-        // across multiple frames. to do this, we should maintain an accumulation
-        // buffer where we can store pixel data while a new frame is being generated.
-        // this can be a collection of two buffers, where one is written to and the 
-        // other is read from. to render a new frame, the newly sampled raw frame can
-        // be combined with data from the previous frame. by repeatedly taking the
-        // average over a large number of frames a cleaner image can be generated.
-        // this is most effective for static images, and will produce motion blur on
-        // a moving image.
-
         // normal drawing pass
         shader = shaders.quad;
         glUseProgram(shader);
@@ -113,7 +87,7 @@ int main(int argc, char* argv[])
     return 0;
 }
 
-void updateUniforms(GLuint shaderProgram, float t)
+void updateUniforms(GLuint shaderProgram)
 {
     // update random values
     vec4 seed = 
@@ -126,17 +100,19 @@ void updateUniforms(GLuint shaderProgram, float t)
     int loc = glGetUniformLocation(shaderProgram, "_seed");
     glUniform4fv(loc, 1, seed);
 
+    // update time
+    float t = now();
     float sin_t = sin(t);
     loc = glGetUniformLocation(shaderProgram, "_t");
     glUniform4f(loc, t, sin_t, (1.0 + sin_t)*0.5, 0.0f);
 
     // update camera
     float aspect = (float)WIDTH/(float)HEIGHT;
-    updateCameraUniforms(shaderProgram, aspect, t);
+    updateCameraUniforms(shaderProgram, aspect);
 
     // make and update spheres
-    const int sphereCount = 41;
+    const int sphereCount = 42;
     struct Sphere spheres[sphereCount];
-    makeSpheres(spheres, sphereCount, t);
+    makeSpheres(spheres, sphereCount);
     updateSphereUniforms(shaderProgram, spheres, sphereCount);
 }

src/sphere.c

@@ -1,7 +1,9 @@
 #include "sphere.h"
 
-void makeSpheres(struct Sphere *spheres, int count, float t)
+void makeSpheres(struct Sphere *spheres, int count)
 {
+    //float t = time();
+
     vec3 albedos[] = 
     {
         {0.0,0.0,1.0},
@@ -13,33 +15,21 @@ void makeSpheres(struct Sphere *spheres, int count, float t)
         {1.0,1.0,1.0}
     };
 
-    vec3 sc = {0.0,0.0,1.0};
-
-    int sphereIdx = 0;
-
-    // a shiny sphere in the middle which bobs around
-    int middleSpheres = 1;
-    sc[0] = 2.5*sin(t*2.67);
-    sc[1] = 1.5+sin(t*1.09);
-    sc[2] = 2.5*cos(-t*2.13);
-    float radius = 1.7;
-    vec3 col = {0.0,0.0,0.0};
-    spheres[sphereIdx++] = makeSphere(sc,radius,col,1);
-
-    // rainbow spheres
-    int rainbowSpheres = count - middleSpheres;
     // distance from center
     float d = 6.0;
-    radius = 0.5;
+    float radius = 0.5;
     float x;
-    for (int i = 0; i < rainbowSpheres; i++)
+    vec3 sc = {0.0,0.0,1.0};
+
+    float t = now();
+    for (int i = 0; i < count; i++)
     {
-        x = 2.0*CGLM_PI * (float)i/(float)rainbowSpheres;
+        x = 2.0*CGLM_PI * (float)i/(float)count;
         sc[0] = sin(x)*d;
         sc[1] = sin(x*3.0-5.0*sin(t));
         sc[2] = cos(x)*d;
 
-        float ic = i/(float)rainbowSpheres*CGLM_PI*2.0;
+        float ic = i/(float)count*CGLM_PI*2.0;
         float r = sin(ic);
         float g = sin(ic+CGLM_PI/1.5);
         float b = sin(ic+2.0*CGLM_PI/1.5);
@@ -48,8 +38,8 @@ void makeSpheres(struct Sphere *spheres, int count, float t)
         glm_vec3_scale(col, 0.5, col);
         glm_vec3_adds(col, 0.5, col);
 
-        int material = i % 2;
-        spheres[sphereIdx++] = makeSphere(sc,radius,col,material);
+        int material = i % 2 == 0;
+        spheres[i] = makeSphere(sc,radius,col,material);
     }
 }
 

src/sphere.h

@@ -16,7 +16,7 @@ struct Sphere
     int material;
 };
 
-void makeSpheres(struct Sphere *spheres, int count, float t);
+void makeSpheres(struct Sphere *spheres, int count);
 struct Sphere makeSphere(vec3 center, float radius, vec3 albedo, int material);
 //void updateSphereUniform(GLuint shaderProgram, struct Sphere sphere);
 void updateSphereUniforms(GLuint shaderProgram, struct Sphere *spheres, int count);