HOME | DD

FlyingMatthew — Spiralworks

#fractalart #spiral #fractview
Published: 2021-03-01 13:39:54 +0000 UTC; Views: 783; Favourites: 25; Downloads: 1
Redirect to original
Description ((z^4.5+z)^2)+(z-c)

FractView for Android 

{"source":["// Min/Max Orbit trap: determines the minimum and maximum distance","// to the given trap.","var x int, y int, color int;","","extern maxdepth int \u003d 120;","","// some further arguments","extern juliaset bool \u003d false;","extern juliapoint cplx \u003d -0.8:0.16;","","func escapetime(c, breakcondition) {"," var i int \u003d 0,"," p cplx \u003d juliapoint if juliaset else c,"," zlast cplx \u003d 0,"," z cplx,"," znext cplx \u003d 0;",""," extern mandelinit expr \u003d \"0\";",""," z \u003d c if juliaset else mandelinit;",""," extern function expr \u003d \"mandelbrot(z, p)\";",""," var color quat;",""," while {"," znext \u003d function;"," not breakcondition(i, znext, z, zlast, c, p, color)"," } do {"," // advance to next values"," zlast \u003d z;"," z \u003d znext;"," }",""," // return color"," color","}","","func pipe(z, color) {"," var ret quat \u003d color;"," var v \u003d circlefn (1 - z);"," var a \u003d (200 - 4 ret.a);"," ret.a \u003d a sqr v + (100 - a) v;"," ret.d \u003d 1 - v (1 - opacity);"," ret","}","","func get_color(c, value) { "," var trapcolorfront quat \u003d 0:0:0:0; // later traps in front"," var trapcolorback quat \u003d 0:0:0:0; // later traps behind",""," var trapfirst real, trapindex \u003d -1;"," var traplast real;",""," // draw single pixel"," func breakcondition(i, znext, z, zlast, c, p, color) {"," // if true, color contains the lab color (quat)"," "," func bailoutcolor() {"," extern bailout real \u003d 128;"," extern max_power real \u003d 2;"," var smooth_i \u003d smoothen(znext, bailout, max_power) ;",""," // the next ones are only used in 3d-fractals"," extern bailoutvalue expr \u003d \"log(20 + i + smooth_i)\";"," value \u003d bailoutvalue ;",""," extern bailouttransfer expr \u003d \"value\";",""," extern bailoutpalette palette \u003d ["," [#006, #26c , #fff , #fa0 , #303]];"," "," color \u003d over(trapcolorfront, bailoutpalette bailouttransfer);"," if trapindex \u003e\u003d 0 then value \u003d traplast;"," }",""," func lakecolor() {"," extern epsilon real \u003d 1e-9;"," "," // the next ones are only used in 3d-fractals"," extern lakevalue expr \u003d \"log(1 + rad znext)\";"," value \u003d lakevalue;",""," extern laketransfer expr \u003d "," \"arcnorm znext : value\";",""," extern lakepalette palette \u003d ["," [#000, #000, #000, #000],"," [#f00 , #ff0 , #0f8 , #00f ],"," [#f88 , #ff8 , #afc , #88f ]];",""," color \u003d over(trapcolorback, lakepalette laketransfer);"," if trapindex \u003e\u003d 0 then value \u003d trapfirst;"," }",""," // update trap"," func updatetrap() {"," extern trapwidth real \u003d 0.1;",""," extern trapfn expr \u003d \"line(0:0, 1:0, znext)\";"," "," var trapdist real \u003d 1 - trapfn / trapwidth;",""," if trapdist \u003e 0 then {"," // we are inside the trap"," extern trapvalue expr \u003d \"trapdist\";"," "," extern traptransfer expr \u003d \"log(i + 1)\";"," "," extern trappalette palette \u003d ["," [#4b0082 , #800000, #ffd700 , #000080, #f5deb3 ]];",""," extern opacity real \u003d 0.67;"," "," var fg \u003d pipe(trapdist, trappalette traptransfer);"," "," trapcolorfront \u003d over(fg, trapcolorfront);"," trapcolorback \u003d over(trapcolorback, fg);"," "," if trapindex \u003c 0 then {"," trapfirst \u003d trapvalue;"," traplast \u003d trapfirst;"," trapindex \u003d i"," } else {"," traplast \u003d trapvalue;"," };"," "," trapindex \u003d i;"," } "," }",""," updatetrap();",""," { lakecolor() ; true } if not next(i, maxdepth) else"," radrange(znext, z, bailout, epsilon, bailoutcolor(), lakecolor())",""," }",""," escapetime(c, breakcondition)","}","","// ******************************************","// * Next are just drawing procedures. They *","// * should be the same for all drawings. * ","// ******************************************","","extern supersampling bool \u003d false;","extern light bool \u003d false;","","// drawpixel for 2D","func drawpixel_2d(x, y) { "," var c cplx \u003d map(x, y);"," var value real;"," get_color(c, value) // value is not used","}","","// drawpixel for 3D","func drawpixel_3d(x, y) {"," var c00 cplx \u003d map(x, y),"," c10 cplx \u003d map(x + 1, y + 0.5),"," c01 cplx \u003d map(x + 0.5, y + 1);"," "," var h00 real, h10 real, h01 real; // heights"," "," // color is already kinda super-sampled"," var color \u003d (get_color(c00, h00) + get_color(c10, h10) + get_color(c01, h01)) / 3;",""," // get height out of value"," func height(value) {"," extern valuetransfer expr \u003d \"value\";"," valuetransfer"," }"," "," h00 \u003d height h00; h01 \u003d height h01; h10 \u003d height h10;",""," // get the normal vector (cross product)"," var xp \u003d c10 - c00, xz \u003d h10 - h00;"," var yp \u003d c01 - c00, yz \u003d h01 - h00;"," "," var np cplx \u003d (xp.y yz - xz yp.y) : (xz yp.x - xp.x yz);"," var nz real \u003d xp.x yp.y - xp.y yp.x;"," "," // normalize np and nz"," var nlen \u003d sqrt(rad2 np + sqr nz);"," np \u003d np / nlen; nz \u003d nz / nlen;"," "," // get light direction"," extern lightvector cplx \u003d -0.667 : -0.667; // direction from which the light is coming"," def lz \u003d sqrt(1 - sqr re lightvector - sqr im lightvector); // this is inlined",""," // Lambert\u0027s law."," var cos_a real \u003d dot(lightvector, np) + lz nz;",""," // diffuse reflexion with ambient factor"," extern lightintensity real \u003d 1;"," extern ambientlight real \u003d 0.5;",""," // if lumen is negative it is behind, "," // but I tweak it a bit for the sake of the looks:"," // cos_a \u003d -1 (which is super-behind) \u003d\u003d\u003e 0"," // cos_a \u003d 0 \u003d\u003d\u003e ambientlight"," // cos_a \u003d 1 \u003d\u003d\u003e lightintensity",""," // for a mathematically correct look use the following:"," // if cos_a \u003c 0 then cos_a \u003d 0;"," // color.a \u003d color.a * (ambientlight + lightintensity lumen);"," "," def d \u003d lightintensity / 2; // will be inlined later",""," // Change L in Lab-Color"," color.a \u003d color.a (((d - ambientlight) cos_a + d) cos_a + ambientlight);",""," // Next, specular reflection. Viewer is always assumed to be in direction (0,0,1)"," extern specularintensity real \u003d 1;",""," extern shininess real \u003d 8;",""," // r \u003d 2 n l - l; v \u003d 0:0:1"," var spec_refl \u003d 2 cos_a nz - lz;"," "," // 100 because L in the Lab-Model is between 0 and 100"," if spec_refl \u003e 0 then"," color.a \u003d color.a + 100 * specularintensity * spec_refl ^ shininess;",""," color","}","","func do_pixel(x, y) {"," // two or three dimensions?"," def drawpixel \u003d drawpixel_3d if light else drawpixel_2d;"," "," func drawaapixel(x, y) {"," 0.25 ("," drawpixel(x - 0.375, y - 0.125) + "," drawpixel(x + 0.125, y - 0.375) + "," drawpixel(x + 0.375, y + 0.125) +"," drawpixel(x - 0.125, y + 0.375) "," );"," }",""," // which function to apply?"," def fn \u003d drawpixel if not supersampling else drawaapixel;",""," color \u003d lab2int fn(x, y)","}","","// and finally call the draing procedure","do_pixel(x, y)"],"arguments":{"ints":{"maxdepth":1200},"reals":{"opacity":0.25},"bools":{"light":false},"exprs":{"trapfn":"min(2box(-3:-2, 0.5:0.5, znext), 2box(2:2, -0.5:-0.5, znext))","function":"((z^4.5+z)^2)+(z-p)","traptransfer":"log(i * p)","bailouttransfer":"log(value^5 / p^-value)"},"palettes":{"bailoutpalette":{"width":9,"height":2,"colors":[-16777216,-1,-16777216,-95998,-1114368,-95998,-16777216,-1,-16777216,-16777216,-1,-16777216,-7470849,-459008,-7536385,-16777216,-1,-16777216]},"trappalette":{"width":6,"height":4,"colors":[-520726,-16776961,-16711930,-327935,-559357,-851968,-16777216,-1,-16777216,-1,-16777216,-1,-1,-16777216,-1,-16777216,-1,-16777216,-327935,-559357,-851968,-520726,-16776961,-16711930]}},"scales":{"Scale":[-4.0058129161885135E-5,6.779089759156485E-5,-6.779089759156485E-5,-4.0058129161885135E-5,-0.0402460139780857,-0.00966585373588781]}}}
Related content
Comments: 3

GLO-HE [2021-03-03 12:15:25 +0000 UTC]

👍: 0 ⏩: 1

FlyingMatthew In reply to GLO-HE [2021-03-06 13:59:03 +0000 UTC]

👍: 0 ⏩: 1

GLO-HE In reply to FlyingMatthew [2021-03-06 18:55:00 +0000 UTC]

👍: 1 ⏩: 0