Fraktalism/README.md

@@ -9,8 +9,8 @@ qui montre ma première expérience dans ce domaine.
 ## La technique
 
 Le gros des calculs de fractales est fait dans `mods/fraktals.f90`,
-et la gestion des pixels 'physiques' est faite par les
+et la gestion des pixels 'physiques' est faite par le
-modules externes `spitpgm` et `pixrgb`.
+module externe `spitpgm`.
 
 Les fonctions d'usage général sont dans
 [mods/](répertoire mods/) ave trop peu
@@ -38,6 +38,3 @@ Generally writen as a *sequencial unformated* file.
 - Rajouter des formules
 - Ne pas procastiner sur le reste
 
-## See also
-
-- https://www.maths.town/fractal-articles/mandelbulb/mandelbulb-all-powers/

Fraktalism/fraktals.f90

@@ -153,4 +153,5 @@ function modulus2(pt)
   modulus2 = real(pt)*real(pt) + imag(pt)*imag(pt)
 end
 !-----------------------------------------------------
+
 end module fraktals

Fraktalism/henon.f90

@@ -8,7 +8,7 @@ program henon
   integer                      :: w, h
   integer                      :: foo, bar
   double precision             :: px, py
-  w = 2000   ;   h = 1600
+  w = 2000   ;   h=1600
 
   write(0, *) "###### Mapping of Henon "
 
@@ -25,6 +25,7 @@ program henon
 !-----------------------------------------------------
 contains
 !-----------------------------------------------------
+
 !-----------------------------------------------------
 subroutine compute_pixel_henon(a, b, maxpasse, passe, limit, rx, ry)
   implicit none

Fraktalism/julia.f90

@@ -0,0 +1,33 @@
+!-----------------------------------------------------
+!           JULIA
+!           =====
+!    this is the main program
+!-----------------------------------------------------
+
+program julia
+
+    use spitpgm
+    use fraktals
+
+    implicit none
+
+    integer, dimension(512, 342)    ::    picz
+    integer                         ::    argc
+    character(200)                  ::    filename, string
+    real                            ::    cx, cy
+
+    argc = IARGC()
+    if (3 .NE. argc) then
+      STOP ": JULIA PROGGY NEED PARAMETERS !"
+    endif
+
+    call getarg(1, filename)
+    call getarg(2, string)       ;   read (string, *) cx
+    call getarg(3, string)       ;   read (string, *) cy
+
+    call simple_julia(picz, cx, cy, 2500)
+    call spit_as_pgm_8(picz, trim(filename))
+
+end program
+
+!-----------------------------------------------------

Fraktalism/lorentz.f90

@@ -0,0 +1,31 @@
+! 
+!    this is the main programm
+! 
+!-----------------------------------------------------
+
+program lorentz
+
+    use spitpgm
+    use fraktals
+
+    implicit none
+
+    integer, dimension(800, 600)    ::    picz
+    integer                         ::    argc
+    character(200)                  ::    filename
+
+    argc = IARGC()
+    if (1 .NE. argc) then
+      STOP ": LORENTZ NEED A FILENAME !"
+    endif
+
+    call getarg(1, filename)
+
+    write (0, "(A)") "Lorentz -> "//trim(filename)
+
+    call lorentz_0(picz, 50000)
+    call spit_as_pgm_8(picz, trim(filename))
+
+end program
+
+!-----------------------------------------------------

Fraktalism/mods/xperiment.f90

@@ -8,6 +8,8 @@ module xperiment
 
 subroutine parasites_1(pic, cx, cy, maxiter)
     implicit none
+
+    ! here is the wtf
     integer, intent(inout), dimension (:,:)  ::  pic
 
     real, intent(in)                         ::  cx, cy
@@ -21,18 +23,22 @@ subroutine parasites_1(pic, cx, cy, maxiter)
  
     width  = ubound(pic, 1)    ;    height = ubound(pic, 2)
     coef   = float(maxiter)
+
     do ix = 1, width
       fx = cx + (float(ix) / (float(width)/4.0) - 2.0)
       burps = (RAND() .lt. 0.01)
       do iy = 1, height
         fy = cy + (float(iy) / (float(height)/4.0) - 2.0)
+
         if (burps) then
           pic(ix, iy) = int(fx * fy * coef * 1.005)
         else
           pic(ix, iy) = int(fx * fy * coef)
         endif
+
       enddo
     enddo
+
 end subroutine parasites_1
 
 !---------------------------------------------------------------
@@ -46,9 +52,13 @@ subroutine loop_of_parasites_1(nbre, mode)
     integer                   :: idx
 
     if (mode .NE. 0)  STOP "BAD MODE"
+
     do idx = 0, nbre
+
       write(0, *) "popcorn loop ", idx
+
     enddo
+
 end subroutine loop_of_parasites_1
 
 !===============================================================

GrafAnim/README.md

