Trustworthy Julia set source code

  • 5 Replies
  • 133 Views

0 Members and 1 Guest are viewing this topic.

Offline marcm200

  • *
  • Fractal Furball
  • ***
  • Posts: 294
« on: July 12, 2019, 12:00:03 PM »
I uploaded a first version of the source code for my implementation of the "trustworthy Julia sets algorithm" by Figueiredo et al.

https://github.com/marcm200/julia-tsa-core

The functionality is pretty simple for now, but images up to 2^18 pixel width have been generated by it:
  • command line version
  • generating a Julia set from scratch for the quadratic, cubic, quartic or pentic case.
  • resuming computation from temporary raw data
  • refinement simulation

Offline gerson

  • *
  • Fractal Feline
  • **
  • Posts: 167
« Reply #1 on: July 12, 2019, 06:31:02 PM »
Could you release the .exe file?

Offline marcm200

  • *
  • Fractal Furball
  • ***
  • Posts: 294
« Reply #2 on: July 12, 2019, 09:21:44 PM »
Yes, I will in a couple of days. I'm trying to understand the licensing of the tdm64 compiler I'm using and if I can just compile and upload the exe of my code or if there's more legal things to do beforehand.

Offline marcm200

  • *
  • Fractal Furball
  • ***
  • Posts: 294
« Reply #3 on: July 14, 2019, 11:48:18 AM »
I uploaded 64bit (tested on Windows 10 and Wine for ubuntu 12.04) and 32bit (tested on Windows 10, 8.1, XP and Wine for ubuntu 12.04) executables. Also included are some DOS-batch files to demonstrate automatic computation of increasingly large images.


Offline gerson

  • *
  • Fractal Feline
  • **
  • Posts: 167
« Reply #4 on: July 15, 2019, 04:52:18 PM »
I edited BATCH CREATION file to render an image, see the result:

@echo off
juliatsacore_quadratic_win32.exe 4096 -0.78 -0.05 2 6
copy _tsa_juliaset_Y00X00.bmp julia-tsa-core_19071501.bmp
pause

Offline marcm200

  • *
  • Fractal Furball
  • ***
  • Posts: 294
« Reply #5 on: July 15, 2019, 05:48:51 PM »
Nice, and a relief for me that the binaries actually work outside my house.

The underlying double floating type in my code however cannot accurately represent your value -0.78 (To my knowledge, decimal fractions can only be represented if they end in 5 or zero (mandatory condition but not sufficient alone).

So the algorithm actually computes a slightly different set (whatever the compiler makes out of it), I doubt it makes any visible difference, but to still be on the side of mathematical guarantee, I usually correct any value I like by:

Code: [Select]
floor or ceil(DECIMAL_SEED_VALUE * 2^25) / 2^25


clip
Minuscule source code for Lyapunov images (ASCII style)

Started by marcm200 on Programming

8 Replies
356 Views
Last post March 20, 2019, 07:37:58 PM
by marcm200
question
GitHub source deployment

Started by ygidaslm on Kalles Fraktaler

5 Replies
242 Views
Last post May 20, 2018, 07:21:53 PM
by ygidaslm
xx
Code: (Select) not working?

Started by Anon on Forum Help And Support

3 Replies
154 Views
Last post February 19, 2018, 03:02:31 PM
by 3DickUlus
xx
Code Snippets for Fragmentarium

Started by Sabine62 on Code Snippets (fragments)

0 Replies
105 Views
Last post November 02, 2018, 01:23:11 PM
by Sabine62
xx
Octave/MATLAB Mandelbrot code

Started by gerrit on Programming

11 Replies
944 Views
Last post October 19, 2017, 05:45:09 AM
by gerrit