#ifdef GL_ES
precision highp float;
#endif

#extension GL_OES_standard_derivatives : enable

uniform float time;
uniform vec2 mouse;
uniform vec2 resolution;

vec2 rand2n(vec2 co, float sampleIndex) {
   vec2 seed = co * (sampleIndex + 1.0);
       seed+=vec2(-1,1);
   return vec2(fract(sin(dot(seed.xy ,vec2(12.9898,78.233))) * 43758.5453),
               fract(cos(dot(seed.xy ,vec2(4.898,7.23))) * 23421.631));
}

vec2 circleInvert(vec2 pos, vec3 circle){
       return ((pos - circle.xy) * circle.z * circle.z)/(length(pos - circle.xy) * length(pos - circle.xy) ) + circle.xy;
}

mat2 getRotationMat2(float angleRadians){
       return mat2(cos(angleRadians), -sin(angleRadians),
               sin(angleRadians), cos(angleRadians));
}
const mat2 UNI_MAT = mat2(1, 0, 0, 1);

const float PI = 3.1415926535;
const float PI_2 = 3.1415926535/2.;
const float PI_4 = 3.1415926535/4.;


const float R = 100. * sqrt(2.) / 2.;
vec3 c1 = vec3(100, 0, 0);
vec3 c3 = vec3(-100, 0, 0);
vec3 c2 = vec3(0, 100, 0);
vec3 c4 = vec3(0, -100, 0);
vec3 c5 = vec3(0);

// control point and rotation matrix of l1
float l1Angle = PI_4;
vec2 l1p = vec2(500, 0);
mat2 l1m;
mat2 l1mInv;
// control point and rotation matrix of l1
float l2Angle = PI + PI_4;
vec2 l2p = vec2(-800, 0);
mat2 l2m;
mat2 l2mInv;

float l3Angle = PI + PI_4;
vec2 l3p = vec2(-1500, 0);
mat2 l3m;
mat2 l3mInv;

float l4Angle = PI_2 + PI_4;
vec2 l4p = vec2(0, -1500);
mat2 l4m;
mat2 l4mInv;

bool enableL1 = false;
bool enableL2 = false;
bool enableL3 = false;
bool enableL4 = false;

const int MAX_ITERATIONS = 35;
float IIS(vec2 pos){
   float loopNum = 0.;
       bool cont = false;
       for(int i = 0 ; i < MAX_ITERATIONS ; i++){
               cont = false;

       if(length(pos - c1.xy) < c1.z){
                       pos = circleInvert(pos, c1);
                       cont = true;
           loopNum++;
               }else if(length(pos - c2.xy) < c2.z){
                       pos = circleInvert(pos, c2);
                       cont = true;
           loopNum++;
               }else if(length(pos - c3.xy) < c3.z){
                       pos = circleInvert(pos, c3);
                       cont = true;
           loopNum++;
               }else if(length(pos - c4.xy) < c4.z){
                       pos = circleInvert(pos, c4);
                       cont = true;
           loopNum++;
               }else if(length(pos - c5.xy) < c5.z){
                       pos = circleInvert(pos, c5);
                       cont = true;
           loopNum++;
               }

       if(enableL1){
               pos -= l1p;
               pos = l1mInv * pos;
               if(pos.x > 0.){
                       pos.x *= -1.;
               loopNum++;
               cont = true;
               }
               pos = l1m * pos;
               pos += l1p;
       }

       if(enableL2){
       pos -= l2p;
       pos = l2mInv * pos;
       if(pos.x > 0.){
               pos.x *= -1.;
           loopNum++;
           cont = true;
       }
       pos = l2m * pos;
       pos += l2p;
       }

       if(enableL3){
       pos -= l3p;
       pos = l3mInv * pos;
       if(pos.x > 0.){
               pos.x *= -1.;
           loopNum++;
           cont = true;
       }
       pos = l3m * pos;
       pos += l3p;
       }

       if(enableL4){
               pos -= l4p;
               pos = l4mInv * pos;
               if(pos.x > 0.){
                       pos.x *= -1.;
               loopNum++;
                               cont = true;
               }
               pos = l4m * pos;
               pos += l4p;
       }

               if(cont == false) break;
       }

       return loopNum;
}