@@ -13,8 +13,6 @@ un peu foireux sur les tracés de ligne...
 
 Distorsions approximatives de la courbe de Lissajous.
 
-Expériences inspirées par  https://bleuje.com/tutorial1/ que c'est d'la balle !
-
 ## doubledice
 
 Ou comment dessiner des gaussiennes.

GrafAnim/doubledice.f90

@@ -1,61 +1,9 @@
 program doubledice
   use usegenplot
-  use utils_ga
-
   implicit none
 
-  integer                   :: nbarg, numframe
+  call init_genplot("essai.genplot")
-  character(len=256)        :: arg
-  integer                   :: idx, foo, bar, xpos
-  integer                   :: buckets(12)
 
-  nbarg = IARGC()
+  call end_genplot("OK boomer")
-  if (nbarg .GT. 0) then
-       call GETARG(1, arg)
-       ! write (0, '(A40, A5)') "argument = ",  arg
-       read (arg, *) numframe
-  endif
-
-  write(0, '(" ------------- DOUBLE DICE ----------")')
-  write(0, '(" --> frame number:", I5)') numframe
-  call init_genplot("dummy.genplot")
-  call gplt_rect(0, 0, 800, 600)
-  call gplt_setcol(3)
-  call gplt_rect(8, 8, 48, 38)
-  call gplt_setcol(2)
-
-  buckets = 0
-  do idx=1, numframe
-    call jouer_un_tour(buckets, idx)
-    foo = maxval(buckets)
-    ! write(0, *) "max value", foo
-  enddo
-
-  call gplt_line(5, 40, 795, 40)
-  do bar=1, 12
-    xpos = bar*50 + 25
-    ! write(0, *) bar, xpos, buckets(bar)
-    call gplt_rect(xpos, 41, xpos+50, 41+15*buckets(bar))
-  enddo
-
-  call end_genplot("OK boomer...")
-
-contains
-
-! ---------------------------------------------------------
-
-subroutine jouer_un_tour(table, compte)
-  integer,intent(inout)   :: table(12)
-  integer, intent(in)     :: compte
-
-  integer                 :: DA, DB, valeur
-
-  DA = fair_random_dice()
-  DB = fair_random_dice()
-  valeur = DA + DB
-  ! write(0, *) "pass ", compte, "  =  ", DA, DB, valeur
-  table(valeur) = table(valeur) + 1
-
-end subroutine
 
 end program

GrafAnim/doublegauss.f90

@@ -1,33 +0,0 @@
-program doublegauss
-  use pixrgb
-  use utils_ga
-  implicit none
-
-  type(t_pixrgb), allocatable     ::    pic(:,:)
-  character (len=80)              ::    filename
-  integer                         ::    pass, iter
-  integer                         ::    xrnd, yrnd
-
-  write(0, *) "----- making a doublegauss picture ----"
-
-  allocate(pic(320, 240))
-  call rgbpix_set_to_zero(pic)
-
-  do pass=0, 99
-    do iter=1, 15000
-      xrnd = fair_random_gauss(320)
-      yrnd = fair_random_gauss(240)
-      ! print *, xrnd, yrnd
-      pic(xrnd,yrnd)%r = pic(xrnd,yrnd)%r + 1
-      pic(xrnd,yrnd)%g = pic(xrnd,yrnd)%g + 2
-      pic(xrnd,yrnd)%b = pic(xrnd,yrnd)%b + 3
-    end do
-
-    write (filename, "(a, i5.5, a)") "F/DBG/", pass, ".pnm"
-    print *, trim(filename)
-    call rgbpix_spit_as_pnm_8 (pic, trim(filename))
-
-  end do
-
-end program
-

GrafAnim/runme.sh

@@ -1,19 +0,0 @@
-#!/bin/bash
-
-set -e
-
-SCRATCH="dummy.genplot"
-
-make doubledice
-
-rm -f F/????.tga
-
-for foo in $(seq 1 250)
-do
-	rm -f dummy.genplot
-	./doubledice $foo 
-	filename=$(printf "F/%04d.tga" $foo)
-	genplot2 -s 640x480 ${SCRATCH} $filename
-done
-
-convert -delay 10 F/????.tga foo.gif

GrafAnim/usegenplot.f90

@@ -3,40 +3,26 @@ module usegenplot
   implicit none
 
     integer, private           :: color = 4
-    integer, private           :: iochannel = -1
     logical, private           :: initialised = .FALSE.
-
   contains
 
 ! -------------------------------------------------------------------
 
 subroutine init_genplot(filename)
+
   character(*), intent(in)    :: filename
-
-  integer                     :: errcode
-
   write(0, *)   '--> init genplot "', filename, '"'
-  open(newunit=iochannel, file=filename, iostat=errcode)
-  write(0, *)   'iochannel', iochannel, 'errcode', errcode
-
   initialised = .TRUE.
   color = 4
 
