[原创代码] [Perl]使用 POGL Shader 实现高效渲染 Julia集 动画

523066680

推荐 Perl 环境：5.24 Portable Edition
官方链接：Strawberry Perl Releases

Portable 版环境配置说明：Strawberry Perl 环境配置 以及 版本推荐

需要安装的模块：
https://metacpan.org/release/OpenGL
https://metacpan.org/pod/OpenGL::Shader

PDL版本的Strawberry Perl 虽然自带OPENGL模块，但是加载 Shader 有问题，最好下载源码包重新编译安装。

=info     Auth:523066680     Date:2017-07     https://www.shadertoy.com/view/ld2fzw     按 '[' 和 ']' 调整阀值，按 - 和 = 调整迭代深度     空格键 - 暂停/继续     q or Q - 退出 =cut use OpenGL qw/ :all /; use OpenGL::Config; use OpenGL::Shader; use Time::HiRes 'sleep'; use feature 'state'; use IO::Handle; STDOUT->autoflush(1); our $ang = 0.0; our $iter = 1000.0; our $test = 80.0; our ($cx, $cy); our $vTest; our $vIter; our $vC; our $PAUSE = 0; &Main(); sub display {     glClear(GL_COLOR_BUFFER_BIT);     glDrawArrays(GL_POINTS, 0, 1);     glutSwapBuffers(); } sub idle {     sleep 0.03;     $ang += 0.05 if ($PAUSE == 0);     $cx = cos($ang);     $cy = sin($ang);     glVertexAttrib2fARB($vC, $cx, $cy);     glutPostRedisplay(); } sub init {     glClearColor(0.0, 0.0, 0.0, 1.0);     my $shdr = new OpenGL::Shader('GLSL');     my $ver = $shdr->GetVersion();     my $stat = $shdr->Load( frag_Shader(), vert_Shader() );     $shdr->Enable();     #如果 shader 编译不成功，错误信息会返回到 $stat     print "$stat $ver\n";     glPointSize(480.0);         $vTest = $shdr->MapAttr('vTest');     $vIter = $shdr->MapAttr('vIter');     $vC = $shdr->MapAttr('vC');     glVertexAttrib1fARB($vTest, $test);     glVertexAttrib1fARB($vIter, $iter);     glVertexAttrib2fARB($vC, $cx, $cy); } sub hitkey {     my $key = shift;     my $char = chr($key);     if ( lc($char) eq 'q' ) { glutDestroyWindow($WinID) }     elsif ($char eq '[')  { $test -= 1.0 if ( $test > 1.0 ) }     elsif ($char eq ']')  { $test += 1.0 }     elsif ($char eq '-')  { $iter -= $iter*0.2 if ( $iter > 1.0 ) }     elsif ($char eq '=')  { $iter += $iter*0.2  }     elsif ($char eq ' ')  { $PAUSE = 1 - $PAUSE }     printf("test: %.2f, iter: %.2f\n", $test, $iter );     glVertexAttrib1fARB($vTest, $test);     glVertexAttrib1fARB($vIter, $iter); } sub Main {     glutInit();     glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE );     glutInitWindowSize(500, 500);     glutInitWindowPosition(1,1);     our $WinID = glutCreateWindow("title");     &init();     glutDisplayFunc(\&display);     glutKeyboardFunc(\&hitkey);     glutIdleFunc(\&idle);     glutMainLoop(); } sub frag_Shader {     return '     varying float fTest;     varying float fIter;     varying vec2 fC;     void main(void)     {         vec2 coord = (gl_FragCoord.xy - 250.0) / 100.0;         vec4 color;         vec2 C = fC;         vec2 Z = coord.xy;         int iterations = 0;         int max_iterations = int(fIter);         float rate;         float threshold_squared = fTest;         while ( (iterations < max_iterations) && (dot(Z, Z) < threshold_squared) )         {             vec2 tZ;             tZ.x = Z.x * Z.x - Z.y * Z.y + C.x;             tZ.y = Z.x * Z.y * 2.0 + C.y;             Z = tZ;             iterations++;         }         if (iterations == max_iterations)         {             color = vec4(0.3, 0.1, 0.1, 1.0);         }         else         {             rate =  float(iterations)/float(max_iterations)*100.0;             //color = vec4(rate, rate, rate, 1.0);  //default color             color = vec4(rate, smoothstep(exp(Z.x), 0.0, 0.5), rate, 1.0);         }         gl_FragColor = color;     }     ' } sub vert_Shader {     return '     attribute float vTest;     attribute float vIter;     attribute vec2 vC;     varying float fTest;     varying float fIter;     varying vec2 fC;     void main(void)     {         fTest = vTest;         fIter = vIter;         fC    = vC;         gl_Position = vec4(0.0, 0.0, 0.0, 1.0);     }     '; }

[Finished in 0.3s]

523066680

Sharp
  
迭代公式：
uZ.x = exp(tZ.x) * cos(tZ.y) * C.x;
uZ.y = exp(tZ.x) * sin(tZ.y) * C.y;

Julia之眼

迭代公式：
tZ.x = Z.x * Z.x - Z.y * Z.y + C.x;
tZ.y = Z.x * Z.y * 2.0 + C.y;
着色：
color = vec4(rate, exp(Z.x), rate, 1.0);

机械臂

tZ.x = Z.x * Z.x - Z.y * Z.y + C.x;
tZ.y = Z.x * Z.y * 1.0 + C.y;
着色：
color = vec4(rate, exp(Z.x), rate, 1.0);

图腾


Z=Z^3+C

浩瀚



旋转的阶梯

523066680

Colorful


  

Mandbrot: Z=Z^4+C

happy886rr

回复 3# 523066680
应该整理成屏保发布出来。

老刘1号

回复 4# happy886rr


    如果搞成屏保，那绝对霸气侧漏啊

523066680

回复 4# happy886rr

已上传 Win32 版本，支持全屏
http://523066680.ys168.com/
目录：OpenGL/2017年作品

不知道能不能正常运行，最好是有 Nvidia 显卡。 GT630测试可以
-
按键

    [ or ]  -> threshold_squared
    9 or 0  -> max_iteration
    - or =  -> scale

    w a s d -> move
    f       -> fullscreen/leave fullscreen
   空格键 -> 暂停
复制代码