tTh/Fortraneries
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

finetuning the mandelbrot

Browse Source
This commit is contained in:
tTh 2023-05-03 02:46:52 +02:00
parent 89d1cbda85
commit 6066dee701
2 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Fraktalism/mandelbrots.f90
View File

@ -28,7 +28,7 @@ subroutine mandelbrot_one(pic, start)
integer :: iter, maxiter integer :: iter, maxiter
logical :: escape logical :: escape
write(0,*) "> plotsomething" ! write(0,*) "> plotsomething"
width = ubound(pic, 1) width = ubound(pic, 1)
height = ubound(pic, 2) height = ubound(pic, 2)
@ -37,7 +37,7 @@ subroutine mandelbrot_one(pic, start)
! initialise constants ! initialise constants
! !
maxiter = 1984 maxiter = 3456
! enter megaloop ! enter megaloop
! !

10
Fraktalism/mkmandel.f90
View File

@ -23,9 +23,9 @@ program mkmandel
write(0, *) "-------- making some mandelbrot -------" write(0, *) "-------- making some mandelbrot -------"
allocate(pic(1024, 768)) allocate(pic(640, 480))
do angle = 0, 1500 do angle = 0, 2000
call rgbpix_set_to_zero(pic) call rgbpix_set_to_zero(pic)
@ -33,10 +33,10 @@ program mkmandel
radius = float(angle) / 2000.0 radius = float(angle) / 2000.0
write (filename, "(a, i5.5, a)") "frames/mandel/", angle, ".pnm" write (filename, "(a, i5.5, a)") "frames/mandel/", angle, ".pnm"
! filename = trim(filename) ! filename = trim(filename)
write(0,*) "#### passe ", angle, radangle, trim(filename) write(0,*) "passe ", angle, radangle, trim(filename)
stx = radius * (sin(radangle*2.07) + 0.2131*sin(radangle*7.36)) stx = radius * (sin(radangle*4.07) + 0.2131*sin(radangle*11.36))
sty = radius * (cos(radangle*3.21) + 0.2725*cos(radangle*9.99)) sty = radius * (cos(radangle*6.21) + 0.2725*cos(radangle*10.99))
call mandelbrot_one(pic, complex(stx, sty)) call mandelbrot_one(pic, complex(stx, sty))
call rgbpix_spit_as_pnm_8 (pic, trim(filename)) call rgbpix_spit_as_pnm_8 (pic, trim(filename))