水面波动

[LeeMcQueen]

uniform mat4 viewProjectMatrix; uniform mat4 viewMatrix; uniform float time; uniform vec3 eyePos;

layout (location = 0) in vec3 vertex;

out vec3 vNormal; out vec3 vViewNormal; out vec3 vEye2Water; out vec4 vProjPos;

vec3 CalculateWavePosition(float q, float a, float w, vec3 dir, vec3 meshVert, float ph, float t) { float qa = q a; float theta = dot(w dir.xz, meshVert.xz) + ph t; float cosTh = cos(theta); float sinTh = sin(theta); float x = qa dir.x cosTh; float z = qa dir.z cosTh; float y = a sinTh; return vec3(x, y, z); }

vec3 CalculateWaveNormal(float q, float a, float w, vec3 dir, vec3 waveVert, float ph, float t) { float wa = w a; float theta = dot(w dir, waveVert) + ph t; float cosTh = cos(theta); float sinTh = sin(theta); float x = dir.x wa cosTh; float z = dir.z wa cosTh; float y = q wa * sinTh; return vec3(x, y, z); }

void main() { float l0 = 31.25, a0 = 0.16, s0 = 2.56; vec3 dir0 = vec3(0.58, 0.0, 0.42); float w0 = 2.0 / l0, ph0 = s0 * w0, q0 = 1.28;

float l1 = 25.0, a1 = 0.22, s1 = 5.12;
vec3 dir1 = vec3(0.31, 0.0, 0.69);
float w1 = 2.0 / l1, ph1 = s1 * w1, q1 = 2.56;

float l2 = 25.6, a2 = 0.22, s2 = 1.28;
vec3 dir2 = vec3(0.33, 0.0, 0.67);
float w2 = 2.0 / l2, ph2 = s2 * w2, q2 = 2.56;

float l3 = 52.5, a3 = 0.34, s3 = 0.64;
vec3 dir3 = vec3(0.26, 0.0, 0.74);
float w3 = 2.0 / l3, ph3 = s3 * w3, q3 = 5.12;

float l4 = 25.6, a4 = 0.34, s4 = 2.56;
vec3 dir4 = vec3(0.3, 0.0, -0.7);
float w4 = 2.0 / l4, ph4 = s4 * w4, q4 = 5.12;

float l5 = 42.5, a5 = 0.46, s5 = 5.12;
vec3 dir5 = vec3(0.4, 0.0, -0.6);
float w5 = 2.0 / l5, ph5 = s5 * w5, q5 = 2.56;

float l6 = 25.0, a6 = 0.22, s6 = 1.28;
vec3 dir6 = vec3(0.1, 0.0, -0.9);
float w6 = 2.0 / l6, ph6 = s6 * w6, q6 = 2.56;

float l7 = 31.25, a7 = 0.16, s7 = 0.64;
vec3 dir7 = vec3(0.43, 0.0, -0.57);
float w7 = 2.0 / l7, ph7 = s7 * w7, q7 = 1.28;

vec3 worldVertex = vertex + vec3(eyePos.x, 0.0, eyePos.z);

vec3 pos0 = CalculateWavePosition(q0, a0, w0, dir0, worldVertex, ph0, time);
vec3 pos1 = CalculateWavePosition(q1, a1, w1, dir1, worldVertex, ph1, time);
vec3 pos2 = CalculateWavePosition(q2, a2, w2, dir2, worldVertex, ph2, time); 
vec3 pos3 = CalculateWavePosition(q3, a3, w3, dir3, worldVertex, ph3, time); 
vec3 pos4 = CalculateWavePosition(q4, a4, w4, dir4, worldVertex, ph4, time); 
vec3 pos5 = CalculateWavePosition(q5, a5, w5, dir5, worldVertex, ph5, time);
vec3 pos6 = CalculateWavePosition(q6, a6, w6, dir6, worldVertex, ph6, time);
vec3 pos7 = CalculateWavePosition(q7, a7, w7, dir7, worldVertex, ph7, time);

vec3 position = pos0 + pos1 + pos2 + pos3 + pos4 + pos5 + pos6 + pos7;
position.xz += worldVertex.xz;

vec3 nor0 = CalculateWaveNormal(q0, a0, w0, dir0, position, ph0, time);
vec3 nor1 = CalculateWaveNormal(q1, a1, w1, dir1, position, ph1, time);
vec3 nor2 = CalculateWaveNormal(q2, a2, w2, dir2, position, ph2, time);
vec3 nor3 = CalculateWaveNormal(q3, a3, w3, dir3, position, ph3, time);
vec3 nor4 = CalculateWaveNormal(q4, a4, w4, dir4, position, ph4, time);
vec3 nor5 = CalculateWaveNormal(q5, a5, w5, dir5, position, ph5, time);
vec3 nor6 = CalculateWaveNormal(q6, a6, w6, dir6, position, ph6, time);
vec3 nor7 = CalculateWaveNormal(q7, a7, w7, dir7, position, ph7, time);

vec3 normal = nor0 + nor1 + nor2 + nor3 + nor4 + nor5 + nor6 + nor7;
normal = vec3(0.0, 1.0, 0.0) - normal;

vNormal = normal;
vViewNormal = mat3(viewMatrix) * normal;
vEye2Water = position - eyePos;

vProjPos = viewProjectMatrix * vec4(position, 1.0);
gl_Position = vProjPos;

}