vec3 hsv2rgb(vec3 c)
{
   vec4 K = vec4(1.0, 2.0 / 3.0, 1.0 / 3.0, 3.0);
   vec3 p = abs(fract(c.xxx + K.xyz) * 6.0 - K.www);
   return c.z * mix(K.xxx, clamp(p - K.xxx, 0.0, 1.0), c.y);
}

//w: start time
//s: duration
float scene(in float t, in float w, in float s){
   return clamp(t - w, 0.0, s) / s;
}

float expEzingIn(float t){
   return pow( 2., 13. * (t - 1.) );
}
float expEzingOut(float t) {
       return -pow( 2., -10. * t) + 1.;
}

float circEzingInOut(float t){
       t /= .5;
       if (t < 1.) return -.5 * (sqrt(1. - t*t) - 1.);
       t -= 2.;
       return .5 * (sqrt(1. - t*t) + 1.);
}

float circEzingIn(float t){
       return -  (sqrt(1. - t*t) - 1.);
}

float ezing(in float t){
   //return pow( 2., 13. * (t - 1.) );
   //return t;
   //return t * t * t;
   //return t * t * (3.0 - 2.0 * t);
   //return sin(t * radians(90.0));
   //return 1.0 - pow(cos(t * radians(90.0)), 0.5);
   return pow(sin(t * radians(90.0)), 0.2);
}

const float DISPLAY_GAMMA_COEFF = 1. / 2.2;
vec3 gammaCorrect(vec3 rgb) {
 return vec3((min(pow(rgb.r, DISPLAY_GAMMA_COEFF), 1.)),
             (min(pow(rgb.g, DISPLAY_GAMMA_COEFF), 1.)),
             (min(pow(rgb.b, DISPLAY_GAMMA_COEFF), 1.)));
}

