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

first public release

Browse Source
This commit is contained in:
tTh 2022-12-27 19:59:07 +01:00
parent c6f6ed48a4
commit e099b398f3
1 changed files with 13 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
Fraktalism/mkmandel.f90
View File

@ -1,8 +1,7 @@
!----------------------------------------------------- !-----------------------------------------------------
! IMAGE PROCESSING ! IMAGE PROCESSING
!----------------------------------------------------- !-----------------------------------------------------
!----------------------------------------------------- !-
subroutine plotsomething(pic, start) subroutine plotsomething(pic, start)
! use cmplxmath ! use cmplxmath
! use imagetools ! use imagetools
@ -24,13 +23,12 @@ subroutine plotsomething(pic, start)
width = ubound(pic, 1) width = ubound(pic, 1)
height = ubound(pic, 2) height = ubound(pic, 2)
print *, " pic size ", width, height ! print *, " pic size ", width, height
print *, " start ", start print *, " start ", start
! call print_centermag(cz)
! initialise constants ! initialise constants
! !
maxiter = 999; maxiter = 2500;
! enter megaloop ! enter megaloop
! !
@ -60,17 +58,16 @@ subroutine plotsomething(pic, start)
iter = iter + 1 iter = iter + 1
!! print *, "ZA ITER ESCAPE", za, iter, escape !! print *, "ZA ITER ESCAPE", za, iter, escape
enddo enddo
if (escape) then if (escape) then
pic(ix, iy)%r = mod(iter, 255) pic(ix, iy)%r = mod(iter*22, 255)
pic(ix, iy)%b = mod(iter*7, 255) pic(ix, iy)%b = mod(iter*7, 255)
else else
pic(ix, iy)%g = mod(iter, 255) pic(ix, iy)%g = mod(int(mod2 * 555), 200)
pic(ix, iy)%b = mod(iter*11, 255) ! pic(ix, iy)%g = mod(iter, 255)
! pic(ix, iy)%b = mod(iter*11, 255)
endif endif
!------------------------------------- !-------------------------------------
end do ! fin boucle sur X end do ! fin boucle sur X
end do end do
end end
@ -100,9 +97,10 @@ program mkmandel
print *, "-------- making some mandelbrot -------" print *, "-------- making some mandelbrot -------"
allocate(pic(512, 342)) allocate(pic(800, 600))
do angle = 0, 1200
do angle = 0, 600
call rgbpix_set_to_zero(pic) call rgbpix_set_to_zero(pic)
radangle = float(angle) * 0.017453292522222 radangle = float(angle) * 0.017453292522222
@ -111,8 +109,8 @@ program mkmandel
! filename = trim(filename) ! filename = trim(filename)
print *, "#### passe ", angle, radangle, trim(filename) print *, "#### passe ", angle, radangle, trim(filename)
stx = radius * sin(radangle*4.0) stx = radius * sin(radangle*3.9)
sty = radius * cos(radangle*3.0) sty = radius * cos(radangle*3.3)
call plotsomething (pic, complex(stx, sty)) call plotsomething (pic, complex(stx, sty))
call rgbpix_spit_as_pnm_8 (pic, trim(filename)) call rgbpix_spit_as_pnm_8 (pic, trim(filename))