diff --git a/src/shaders/fragment.frag b/src/shaders/fragment.frag
index adf3e97..0c4596b 100644
--- a/src/shaders/fragment.frag
+++ b/src/shaders/fragment.frag
@@ -57,18 +57,69 @@ float noise(vec2 p, float scale, int a)
                               hash( i + vec2(1.), a), u.x), u.y);
 }
 
-// Hexagonal prism, circumcircle variant
+float sdHexPrism( vec3 p, vec2 h )
+{
+  const vec3 k = vec3(-0.8660254, 0.5, 0.57735);
+  p = abs(p);
+  p.xy -= 2.0*min(dot(k.xy, p.xy), 0.0)*k.xy;
+  vec2 d = vec2(
+       length(p.xy-vec2(clamp(p.x,-k.z*h.x,k.z*h.x), h.x))*sign(p.y-h.x),
+       p.z-h.y );
+  return min(max(d.x,d.y),0.0) + length(max(d,0.0));
+}
+
 float fHexagonCircumcircle(vec3 p, vec2 h) {
     vec3 q = abs(p);
     return max(q.y - h.y, max(q.x * sqrt(3.) * 0.5 + q.z * 0.5, q.z) - h.x);
-	//this is mathematically equivalent to this line, but less efficient:
-	//return max(q.y - h.y, max(dot(vec2(cos(PI/3), sin(PI/3)), q.zx), q.z) - h.x);
 }
 