const float SAMPLE_NUM = 20.;
void main(){
   vec3 sum = vec3(0);
   float t = mod(time, 22.);
       float ratio = resolution.x / resolution.y / 2.0;

   enableL1 = enableL2 = enableL3 = enableL4 = false;

       float start = 1.5;
  c1.z = mix(0., R, expEzingIn(scene(t, start, 1.)));
   c2.z = mix(0., R, expEzingIn(scene(t, start + .1, 1.)));
       c3.z = mix(0., R, expEzingIn(scene(t, start + .3, 1.)));
       c4.z = mix(0., R, expEzingIn(scene(t, start + .5, 1.)));

   float rotationStart = start + 1.7;
   float scaleFactor = 150.;
   float theta = mix(0., 2. * PI + PI_4, circEzingInOut(scene(t, rotationStart, 1.)));
   scaleFactor += mix(0., 800., circEzingInOut(scene(t, rotationStart, 1.)));

   //1.8---3.
       l1p.x = mix(1500., 200., expEzingIn(scene(t, rotationStart, 1.)));
       l2p.x = mix(-1500., -200., expEzingIn(scene(t, rotationStart, 1.5)));
   enableL1 = t > rotationStart + .5;
   enableL2 = t > rotationStart + .5;

       float lineStart = rotationStart + 1.5;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, lineStart, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, lineStart, 1.)));
   float lpRotateAngle1 = mix(0., PI_2, circEzingInOut(scene(t, lineStart, 1.)));
       l3p.x = mix(-1500., -200., expEzingIn(scene(t, lineStart, 1.)));
   enableL3 = t > lineStart;

   float rotateStart2 = lineStart + 1.1;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart2, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart2, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart2, 1.)));
       lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart2, 1.)));
   l4p.y = mix(-1500., -200., expEzingIn(scene(t, rotateStart2, 1.)));
   enableL4 = t > rotateStart2;

   float rotateStart3 = rotateStart2 + 1.1;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart3, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart3, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart3, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart3, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart3, 1.)));
       float lpRotateAngle2 = mix(0., PI_2, circEzingInOut(scene(t, rotateStart3, 1.)));

   float rotateStart4 = rotateStart3 + 1.3;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart4, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart4, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart4, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart4, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart4, 1.)));
       lpRotateAngle2 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart4, 1.)));
   scaleFactor += mix(0., -600., circEzingInOut(scene(t, rotateStart4, 1.)));

   float rotateStart5 = rotateStart4 + 1.1;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart5, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart5, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart5, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart5, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart5, 1.)));
       lpRotateAngle2 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart5, 1.)));

   float rotateStart6 = rotateStart5 + 1.3;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart6, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart6, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart6, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart6, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart6, 1.)));
       lpRotateAngle2 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart6, 1.)));

   float rotateStart7 = rotateStart6 + 1.3;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart7, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart7, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart7, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart7, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart7, 1.)));
       lpRotateAngle2 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart7, 1.)));

   float rotateStart8 = rotateStart7 + 1.3;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart8, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart8, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart8, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart8, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart8, 1.)));
       lpRotateAngle2 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart8, 1.)));

   float rotateStart9 = rotateStart8 + 1.3;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart9, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart9, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart9, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart9, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart9, 1.)));
       lpRotateAngle2 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart9, 1.)));

   float rotateStart10 = rotateStart9 + 1.2;
   l1Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart10, 1.)));
   l2Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart10, 1.)));
       l3Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart10, 1.)));
       l4Angle += mix(0., PI_2 , circEzingInOut(scene(t, rotateStart10, 1.)));
   lpRotateAngle1 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart10, 1.)));
       lpRotateAngle2 += mix(0., PI_2, circEzingInOut(scene(t, rotateStart10, 1.)));


   float circleStart = rotateStart5;
   float deformR = 100.;
   c1.z += mix(0., deformR, expEzingOut(scene(t, rotateStart5, 1.)));
   c1.z += mix(0., -deformR, expEzingOut(scene(t, rotateStart6, 1.)));

   //c2.z += mix(0., deformR, expEzingOut(scene(t, circleStart+.5, 1.)));
       //c2.z += mix(0., -deformR, expEzingOut(scene(t, circleStart+.5 +1., 1.)));

       c3.z += mix(0., deformR, expEzingOut(scene(t, rotateStart6, 1.)));
       c3.z += mix(0., -deformR, expEzingOut(scene(t, rotateStart7, 1.)));

       //c4.z += mix(0., deformR, expEzingOut(scene(t, circleStart+1.5, 1.)));
       //c4.z += mix(0., -deformR, expEzingOut(scene(t, circleStart+1.5+1., 1.)));

   c5.z += mix(0., 50., circEzingInOut(scene(t, rotateStart7, 1.)));

   scaleFactor += mix(0., -300., circEzingInOut(scene(t, rotateStart7, 1.)));
   scaleFactor /= mix(1., 4., circEzingInOut(scene(t, rotateStart8, 1.)));
   scaleFactor /= mix(1., 4., circEzingInOut(scene(t, rotateStart9, 1.)));
   scaleFactor /= mix(1., 4., circEzingInOut(scene(t, rotateStart10, 1.)));

   float endingStart = rotateStart10 + 1.3;
   if(t >= endingStart){

       float rotateStart11 = endingStart;
               l1Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart11, 1.)));
       l2Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart11, 1.)));
               l3Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart11, 1.)));
               l4Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart11, 1.)));
       lpRotateAngle1 -= mix(0., PI_2, circEzingInOut(scene(t, rotateStart11, 1.)));
               lpRotateAngle2 -= mix(0., PI_2, circEzingInOut(scene(t, rotateStart11, 1.)));

       float rotateStart12 = rotateStart11 + 1.3;
               l1Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart12, 1.)));
       l2Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart12, 1.)));
               l3Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart12, 1.)));
               l4Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart12, 1.)));
       lpRotateAngle1 -= mix(0., PI_2, circEzingInOut(scene(t, rotateStart12, 1.)));
               lpRotateAngle2 -= mix(0., PI_2, circEzingInOut(scene(t, rotateStart12, 1.)));

       float rotateStart13 = rotateStart12+1.3;
               l1Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart13, 1.)));
       l2Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart13, 1.)));
               l3Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart13, 1.)));
               l4Angle -= mix(0., PI_2 , circEzingInOut(scene(t, rotateStart13, 1.)));
       lpRotateAngle1 -= mix(0., PI_2, circEzingInOut(scene(t, rotateStart13, 1.)));
               lpRotateAngle2 -= mix(0., PI_2, circEzingInOut(scene(t, rotateStart13, 1.)));


       float endingTime = 1.;

               c5.z  = mix(50., 0., circEzingInOut(scene(t, rotateStart12, endingTime)));
       l1p.x = mix(200., 5000., circEzingIn(scene(t, rotateStart13, endingTime)));
               l2p.x = mix(-200., -5000., circEzingIn(scene(t, rotateStart13, endingTime)));
       l3p.y = mix(0., -5000., circEzingIn(scene(t, rotateStart13, endingTime)));
               l4p.y = mix(-200., -5000., circEzingIn(scene(t, rotateStart13, endingTime)));
       c1.z = mix(R, 0., circEzingInOut(scene(t, rotateStart13, endingTime)));
       c2.z = mix(R, 0., circEzingInOut(scene(t, rotateStart13, endingTime)));
               c3.z = mix(R, 0., circEzingInOut(scene(t, rotateStart13, endingTime)));
               c4.z = mix(R, 0., circEzingInOut(scene(t, rotateStart13, endingTime)));

       theta += mix(0., -PI_2, circEzingInOut(scene(t, rotateStart13, 1.)));

       scaleFactor *= mix(1., 650., circEzingInOut(scene(t, endingStart, 1.)));

   }

   mat2 lpRotate1 = getRotationMat2(lpRotateAngle1);
   l1m = getRotationMat2(l1Angle);
   l1mInv = getRotationMat2(-l1Angle);
   l2m = getRotationMat2(l2Angle);
   l2mInv = getRotationMat2(-l2Angle);
   l3m = getRotationMat2(l3Angle);
   l3mInv = getRotationMat2(-l3Angle);
       l4m = getRotationMat2(l4Angle);
   l4mInv = getRotationMat2(-l4Angle);

   l1p = lpRotate1 * l1p;
       l2p = lpRotate1 * l2p;
       l3p = lpRotate1 * l3p;

   mat2 lpRotate2 = getRotationMat2(lpRotateAngle2);
       l4p = lpRotate2 * l4p;

   mat2 m = getRotationMat2(theta);

   for(float i = 0. ; i < SAMPLE_NUM ; i++){
       vec2 position = ( (gl_FragCoord.xy + rand2n(gl_FragCoord.xy, i)) / resolution.yy ) - vec2(ratio, 0.5);

       position = position * scaleFactor;
       position = m * position;


       float loopNum = IIS(position);

       if(loopNum >  0.){
           sum += hsv2rgb(vec3(0.01 + 0.05 * (loopNum-1.),1.0,1.0));
       }
   }
   gl_FragColor = vec4(gammaCorrect(sum/SAMPLE_NUM), 1.);
}

You are viewing proxied material from ftp.icm.edu.pl. The copyright of proxied material belongs to its original authors. Any comments or complaints in relation to proxied material should be directed to the original authors of the content concerned. Please see the disclaimer for more details.