// Generated by Haxe 3.4.7 (function () { "use strict"; function $extend(from, fields) { function Inherit() {} Inherit.prototype = from; var proto = new Inherit(); for (var name in fields) proto[name] = fields[name]; if( fields.toString !== Object.prototype.toString ) proto.toString = fields.toString; return proto; } var Main = function() { this._amp = 300; this._mouseY = 0; this._mouseX = 0; this._count = 0; window.onload = $bind(this,this.initialize); }; Main.main = function() { new Main(); }; Main.prototype = { initialize: function(e) { var _gthis = this; this._scene = new THREE.Scene(); this._camera = new THREE.PerspectiveCamera(48,window.innerWidth / window.innerHeight,1,5000); this._camera.position.z = 100; this._camera.lookAt(new THREE.Vector3()); var s = 3; this._camera.lookAt(new THREE.Vector3()); this._renderer = new THREE.WebGLRenderer({ devicePixelRatio : 1}); this._renderer.setClearColor(new THREE.Color(16772812)); this._renderer.shadowMap.enabled = true; this._renderer.setSize(window.innerWidth,window.innerHeight); window.document.body.appendChild(this._renderer.domElement); this.directionalLight = new THREE.DirectionalLight(10070783,1); this.directionalLight.position.x = 50; this.directionalLight.position.y = 100; this.directionalLight.position.z = 30; this.directionalLight.castShadow = true; this.directionalLight.shadow.camera.near = 1; this.directionalLight.shadow.camera.far = 300; this.directionalLight.shadow.camera.right = 150; this.directionalLight.shadow.camera.left = -150; this.directionalLight.shadow.camera.top = 150; this.directionalLight.shadow.camera.bottom = -150; this.directionalLight.shadow.mapSize.width = 1024; this.directionalLight.shadow.mapSize.height = 1024; this._scene.add(this.directionalLight); var ambient = new THREE.AmbientLight(5592405); this._scene.add(ambient); this._wave = new Wave(); this._scene.add(this._wave); var mm = new THREE.MeshPhongMaterial({ color : 16777215}); mm.side = 2; this.p = new THREE.Mesh(new THREE.PlaneGeometry(200,200,1,1),mm); this.p.receiveShadow = true; var helper = new THREE.CameraHelper(this.directionalLight.shadow.camera); this._run(); this._onResize(); window.onresize = $bind(this,this._onResize); window.onwheel = function(e1) { _gthis._amp += e1.wheelDelta * 0.1; if(_gthis._amp < 1) { _gthis._amp = 1; } if(_gthis._amp > 1000) { _gthis._amp = 1000; } }; } ,_run: function() { if(this._wave != null) { this._wave.update(); } this.p.rotation.x += 0.005; this.p.rotation.y += 0.008; this.p.rotation.z += 0.002; this._camera.position.z = this._amp; this._renderer.render(this._scene,this._camera); window.requestAnimationFrame($bind(this,this._run)); } ,_onResize: function() { this._camera.aspect = window.innerWidth / window.innerHeight; this._camera.updateProjectionMatrix(); this._renderer.setSize(window.innerWidth,window.innerHeight); } }; var Three = function() { }; Three.requestAnimationFrame = function(f) { return window.requestAnimationFrame(f); }; Three.cancelAnimationFrame = function(f) { window.cancelAnimationFrame(id); }; var Wave = function() { this._count = 0; this._attenuation = 0.95; this._k = 0.08; this._mat = new mat_MyShaderMat(); this._mat.side = 0; this._geo = new THREE.PlaneGeometry(200,200.,99,99); this._velocity = []; var _g = 0; while(_g < 100) { var i = _g++; this._velocity[i] = []; var _g1 = 0; while(_g1 < 100) { var j = _g1++; this._velocity[i][j] = 0; this._geo.vertices[this._getIndex(i,j)].z = 2.0; } } this._geo.verticesNeedUpdate = true; THREE.Mesh.call(this,this._geo,this._mat); this.scale.z = 10; this.receiveShadow = true; this.castShadow = true; }; Wave.__super__ = THREE.Mesh; Wave.prototype = $extend(THREE.Mesh.prototype,{ update: function() { var pos = this._geo.vertices; var _g1 = 1; var _g = 99; while(_g1 < _g) { var i = _g1++; var _g3 = 1; var _g2 = 99; while(_g3 < _g2) { var j = _g3++; var idxX0Yp = this._getIndex(i,j - 1); var idxX0Yn = this._getIndex(i,j + 1); var idxXpY0 = this._getIndex(i - 1,j); var idxXnY0 = this._getIndex(i + 1,j); var idxX0Y0 = this._getIndex(i,j); var accel = pos[idxX0Yp].z + pos[idxX0Yn].z + pos[idxXpY0].z + pos[idxXnY0].z - pos[idxX0Y0].z * 4; accel *= this._k; this._velocity[i][j] = (this._velocity[i][j] + accel) * this._attenuation; } } var _g11 = 1; var _g4 = 99; while(_g11 < _g4) { var i1 = _g11++; var _g31 = 1; var _g21 = 99; while(_g31 < _g21) { var j1 = _g31++; pos[this._getIndex(i1,j1)].z += this._velocity[i1][j1]; } } this._count += 0.1; if(Math.random() < 0.11) { var xx = Math.floor(2 + 96 * Math.random()); var yy = Math.floor(2 + 96 * Math.random()); var amp = 30 * Math.random(); pos[this._getIndex(xx,yy)].z = amp; pos[this._getIndex(xx - 1,yy)].z = amp * 0.3; pos[this._getIndex(xx + 1,yy)].z = amp * 0.3; pos[this._getIndex(xx,yy - 1)].z = amp * 0.3; pos[this._getIndex(xx,yy + 1)].z = amp * 0.3; } this._geo.computeVertexNormals(); this._geo.verticesNeedUpdate = true; this._geo.normalsNeedUpdate = true; this._mat.update(); } ,_getIndex: function(i,j) { return j * 100 + i % 100; } ,getNoise: function(xx,yy) { var f = noise.perlin3; return f(xx,yy,0); } }); var mat_MyShaderMat = function(deform) { if(deform == null) { deform = 1; } this.ff = "\r\n//\r\n// Description : Array and textureless GLSL 2D/3D/4D simplex \r\n// noise functions.\r\n// Author : Ian McEwan, Ashima Arts.\r\n// Maintainer : ijm\r\n// Lastmod : 20110822 (ijm)\r\n// License : Copyright (C) 2011 Ashima Arts. All rights reserved.\r\n// Distributed under the MIT License. See LICENSE file.\r\n// https://github.com/ashima/webgl-noise\r\n// \r\n\r\nvec3 mod289(vec3 x) {\r\n\treturn x - floor(x * (1.0 / 289.0)) * 289.0;\r\n}\r\n\r\nvec4 mod289(vec4 x) {\r\n\treturn x - floor(x * (1.0 / 289.0)) * 289.0;\r\n}\r\n\r\nvec4 permute(vec4 x) {\r\n\treturn mod289(((x*34.0)+1.0)*x);\r\n}\r\n\r\nvec4 taylorInvSqrt(vec4 r){\r\n\treturn 1.79284291400159 - 0.85373472095314 * r;\r\n}\r\n\r\nfloat snoise(vec3 v) { \r\n\r\n\tconst vec2 C = vec2(1.0/6.0, 1.0/3.0) ;\r\n\tconst vec4 D = vec4(0.0, 0.5, 1.0, 2.0);\r\n\r\n\t// First corner\r\n\tvec3 i = floor(v + dot(v, C.yyy) );\r\n\tvec3 x0 = v - i + dot(i, C.xxx) ;\r\n\r\n\t// Other corners\r\n\tvec3 g = step(x0.yzx, x0.xyz);\r\n\tvec3 l = 1.0 - g;\r\n\tvec3 i1 = min( g.xyz, l.zxy );\r\n\tvec3 i2 = max( g.xyz, l.zxy );\r\n\r\n\t// x0 = x0 - 0.0 + 0.0 * C.xxx;\r\n\t// x1 = x0 - i1 + 1.0 * C.xxx;\r\n\t// x2 = x0 - i2 + 2.0 * C.xxx;\r\n\t// x3 = x0 - 1.0 + 3.0 * C.xxx;\r\n\tvec3 x1 = x0 - i1 + C.xxx;\r\n\tvec3 x2 = x0 - i2 + C.yyy; // 2.0*C.x = 1/3 = C.y\r\n\tvec3 x3 = x0 - D.yyy; // -1.0+3.0*C.x = -0.5 = -D.y\r\n\r\n\t// Permutations\r\n\ti = mod289(i); \r\n\tvec4 p = permute( permute( permute( \r\n\t\t i.z + vec4(0.0, i1.z, i2.z, 1.0 ))\r\n\t\t+ i.y + vec4(0.0, i1.y, i2.y, 1.0 )) \r\n\t\t+ i.x + vec4(0.0, i1.x, i2.x, 1.0 ));\r\n\r\n\t// Gradients: 7x7 points over a square, mapped onto an octahedron.\r\n\t// The ring size 17*17 = 289 is close to a multiple of 49 (49*6 = 294)\r\n\tfloat n_ = 0.142857142857; // 1.0/7.0\r\n\tvec3 ns = n_ * D.wyz - D.xzx;\r\n\r\n\tvec4 j = p - 49.0 * floor(p * ns.z * ns.z); // mod(p,7*7)\r\n\r\n\tvec4 x_ = floor(j * ns.z);\r\n\tvec4 y_ = floor(j - 7.0 * x_ ); // mod(j,N)\r\n\r\n\tvec4 x = x_ *ns.x + ns.yyyy;\r\n\tvec4 y = y_ *ns.x + ns.yyyy;\r\n\tvec4 h = 1.0 - abs(x) - abs(y);\r\n\r\n\tvec4 b0 = vec4( x.xy, y.xy );\r\n\tvec4 b1 = vec4( x.zw, y.zw );\r\n\r\n\t//vec4 s0 = vec4(lessThan(b0,0.0))*2.0 - 1.0;\r\n\t//vec4 s1 = vec4(lessThan(b1,0.0))*2.0 - 1.0;\r\n\tvec4 s0 = floor(b0)*2.0 + 1.0;\r\n\tvec4 s1 = floor(b1)*2.0 + 1.0;\r\n\tvec4 sh = -step(h, vec4(0.0));\r\n\r\n\tvec4 a0 = b0.xzyw + s0.xzyw*sh.xxyy ;\r\n\tvec4 a1 = b1.xzyw + s1.xzyw*sh.zzww ;\r\n\r\n\tvec3 p0 = vec3(a0.xy,h.x);\r\n\tvec3 p1 = vec3(a0.zw,h.y);\r\n\tvec3 p2 = vec3(a1.xy,h.z);\r\n\tvec3 p3 = vec3(a1.zw,h.w);\r\n\r\n\t//Normalise gradients\r\n\tvec4 norm = taylorInvSqrt(vec4(dot(p0,p0), dot(p1,p1), dot(p2, p2), dot(p3,p3)));\r\n\tp0 *= norm.x;\r\n\tp1 *= norm.y;\r\n\tp2 *= norm.z;\r\n\tp3 *= norm.w;\r\n\r\n\t// Mix final noise value\r\n\tvec4 m = max(0.6 - vec4(dot(x0,x0), dot(x1,x1), dot(x2,x2), dot(x3,x3)), 0.0);\r\n\tm = m * m;\r\n\treturn 42.0 * dot( m*m, vec4( dot(p0,x0), dot(p1,x1), dot(p2,x2), dot(p3,x3) ) );\r\n\r\n}\r\n\r\nvec3 snoiseVec3( vec3 x ){\r\n\r\n\tfloat s = snoise(vec3( x ));\r\n\tfloat s1 = snoise(vec3( x.y - 19.1 , x.z + 33.4 , x.x + 47.2 ));\r\n\tfloat s2 = snoise(vec3( x.z + 74.2 , x.x - 124.5 , x.y + 99.4 ));\r\n\tvec3 c = vec3( s , s1 , s2 );\r\n\treturn c;\r\n\r\n}\r\n\r\nvec3 snoiseVec3Abs( vec3 x ){\r\n\r\n\tfloat s = snoise(vec3( x ));\r\n\tfloat s1 = snoise(vec3( x.y - 19.1 , x.z + 33.4 , x.x + 47.2 ));\r\n\tfloat s2 = snoise(vec3( x.z + 74.2 , x.x - 124.5 , x.y + 99.4 ));\r\n\tvec3 c = vec3( abs(s) , abs(s1) , abs(s2) );\r\n\treturn c;\r\n\r\n}\r\n\r\n\r\nvec3 curlNoise( vec3 p ){\r\n \r\n\tconst float e = .1;\r\n\tvec3 dx = vec3( e , 0.0 , 0.0 );\r\n\tvec3 dy = vec3( 0.0 , e , 0.0 );\r\n\tvec3 dz = vec3( 0.0 , 0.0 , e );\r\n\r\n\tvec3 p_x0 = snoiseVec3( p - dx );\r\n\tvec3 p_x1 = snoiseVec3( p + dx );\r\n\tvec3 p_y0 = snoiseVec3( p - dy );\r\n\tvec3 p_y1 = snoiseVec3( p + dy );\r\n\tvec3 p_z0 = snoiseVec3( p - dz );\r\n\tvec3 p_z1 = snoiseVec3( p + dz );\r\n\r\n\tfloat x = p_y1.z - p_y0.z - p_z1.y + p_z0.y;\r\n\tfloat y = p_z1.x - p_z0.x - p_x1.z + p_x0.z;\r\n\tfloat z = p_x1.y - p_x0.y - p_y1.x + p_y0.x;\r\n\r\n\tconst float divisor = 1.0 / ( 2.0 * e );\r\n\treturn normalize( vec3( x , y , z ) * divisor );\r\n\r\n}\r\n\r\nvec3 curlNoise2( vec3 p ) {\r\n\r\n\tconst float e = .1;\r\n\r\n\tvec3 xNoisePotentialDerivatives = snoiseVec3( p );\r\n\tvec3 yNoisePotentialDerivatives = snoiseVec3( p + e * vec3( 3., -3., 1. ) );\r\n\tvec3 zNoisePotentialDerivatives = snoiseVec3( p + e * vec3( 2., 4., -3. ) );\r\n\r\n\tvec3 noiseVelocity = vec3(\r\n\t\tzNoisePotentialDerivatives.y - yNoisePotentialDerivatives.z,\r\n\t\txNoisePotentialDerivatives.z - zNoisePotentialDerivatives.x,\r\n\t\tyNoisePotentialDerivatives.x - xNoisePotentialDerivatives.y\r\n\t);\r\n\r\n\treturn normalize( noiseVelocity );\r\n\r\n}\r\n\r\nvec4 snoiseD(vec3 v) { //returns vec4(value, dx, dy, dz)\r\n const vec2 C = vec2(1.0/6.0, 1.0/3.0) ;\r\n const vec4 D = vec4(0.0, 0.5, 1.0, 2.0);\r\n \r\n vec3 i = floor(v + dot(v, C.yyy) );\r\n vec3 x0 = v - i + dot(i, C.xxx) ;\r\n \r\n vec3 g = step(x0.yzx, x0.xyz);\r\n vec3 l = 1.0 - g;\r\n vec3 i1 = min( g.xyz, l.zxy );\r\n vec3 i2 = max( g.xyz, l.zxy );\r\n \r\n vec3 x1 = x0 - i1 + C.xxx;\r\n vec3 x2 = x0 - i2 + C.yyy;\r\n vec3 x3 = x0 - D.yyy;\r\n \r\n i = mod289(i);\r\n vec4 p = permute( permute( permute(\r\n i.z + vec4(0.0, i1.z, i2.z, 1.0 ))\r\n + i.y + vec4(0.0, i1.y, i2.y, 1.0 ))\r\n + i.x + vec4(0.0, i1.x, i2.x, 1.0 ));\r\n \r\n float n_ = 0.142857142857; // 1.0/7.0\r\n vec3 ns = n_ * D.wyz - D.xzx;\r\n \r\n vec4 j = p - 49.0 * floor(p * ns.z * ns.z);\r\n \r\n vec4 x_ = floor(j * ns.z);\r\n vec4 y_ = floor(j - 7.0 * x_ );\r\n \r\n vec4 x = x_ *ns.x + ns.yyyy;\r\n vec4 y = y_ *ns.x + ns.yyyy;\r\n vec4 h = 1.0 - abs(x) - abs(y);\r\n \r\n vec4 b0 = vec4( x.xy, y.xy );\r\n vec4 b1 = vec4( x.zw, y.zw );\r\n \r\n vec4 s0 = floor(b0)*2.0 + 1.0;\r\n vec4 s1 = floor(b1)*2.0 + 1.0;\r\n vec4 sh = -step(h, vec4(0.0));\r\n \r\n vec4 a0 = b0.xzyw + s0.xzyw*sh.xxyy ;\r\n vec4 a1 = b1.xzyw + s1.xzyw*sh.zzww ;\r\n \r\n vec3 p0 = vec3(a0.xy,h.x);\r\n vec3 p1 = vec3(a0.zw,h.y);\r\n vec3 p2 = vec3(a1.xy,h.z);\r\n vec3 p3 = vec3(a1.zw,h.w);\r\n \r\n vec4 norm = taylorInvSqrt(vec4(dot(p0,p0), dot(p1,p1), dot(p2, p2), dot(p3,p3)));\r\n p0 *= norm.x;\r\n p1 *= norm.y;\r\n p2 *= norm.z;\r\n p3 *= norm.w;\r\n \r\n vec4 values = vec4( dot(p0,x0), dot(p1,x1), dot(p2,x2), dot(p3,x3) ); //value of contributions from each corner (extrapolate the gradient)\r\n \r\n vec4 m = max(0.5 - vec4(dot(x0,x0), dot(x1,x1), dot(x2,x2), dot(x3,x3)), 0.0); //kernel function from each corner\r\n \r\n vec4 m2 = m * m;\r\n vec4 m3 = m * m * m;\r\n \r\n vec4 temp = -6.0 * m2 * values;\r\n float dx = temp[0] * x0.x + temp[1] * x1.x + temp[2] * x2.x + temp[3] * x3.x + m3[0] * p0.x + m3[1] * p1.x + m3[2] * p2.x + m3[3] * p3.x;\r\n float dy = temp[0] * x0.y + temp[1] * x1.y + temp[2] * x2.y + temp[3] * x3.y + m3[0] * p0.y + m3[1] * p1.y + m3[2] * p2.y + m3[3] * p3.y;\r\n float dz = temp[0] * x0.z + temp[1] * x1.z + temp[2] * x2.z + temp[3] * x3.z + m3[0] * p0.z + m3[1] * p1.z + m3[2] * p2.z + m3[3] * p3.z;\r\n \r\n return vec4(dot(m3, values), dx, dy, dz) * 42.0;\r\n}\r\n\r\n\r\nvec3 curlNoise3 (vec3 p) {\r\n\r\n vec3 xNoisePotentialDerivatives = snoiseD( p ).yzw; //yzw are the xyz derivatives\r\n vec3 yNoisePotentialDerivatives = snoiseD(vec3( p.y - 19.1 , p.z + 33.4 , p.x + 47.2 )).zwy;\r\n vec3 zNoisePotentialDerivatives = snoiseD(vec3( p.z + 74.2 , p.x - 124.5 , p.y + 99.4 )).wyz;\r\n vec3 noiseVelocity = vec3(\r\n zNoisePotentialDerivatives.y - yNoisePotentialDerivatives.z,\r\n xNoisePotentialDerivatives.z - zNoisePotentialDerivatives.x,\r\n yNoisePotentialDerivatives.x - xNoisePotentialDerivatives.y\r\n );\r\n\t\r\n\tconst float e = .1;\r\n\tconst float divisor = 1.0 / ( 2.0 * e );\r\n\treturn normalize( noiseVelocity * divisor );\r\n\r\n}\r\n\t\r\n\t" + "\r\n\t\r\n\t\t#define PHONG\r\n\r\n\t\tuniform vec3 diffuse;\r\n\t\tuniform vec3 emissive;\r\n\t\tuniform vec3 specular;\r\n\t\tuniform float shininess;\r\n\t\tuniform float opacity;\r\n\r\n\t\tuniform float counter;\r\n\t\tvarying vec3 vPos;\r\n\t\t\r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\r\n\t\tfloat rand1(vec2 co){\r\n\t\t\treturn fract(sin(dot(co.xy ,vec2(12.9898,78.233))) * 43758.5453);\r\n\t\t}\t\t\r\n\t\t\r\n\t\tvoid main() {\r\n\r\n\t\t\t#include \r\n\r\n\t\t\tvec4 diffuseColor = vec4( diffuse, opacity );\r\n\t\t\tReflectedLight reflectedLight = ReflectedLight( vec3( 0.0 ), vec3( 0.0 ), vec3( 0.0 ), vec3( 0.0 ) );\r\n\t\t\tvec3 totalEmissiveRadiance = emissive;\r\n\t\t\t\r\n\t\t\tdiffuseColor.xyz = snoiseVec3( vec3(vPos.x * 0.01, vPos.y * 0.01, vPos.z* 30.2 + counter) );\r\n\t\t\tdiffuseColor.x = abs( 1.0 * diffuseColor.z + 0.7 ) + sin( vPos.z * 0.2 );\r\n\t\t\tdiffuseColor.y = abs( 0.8 * diffuseColor.y + 0.5 ) + sin( vPos.z * 0.1 );\r\n\t\t\tdiffuseColor.z = abs( 0.2 * diffuseColor.z + 0.7 ) + sin( vPos.z * 0.05 );\r\n\t\t\tdiffuseColor.x = mix(1.00,diffuseColor.x, smoothstep(0.0, 0.1, abs(vPos.z)) );\r\n\t\t\tdiffuseColor.y = mix(0.97,diffuseColor.y, smoothstep(0.0, 0.1, abs(vPos.z)) );\r\n\t\t\tdiffuseColor.z = mix(0.85,diffuseColor.z, smoothstep(0.0, 0.1, abs(vPos.z)) );\r\n\t\t\t\r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\r\n\t\t\t// accumulation\r\n\t\t\t#include \r\n\t\t\t#include \r\n\r\n\t\t\t// modulation\r\n\t\t\t#include \r\n\r\n\t\t\tvec3 outgoingLight = reflectedLight.directDiffuse + reflectedLight.indirectDiffuse + reflectedLight.directSpecular + reflectedLight.indirectSpecular + totalEmissiveRadiance;\r\n\r\n\t\t\t#include \r\n\r\n\t\t\tvec3 vv = vec3(0.15 * rand1(diffuseColor.xy * 100.0 + counter));\r\n\t\t\toutgoingLight.xyz += vv;\t\t\t\r\n\t\t\t\r\n\t\t\tgl_FragColor = vec4( outgoingLight, diffuseColor.a );\r\n\r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\r\n\t\t}\r\n\t\r\n\t"; this.vv = "\r\n//\r\n// Description : Array and textureless GLSL 2D/3D/4D simplex \r\n// noise functions.\r\n// Author : Ian McEwan, Ashima Arts.\r\n// Maintainer : ijm\r\n// Lastmod : 20110822 (ijm)\r\n// License : Copyright (C) 2011 Ashima Arts. All rights reserved.\r\n// Distributed under the MIT License. See LICENSE file.\r\n// https://github.com/ashima/webgl-noise\r\n// \r\n\r\nvec3 mod289(vec3 x) {\r\n\treturn x - floor(x * (1.0 / 289.0)) * 289.0;\r\n}\r\n\r\nvec4 mod289(vec4 x) {\r\n\treturn x - floor(x * (1.0 / 289.0)) * 289.0;\r\n}\r\n\r\nvec4 permute(vec4 x) {\r\n\treturn mod289(((x*34.0)+1.0)*x);\r\n}\r\n\r\nvec4 taylorInvSqrt(vec4 r){\r\n\treturn 1.79284291400159 - 0.85373472095314 * r;\r\n}\r\n\r\nfloat snoise(vec3 v) { \r\n\r\n\tconst vec2 C = vec2(1.0/6.0, 1.0/3.0) ;\r\n\tconst vec4 D = vec4(0.0, 0.5, 1.0, 2.0);\r\n\r\n\t// First corner\r\n\tvec3 i = floor(v + dot(v, C.yyy) );\r\n\tvec3 x0 = v - i + dot(i, C.xxx) ;\r\n\r\n\t// Other corners\r\n\tvec3 g = step(x0.yzx, x0.xyz);\r\n\tvec3 l = 1.0 - g;\r\n\tvec3 i1 = min( g.xyz, l.zxy );\r\n\tvec3 i2 = max( g.xyz, l.zxy );\r\n\r\n\t// x0 = x0 - 0.0 + 0.0 * C.xxx;\r\n\t// x1 = x0 - i1 + 1.0 * C.xxx;\r\n\t// x2 = x0 - i2 + 2.0 * C.xxx;\r\n\t// x3 = x0 - 1.0 + 3.0 * C.xxx;\r\n\tvec3 x1 = x0 - i1 + C.xxx;\r\n\tvec3 x2 = x0 - i2 + C.yyy; // 2.0*C.x = 1/3 = C.y\r\n\tvec3 x3 = x0 - D.yyy; // -1.0+3.0*C.x = -0.5 = -D.y\r\n\r\n\t// Permutations\r\n\ti = mod289(i); \r\n\tvec4 p = permute( permute( permute( \r\n\t\t i.z + vec4(0.0, i1.z, i2.z, 1.0 ))\r\n\t\t+ i.y + vec4(0.0, i1.y, i2.y, 1.0 )) \r\n\t\t+ i.x + vec4(0.0, i1.x, i2.x, 1.0 ));\r\n\r\n\t// Gradients: 7x7 points over a square, mapped onto an octahedron.\r\n\t// The ring size 17*17 = 289 is close to a multiple of 49 (49*6 = 294)\r\n\tfloat n_ = 0.142857142857; // 1.0/7.0\r\n\tvec3 ns = n_ * D.wyz - D.xzx;\r\n\r\n\tvec4 j = p - 49.0 * floor(p * ns.z * ns.z); // mod(p,7*7)\r\n\r\n\tvec4 x_ = floor(j * ns.z);\r\n\tvec4 y_ = floor(j - 7.0 * x_ ); // mod(j,N)\r\n\r\n\tvec4 x = x_ *ns.x + ns.yyyy;\r\n\tvec4 y = y_ *ns.x + ns.yyyy;\r\n\tvec4 h = 1.0 - abs(x) - abs(y);\r\n\r\n\tvec4 b0 = vec4( x.xy, y.xy );\r\n\tvec4 b1 = vec4( x.zw, y.zw );\r\n\r\n\t//vec4 s0 = vec4(lessThan(b0,0.0))*2.0 - 1.0;\r\n\t//vec4 s1 = vec4(lessThan(b1,0.0))*2.0 - 1.0;\r\n\tvec4 s0 = floor(b0)*2.0 + 1.0;\r\n\tvec4 s1 = floor(b1)*2.0 + 1.0;\r\n\tvec4 sh = -step(h, vec4(0.0));\r\n\r\n\tvec4 a0 = b0.xzyw + s0.xzyw*sh.xxyy ;\r\n\tvec4 a1 = b1.xzyw + s1.xzyw*sh.zzww ;\r\n\r\n\tvec3 p0 = vec3(a0.xy,h.x);\r\n\tvec3 p1 = vec3(a0.zw,h.y);\r\n\tvec3 p2 = vec3(a1.xy,h.z);\r\n\tvec3 p3 = vec3(a1.zw,h.w);\r\n\r\n\t//Normalise gradients\r\n\tvec4 norm = taylorInvSqrt(vec4(dot(p0,p0), dot(p1,p1), dot(p2, p2), dot(p3,p3)));\r\n\tp0 *= norm.x;\r\n\tp1 *= norm.y;\r\n\tp2 *= norm.z;\r\n\tp3 *= norm.w;\r\n\r\n\t// Mix final noise value\r\n\tvec4 m = max(0.6 - vec4(dot(x0,x0), dot(x1,x1), dot(x2,x2), dot(x3,x3)), 0.0);\r\n\tm = m * m;\r\n\treturn 42.0 * dot( m*m, vec4( dot(p0,x0), dot(p1,x1), dot(p2,x2), dot(p3,x3) ) );\r\n\r\n}\r\n\r\nvec3 snoiseVec3( vec3 x ){\r\n\r\n\tfloat s = snoise(vec3( x ));\r\n\tfloat s1 = snoise(vec3( x.y - 19.1 , x.z + 33.4 , x.x + 47.2 ));\r\n\tfloat s2 = snoise(vec3( x.z + 74.2 , x.x - 124.5 , x.y + 99.4 ));\r\n\tvec3 c = vec3( s , s1 , s2 );\r\n\treturn c;\r\n\r\n}\r\n\r\nvec3 snoiseVec3Abs( vec3 x ){\r\n\r\n\tfloat s = snoise(vec3( x ));\r\n\tfloat s1 = snoise(vec3( x.y - 19.1 , x.z + 33.4 , x.x + 47.2 ));\r\n\tfloat s2 = snoise(vec3( x.z + 74.2 , x.x - 124.5 , x.y + 99.4 ));\r\n\tvec3 c = vec3( abs(s) , abs(s1) , abs(s2) );\r\n\treturn c;\r\n\r\n}\r\n\r\n\r\nvec3 curlNoise( vec3 p ){\r\n \r\n\tconst float e = .1;\r\n\tvec3 dx = vec3( e , 0.0 , 0.0 );\r\n\tvec3 dy = vec3( 0.0 , e , 0.0 );\r\n\tvec3 dz = vec3( 0.0 , 0.0 , e );\r\n\r\n\tvec3 p_x0 = snoiseVec3( p - dx );\r\n\tvec3 p_x1 = snoiseVec3( p + dx );\r\n\tvec3 p_y0 = snoiseVec3( p - dy );\r\n\tvec3 p_y1 = snoiseVec3( p + dy );\r\n\tvec3 p_z0 = snoiseVec3( p - dz );\r\n\tvec3 p_z1 = snoiseVec3( p + dz );\r\n\r\n\tfloat x = p_y1.z - p_y0.z - p_z1.y + p_z0.y;\r\n\tfloat y = p_z1.x - p_z0.x - p_x1.z + p_x0.z;\r\n\tfloat z = p_x1.y - p_x0.y - p_y1.x + p_y0.x;\r\n\r\n\tconst float divisor = 1.0 / ( 2.0 * e );\r\n\treturn normalize( vec3( x , y , z ) * divisor );\r\n\r\n}\r\n\r\nvec3 curlNoise2( vec3 p ) {\r\n\r\n\tconst float e = .1;\r\n\r\n\tvec3 xNoisePotentialDerivatives = snoiseVec3( p );\r\n\tvec3 yNoisePotentialDerivatives = snoiseVec3( p + e * vec3( 3., -3., 1. ) );\r\n\tvec3 zNoisePotentialDerivatives = snoiseVec3( p + e * vec3( 2., 4., -3. ) );\r\n\r\n\tvec3 noiseVelocity = vec3(\r\n\t\tzNoisePotentialDerivatives.y - yNoisePotentialDerivatives.z,\r\n\t\txNoisePotentialDerivatives.z - zNoisePotentialDerivatives.x,\r\n\t\tyNoisePotentialDerivatives.x - xNoisePotentialDerivatives.y\r\n\t);\r\n\r\n\treturn normalize( noiseVelocity );\r\n\r\n}\r\n\r\nvec4 snoiseD(vec3 v) { //returns vec4(value, dx, dy, dz)\r\n const vec2 C = vec2(1.0/6.0, 1.0/3.0) ;\r\n const vec4 D = vec4(0.0, 0.5, 1.0, 2.0);\r\n \r\n vec3 i = floor(v + dot(v, C.yyy) );\r\n vec3 x0 = v - i + dot(i, C.xxx) ;\r\n \r\n vec3 g = step(x0.yzx, x0.xyz);\r\n vec3 l = 1.0 - g;\r\n vec3 i1 = min( g.xyz, l.zxy );\r\n vec3 i2 = max( g.xyz, l.zxy );\r\n \r\n vec3 x1 = x0 - i1 + C.xxx;\r\n vec3 x2 = x0 - i2 + C.yyy;\r\n vec3 x3 = x0 - D.yyy;\r\n \r\n i = mod289(i);\r\n vec4 p = permute( permute( permute(\r\n i.z + vec4(0.0, i1.z, i2.z, 1.0 ))\r\n + i.y + vec4(0.0, i1.y, i2.y, 1.0 ))\r\n + i.x + vec4(0.0, i1.x, i2.x, 1.0 ));\r\n \r\n float n_ = 0.142857142857; // 1.0/7.0\r\n vec3 ns = n_ * D.wyz - D.xzx;\r\n \r\n vec4 j = p - 49.0 * floor(p * ns.z * ns.z);\r\n \r\n vec4 x_ = floor(j * ns.z);\r\n vec4 y_ = floor(j - 7.0 * x_ );\r\n \r\n vec4 x = x_ *ns.x + ns.yyyy;\r\n vec4 y = y_ *ns.x + ns.yyyy;\r\n vec4 h = 1.0 - abs(x) - abs(y);\r\n \r\n vec4 b0 = vec4( x.xy, y.xy );\r\n vec4 b1 = vec4( x.zw, y.zw );\r\n \r\n vec4 s0 = floor(b0)*2.0 + 1.0;\r\n vec4 s1 = floor(b1)*2.0 + 1.0;\r\n vec4 sh = -step(h, vec4(0.0));\r\n \r\n vec4 a0 = b0.xzyw + s0.xzyw*sh.xxyy ;\r\n vec4 a1 = b1.xzyw + s1.xzyw*sh.zzww ;\r\n \r\n vec3 p0 = vec3(a0.xy,h.x);\r\n vec3 p1 = vec3(a0.zw,h.y);\r\n vec3 p2 = vec3(a1.xy,h.z);\r\n vec3 p3 = vec3(a1.zw,h.w);\r\n \r\n vec4 norm = taylorInvSqrt(vec4(dot(p0,p0), dot(p1,p1), dot(p2, p2), dot(p3,p3)));\r\n p0 *= norm.x;\r\n p1 *= norm.y;\r\n p2 *= norm.z;\r\n p3 *= norm.w;\r\n \r\n vec4 values = vec4( dot(p0,x0), dot(p1,x1), dot(p2,x2), dot(p3,x3) ); //value of contributions from each corner (extrapolate the gradient)\r\n \r\n vec4 m = max(0.5 - vec4(dot(x0,x0), dot(x1,x1), dot(x2,x2), dot(x3,x3)), 0.0); //kernel function from each corner\r\n \r\n vec4 m2 = m * m;\r\n vec4 m3 = m * m * m;\r\n \r\n vec4 temp = -6.0 * m2 * values;\r\n float dx = temp[0] * x0.x + temp[1] * x1.x + temp[2] * x2.x + temp[3] * x3.x + m3[0] * p0.x + m3[1] * p1.x + m3[2] * p2.x + m3[3] * p3.x;\r\n float dy = temp[0] * x0.y + temp[1] * x1.y + temp[2] * x2.y + temp[3] * x3.y + m3[0] * p0.y + m3[1] * p1.y + m3[2] * p2.y + m3[3] * p3.y;\r\n float dz = temp[0] * x0.z + temp[1] * x1.z + temp[2] * x2.z + temp[3] * x3.z + m3[0] * p0.z + m3[1] * p1.z + m3[2] * p2.z + m3[3] * p3.z;\r\n \r\n return vec4(dot(m3, values), dx, dy, dz) * 42.0;\r\n}\r\n\r\n\r\nvec3 curlNoise3 (vec3 p) {\r\n\r\n vec3 xNoisePotentialDerivatives = snoiseD( p ).yzw; //yzw are the xyz derivatives\r\n vec3 yNoisePotentialDerivatives = snoiseD(vec3( p.y - 19.1 , p.z + 33.4 , p.x + 47.2 )).zwy;\r\n vec3 zNoisePotentialDerivatives = snoiseD(vec3( p.z + 74.2 , p.x - 124.5 , p.y + 99.4 )).wyz;\r\n vec3 noiseVelocity = vec3(\r\n zNoisePotentialDerivatives.y - yNoisePotentialDerivatives.z,\r\n xNoisePotentialDerivatives.z - zNoisePotentialDerivatives.x,\r\n yNoisePotentialDerivatives.x - xNoisePotentialDerivatives.y\r\n );\r\n\t\r\n\tconst float e = .1;\r\n\tconst float divisor = 1.0 / ( 2.0 * e );\r\n\treturn normalize( noiseVelocity * divisor );\r\n\r\n}\r\n\t\r\n\t" + "\r\n\t\t#define PHONG\r\n\r\n\t\tvarying vec3 vViewPosition;\r\n\t\t\r\n\t\tvarying vec3 vPos;\r\n\t\t\r\n\t\t#ifndef FLAT_SHADED\r\n\r\n\t\t\tvarying vec3 vNormal;\r\n\r\n\t\t#endif\r\n\r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\t\t#include \r\n\r\n\t\tvoid main() {\r\n\r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\r\n\t\t#ifndef FLAT_SHADED // Normal computed with derivatives when FLAT_SHADED\r\n\r\n\t\t\tvNormal = normalize( transformedNormal );\r\n\r\n\t\t#endif\r\n\r\n\t\t\t#include \r\n\t\t\tvPos = position;\r\n\t\t\t\r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\r\n\t\t\tvViewPosition = - mvPosition.xyz;\r\n\r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\t\t\t#include \r\n\r\n\t\t}\t\r\n\t"; var phongShader = THREE.ShaderLib.phong; var unif = THREE.UniformsUtils .clone(phongShader.uniforms); unif.counter = { type : "f", value : 0}; unif.specular.value = new THREE.Color(8947848); unif.emissive.value = new THREE.Color(2236962); THREE.ShaderMaterial.call(this,{ uniforms : unif, vertexShader : this.vv, fragmentShader : this.ff, lights : true}); this.vertexColors = false; this.side = 2; }; mat_MyShaderMat.__super__ = THREE.ShaderMaterial; mat_MyShaderMat.prototype = $extend(THREE.ShaderMaterial.prototype,{ update: function() { this.uniforms.counter.value += 0.005; } }); var shader_SimplexNoise = function() { }; var three_Face = function() { }; var three_IFog = function() { }; var three_Mapping = function() { }; var three_Renderer = function() { }; var utils_VectorUtils = function() { }; utils_VectorUtils.limit = function(tgt,lim) { if(tgt.length() > lim) { tgt = tgt.normalize().multiplyScalar(lim); } return tgt; }; var $_, $fid = 0; function $bind(o,m) { if( m == null ) return null; if( m.__id__ == null ) m.__id__ = $fid++; var f; if( o.hx__closures__ == null ) o.hx__closures__ = {}; else f = o.hx__closures__[m.__id__]; if( f == null ) { f = function(){ return f.method.apply(f.scope, arguments); }; f.scope = o; f.method = m; o.hx__closures__[m.__id__] = f; } return f; } var q = window.jQuery; var js = js || {} js.JQuery = q; Three.CullFaceNone = 0; Three.CullFaceBack = 1; Three.CullFaceFront = 2; Three.CullFaceFrontBack = 3; Three.FrontFaceDirectionCW = 0; Three.FrontFaceDirectionCCW = 1; Three.BasicShadowMap = 0; Three.PCFShadowMap = 1; Three.PCFSoftShadowMap = 2; Three.FrontSide = 0; Three.BackSide = 1; Three.DoubleSide = 2; Three.NoShading = 0; Three.FlatShading = 1; Three.SmoothShading = 2; Three.NoColors = 0; Three.FaceColors = 1; Three.VertexColors = 2; Three.NoBlending = 0; Three.NormalBlending = 1; Three.AdditiveBlending = 2; Three.SubtractiveBlending = 3; Three.MultiplyBlending = 4; Three.CustomBlending = 5; Three.AddEquation = 100; Three.SubtractEquation = 101; Three.ReverseSubtractEquation = 102; Three.ZeroFactor = 200; Three.OneFactor = 201; Three.SrcColorFactor = 202; Three.OneMinusSrcColorFactor = 203; Three.SrcAlphaFactor = 204; Three.OneMinusSrcAlphaFactor = 205; Three.DstAlphaFactor = 206; Three.OneMinusDstAlphaFactor = 207; Three.MultiplyOperation = 0; Three.MixOperation = 1; Three.AddOperation = 2; Three.RepeatWrapping = 1000; Three.ClampToEdgeWrapping = 1001; Three.MirroredRepeatWrapping = 1002; Three.NearestFilter = 1003; Three.NearestMipMapNearestFilter = 1004; Three.NearestMipMapLinearFilter = 1005; Three.LinearFilter = 1006; Three.LinearMipMapNearestFilter = 1007; Three.LinearMipMapLinearFilter = 1008; Three.UnsignedByteType = 1009; Three.ByteType = 1010; Three.ShortType = 1011; Three.UnsignedShortType = 1012; Three.IntType = 1013; Three.UnsignedIntType = 1014; Three.FloatType = 1015; Three.UnsignedShort4444Type = 1016; Three.UnsignedShort5551Type = 1017; Three.UnsignedShort565Type = 1018; Three.AlphaFormat = 1019; Three.RGBFormat = 1020; Three.RGBAFormat = 1021; Three.LuminanceFormat = 1022; Three.LuminanceAlphaFormat = 1023; Three.RGB_S3TC_DXT1_Format = 2001; Three.RGBA_S3TC_DXT1_Format = 2002; Three.RGBA_S3TC_DXT3_Format = 2003; Three.RGBA_S3TC_DXT5_Format = 2004; Three.LineStrip = 0; Three.LinePieces = 1; Wave.MAX_X = 100; Wave.MAX_Y = 100; shader_SimplexNoise.glsl = "\r\n//\r\n// Description : Array and textureless GLSL 2D/3D/4D simplex \r\n// noise functions.\r\n// Author : Ian McEwan, Ashima Arts.\r\n// Maintainer : ijm\r\n// Lastmod : 20110822 (ijm)\r\n// License : Copyright (C) 2011 Ashima Arts. All rights reserved.\r\n// Distributed under the MIT License. See LICENSE file.\r\n// https://github.com/ashima/webgl-noise\r\n// \r\n\r\nvec3 mod289(vec3 x) {\r\n\treturn x - floor(x * (1.0 / 289.0)) * 289.0;\r\n}\r\n\r\nvec4 mod289(vec4 x) {\r\n\treturn x - floor(x * (1.0 / 289.0)) * 289.0;\r\n}\r\n\r\nvec4 permute(vec4 x) {\r\n\treturn mod289(((x*34.0)+1.0)*x);\r\n}\r\n\r\nvec4 taylorInvSqrt(vec4 r){\r\n\treturn 1.79284291400159 - 0.85373472095314 * r;\r\n}\r\n\r\nfloat snoise(vec3 v) { \r\n\r\n\tconst vec2 C = vec2(1.0/6.0, 1.0/3.0) ;\r\n\tconst vec4 D = vec4(0.0, 0.5, 1.0, 2.0);\r\n\r\n\t// First corner\r\n\tvec3 i = floor(v + dot(v, C.yyy) );\r\n\tvec3 x0 = v - i + dot(i, C.xxx) ;\r\n\r\n\t// Other corners\r\n\tvec3 g = step(x0.yzx, x0.xyz);\r\n\tvec3 l = 1.0 - g;\r\n\tvec3 i1 = min( g.xyz, l.zxy );\r\n\tvec3 i2 = max( g.xyz, l.zxy );\r\n\r\n\t// x0 = x0 - 0.0 + 0.0 * C.xxx;\r\n\t// x1 = x0 - i1 + 1.0 * C.xxx;\r\n\t// x2 = x0 - i2 + 2.0 * C.xxx;\r\n\t// x3 = x0 - 1.0 + 3.0 * C.xxx;\r\n\tvec3 x1 = x0 - i1 + C.xxx;\r\n\tvec3 x2 = x0 - i2 + C.yyy; // 2.0*C.x = 1/3 = C.y\r\n\tvec3 x3 = x0 - D.yyy; // -1.0+3.0*C.x = -0.5 = -D.y\r\n\r\n\t// Permutations\r\n\ti = mod289(i); \r\n\tvec4 p = permute( permute( permute( \r\n\t\t i.z + vec4(0.0, i1.z, i2.z, 1.0 ))\r\n\t\t+ i.y + vec4(0.0, i1.y, i2.y, 1.0 )) \r\n\t\t+ i.x + vec4(0.0, i1.x, i2.x, 1.0 ));\r\n\r\n\t// Gradients: 7x7 points over a square, mapped onto an octahedron.\r\n\t// The ring size 17*17 = 289 is close to a multiple of 49 (49*6 = 294)\r\n\tfloat n_ = 0.142857142857; // 1.0/7.0\r\n\tvec3 ns = n_ * D.wyz - D.xzx;\r\n\r\n\tvec4 j = p - 49.0 * floor(p * ns.z * ns.z); // mod(p,7*7)\r\n\r\n\tvec4 x_ = floor(j * ns.z);\r\n\tvec4 y_ = floor(j - 7.0 * x_ ); // mod(j,N)\r\n\r\n\tvec4 x = x_ *ns.x + ns.yyyy;\r\n\tvec4 y = y_ *ns.x + ns.yyyy;\r\n\tvec4 h = 1.0 - abs(x) - abs(y);\r\n\r\n\tvec4 b0 = vec4( x.xy, y.xy );\r\n\tvec4 b1 = vec4( x.zw, y.zw );\r\n\r\n\t//vec4 s0 = vec4(lessThan(b0,0.0))*2.0 - 1.0;\r\n\t//vec4 s1 = vec4(lessThan(b1,0.0))*2.0 - 1.0;\r\n\tvec4 s0 = floor(b0)*2.0 + 1.0;\r\n\tvec4 s1 = floor(b1)*2.0 + 1.0;\r\n\tvec4 sh = -step(h, vec4(0.0));\r\n\r\n\tvec4 a0 = b0.xzyw + s0.xzyw*sh.xxyy ;\r\n\tvec4 a1 = b1.xzyw + s1.xzyw*sh.zzww ;\r\n\r\n\tvec3 p0 = vec3(a0.xy,h.x);\r\n\tvec3 p1 = vec3(a0.zw,h.y);\r\n\tvec3 p2 = vec3(a1.xy,h.z);\r\n\tvec3 p3 = vec3(a1.zw,h.w);\r\n\r\n\t//Normalise gradients\r\n\tvec4 norm = taylorInvSqrt(vec4(dot(p0,p0), dot(p1,p1), dot(p2, p2), dot(p3,p3)));\r\n\tp0 *= norm.x;\r\n\tp1 *= norm.y;\r\n\tp2 *= norm.z;\r\n\tp3 *= norm.w;\r\n\r\n\t// Mix final noise value\r\n\tvec4 m = max(0.6 - vec4(dot(x0,x0), dot(x1,x1), dot(x2,x2), dot(x3,x3)), 0.0);\r\n\tm = m * m;\r\n\treturn 42.0 * dot( m*m, vec4( dot(p0,x0), dot(p1,x1), dot(p2,x2), dot(p3,x3) ) );\r\n\r\n}\r\n\r\nvec3 snoiseVec3( vec3 x ){\r\n\r\n\tfloat s = snoise(vec3( x ));\r\n\tfloat s1 = snoise(vec3( x.y - 19.1 , x.z + 33.4 , x.x + 47.2 ));\r\n\tfloat s2 = snoise(vec3( x.z + 74.2 , x.x - 124.5 , x.y + 99.4 ));\r\n\tvec3 c = vec3( s , s1 , s2 );\r\n\treturn c;\r\n\r\n}\r\n\r\nvec3 snoiseVec3Abs( vec3 x ){\r\n\r\n\tfloat s = snoise(vec3( x ));\r\n\tfloat s1 = snoise(vec3( x.y - 19.1 , x.z + 33.4 , x.x + 47.2 ));\r\n\tfloat s2 = snoise(vec3( x.z + 74.2 , x.x - 124.5 , x.y + 99.4 ));\r\n\tvec3 c = vec3( abs(s) , abs(s1) , abs(s2) );\r\n\treturn c;\r\n\r\n}\r\n\r\n\r\nvec3 curlNoise( vec3 p ){\r\n \r\n\tconst float e = .1;\r\n\tvec3 dx = vec3( e , 0.0 , 0.0 );\r\n\tvec3 dy = vec3( 0.0 , e , 0.0 );\r\n\tvec3 dz = vec3( 0.0 , 0.0 , e );\r\n\r\n\tvec3 p_x0 = snoiseVec3( p - dx );\r\n\tvec3 p_x1 = snoiseVec3( p + dx );\r\n\tvec3 p_y0 = snoiseVec3( p - dy );\r\n\tvec3 p_y1 = snoiseVec3( p + dy );\r\n\tvec3 p_z0 = snoiseVec3( p - dz );\r\n\tvec3 p_z1 = snoiseVec3( p + dz );\r\n\r\n\tfloat x = p_y1.z - p_y0.z - p_z1.y + p_z0.y;\r\n\tfloat y = p_z1.x - p_z0.x - p_x1.z + p_x0.z;\r\n\tfloat z = p_x1.y - p_x0.y - p_y1.x + p_y0.x;\r\n\r\n\tconst float divisor = 1.0 / ( 2.0 * e );\r\n\treturn normalize( vec3( x , y , z ) * divisor );\r\n\r\n}\r\n\r\nvec3 curlNoise2( vec3 p ) {\r\n\r\n\tconst float e = .1;\r\n\r\n\tvec3 xNoisePotentialDerivatives = snoiseVec3( p );\r\n\tvec3 yNoisePotentialDerivatives = snoiseVec3( p + e * vec3( 3., -3., 1. ) );\r\n\tvec3 zNoisePotentialDerivatives = snoiseVec3( p + e * vec3( 2., 4., -3. ) );\r\n\r\n\tvec3 noiseVelocity = vec3(\r\n\t\tzNoisePotentialDerivatives.y - yNoisePotentialDerivatives.z,\r\n\t\txNoisePotentialDerivatives.z - zNoisePotentialDerivatives.x,\r\n\t\tyNoisePotentialDerivatives.x - xNoisePotentialDerivatives.y\r\n\t);\r\n\r\n\treturn normalize( noiseVelocity );\r\n\r\n}\r\n\r\nvec4 snoiseD(vec3 v) { //returns vec4(value, dx, dy, dz)\r\n const vec2 C = vec2(1.0/6.0, 1.0/3.0) ;\r\n const vec4 D = vec4(0.0, 0.5, 1.0, 2.0);\r\n \r\n vec3 i = floor(v + dot(v, C.yyy) );\r\n vec3 x0 = v - i + dot(i, C.xxx) ;\r\n \r\n vec3 g = step(x0.yzx, x0.xyz);\r\n vec3 l = 1.0 - g;\r\n vec3 i1 = min( g.xyz, l.zxy );\r\n vec3 i2 = max( g.xyz, l.zxy );\r\n \r\n vec3 x1 = x0 - i1 + C.xxx;\r\n vec3 x2 = x0 - i2 + C.yyy;\r\n vec3 x3 = x0 - D.yyy;\r\n \r\n i = mod289(i);\r\n vec4 p = permute( permute( permute(\r\n i.z + vec4(0.0, i1.z, i2.z, 1.0 ))\r\n + i.y + vec4(0.0, i1.y, i2.y, 1.0 ))\r\n + i.x + vec4(0.0, i1.x, i2.x, 1.0 ));\r\n \r\n float n_ = 0.142857142857; // 1.0/7.0\r\n vec3 ns = n_ * D.wyz - D.xzx;\r\n \r\n vec4 j = p - 49.0 * floor(p * ns.z * ns.z);\r\n \r\n vec4 x_ = floor(j * ns.z);\r\n vec4 y_ = floor(j - 7.0 * x_ );\r\n \r\n vec4 x = x_ *ns.x + ns.yyyy;\r\n vec4 y = y_ *ns.x + ns.yyyy;\r\n vec4 h = 1.0 - abs(x) - abs(y);\r\n \r\n vec4 b0 = vec4( x.xy, y.xy );\r\n vec4 b1 = vec4( x.zw, y.zw );\r\n \r\n vec4 s0 = floor(b0)*2.0 + 1.0;\r\n vec4 s1 = floor(b1)*2.0 + 1.0;\r\n vec4 sh = -step(h, vec4(0.0));\r\n \r\n vec4 a0 = b0.xzyw + s0.xzyw*sh.xxyy ;\r\n vec4 a1 = b1.xzyw + s1.xzyw*sh.zzww ;\r\n \r\n vec3 p0 = vec3(a0.xy,h.x);\r\n vec3 p1 = vec3(a0.zw,h.y);\r\n vec3 p2 = vec3(a1.xy,h.z);\r\n vec3 p3 = vec3(a1.zw,h.w);\r\n \r\n vec4 norm = taylorInvSqrt(vec4(dot(p0,p0), dot(p1,p1), dot(p2, p2), dot(p3,p3)));\r\n p0 *= norm.x;\r\n p1 *= norm.y;\r\n p2 *= norm.z;\r\n p3 *= norm.w;\r\n \r\n vec4 values = vec4( dot(p0,x0), dot(p1,x1), dot(p2,x2), dot(p3,x3) ); //value of contributions from each corner (extrapolate the gradient)\r\n \r\n vec4 m = max(0.5 - vec4(dot(x0,x0), dot(x1,x1), dot(x2,x2), dot(x3,x3)), 0.0); //kernel function from each corner\r\n \r\n vec4 m2 = m * m;\r\n vec4 m3 = m * m * m;\r\n \r\n vec4 temp = -6.0 * m2 * values;\r\n float dx = temp[0] * x0.x + temp[1] * x1.x + temp[2] * x2.x + temp[3] * x3.x + m3[0] * p0.x + m3[1] * p1.x + m3[2] * p2.x + m3[3] * p3.x;\r\n float dy = temp[0] * x0.y + temp[1] * x1.y + temp[2] * x2.y + temp[3] * x3.y + m3[0] * p0.y + m3[1] * p1.y + m3[2] * p2.y + m3[3] * p3.y;\r\n float dz = temp[0] * x0.z + temp[1] * x1.z + temp[2] * x2.z + temp[3] * x3.z + m3[0] * p0.z + m3[1] * p1.z + m3[2] * p2.z + m3[3] * p3.z;\r\n \r\n return vec4(dot(m3, values), dx, dy, dz) * 42.0;\r\n}\r\n\r\n\r\nvec3 curlNoise3 (vec3 p) {\r\n\r\n vec3 xNoisePotentialDerivatives = snoiseD( p ).yzw; //yzw are the xyz derivatives\r\n vec3 yNoisePotentialDerivatives = snoiseD(vec3( p.y - 19.1 , p.z + 33.4 , p.x + 47.2 )).zwy;\r\n vec3 zNoisePotentialDerivatives = snoiseD(vec3( p.z + 74.2 , p.x - 124.5 , p.y + 99.4 )).wyz;\r\n vec3 noiseVelocity = vec3(\r\n zNoisePotentialDerivatives.y - yNoisePotentialDerivatives.z,\r\n xNoisePotentialDerivatives.z - zNoisePotentialDerivatives.x,\r\n yNoisePotentialDerivatives.x - xNoisePotentialDerivatives.y\r\n );\r\n\t\r\n\tconst float e = .1;\r\n\tconst float divisor = 1.0 / ( 2.0 * e );\r\n\treturn normalize( noiseVelocity * divisor );\r\n\r\n}\r\n\t\r\n\t"; Main.main(); })();