GLSL Sandboxで波動方程式を解いてみた

きっかけ

先に言っておきますが、私は今日は間に合いませんので誰かお願いします

— h_doxas (@h_doxas) 2016年12月6日

だれか！ お客様のなかに、シェーダ書きたくてしょうがないマンはいらっしゃいませんか！？！？！？！？

— h_doxas (@h_doxas) 2016年12月6日

しゃあない雑に書くか

参考

動作

GLSL Sandbox

コード

#ifdef GL_ES
precision mediump float;
#endif

#extension GL_OES_standard_derivatives : enable

uniform float time;
uniform vec2 mouse;
uniform vec2 resolution;
uniform sampler2D backbuffer;

vec2 tex(vec2 uv)
{
	return texture2D(backbuffer, uv).xy - 0.5;
}

void main( void ) {

	vec2 pos = ( gl_FragCoord.xy - resolution.xy / 2.0) / resolution.y - mouse + 0.5;
	vec2 uv =  ( gl_FragCoord.xy / resolution.xy );
	vec2 prev = tex(uv);
	vec2 pixel = 1./resolution;

	// ラプラシアンフィルタで加速度を計算
	float accel =
		tex(uv + pixel * vec2(1, 0)).x +
		tex(uv - pixel * vec2(1, 0)).x +
		tex(uv + pixel * vec2(0, 1)).x +
		tex(uv - pixel * vec2(0, 1)).x -
		prev.x * 4.0;

	// 伝播速度を掛ける
	accel *= 0.2;

	// 現在の速度に加速度を足し、さらに減衰率を掛ける
	float velocity = (prev.y + accel) * 0.95;

	// 高さを更新
	float height = prev.x + velocity;

	// マウス位置に波紋を出す
	height += max(0.0, 1.0 - length(pos) * 30.0);
	
	gl_FragColor = vec4(height + 0.5, velocity + 0.5, 0, 1);

}