Another Shader Toy sketch translated to JRubyArt

From time to time I find a nice Shadertoy sketch to convert to JRubyArt, I usually start using a default wrapper (simply inserting the shader code into the default wrapper see here. But often is quite easy to unwrap the shaderoy.glsl to conform as a processing shader, as I have done here with a sketch by Victor Shepardson (aka @victorrileys twitter).

attr_reader :previous_time, :wrapper, :start
# Shader by Victor Shepardson

def settings
  size(800, 450, P2D)
end

def setup
  sketch_title 'Sigmoids n sine (by Victor Shepardson)'
  @previous_time = 0.0
  @mouse_dragged = false
  @mouse_click_state = 0.0
  @wrapper = load_shader(data_path('sigmoids.glsl'))
  # Assume the dimension of the window will not change over time
  wrapper.set('iResolution', width.to_f, height.to_f, 0.0)
  @last_mouse_position = Vec2D.new(mouse_x, mouse_y)
  @start = java.lang.System.current_time_millis
end

def playback_time_seconds
  (java.lang.System.current_time_millis - start) / 1000.0
end

def render_time
  playback_time_seconds - previous_time
end

def draw
  wrapper.set('iTime', playback_time_seconds)
  @previous_time = playback_time_seconds
  shader(wrapper)
  # Draw the output of the shader onto a rectangle that covers the whole viewport.
  rect(0, 0, width, height)
end

Shader

#ifdef GL_ES
precision mediump float;
precision mediump int;
#endif

// From a shadertoy sketch by Victor Shepardson
// https://www.shadertoy.com/view/ltfXzj

uniform vec3      iResolution;           // viewport resolution (in pixels)
uniform float     iTime;                 // shader playback time (in seconds)

const float pi = 3.14159;

float sigmoid(float x){
 	return x/(1.+abs(x));
}

float iter(vec2 p, vec4 a, vec4 wt, vec4 ws, float t, float m, float stereo){
    float wp = .2;
    vec4 phase = vec4(mod(t, wp), mod(t+wp*.25, wp), mod(t+wp*.5, wp), mod(t+wp*.75, wp))/wp;
    float zoom = 1./(1.+.5*(p.x*p.x+p.y*p.y));
    vec4 scale = zoom*pow(vec4(2.), -4.*phase);
    vec4 ms = .5-.5*cos(2.*pi*phase);
    vec4 pan = stereo/scale*(1.-phase)*(1.-phase);
    vec4 v = ms*sin( wt*(t+m) + (m+ws*scale)*((p.x+pan) * cos((t+m)*a) + p.y * sin((t+m)*a)));
    return sigmoid(v.x+v.y+v.z+v.w+m);
}

vec3 scene(float gt, vec2 uv, vec4 a0, vec4 wt0, vec4 ws0, float blur){
    //time modulation
    float tm = mod(.0411*gt, 1.);
    tm = sin(2.*pi*tm*tm);
    float t = (.04*gt + .05*tm);

    float stereo = 1.*(sigmoid(2.*(sin(1.325*t*cos(.5*t))+sin(-.7*t*sin(.77*t)))));//+sin(-17.*t)+sin(10.*t))));
    //t = 0.;
    //also apply spatial offset
    uv+= .5*sin(.33*t)*vec2(cos(t), sin(t));

    //wildly iterate and divide
    float p0 = iter(uv, a0, wt0, ws0, t, 0., stereo);

   	float p1 = iter(uv, a0, wt0, ws0, t, p0, stereo);

    float p2 = sigmoid(p0/(p1+blur));

    float p3 = iter(uv, a0, wt0, ws0, t, p2, stereo);

    float p4 = sigmoid(p3/(p2+blur));

    float p5 = iter(uv, a0, wt0, ws0, t, p4, stereo);

    float p6 = sigmoid(p4/(p5+blur));

    float p7 = iter(uv, a0, wt0, ws0, t, p6, stereo);

    float p8 = sigmoid(p4/(p2+blur));

    float p9 = sigmoid(p8/(p7+blur));

    float p10 = iter(uv, a0, wt0, ws0, t, p8, stereo);

    float p11 = iter(uv, a0, wt0, ws0, t, p9, stereo);

    float p12 = sigmoid(p11/(p10+blur));

    float p13 = iter(uv, a0, wt0, ws0, t, p12, stereo);

    //colors
    vec3 accent_color = vec3(1.,0.2,0.);//vec3(0.99,0.5,0.2);
    /*float r = sigmoid(-1.+2.*p0+p1-max(1.*p3,0.)+p5+p7+p10+p11+p13);
    float g = sigmoid(-1.+2.*p0-max(1.*p1,0.)-max(2.*p3,0.)-max(2.*p5,0.)+p7+p10+p11+p13);
    float b = sigmoid(0.+1.5*p0+p1+p3+-max(2.*p5,0.)+p7+p10+p11+p13);
    */
    float r = sigmoid(p0+p1+p5+p7+p10+p11+p13);
    float g = sigmoid(p0-p1+p3+p7+p10+p11);
    float b = sigmoid(p0+p1+p3+p5+p11+p13);


    vec3 c = max(vec3(0.), .4+.6*vec3(r,g,b));

    float eps = .4;
    float canary = min(abs(p1), abs(p2));
    canary = min(canary, abs(p5));
    //canary = min(canary, abs(p6));
    canary = min(canary, abs(p7));
    canary = min(canary, abs(p10));
    float m = max(0.,eps-canary)/eps;
    m = sigmoid((m-.5)*700./(1.+10.*blur))*.5+.5;
    //m = m*m*m*m*m*m*m*m*m*m;
    vec3 m3 = m*(1.-accent_color);
    c *= .8*(1.-m3)+.3;//mix(c, vec3(0.), m);

    return c;
}

void main( void )
{
    float s = min(iResolution.x, iResolution.y);
   	vec2 uv = (2.*gl_FragCoord.xy - vec2(iResolution.xy)) / s;

    float blur = .5*(uv.x*uv.x+uv.y*uv.y);

    //angular, spatial and temporal frequencies
    vec4 a0 = pi*vec4(.1, -.11, .111, -.1111);
    vec4 wt0 = 2.*pi*vec4(.3);//.3333, .333, .33, .3);
    vec4 ws0 = 2.5*vec4(11., 13., 11., 5.);

    //aa and motion blur
    float mb = 1.;
    float t = 1100.+iTime;
    vec3 c = scene(t, uv, a0, wt0, ws0, blur)
        + scene(t-mb*.00185, uv+(1.+blur)*vec2(.66/s, 0.), a0, wt0, ws0, blur)
        + scene(t-mb*.00370, uv+(1.+blur)*vec2(-.66/s, 0.), a0, wt0, ws0, blur)
        + scene(t-mb*.00555, uv+(1.+blur)*vec2(0., .66/s), a0, wt0, ws0, blur)
        + scene(t-mb*.00741, uv+(1.+blur)*vec2(0., -.66/s), a0, wt0, ws0, blur)
        + scene(t-mb*.00926, uv+(1.+blur)*vec2(.5/s, .5/s), a0, wt0, ws0, blur)
        + scene(t-mb*.01111, uv+(1.+blur)*vec2(-.5/s, .5/s), a0, wt0, ws0, blur)
        + scene(t-mb*.01296, uv+(1.+blur)*vec2(-.5/s, -.5/s), a0, wt0, ws0, blur)
        + scene(t-mb*.01481, uv+(1.+blur)*vec2(.5/s, -.5/s), a0, wt0, ws0, blur)

        ;
    c/=9.;

    gl_FragColor = vec4(c,1.0);
}

See shadertoy sketch. PS: my converted shader should work for vanilla processing and py5.