-  ! XXX STOP 'ABEND'
-
 end subroutine
 
 subroutine end_genplot(message)
+
   character(*), intent(in)    :: message
-
-  integer                     :: errcode
-
   write(0, *)   '--> end genplot "', message, '"'
   initialised = .FALSE.
 
-  close(unit=iochannel, iostat=errcode)
-  write(0, *)   'close errcode', errcode
-
 end subroutine
 
 ! -------------------------------------------------------------------
@@ -52,7 +38,7 @@ subroutine do_initialise_once(motif)
 end subroutine
 
 ! -------------------------------------------------------------------
-!-    getter, setter, wtf ?
+
 subroutine gplt_setcol(col)
   integer, intent(in)           :: col
   color = col
@@ -69,7 +55,8 @@ subroutine gplt_move(x, y)
   if (.NOT. initialised) then
     call do_initialise_once('in gplt_move')
   endif
-  write(iochannel, '(3I8)') x, y, 0
+
+  print *, x, y, 0
 end subroutine
 
 ! -------------------------------------------------------------------
@@ -79,7 +66,7 @@ subroutine gplt_draw(x, y)
   if (.NOT. initialised) then
     call do_initialise_once('in gplt_draw')
   endif
-  write(iochannel, '(3I8)') x, y, color
+  print *, x, y, color
 end subroutine
 
 ! -------------------------------------------------------------------
@@ -91,23 +78,10 @@ subroutine gplt_line(x1, y1, x2, y2)
 end subroutine
 
 ! -------------------------------------------------------------------
-!-
-!                           sx, sy
-!        +-------------------+
-!        |             x2,y2 |
-!        |                   |
-!        |                   |
-!        | x1,y1             |
-!        +-------------------+
-!       0,0
 
 subroutine gplt_rect(x1, y1, x2, y2)
    integer, intent(in)           :: x1, y1, x2, y2
 
-  if (.NOT. initialised) then
-    call do_initialise_once('in gplt_rect')
-  endif
-
    call gplt_move(x1, y1)
    call gplt_draw(x2, y1)
    call gplt_draw(x2, y2)

GrafAnim/utils_ga.f90

@@ -1,41 +0,0 @@
-! -------------------------------------------------------------------
-!-      fonctions diverses pour faire des images bizarres
-! -------------------------------------------------------------------
-
-module utils_ga
-  use pixrgb
-  implicit none
-
-  contains
-
-! -------------------------------------------------------------------
-function fair_random_dice()
-  integer            :: fair_random_dice
-
-  fair_random_dice = 1 + int(rand()*6.0)
-
-end function
-! -------------------------------------------------------------------
-!     usage  -->  see doublegauss.f90
-function fair_random_gauss(hilevel)
-  integer, intent(in) :: hilevel
-  integer             :: fair_random_gauss
-  integer             :: foo, bar
-
-  foo = int(rand()*hilevel/2)
-  bar = int(rand()*hilevel/2)
-  fair_random_gauss = 1 + foo + bar
-
-end function
-! -------------------------------------------------------------------
-!     usage  -->  see doublegauss.f90
-subroutine increment_pixel(pix, k)
-  type(t_pixrgb), intent(inout)  :: pix
-  integer                        :: k
-
-
-
-end subroutine
-! -------------------------------------------------------------------
-
-end module

GrafAnim/vue3axes.f90

@@ -1,12 +0,0 @@
-module vue3axes
-  use pixrgb
-  implicit none
-!-------------------------------------------------------------
-
-
-!-------------------------------------------------------------
-!-
-!-	ouf, c'est fini
-!-
-end module
-

README.md

@@ -16,10 +16,6 @@ l'existence des _pointeurs_ compense largement.
 - [RandomStuff](RandomStuff/) : on a tous droit à notre jardin secret
 - [GrafAnim](GrafAnim/) : Ah, enfin de la gif89a en vue !
 
-## Prérequis
-
-- GNUtrucs : bash, make, awk... 
-
 ## hotline
 
 - Le canal `#tetalab` sur le réseau IRC de

SoundBrotching/c-tools/text2wav.c

@@ -65,12 +65,6 @@ while (2==fscanf(input, "%d %d", &left, &right)) {
 	buffer[idx++] = left;
 	buffer[idx++] = right;
 
-	/* SANITY CONTROL */
-	if ( (abs(left) > 32767) || (abs(right) > 32767) ) {
-		fprintf(stderr, "OVERFLOW sample %9d   values %7d %7d\n",
-						idx, left, right);
-		}
-
 	if (idx >= SMPL_COUNT) {
 		/* flush buffer to file */
 		sf_write_short(sndf, buffer, idx);

SoundBrotching/soundbrotch.f90

@@ -6,10 +6,6 @@ module soundbrotch
     integer             :: left
     integer             :: right
   end type
-  type t_sample2r
-    real                :: left
-    real                :: right
-  end type
   ! ---------------------------------------------------------
   !      some private variables