-float sdHex(vec3 pos, vec2 id, float angle) {
-    pos.xz *= rot2D(0.1);
-    float d1 = fHexagonCircumcircle(pos, vec2(0.86, 1.));
-    return d1;
+float hexPylon(vec3 p, vec2 h) {//float r, float ht){
+
+    //vec3 p = vec3(p.x, p.z, p2.y);
+    vec3 b = vec3(h.x, h.y, h.x);
+    
+    // Hexagon.
+    p.xz = abs(p.xz);
+    p.xz = vec2(p.x*.866025 + p.z*.5, p.z);
+    // The ".015" is a subtle rounding factor. Zero gives sharp edges,
+    // and larger numbers give a more rounded look.
+  	return length(max(abs(p) - b + .015, 0.)) - .015;
+}
+
+#define zclamp(a) max(a,0.0) //Clamp negative values at zero
+float DF_RoundedHex( vec3 p, vec2 h) //float width, float height)
+{
+    float width = h.x;
+    float height = h.y;
+    //Modified version (smooth edges) of the exagon prism found here:
+    //https://iquilezles.org/articles/distfunctions
+    float smoothRadius = 0.05;
+    width -= smoothRadius*2.0;
+    
+    //Hexagon prism constructed using X,Y,Z symmetry.
+    //Only quadrant 1 needs to be solved, but the joining diagonal to quadrant IV is also
+    //required for distance blending (see db).
+    p = abs(p);
+    
+    //Hexagonal edge distances :
+    //Note : [.8666,0.5] = [sin(PI/3,cos(PI/3)] -> Hexagon edges rotation coeff (60 degrees).
+    float da = (p.x*0.866025+p.z*0.5)-width; //quadrant I diagonal edge distance
+    float db = (p.x*0.866025-p.z*0.5)-width; //quadrant IV diagonal edge distance (needed for blending)
+    float dc = p.z-width; //upper distance
+    
+    vec3 d = zclamp(vec3(da,db,dc));
+    //Note: this is not an euclidian length, therefore this operation slightly distorts our distance field.
+    //Yet, it is harmless to convergence, and does the smoothing job quite well.
+    float dw = length(d)-smoothRadius; //hexagonal part smoothness (blending at 60 deg)
+    float dh = p.y-height;
+    
+    //Now that we have xz distance(dw) and y distance (dh), we can compute the distance 
+    //for the given isovalue (the smoothing radius).
+    //Note : internal distance (maxX,maxY,maxZ) is also used to genereate internal signed dist,
+    //       helping convergence when overstepping (very frequent with domain repetition).
+    float externalDistance = length(zclamp(vec2(dh,dw)))-smoothRadius; //Smoothed, unsigned
+	float internalDistance = max(max(da,dc),dh); //Sharp, signed.
+    return min(externalDistance,internalDistance);
 }
 
 float getScaledFFT(int index, float scale, float offset) {
@@ -135,16 +186,17 @@ vec3 mapScene(vec3 p) {
 
     // Use axial coordinates as a stable hex ID
     float distFromCenter = hexDistance(hex.axial);
-    int fftIndex = int(clamp(distFromCenter * 1.0, 0.0, 511.0)); // tweak 15.0 to taste
+    int fftIndex = int(clamp(distFromCenter +1.0, 0.0, 511.0)); // tweak 15.0 to taste
 
     float fftVal = fft_output[fftIndex];
     float hexHeight = 1.0 + fftVal * 1.0;
 
     // Rotate individual hex tiles if needed
     vec3 r = hex.local;
+    //r.yz *= rot2D(1.0);
     r.xz *= rot2D(0.5);
 
-    float d1 = fHexagonCircumcircle(r, vec2(hexRadius, hexHeight));
+    float d1 = fHexagonCircumcircle(vec3(r.x,(r.y-hexHeight/2),r.z), vec2(hexRadius, hexHeight/2));
     d = min(d,d1);
 
     return vec3(d, 0.0, 0.0);
@@ -156,13 +208,27 @@ vec3 mapScene(vec3 p) {
 
 float rayMarch(vec3 ro, vec3 rd, int a) {
     vec3 d;
-    float t = 0.; // total distance travelled
+    float t = 0.,ad,tmax=100.; // total distance travelled
+    const float tolerance = 0.00001;
+    const float Z_REPEAT_DIST = 1.;
+
     // Raymarching
-    for (int i = 0; i < 100; i++) {
+    for (int i = 0; i < 80; i++) {
         d = mapScene(ro + rd * t); // Get distance to objects
+        ad = abs(d.x);
+        if (ad < tolerance*(t*0.125 + 1.0) || t > tmax) break;
         t += d.x; // "march" the ray
-        if (d.x < 1e-3 || t > 500.) break;
     }
+    t -= Z_REPEAT_DIST*15.;
+    
+    for( int i=0; i<80; i++ )
+    {
+        d = mapScene(ro + rd * t); // get distance to objects
+        ad = abs(d.x);
+        if (ad < tolerance*(t*0.00125) || t > tmax) break;
+        t += min(d.x, Z_REPEAT_DIST/5.0); // "march" the ray
+    }
+    if (ad >= tmax) t= - 1.0;
     return t;
 }
 
@@ -179,7 +245,7 @@ float getLight(vec3 p, vec3 lightPos, float intensity, float shadow, vec3 n, flo
     vec3 l = no(lightPos - p);
     float len = length( lightPos - p ); // Distance from the light to the surface point.
     float dif = cl(dot(n, l)*intensity, 0., intensity) * 1.0 / (1.0 + atte*len),
-    d = rayMarch(p+n*.0025, l, 1);
+    d = rayMarch(p+n*.025, l, 1);
     if(d<length(lightPos-p)) dif *= shadow;
     return dif;
 }
@@ -194,14 +260,6 @@ float specular(vec3 normal,vec3 lightPos,vec3 rayOrigin,float specular) {
     return pow(max(dot(reflect(rayOrigin,normal),lightPos),0.0),specular) * nrm;
 }
 
-vec3 getSeaColor(vec3 p, vec3 n, vec3 l, vec3 eye) {
-
-    vec3 color = vec3(0.0, 0.1, 0.3) + diffuse(n,l,60.0) * vec3(0.11, 0.16, 0.18) * 0.3;
-    color -= vec3(0.73, 0.15, 0.66) * pow(cl(1.-dot(n, -eye), 0., 1.), .7);
-    color += vec3(specular(n,l,eye,60.0))*0.2;
-    return color;
-}
-
 
 vec3 applyFog(vec3 col, float t, vec3 rd, vec3  lightDir, float b ) {
     vec3 fogColor = mix( vec3(0.34, 0.11, 0.34), // blue
@@ -280,37 +338,10 @@ vec3 sceneGate(vec2 uv)
 
         if(mat==0.)
             col *= vec3(0.2, 0.3, 0.3) + addSpecular(n,rd,.5, 2.);
-        if(mat==1.)
-            col *= getSeaColor(p, n, no(vec3(0.0,0.3,0.8)),no(rd));
-        if(mat==2.)
-            col *= vec3(0.7, 0.7, 0.4) + noise (p.xz*3.+1.5, 0.1,0);
-        if(mat==3.)
-            col *= vec3(.8, .8, .5) + noise(p.xz*500.+1e5, .3,0);
-        if(mat==4.)
-            col *= vec3(0.0, 0.08, 0.11) + addSpecular(n,rd,0.3,0.7);
-        if(mat ==5.)
-           col = vec3(0.0, 0.08, 0.11);// * pow(cl(1. -dot(n, -rd), 0., 1.), .3);
-        if(mat == 6.)
-            col *= vec3(0.01, 0.04, 0.06) + addSpecular(n,rd,.1, 2.5);
    }
    return postProcess(applyFog(col, d, rd, vec3(0., -.1, -1.), .01));
 }
 
 void main() {
-    // Determine FFT bin index for current x position
-//    vec2 uv = (gl_FragCoord.xy * 2. - u_resolution.xy) / u_resolution.y;
-//    int index = int(floor((1.0+uv.x)*128.0));
-//    index = clamp(index, 0, 255);
-//
-//    // Get the FFT energy (clamped to avoid NaNs or overflow)
-//    float energy = clamp(texture(u_fft_texture, vec2(index/255., 0.0)).r, 0.0, 1.0);
-//
-//    float bar_height = energy;
-//    float fade = smoothstep(bar_height, bar_height + 0.02, 1.0 - uv.y);
-//
-//    vec3 color = vec3(fade);
-//
-//    o = vec4(color, 1.0);
-
     o = vec4(sceneGate( (gl_FragCoord.xy * 2. - u_resolution.xy) / u_resolution.y), 1.);
 }
\ No newline at end of file
diff --git a/src/shaders/fragment.inl b/src/shaders/fragment.inl
index 5130d24..97eb118 100644
--- a/src/shaders/fragment.inl
+++ b/src/shaders/fragment.inl
@@ -2,15 +2,15 @@
 #ifndef FRAGMENT_INL_
 # define FRAGMENT_INL_
 # define VAR_fft_output "H"
-# define VAR_o "v"
+# define VAR_o "f"
 # define VAR_syncs "a"
 # define VAR_u_hexGridTex "l"
 
 const char *fragment_frag =
  "#version 460\n"
  "precision mediump float;"
- "out vec4 v;"
- "const float f=2.*acos(-1.),m=sqrt(5.)*.5+.5;"
+ "out vec4 f;"
+ "const float m=2.*acos(-1.),v=sqrt(5.)*.5+.5;"
  "layout(location=0)uniform float a[7];"
  "layout(location=8)uniform float H[512];"
  "uniform sampler2D l;"
@@ -29,128 +29,111 @@ const char *fragment_frag =
    "float v=sin(.5),f=cos(.5);"
    "return mat2(f,-v,v,f);"
  "}"
- "float s(vec2 v)"
- "{"
-   "v=50.*fract(v*.3183099);"
-   "return fract(v.x*v.y*(v.x+v.y));"
- "}"
- "float s(vec2 v,float x)"
- "{"
-   "vec2 f=floor(v);"
-   "v=fract(v);"
-   "v=v*v*(3.-2.*v);"
-   "return-x+x*mix(mix(s(f+vec2(0)),s(f+vec2(1,0)),v.x),mix(s(f+vec2(0,1)),s(f+vec2(1)),v.x),v.y);"
- "}"
- "float s(vec3 v,vec2 f)"
+ "float s(vec3 v,vec2 m)"
  "{"
    "v=abs(v);"
-   "return max(v.y-f.y,max(v.x*sqrt(3.)*.5+v.z*.5,v.z)-f.x);"
- "}"
+   "return max(v.y-m.y,max(v.x*sqrt(3.)*.5+v.z*.5,v.z)-m.x);"
+ "}\n"
+ "#define zclamp(a)max(a,0.0)\n"
  "struct HexData{vec3 local;vec2 axial;};"
  "HexData t(vec3 v)"
  "{"
-   "float f=sqrt(3.)/3.*v.x-1./3.*v.z,i=2./3.*v.z,x=round(f),m=round(i),l=round(-f-i),d=abs(x-f),s=abs(m-i);"
-   "f=abs(l+f+i);"
-   "if(d>s&&d>f)"
-     "x=-m-l;"
-   "else if(s>f)"
-     "m=-x-l;"
-   "f=sqrt(3.)*(x+m*.5);"
-   "i=1.5*m;"
+   "float f=sqrt(3.)/3.*v.x-1./3.*v.z,m=2./3.*v.z,x=round(f),l=round(m),a=round(-f-m),d=abs(x-f),p=abs(l-m);"
+   "f=abs(a+f+m);"
+   "if(d>p&&d>f)"
+     "x=-l-a;"
+   "else if(p>f)"
+     "l=-x-a;"
+   "f=sqrt(3.)*(x+l*.5);"
+   "m=1.5*l;"
    "HexData r;"
-   "r.local=v-vec3(f,0,i);"
-   "r.axial=vec2(x,m);"
+   "r.local=v-vec3(f,0,m);"
+   "r.axial=vec2(x,l);"
    "return r;"
  "}"
  "struct HexData{vec3 local;vec2 axial;};"
- "float t(vec2 v)"
+ "float s(vec2 v)"
  "{"
-   "float f=v.x,m=v.y;"
-   "return max(abs(f),max(abs(m),abs(-f-m)));"
+   "float m=v.x,f=v.y;"
+   "return max(abs(m),max(abs(f),abs(-m-f)));"
  "}"
- "vec3 x(vec3 v)"
+ "vec3 p(vec3 v)"
  "{"
    "float f=1e9;"
    "HexData m=t(vec3(v.x,v.y-10.,v.z));"
+   "float l=1.+H[int(clamp(s(m.axial)+1.,0.,511.))];"
    "v=m.local;"
    "v.xz*=s();"
-   "float x=s(v,vec2(.83,1.+H[int(clamp(t(m.axial),0.,511.))]));"
-   "f=min(f,x);"
+   "l=s(vec3(v.x,v.y-l/2,v.z),vec2(.83,l/2));"
+   "f=min(f,l);"
    "return vec3(f,0,0);"
  "}"
- "float s(vec3 v,vec3 f,int m)"
+ "float p(vec3 v,vec3 f,int m)"
  "{"
    "vec3 l;"
-   "float i=0.;"
-   "for(int r=0;r<100;r++)"
+   "float x=0.,r;"
+   "for(int m=0;m<80;m++)"
      "{"
-       "l=x(v+f*i);"
-       "i+=l.x;"
-       "if(l.x<.001||i>5e2)"
+       "l=p(v+f*x);"
+       "r=abs(l.x);"
+       "if(r<1e-5*(x*.125+1.)||x>1e2)"
          "break;"
+       "x+=l.x;"
      "}"
-   "return i;"
+   "x-=15.;"
+   "for(int m=0;m<80;m++)"
+     "{"
+       "l=p(v+f*x);"
+       "r=abs(l.x);"
+       "if(r<x*.00125*1e-5||x>1e2)"
+         "break;"
+       "x+=min(l.x,.2);"
+     "}"
+   "if(r>=1e2)"
+     "x=-1.;"
+   "return x;"
  "}"
- "vec3 w(vec3 v)"
+ "vec3 x(vec3 v)"
  "{"
-   "vec2 f=vec2(.01,0);"
-   "return s(x(v).x-vec3(x(v-f.xyy).x,x(v-f.yxy).x,x(v-f.yyx)));"
+   "vec2 m=vec2(.01,0);"
+   "return s(p(v).x-vec3(p(v-m.xyy).x,p(v-m.yxy).x,p(v-m.yyx)));"
  "}"
- "float s(vec3 v,vec3 f,float m,float l,vec3 x,float i)"
+ "float p(vec3 v,vec3 f,float m,float l,vec3 x,float y)"
  "{"
-   "vec3 y=s(f-v);"
-   "m=s(dot(x,y)*m,0.,m)/(1.+i*length(f-v));"
-   "if(s(v+x*.0025,y,1)<length(f-v))"
+   "vec3 a=s(f-v);"
+   "m=s(dot(x,a)*m,0.,m)/(1.+y*length(f-v));"
+   "if(p(v+x*.025,a,1)<length(f-v))"
      "m*=l;"
    "return m;"
  "}"
- "float s(vec3 v,vec3 f,vec3 m,float x)"
+ "float p(vec3 v,vec3 m,vec3 f)"
  "{"
-   "return pow(max(dot(reflect(m,v),f),0.),x)*((x+8.)/(acos(-1.)*8.));"
+   "float l=pow(10.,2.);"
+   "return pow(max(dot(reflect(f,v),m),0.),l)*((l+8.)/(acos(-1.)*8.));"
  "}"
- "vec3 s(vec3 v,vec3 f,vec3 m,vec3 x)"
- "{"
-   "return vec3(0,.1,.3)+pow(dot(f,m)*.4+.6,60.)*vec3(.11,.16,.18)*.3-vec3(.73,.15,.66)*pow(s(1.-dot(f,-x),0.,1.),.7)+vec3(s(f,m,x,60.))*.2;"
- "}"
- "vec3 s(vec2 v,vec3 f,float m)"
+ "vec3 p(vec2 v,vec3 f,float m)"
  "{"
    "f=s(vec3(0,0,-5)-f);"
-   "vec3 x=s(cross(vec3(0,1,0),f));"
-   "return s(f*m+v.x*x+v.y*cross(f,x));"
+   "vec3 l=s(cross(vec3(0,1,0),f));"
+   "return s(f*m+v.x*l+v.y*cross(f,l));"
  "}"
- "vec3 s(vec3 v,vec3 f,float m,float x)"
+ "vec3 p(vec2 v)"
  "{"
-   "return vec3(s(v,s(vec3(0,.3,.8)),s(f),pow(10.,x)))*m;"
- "}"
- "vec3 w(vec2 v)"
- "{"
-   "vec3 f=vec3(0,25,25),m=s(v,f,s(1.,2.,25.)),i=vec3(0);"
-   "float l=s(f,m,0),r=0.;"
-   "if(l<5e2)"
+   "vec3 f=vec3(0,25,25),m=p(v,f,s(1.,2.,25.)),l=vec3(0);"
+   "float a=p(f,m,0),r=0.;"
+   "if(a<5e2)"
      "{"
-       "vec3 v=f+m*l,y=w(v);"
-       "r=x(v).y;"
-       "i=i+vec3(.82,.5,.9)*s(v,vec3(10,15,25),1.,.2,y,1e-10)+vec3(.79,.66,.43)*s(v,vec3(4,2,-15),1.,1.,y,1e-10)+vec3(0,.06,.7)*s(v,vec3(0,0,5),s((d-29.)*1e2,0.,50.),0.,y,3.1)+vec3(.29,.28,.33)*s(dot(y,s(vec3(0,1,10))),0.,1.);"
+       "vec3 v=f+m*a,i=x(v);"
+       "r=p(v).y;"
+       "l=l+vec3(.82,.5,.9)*p(v,vec3(10,15,25),1.,.2,i,1e-10)+vec3(.79,.66,.43)*p(v,vec3(4,2,-15),1.,1.,i,1e-10)+vec3(0,.06,.7)*p(v,vec3(0,0,5),s((d-29.)*1e2,0.,50.),0.,i,3.1)+vec3(.29,.28,.33)*s(dot(i,s(vec3(0,1,10))),0.,1.);"
        "if(r==0.)"
-         "i*=vec3(.2,.3,.3)+s(y,m,.5,2.);"
-       "if(r==1.)"
-         "i*=s(v,y,s(vec3(0,.3,.8)),s(m));"
-       "if(r==2.)"
-         "i*=vec3(.7,.7,.4)+s(v.xz*3.+1.5,.1);"
-       "if(r==3.)"
-         "i*=vec3(.8,.8,.5)+s(v.xz*5e2+1e5,.3);"
-       "if(r==4.)"
-         "i*=vec3(0,.08,.11)+s(y,m,.3,.7);"
-       "if(r==5.)"
-         "i=vec3(0,.08,.11);"
-       "if(r==6.)"
-         "i*=vec3(.01,.04,.06)+s(y,m,.1,2.5);"
+         "l*=vec3(.2,.3,.3)+vec3(p(i,s(vec3(0,.3,.8)),s(m)))*.5;"
      "}"
-   "return smoothstep(0.,1.,pow(mix(i,mix(vec3(.34,.11,.34),vec3(.93,.37,.16),pow(max(dot(m,vec3(0,-.1,-1)),0.),8.)),1.-exp(-l*.01))*vec3(.9,.8,.7),vec3(.45)));"
+   "return smoothstep(0.,1.,pow(mix(l,mix(vec3(.34,.11,.34),vec3(.93,.37,.16),pow(max(dot(m,vec3(0,-.1,-1)),0.),8.)),1.-exp(-a*.01))*vec3(.9,.8,.7),vec3(.45)));"
  "}"
  "void main()"
  "{"
-   "v=vec4(w((gl_FragCoord.xy*2.-n.xy)/n.y),1);"
+   "f=vec4(p((gl_FragCoord.xy*2.-n.xy)/n.y),1);"
  "}";
 
 #endif // FRAGMENT_INL_