如何表达复数的四聚函数
How to express tetration function, for complex numbers
确实存在所谓的hyperoperation sequence. It works like you construct multiplication a*b=a+a+a+a...+a with many additions of a repeated b times. Then there goes exponentiation a^b = a*a*a*a*...*a with many multiplicaitions of a repeated b times. Then, there goes tetration，表示为指数塔，与a^^b == a^a^a^...^a相同，重复b次。
我有兴趣如何编写此函数，用于浮点数和复数？
我已经用 glsl 编写了乘法和求幂函数:
// complex multiplication:
vec2 cmul(in vec2 a, in vec2 b) {
return vec2(a.x*b.x-a.y*b.y, a.x*b.y+a.y*b.x);
}
// complex exponent e^a
vec2 cexp(in vec2 a) {
float ea = exp(a.x);
float vl = a.y;
return ea * vec2( cos(vl), sin(vl) );
}
// complex natural logarithm ln(a)
vec2 cln(in vec2 a) {
float ql = length(a);
return vec2( log(ql), atan(a.y, a.x));
}
// complex power function a^b
vec2 cpow(in vec2 a, in vec2 b) {
return cexp(cmul(cln(a), b));
}
但后来我卡住了！我们如何编写 ctet(in vec2 a, in vec2 b) tetration 函数，不仅针对浮点数，而且针对整个复平面本身？
好吧，让我们从 真实域 和 整数 b 开始：
a^^b = a^a^a^a^a...^a // a is there b times
在C++中可以这样计算：
double tetration(double a,int b) // a^^b = a^a^a^a... b times
{
double c;
if (b<=0) return 0;
for (c=a;b>1;b--) c=pow(a,c);
return c;
}
因为您已经获得了复杂域的 pow，您也可以在那里做同样的事情...为了简单起见，我暂时不会触及它...
这里是一些结果：
a\b| 1| 2| 3| 4
-------------------------------
1 | 1| 1| 1| 1
2 | 2| 4| 16|65536
3 | 3| 27|7625597484987|
4 | 4| 256| |
5 | 5|3125| |
顺便说一句。 所有这些 超级操作 都与 Ackermann function 有关，您可以在 C++ 中找到我的迭代实现这里：
然而，由于极快的增长，甚至 double 也将很快超出范围（因此缺少值）...
现在如何将b移动到真实域？对此没有代数方法的线索，但几何方法是可能的。
简单地“绘制”a^^b 作为变量 b 和常量 a 的函数，用于围绕您想要的实数 b 的整数值 b 和然后使用整数域 b作为控制点插入实域b。类似于函数的非整数阶求导
所以 (X,Y) 将是您的 (a^^b,b)。现在使用任意插值构造实域函数。
线性插值将如下所示：
y0 = a^^(int(b))
y1 = a^^(int(b)+1)
a^^b = y0 + (b-int(b))*(y1-y0)
然而，需要更高阶的插值，而且插值参数应该缩放到非线性度量。有关详细信息，请参阅：
经过一些详细说明，立方体 (t^3) 和 log^2 比例被证明是足够的（C++ 示例使用我的 128 位浮点数 f128 class 只需将其重命名为double):
f128 tetration_fi(f128 a,int b) // a^^b = a^a^a^a... b times
{
f128 c;
if (b==-1) return 0.0; // first singularity
if (b== 0) return 1.0; // second singularity
if (b< -1) return 0.0; // uncomputed
for (c=a;b>1;b--) c=pow(a,c);
return c;
}
//---------------------------------------------------------------------------
f128 tetration_ff(f128 a,f128 b) // a^^b = a^a^a^a... b times
{
int bi;
f128 z0,z1,z2,z3,a0,a1,a2,a3,t,tt,ttt,o=2.0;
if (b==-1) return 0.0; // first singularity
if (b== 0) return 1.0; // second singularity
if (b< -1) return 0.0; // uncomputed
bi=b.ToInt(); b-=bi;
if (b.iszero()) return tetration_fi(a,bi);
z0=tetration_fi(a,bi-1); // known points around a^^b
z1=pow(a,z0);
z2=pow(a,z1);
z3=pow(a,z2);
z0=log2(log2(z0+o)+o); // log^2 scale
z1=log2(log2(z1+o)+o);
z2=log2(log2(z2+o)+o);
z3=log2(log2(z2+o)+o);
t =0.5*(z2-z0); // cubic interpolation coeff.
tt=0.5*(z3-z1);
a0=z1;
a1=t;
a2=(3.0*(z2-z1))-(2.0*t)-tt;
a3=t+tt+(2.0*(z1-z2));
t=b-bi; // cubic interpolation
tt=t*t;
ttt=tt*t;
z0=a0+(a1*t)+(a2*t*t)+(a3*t*t*t);
z0=exp2(exp2(z0)-o)-o; // linear scale
return z0;
}
//---------------------------------------------------------------------------
这是我比较的：
我 select 相同的图表基于 a^^b 的 a，正如您所看到的，它非常匹配，只有低于 1.0 的范围略有偏差。
让我们开始复杂域分形
现在当你想进入复杂域时，你不能像在Real中那样做，因为结果太混乱而无法插值。所以只能坚持整数b或者用Kneser算法来计算。
幸运的是，我们有更多的方法来显示分形...例如，我们可以从 a^^b 计算整数 b，其中只有 a 是复数，并使用为输出着色的结果。这里 GLSL 示例（基于我的 和你的复杂数学）：
片段：
// Fragment
#version 450 core
uniform dvec2 p0=dvec2(0.0,0.0); // mouse position <-1,+1>
uniform double zoom=1.000; // zoom [-]
in smooth vec2 p32;
out vec4 col;
//---------------------------------------------------------------------------
// All components are in the range [0…1], including hue.
vec3 rgb2hsv(vec3 c)
{
vec4 K = vec4(0.0, -1.0 / 3.0, 2.0 / 3.0, -1.0);
vec4 p = mix(vec4(c.bg, K.wz), vec4(c.gb, K.xy), step(c.b, c.g));
vec4 q = mix(vec4(p.xyw, c.r), vec4(c.r, p.yzx), step(p.x, c.r));
float d = q.x - min(q.w, q.y);
float e = 1.0e-10;
return vec3(abs(q.z + (q.w - q.y) / (6.0 * d + e)), d / (q.x + e), q.x);
}
//---------------------------------------------------------------------------
// All components are in the range [0…1], including hue.
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);
}
//---------------------------------------------------------------------------
vec3 spectral_color(float l) // RGB <0,1> <- lambda l <400,700> [nm]
{
float t; vec3 c=vec3(0.0,0.0,0.0);
if ((l>=400.0)&&(l<410.0)) { t=(l-400.0)/(410.0-400.0); c.r= +(0.33*t)-(0.20*t*t); }
else if ((l>=410.0)&&(l<475.0)) { t=(l-410.0)/(475.0-410.0); c.r=0.14 -(0.13*t*t); }
else if ((l>=545.0)&&(l<595.0)) { t=(l-545.0)/(595.0-545.0); c.r= +(1.98*t)-( t*t); }
else if ((l>=595.0)&&(l<650.0)) { t=(l-595.0)/(650.0-595.0); c.r=0.98+(0.06*t)-(0.40*t*t); }
else if ((l>=650.0)&&(l<700.0)) { t=(l-650.0)/(700.0-650.0); c.r=0.65-(0.84*t)+(0.20*t*t); }
if ((l>=415.0)&&(l<475.0)) { t=(l-415.0)/(475.0-415.0); c.g= +(0.80*t*t); }
else if ((l>=475.0)&&(l<590.0)) { t=(l-475.0)/(590.0-475.0); c.g=0.8 +(0.76*t)-(0.80*t*t); }
else if ((l>=585.0)&&(l<639.0)) { t=(l-585.0)/(639.0-585.0); c.g=0.84-(0.84*t) ; }
if ((l>=400.0)&&(l<475.0)) { t=(l-400.0)/(475.0-400.0); c.b= +(2.20*t)-(1.50*t*t); }
else if ((l>=475.0)&&(l<560.0)) { t=(l-475.0)/(560.0-475.0); c.b=0.7 -( t)+(0.30*t*t); }
return c;
}
//---------------------------------------------------------------------------
// complex domain math
vec3 color_wheel(vec2 a) // complex -> polar -> HSV -> RGB
{
float an=(atan(-a.y,-a.x)*0.15915494309189533576888376337251)+0.5;
float r=length(a); r-=floor(r); r*=0.75; r+=0.25;
return hsv2rgb(vec3(an,1.0,r));
}
vec3 color_spectral(vec2 a) // complex -> wavelength -> RGB
{
float r=length(a); r-=floor(r);
return spectral_color(400.0+(300.0*r));
}
vec2 cadd(vec2 a,vec2 b) // a+b
{
return a+b;
}
vec2 csub(vec2 a,vec2 b) // a-b
{
return a-b;
}
vec2 cmul(vec2 a,vec2 b) // a*b
{
return vec2((a.x*b.x)-(a.y*b.y),(a.x*b.y)+(a.y*b.x));
}
vec2 cdiv(vec2 a,vec2 b) // a/b
{
float an=atan(-a.y,-a.x)-atan(-b.y,-b.x);
float r=length(a)/length(b);
return r*vec2(cos(an),sin(an));
}
vec2 csqr(vec2 a) // a^2
{
return cmul(a,a);
}
vec2 cexp(vec2 a) // e^a
{
// e^(x+y*i)= e^x * e^(y*i) = e^x * ( cos(y) + i*sin(y) )
return exp(a.x)*vec2(cos(a.y),sin(a.y));
}
vec2 cln(vec2 a) // ln(a)
{
return vec2(log(length(a)),atan(a.y,a.x));
}
vec2 cpow(vec2 a,vec2 b) // a^b
{
return cexp(cmul(cln(a),b));
}
vec2 ctet(vec2 a,int b) // a^^b
{
vec2 c=vec2(1.0,0.0);
for (;b>0;b--) c=cpow(a,c);
return c;
}
//---------------------------------------------------------------------------
void main()
{
// poistion (double)
dvec2 p=dvec2(p32);
p=(p/zoom)-p0; // x,y (-1.0, 1.0)
// position (float)
vec2 pp=vec2(p);
// [chose function]
// complex domain test function 1 (color wheel)
// vec2 a=cdiv(cmul(csub(cmul(pp,pp),vec2(1.0,0.0)),csqr(csub(pp,vec2(2.0,1.0)))),cadd(csqr(pp),vec2(2.0,2.0)));
// complex domain test function 2 (color wheel)
// vec2 a=pp; a=cln(a);
// exponentiation escape fractal 1 (color wheel)
// vec2 a=cpow(pp,vec2(100,0));
// exponentiation escape fractal 2 (color wheel)
// vec2 a=vec2(1.0,1.0); for (int i=0;i<100;i++) a=cpow(a,pp);
// exponentiation escape fractal 3 (color wheel)
// vec2 a=vec2(0.0,0.0),b=vec2(1.0,0.0); float r=0.5,rr=1.0,wt=0.1; for (int i=0;i<20;i++){ a+=rr*cexp(vec2(-b.y,b.x)*wt); b=cmul(b,pp); rr*=r; } a*=(1.0-r);
// tetration escape fractal (grayscale)
// vec2 a=ctet(pp,100);
// pentation escape fractal (grayscale)
vec2 a=pp; for (int i=0;i<20;i++) a=ctet(a,20); a*=100.0;
// [chose coloring method]
// grayscale based on escape
float r=0.2*length(a); r-=floor(r); r=0.25+0.75*r; col=vec4(r,r,r,1.0);
// RGB based on result
// col=vec4(a,a.x+a.y,1.0);
// result -> wavelength+intensity
// col=vec4(color_wheel(a),1.0);
// result -> spectral color
// col=vec4(color_spectral(a),1.0);
}
和四舍五入预览：
这是我对比的：
它与我的结果相匹配 x,y
所以我所做的是计算 a^^100，其中 a 是屏幕上片段的复杂域位置 <-1,+1> 以及一些 panning 和 zooming 并渲染根据结果构建的颜色 ...
我还留下了一个测试函数（不是分形）我用来测试着色方法和复杂的数学从 here 第一个来自维基，第二个是着色器结果（色轮）：
您可以像 mandelbrot 或其他任何算法一样进行逃逸测试来代替显示分形。
这里是 zoom=500.0 pos=-0.188418+0.234466i
的着色选项截图（我喜欢灰度）
最后五字：
确实存在所谓的hyperoperation sequence. It works like you construct multiplication a*b=a+a+a+a...+a with many additions of a repeated b times. Then there goes exponentiation a^b = a*a*a*a*...*a with many multiplicaitions of a repeated b times. Then, there goes tetration，表示为指数塔，与a^^b == a^a^a^...^a相同，重复b次。
我有兴趣如何编写此函数，用于浮点数和复数？
我已经用 glsl 编写了乘法和求幂函数:
// complex multiplication:
vec2 cmul(in vec2 a, in vec2 b) {
return vec2(a.x*b.x-a.y*b.y, a.x*b.y+a.y*b.x);
}
// complex exponent e^a
vec2 cexp(in vec2 a) {
float ea = exp(a.x);
float vl = a.y;
return ea * vec2( cos(vl), sin(vl) );
}
// complex natural logarithm ln(a)
vec2 cln(in vec2 a) {
float ql = length(a);
return vec2( log(ql), atan(a.y, a.x));
}
// complex power function a^b
vec2 cpow(in vec2 a, in vec2 b) {
return cexp(cmul(cln(a), b));
}
但后来我卡住了！我们如何编写 ctet(in vec2 a, in vec2 b) tetration 函数，不仅针对浮点数，而且针对整个复平面本身？
好吧，让我们从 真实域 和 整数 b 开始：
a^^b = a^a^a^a^a...^a // a is there b times
在C++中可以这样计算：
double tetration(double a,int b) // a^^b = a^a^a^a... b times
{
double c;
if (b<=0) return 0;
for (c=a;b>1;b--) c=pow(a,c);
return c;
}
因为您已经获得了复杂域的 pow，您也可以在那里做同样的事情...为了简单起见，我暂时不会触及它...
这里是一些结果：
a\b| 1| 2| 3| 4
-------------------------------
1 | 1| 1| 1| 1
2 | 2| 4| 16|65536
3 | 3| 27|7625597484987|
4 | 4| 256| |
5 | 5|3125| |
顺便说一句。 所有这些 超级操作 都与 Ackermann function 有关，您可以在 C++ 中找到我的迭代实现这里：
然而，由于极快的增长，甚至 double 也将很快超出范围（因此缺少值）...
现在如何将b移动到真实域？对此没有代数方法的线索，但几何方法是可能的。
简单地“绘制”a^^b 作为变量 b 和常量 a 的函数，用于围绕您想要的实数 b 的整数值 b 和然后使用整数域 b作为控制点插入实域b。类似于函数的非整数阶求导
所以 (X,Y) 将是您的 (a^^b,b)。现在使用任意插值构造实域函数。
线性插值将如下所示：
y0 = a^^(int(b))
y1 = a^^(int(b)+1)
a^^b = y0 + (b-int(b))*(y1-y0)
然而，需要更高阶的插值，而且插值参数应该缩放到非线性度量。有关详细信息，请参阅：
经过一些详细说明，立方体 (t^3) 和 log^2 比例被证明是足够的（C++ 示例使用我的 128 位浮点数 f128 class 只需将其重命名为double):
f128 tetration_fi(f128 a,int b) // a^^b = a^a^a^a... b times
{
f128 c;
if (b==-1) return 0.0; // first singularity
if (b== 0) return 1.0; // second singularity
if (b< -1) return 0.0; // uncomputed
for (c=a;b>1;b--) c=pow(a,c);
return c;
}
//---------------------------------------------------------------------------
f128 tetration_ff(f128 a,f128 b) // a^^b = a^a^a^a... b times
{
int bi;
f128 z0,z1,z2,z3,a0,a1,a2,a3,t,tt,ttt,o=2.0;
if (b==-1) return 0.0; // first singularity
if (b== 0) return 1.0; // second singularity
if (b< -1) return 0.0; // uncomputed
bi=b.ToInt(); b-=bi;
if (b.iszero()) return tetration_fi(a,bi);
z0=tetration_fi(a,bi-1); // known points around a^^b
z1=pow(a,z0);
z2=pow(a,z1);
z3=pow(a,z2);
z0=log2(log2(z0+o)+o); // log^2 scale
z1=log2(log2(z1+o)+o);
z2=log2(log2(z2+o)+o);
z3=log2(log2(z2+o)+o);
t =0.5*(z2-z0); // cubic interpolation coeff.
tt=0.5*(z3-z1);
a0=z1;
a1=t;
a2=(3.0*(z2-z1))-(2.0*t)-tt;
a3=t+tt+(2.0*(z1-z2));
t=b-bi; // cubic interpolation
tt=t*t;
ttt=tt*t;
z0=a0+(a1*t)+(a2*t*t)+(a3*t*t*t);
z0=exp2(exp2(z0)-o)-o; // linear scale
return z0;
}
//---------------------------------------------------------------------------
这是我比较的：
我 select 相同的图表基于 a^^b 的 a，正如您所看到的，它非常匹配，只有低于 1.0 的范围略有偏差。
让我们开始复杂域分形
现在当你想进入复杂域时，你不能像在Real中那样做，因为结果太混乱而无法插值。所以只能坚持整数b或者用Kneser算法来计算。
幸运的是，我们有更多的方法来显示分形...例如，我们可以从 a^^b 计算整数 b，其中只有 a 是复数，并使用为输出着色的结果。这里 GLSL 示例（基于我的
片段：
// Fragment
#version 450 core
uniform dvec2 p0=dvec2(0.0,0.0); // mouse position <-1,+1>
uniform double zoom=1.000; // zoom [-]
in smooth vec2 p32;
out vec4 col;
//---------------------------------------------------------------------------
// All components are in the range [0…1], including hue.
vec3 rgb2hsv(vec3 c)
{
vec4 K = vec4(0.0, -1.0 / 3.0, 2.0 / 3.0, -1.0);
vec4 p = mix(vec4(c.bg, K.wz), vec4(c.gb, K.xy), step(c.b, c.g));
vec4 q = mix(vec4(p.xyw, c.r), vec4(c.r, p.yzx), step(p.x, c.r));
float d = q.x - min(q.w, q.y);
float e = 1.0e-10;
return vec3(abs(q.z + (q.w - q.y) / (6.0 * d + e)), d / (q.x + e), q.x);
}
//---------------------------------------------------------------------------
// All components are in the range [0…1], including hue.
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);
}
//---------------------------------------------------------------------------
vec3 spectral_color(float l) // RGB <0,1> <- lambda l <400,700> [nm]
{
float t; vec3 c=vec3(0.0,0.0,0.0);
if ((l>=400.0)&&(l<410.0)) { t=(l-400.0)/(410.0-400.0); c.r= +(0.33*t)-(0.20*t*t); }
else if ((l>=410.0)&&(l<475.0)) { t=(l-410.0)/(475.0-410.0); c.r=0.14 -(0.13*t*t); }
else if ((l>=545.0)&&(l<595.0)) { t=(l-545.0)/(595.0-545.0); c.r= +(1.98*t)-( t*t); }
else if ((l>=595.0)&&(l<650.0)) { t=(l-595.0)/(650.0-595.0); c.r=0.98+(0.06*t)-(0.40*t*t); }
else if ((l>=650.0)&&(l<700.0)) { t=(l-650.0)/(700.0-650.0); c.r=0.65-(0.84*t)+(0.20*t*t); }
if ((l>=415.0)&&(l<475.0)) { t=(l-415.0)/(475.0-415.0); c.g= +(0.80*t*t); }
else if ((l>=475.0)&&(l<590.0)) { t=(l-475.0)/(590.0-475.0); c.g=0.8 +(0.76*t)-(0.80*t*t); }
else if ((l>=585.0)&&(l<639.0)) { t=(l-585.0)/(639.0-585.0); c.g=0.84-(0.84*t) ; }
if ((l>=400.0)&&(l<475.0)) { t=(l-400.0)/(475.0-400.0); c.b= +(2.20*t)-(1.50*t*t); }
else if ((l>=475.0)&&(l<560.0)) { t=(l-475.0)/(560.0-475.0); c.b=0.7 -( t)+(0.30*t*t); }
return c;
}
//---------------------------------------------------------------------------
// complex domain math
vec3 color_wheel(vec2 a) // complex -> polar -> HSV -> RGB
{
float an=(atan(-a.y,-a.x)*0.15915494309189533576888376337251)+0.5;
float r=length(a); r-=floor(r); r*=0.75; r+=0.25;
return hsv2rgb(vec3(an,1.0,r));
}
vec3 color_spectral(vec2 a) // complex -> wavelength -> RGB
{
float r=length(a); r-=floor(r);
return spectral_color(400.0+(300.0*r));
}
vec2 cadd(vec2 a,vec2 b) // a+b
{
return a+b;
}
vec2 csub(vec2 a,vec2 b) // a-b
{
return a-b;
}
vec2 cmul(vec2 a,vec2 b) // a*b
{
return vec2((a.x*b.x)-(a.y*b.y),(a.x*b.y)+(a.y*b.x));
}
vec2 cdiv(vec2 a,vec2 b) // a/b
{
float an=atan(-a.y,-a.x)-atan(-b.y,-b.x);
float r=length(a)/length(b);
return r*vec2(cos(an),sin(an));
}
vec2 csqr(vec2 a) // a^2
{
return cmul(a,a);
}
vec2 cexp(vec2 a) // e^a
{
// e^(x+y*i)= e^x * e^(y*i) = e^x * ( cos(y) + i*sin(y) )
return exp(a.x)*vec2(cos(a.y),sin(a.y));
}
vec2 cln(vec2 a) // ln(a)
{
return vec2(log(length(a)),atan(a.y,a.x));
}
vec2 cpow(vec2 a,vec2 b) // a^b
{
return cexp(cmul(cln(a),b));
}
vec2 ctet(vec2 a,int b) // a^^b
{
vec2 c=vec2(1.0,0.0);
for (;b>0;b--) c=cpow(a,c);
return c;
}
//---------------------------------------------------------------------------
void main()
{
// poistion (double)
dvec2 p=dvec2(p32);
p=(p/zoom)-p0; // x,y (-1.0, 1.0)
// position (float)
vec2 pp=vec2(p);
// [chose function]
// complex domain test function 1 (color wheel)
// vec2 a=cdiv(cmul(csub(cmul(pp,pp),vec2(1.0,0.0)),csqr(csub(pp,vec2(2.0,1.0)))),cadd(csqr(pp),vec2(2.0,2.0)));
// complex domain test function 2 (color wheel)
// vec2 a=pp; a=cln(a);
// exponentiation escape fractal 1 (color wheel)
// vec2 a=cpow(pp,vec2(100,0));
// exponentiation escape fractal 2 (color wheel)
// vec2 a=vec2(1.0,1.0); for (int i=0;i<100;i++) a=cpow(a,pp);
// exponentiation escape fractal 3 (color wheel)
// vec2 a=vec2(0.0,0.0),b=vec2(1.0,0.0); float r=0.5,rr=1.0,wt=0.1; for (int i=0;i<20;i++){ a+=rr*cexp(vec2(-b.y,b.x)*wt); b=cmul(b,pp); rr*=r; } a*=(1.0-r);
// tetration escape fractal (grayscale)
// vec2 a=ctet(pp,100);
// pentation escape fractal (grayscale)
vec2 a=pp; for (int i=0;i<20;i++) a=ctet(a,20); a*=100.0;
// [chose coloring method]
// grayscale based on escape
float r=0.2*length(a); r-=floor(r); r=0.25+0.75*r; col=vec4(r,r,r,1.0);
// RGB based on result
// col=vec4(a,a.x+a.y,1.0);
// result -> wavelength+intensity
// col=vec4(color_wheel(a),1.0);
// result -> spectral color
// col=vec4(color_spectral(a),1.0);
}
和四舍五入预览：
这是我对比的：
它与我的结果相匹配 x,y
所以我所做的是计算 a^^100，其中 a 是屏幕上片段的复杂域位置 <-1,+1> 以及一些 panning 和 zooming 并渲染根据结果构建的颜色 ...
我还留下了一个测试函数（不是分形）我用来测试着色方法和复杂的数学从 here 第一个来自维基，第二个是着色器结果（色轮）：
您可以像 mandelbrot 或其他任何算法一样进行逃逸测试来代替显示分形。
这里是 zoom=500.0 pos=-0.188418+0.234466i
最后五字：