fixed the tunnel precalc to produce a continuous 10 bit depth value
[gbajam21] / tools / tungen.c
index f480758..d1ff400 100644 (file)
@@ -1,10 +1,11 @@
 #include <stdio.h>
 #include <stdlib.h>
+#include <string.h>
 #include <stdint.h>
 #include <math.h>
 
 struct vec2 {
-       float x, y;
+       int x, y;
 };
 
 int main(int argc, char **argv)
@@ -89,6 +90,8 @@ invalopt:     fprintf(stderr, "invalid argument: %s\n", argv[i]);
        for(frm=0; frm<num_frames; frm++) {
                int coffs = num_frames > 1 ? frm * center / (num_frames - 1) : center;
 
+               memset(buf, 0xff, xsz * ysz * sizeof *buf);
+
 #define UDIV   2048
 #define VDIV   32768
                prev_r = 0.0f;
@@ -99,10 +102,10 @@ invalopt:  fprintf(stderr, "invalid argument: %s\n", argv[i]);
                        float z = v * coffs;
 
                        /* don't bother drawing rings < 1 pixel apart */
-                       if(fabs(r - prev_r) < 0.05) continue;
+                       if(r < 0 || fabs(r - prev_r) < 0.05) continue;
 
                        for(j=0; j<UDIV; j++) {
-                               float u = (float)j / (float)UDIV;
+                               float u = (float)j / (float)(UDIV - 1);
                                float theta = 2.0f * u * M_PI;
 
                                int x = (int)(cos(theta) * r - z) + xsz / 2;
@@ -110,8 +113,8 @@ invalopt:   fprintf(stderr, "invalid argument: %s\n", argv[i]);
 
                                if(x >= 0 && x < xsz && y >= 0 && y < ysz) {
                                        ptr = buf + y * xsz + x;
-                                       ptr->x = u;
-                                       ptr->y = v * 8;
+                                       ptr->x = (j << 8) / UDIV;
+                                       ptr->y = ((VDIV - i) << 11) / VDIV;
                                }
                        }
                        prev_r = r;
@@ -120,26 +123,22 @@ invalopt: fprintf(stderr, "invalid argument: %s\n", argv[i]);
                ptr = buf;
                for(i=0; i<out_nlines; i++) {
                        for(j=0; j<xsz; j++) {
-                               float u = ptr->x;
-                               float v = ptr->y;
-                               int r = (int)(u * 8.0 * 255.0f) & 0xff;
-                               int g = (int)(v * 8.0 * 255.0f) & 0xff;
-                               int b = (~(int)(v * 0.5 * 255.0f) & 0xff) + 105;
-                               if(b > 255) b = 255;
-                               if(b < 0) b = 0;
+                               int u = ptr->x;
+                               int v = ptr->y;
+                               int r = (u << 3) & 0xff;
+                               int g = (v >> 3) & 0xff;
 
                                /*if(v > 2.0) r = g = b = 0;*/
 
                                ptr++;
 
-                               uint16_t out = ((uint16_t)(u * 255.0f) & 0xff) |
-                                       (((uint16_t)(v * 255.0f) & 0xff) << 8);
+                               uint16_t out = ((uint16_t)u & 0x3f) | (((uint16_t)v & 0x3ff) << 6);
                                fwrite(&out, sizeof out, 1, stdout);
 
                                if(fp) {
                                        fputc(r, fp);
                                        fputc(g, fp);
-                                       fputc(b, fp);
+                                       fputc(0, fp);
                                }
                        }
                }