OpenGL4.0 Shader 绘制MandelBrot集合

OpenGL4.0 Shader 绘制的MandelBrot集合

OpenGL4.0 Shader 绘制的MandelBrot集合，先上一幅图

这个图形是用OpenGL 的Fragment计算生成的。CSDN不支持OpenGL Shader类型代码，这里用C类型贴出来试试。

#version 400
//#extension GL_ARB_gpu_shader_fp64 enable
uniform vec2 viewport;
uniform dvec2 pfrom;
uniform double psize;

#define GAP   20
#define INDSIZE 32
#define LIMIT GAP * GAP * GAP

in vec2 guv;
void main() {
    int ix, iy;
	ix = int(floor(guv.x * 1024));
	iy = int(floor(guv.y * 1024));
	if (ix >= (512 - 1) && ix <= (512 + 1) && (iy >= 512-1) && (iy <= 512+1) ) {
	 gl_FragColor = vec4(1.0f, 1.0f, 1.0f, 1.0f);
	} else {
	dvec2 z, c;
	int i;

	z = dvec2(0.0lf, 0.0lf);
	c = dvec2(double(guv.x) - 0.5lf, double(guv.y) - 0.5lf);
	c *= psize;
	c += pfrom;

	for (i = 0;i<LIMIT;i++) {
	    double zx = z.x;
		z.x = z.x * z.x - z.y*z.y + c.x;
		z.y = 2 * zx * z.y + c.y;
		if (abs(z.x) >= 2.0lf || abs(z.y) >= 2.0lf)
		break;
	}

	if (i == LIMIT) {
		gl_FragColor = vec4(0.1, 0.1, 0.3, 1.0);
	} else {
		float ir = (i % GAP) / float(GAP);
		float ig = ((i / GAP) % GAP) / float(GAP);
		float ib = ((i / (GAP * GAP)) % GAP) / float(GAP);

		gl_FragColor = vec4(ir, ig, ib, 1.0f);	
	}
	}	       
}

$$\begin{gathered} MandelBrot集合是由迭代方程z=z^2+c生成的， \\ 其中z和c都是复平面上的点,z的初值为0，容易证明 \\ |c|>2时，迭代会发散到无穷，后面证明|c|<1/4时 \\ 迭代收敛，计算时可以依此做一些优化工作，减少计算量。 \end{gathered}$$
下面时c=-0.683920704845814时的迭代结果：

序号	迭代次数	迭代值
1	0	-0.683920704845814
2	1	-0.216173174329018
3	64	-0.468359644296487
4	65	-0.464559948440282
5	66	-0.468104759150976
6	67	-0.464798639306020
7	68	-0.467882929745086
8	69	-0.465006268898969
9	70	-0.467689874730474
10	71	-0.465186885920408
11	72	-0.467521866013488
12	73	-0.465344009645080
13	74	-0.467375657533253
14	75	-0.465480699591173
15	76	-0.467248423153926
16	77	-0.465599615905984
17	78	-0.467137702514014
18	79	-0.465703071735742
19	99984	-0.466395728904994
20	99985	-0.466395728904993
21	99986	-0.466395728904994
22	99987	-0.466395728904993
23	99988	-0.466395728904994
24	99989	-0.466395728904993
25	99990	-0.466395728904994
26	99991	-0.466395728904993
27	99992	-0.466395728904994
28	99993	-0.466395728904993
29	99994	-0.466395728904994
30	99995	-0.466395728904993
31	99996	-0.466395728904994
32	99997	-0.466395728904993
33	99998	-0.466395728904994
34	99999	-0.466395728904993

后面应该时收敛到一个数了，最低位出现跳动应该时计算精度不够导致的。

以下来证明|c|<1/4时，迭代收敛，只需对c时正实数时证明即可：
$$\begin{gathered} \begin{array}{rl}对z& =0,c为正实数而且c<1/4，迭代：\\ z_n& =z_{n-1}^2+c，\\ 求证& ：z_n收敛\\ 证明& ：首先证明z_n<1/2，用数学归纳法，\end{array} \\ \begin{array}{rl}显然：z_0& =0<1/2,\\ z_1& <1/4<1/2\\ 假设：z_{n-1}& <1/2，\\ 记：z_{n-1}& =1/2-x,\\ 则：z_n& =z_{n-1}^2+c\\ & =1/4-x+x^2+c\\ & <1/4-x+x^2+1/4\\ & =1/2+x(x-1)\end{array} \\ 注意到x在（0,1/2）中，所以x(x-1)<0,因此z_n<1/2 \\ 其次证明z_n单调递增。注意到， \\ \begin{array}{rl}{z}_n-z_{n-1}& =(z_{n-1}^2+c)-(z_{n-2}^2+c)\\ & =z_{n-1}^2-z_{n-2}^2，\\ 可知z_n-z_{n-1}& 与z_{n-1}-z_{n-2}同号，注意到z_1>z_0，\end{array} \\ 用数学归纳法同样证明z_n单调递增。 \\ 单调递增数列必收敛。证毕。 \end{gathered}$$

原文有绘制的c语言代码和对应的OpenGL Shader代码，发布在微信公众号上。点此查看微信公众号原文。
以上测试了CSDN对图片显示，源代码代码，表格，数学公式等方面的支持，功能还是不错的。后面发类似的文干脆搬这边来算了。数学公式这边对齐是个问题，估计是由于对公式编辑这部分还不熟悉，后面有机会好好学习一下。