fp64_escaped_geom.glsl

Go to the documentation of this file.

// mightymandel -- GPU-based Mandelbrot Set explorer
// Copyright (C) 2012,2013,2014,2015 Claude Heiland-Allen
// License GPL3+ http://www.gnu.org/licenses/gpl.html

layout(points) in;
layout(points, max_vertices = 1) out;
uniform double loger2;
uniform double pxs;
uniform dvec2 center;
uniform double radius;
uniform double aspect;
in  dvec4 cne1[1];
#ifdef DE
in  dvec4 zdz1[1];
#else
in  dvec2 zdz1[1];
#endif
flat out vec3 ida1;
#ifdef DE
flat out dvec4 zdz;
#else
flat out dvec2 zdz;
#endif
void main() {
  bool escaped = ! (cne1[0].w <= 0.0);
  if (escaped) {
    double z2 = cmag2(zdz1[0].xy);
    double logz2 = double(log(float(z2)));
    float i = float(1.0 + cne1[0].z - clamp(double(log2(float(logz2 / loger2))), 0.00001, 0.99999));
#ifdef DE
    double dz2 = cmag2(zdz1[0].zw);
    float d = float(logz2 * sqrt(z2 / dz2) * pxs);
#else
    float d = 0.0;
#endif
    ida1 = vec3(i, 0.0, d);
    gl_Position = vec4(vec2((cne1[0].xy - center) / radius * dvec2(aspect, 1.0)), 0.0, 1.0);
    EmitVertex();
    EndPrimitive();
  }
}