Compare commits
No commits in common. "a6933029695b4173122a7dccdc6fb4f638f1c112" and "d5ad74cf77471f8855d6ba182c18898bc5555105" have entirely different histories.
a693302969
...
d5ad74cf77
@ -1,21 +0,0 @@
|
||||
# Fraktalism
|
||||
|
||||
## Iterative computing inside !
|
||||
|
||||
Voyons d'abord
|
||||
[une vidéo](http://la.buvette.org/fractales/f90/video.html)
|
||||
qui montre ma première expérience dans ce domaine.
|
||||
|
||||
## La technique
|
||||
|
||||
Le gros des calculs de fractales est fait dans XXX, et la gestion
|
||||
des pixels 'physiques' est fait dans YYY
|
||||
|
||||
Des scripts _shell_ sont utilisés pour construire les vidéos.
|
||||
|
||||
## TODO
|
||||
|
||||
- Voir de près le calcul du cadrage
|
||||
- Rajouter des formules
|
||||
- Ne pas procastiner
|
||||
|
||||
@ -1,130 +0,0 @@
|
||||
!-----------------------------------------------------
|
||||
! IMAGE PROCESSING
|
||||
!-----------------------------------------------------
|
||||
!-----------------------------------------------------
|
||||
|
||||
subroutine plotsomething(pic, start, cz)
|
||||
use cmplxmath
|
||||
use imagetools
|
||||
|
||||
implicit none
|
||||
|
||||
integer, intent(inout), dimension (:,:) :: pic
|
||||
complex, intent(in) :: start
|
||||
type (CenterMag), intent(in) :: cz
|
||||
|
||||
integer :: ix, iy, width, height
|
||||
real :: fx, fy, mod2
|
||||
complex :: za, zb, cste
|
||||
integer :: iter, maxiter
|
||||
logical :: escape
|
||||
|
||||
print *, "> plotsomething"
|
||||
|
||||
width = ubound(pic, 1)
|
||||
height = ubound(pic, 2)
|
||||
print *, " pic size ", width, height
|
||||
print *, " start ", start
|
||||
call print_centermag(cz)
|
||||
|
||||
! initialise constants
|
||||
!
|
||||
maxiter = 999;
|
||||
|
||||
! enter megaloop
|
||||
!
|
||||
do iy = 1, height
|
||||
fy = (float(iy) / float(height/3)) - 1.5
|
||||
!! print *, "line ", iy, fy
|
||||
do ix = 1, width
|
||||
fx = (float(ix) / float(width/3)) - 2.0
|
||||
|
||||
!! print *, "pixel ", ix, iy, " at ", fx, fy
|
||||
!-------------------------------------
|
||||
! working on the current pixel
|
||||
za = start
|
||||
cste = complex ( fx, fy )
|
||||
iter = 0
|
||||
escape = .FALSE.
|
||||
do while (iter .lt. maxiter)
|
||||
zb = (za * za) + cste
|
||||
! if (modulus2(zb) .gt. 4.0) then
|
||||
mod2 = real(zb)*real(zb) + aimag(zb)*aimag(zb)
|
||||
!! print *, "mod2 ", mod2
|
||||
if (mod2 .GT. 4.0) then
|
||||
escape = .TRUE.
|
||||
exit
|
||||
endif
|
||||
za = zb
|
||||
iter = iter + 1
|
||||
!! print *, "ZA ITER ESCAPE", za, iter, escape
|
||||
enddo
|
||||
|
||||
if (escape) then
|
||||
pic(ix, iy) = mod(iter, 333)
|
||||
else
|
||||
! esoteric computation here
|
||||
! pic(ix, iy) = mod(8*floor(mod2*11.11), 24)
|
||||
pic(ix, iy) = mod(iter, 222)
|
||||
endif
|
||||
!-------------------------------------
|
||||
end do ! fin boucle sur X
|
||||
|
||||
end do
|
||||
|
||||
end
|
||||
|
||||
!-----------------------------------------------------
|
||||
!
|
||||
! this is the main programm
|
||||
!
|
||||
program image
|
||||
|
||||
use imagetools
|
||||
|
||||
implicit none
|
||||
|
||||
interface
|
||||
subroutine plotsomething (pic, start, cz)
|
||||
use imagetools
|
||||
integer, intent(inout), dimension (:,:) :: pic
|
||||
complex, intent(in) :: start
|
||||
type (CenterMag), intent(in) :: cz
|
||||
end subroutine plotsomething
|
||||
end interface
|
||||
|
||||
integer, dimension(768, 768) :: picz
|
||||
type (CenterMag) :: cm
|
||||
integer :: angle
|
||||
real :: radangle, radius
|
||||
real :: stx, sty
|
||||
character (len=80) :: filename
|
||||
|
||||
cm%cx = 0.0 ; cm%cy = 0.0 ; cm%mag = 3.0
|
||||
picz = 0 ! clear screen
|
||||
|
||||
print *, "-------- making some mandelbrot -------"
|
||||
|
||||
do angle = 0, 1800
|
||||
|
||||
radangle = float(angle) * 0.017453292522222
|
||||
radius = float(angle) / 2000.0
|
||||
write (filename, "(a, i5.5, a)") "img/", angle, ".pnm"
|
||||
! filename = trim(filename)
|
||||
print *, "#### passe ", angle, radangle, trim(filename)
|
||||
|
||||
stx = radius * sin(radangle*4.0)
|
||||
sty = radius * cos(radangle*3.0)
|
||||
|
||||
call plotsomething (picz, complex(stx, sty), cm)
|
||||
call spitaspnm (picz, trim(filename))
|
||||
|
||||
print *
|
||||
|
||||
enddo
|
||||
|
||||
print *, "[DONE]"
|
||||
|
||||
end
|
||||
|
||||
!-----------------------------------------------------
|
||||
22
README.md
22
README.md
@ -1,23 +1,3 @@
|
||||
# Fortraneries
|
||||
|
||||
Le Fortran moderne, c'est quand même assez cool, alors autant diffuser
|
||||
mes gruikeries.
|
||||
|
||||
J'ai découvert le Fortran moderne lors de mon reclufinement
|
||||
de Janvier 2022, et j'ai bien aimé. Bon, contrairement à la
|
||||
version de 77, les `GOTO`s sont moins agréables à faire, mais
|
||||
l'existence des _pointeurs_ compense largement.
|
||||
|
||||
## content
|
||||
|
||||
- [SoundBrotching](SoundBrotching/) : faire gémir vos tympans
|
||||
- [BloubWorld](BloubWorld/) : la vie des particules
|
||||
- [Fraktalism](Fraktalism/) : du chaos dans les pixels
|
||||
- [RandomStuff](RandomStuff) : on a tous droit à notre jardin secret
|
||||
|
||||
## hotline
|
||||
|
||||
- Le canal `#tetalab` sur le réseau IRC de
|
||||
[Libera](https://libera.chat/)
|
||||
- La [mailing-list publique](https://lists.tetalab.org/mailman/listinfo/tetalab) du Tetalab.
|
||||
|
||||
Le fortran moderne, c'est quand même assez cool, alors autant diffuser mes gruikeries.
|
||||
15
SoundBrotching/.gitignore
vendored
15
SoundBrotching/.gitignore
vendored
@ -1,15 +0,0 @@
|
||||
#
|
||||
# please ignore those files
|
||||
#
|
||||
|
||||
*.wav
|
||||
*.text
|
||||
|
||||
text2wav
|
||||
wav2text
|
||||
text2ao
|
||||
|
||||
panoramix
|
||||
genwaves
|
||||
soundbrotch
|
||||
|
||||
@ -1,46 +0,0 @@
|
||||
#
|
||||
# tth@konrad:~/Devel/Fortraneries/SoundBrotching$
|
||||
#
|
||||
|
||||
COPT = -Wall -Wextra -g -DDEBUG_LEVEL=1
|
||||
|
||||
all: text2wav wav2text text2ao \
|
||||
panoramix genwaves
|
||||
|
||||
# ----------------------------------------------------------
|
||||
#
|
||||
# C tools
|
||||
#
|
||||
CLIBS = -lsndfile support.o
|
||||
|
||||
support.o: support.c Makefile support.h
|
||||
gcc $(COPT) -c $<
|
||||
|
||||
text2wav: text2wav.c Makefile support.h support.o
|
||||
gcc $(COPT) $< $(CLIBS) -o $@
|
||||
|
||||
wav2text: wav2text.c Makefile support.h support.o
|
||||
gcc $(COPT) $< $(CLIBS) -o $@
|
||||
|
||||
text2ao: text2ao.c Makefile support.h support.o
|
||||
gcc $(COPT) $< $(CLIBS) -o $@
|
||||
|
||||
# ----------------------------------------------------------
|
||||
#
|
||||
# real softwares
|
||||
#
|
||||
FLIBS = soundbrotch.o
|
||||
|
||||
soundbrotch.o: soundbrotch.f90 Makefile
|
||||
gfortran $(COPT) -c $<
|
||||
|
||||
panoramix: panoramix.f90 Makefile $(FLIBS)
|
||||
gfortran $(COPT) $< $(FLIBS) -o $@
|
||||
|
||||
genwaves: genwaves.f90 Makefile $(FLIBS)
|
||||
gfortran $(COPT) $< $(FLIBS) -o $@
|
||||
|
||||
# ----------------------------------------------------------
|
||||
|
||||
|
||||
# ----------------------------------------------------------
|
||||
@ -2,10 +2,3 @@
|
||||
|
||||
Stay tuned, film at 11.
|
||||
|
||||
## Cheat Code
|
||||
|
||||
Certains composants de ce sous-projet fortranique sont
|
||||
(pour le moment)
|
||||
ecrits en C, pour un accès facile à des bibliothèques tierces
|
||||
comme `libsndfile`, composant essentiel du SoundBrotching.
|
||||
|
||||
|
||||
@ -1,8 +0,0 @@
|
||||
#!/bin/bash
|
||||
|
||||
INWAV="1637.wav"
|
||||
OUTWAV="foo.wav"
|
||||
|
||||
./wav2text $INWAV | ./panoramix | ./text2wav $OUTWAV
|
||||
|
||||
|
||||
@ -1,15 +0,0 @@
|
||||
program genwaves
|
||||
|
||||
use soundbrotch
|
||||
! -----------------------------
|
||||
|
||||
implicit none
|
||||
|
||||
|
||||
print *, "genwaves is coming soon..."
|
||||
|
||||
call soundbrotch_version()
|
||||
|
||||
|
||||
end program
|
||||
|
||||
@ -1,9 +0,0 @@
|
||||
|
||||
program panoramix
|
||||
|
||||
implicit none
|
||||
|
||||
print *, "Panoramix coming soon..."
|
||||
|
||||
end program
|
||||
|
||||
@ -1,18 +0,0 @@
|
||||
module soundbrotch
|
||||
|
||||
implicit none
|
||||
|
||||
contains
|
||||
|
||||
! ---------------------------------------------------------
|
||||
|
||||
subroutine soundbrotch_version ()
|
||||
|
||||
write(0, '(A)') "*** this is soundbrotch version alpha"
|
||||
|
||||
end subroutine
|
||||
|
||||
! ---------------------------------------------------------
|
||||
|
||||
end module
|
||||
|
||||
@ -1,31 +0,0 @@
|
||||
/*
|
||||
* C SUPPORT FUNCTIONS
|
||||
*/
|
||||
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
|
||||
#include <sndfile.h>
|
||||
|
||||
#include "support.h"
|
||||
|
||||
/* --------------------------------------------------------------- */
|
||||
int display_sf_info(SF_INFO *psf, char *text)
|
||||
{
|
||||
|
||||
fprintf(stderr, " +-- sf info [%s] %p\n", text, psf);
|
||||
|
||||
fprintf(stderr, " | samplerate %d\n", psf->samplerate);
|
||||
fprintf(stderr, " | channels %d\n", psf->channels);
|
||||
fprintf(stderr, " | frames %d\n", psf->frames);
|
||||
|
||||
return 0;
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
|
||||
void print_version(char *msg)
|
||||
{
|
||||
fprintf(stderr, "======== %s === compiled %s, %s\n", \
|
||||
msg, __DATE__, __TIME__);
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
@ -1,13 +0,0 @@
|
||||
/*
|
||||
* C SUPPORT FUNCTIONS
|
||||
*/
|
||||
|
||||
#define BUFFER_SIZE 8192
|
||||
|
||||
/* --------------------------------------------------------- */
|
||||
|
||||
int display_sf_info(SF_INFO *psf, char *text);
|
||||
|
||||
void print_version(char *title);
|
||||
|
||||
/* --------------------------------------------------------- */
|
||||
@ -1,21 +0,0 @@
|
||||
/*
|
||||
* TEXT TO AUDIO OUTPUT
|
||||
*/
|
||||
|
||||
#include <stdio.h>
|
||||
|
||||
#include <sndfile.h>
|
||||
|
||||
#include "support.h"
|
||||
|
||||
/* --------------------------------------------------------------- */
|
||||
|
||||
/* --------------------------------------------------------------- */
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
|
||||
print_version(argv[0]);
|
||||
|
||||
return 0;
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
@ -1,34 +0,0 @@
|
||||
/*
|
||||
* TEXT TO WAV
|
||||
*/
|
||||
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <sndfile.h>
|
||||
|
||||
#include "support.h"
|
||||
|
||||
/* --------------------------------------------------------------- */
|
||||
int convert_text_to_wav(char *outfname)
|
||||
{
|
||||
|
||||
#if DEBUG_LEVEL
|
||||
fprintf(stderr, ">>> %s ( '%s' )\n", __func__, outfname);
|
||||
#endif
|
||||
|
||||
return -1;
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
|
||||
print_version(argv[0]);
|
||||
|
||||
if (2 != argc) {
|
||||
fprintf(stderr, "fubar\n");
|
||||
exit(1);
|
||||
}
|
||||
|
||||
return 0;
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
@ -1,81 +0,0 @@
|
||||
/*
|
||||
* WAV TO TEXT
|
||||
*/
|
||||
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <string.h>
|
||||
|
||||
#include <sndfile.h>
|
||||
|
||||
#include "support.h"
|
||||
|
||||
/* --------------------------------------------------------------- */
|
||||
int convert_wav_to_text(char *infname)
|
||||
{
|
||||
SNDFILE * sndf;
|
||||
SF_INFO sfinfo;
|
||||
int foo, lu;
|
||||
short * samples;
|
||||
|
||||
#if DEBUG_LEVEL
|
||||
fprintf(stderr, ">>> %s ( '%s' )\n", __func__, infname);
|
||||
#endif
|
||||
|
||||
memset(&sfinfo, 0, sizeof(sfinfo)); /* be clean */
|
||||
|
||||
sndf = sf_open(infname, SFM_READ, &sfinfo);
|
||||
if (sndf==NULL)
|
||||
{
|
||||
/*catch the snfile errmsg here XXX */
|
||||
fprintf(stderr, "error sf_opening %s\n", infname);
|
||||
exit(1);
|
||||
}
|
||||
|
||||
foo = display_sf_info(&sfinfo, "why not ?");
|
||||
if (foo) {
|
||||
fprintf(stderr, "%s: corrupted sf_info ?\n", __func__);
|
||||
abort();
|
||||
}
|
||||
|
||||
/* get memory for bufferins read from sound file */
|
||||
if ( NULL == (samples = malloc(BUFFER_SIZE*sizeof(short))) )
|
||||
{
|
||||
perror("\n no memory in converter");
|
||||
abort();
|
||||
}
|
||||
|
||||
while ( (lu=sf_read_short(sndf, samples, BUFFER_SIZE)) > 0 )
|
||||
{
|
||||
fprintf(stderr, " LU = %5lu\n", lu);
|
||||
|
||||
}
|
||||
/*
|
||||
* job done, some cleanup
|
||||
*/
|
||||
free(samples);
|
||||
|
||||
|
||||
return -1;
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
void usage(void)
|
||||
{
|
||||
printf("usage:\n\twav2txt 1337.wav\n");
|
||||
exit(0);
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
int foo;
|
||||
|
||||
print_version(argv[0]);
|
||||
|
||||
if (2 != argc) usage();
|
||||
|
||||
foo = convert_wav_to_text(argv[1]);
|
||||
fprintf(stderr, "got a %d from converter\n", foo);
|
||||
|
||||
return 0;
|
||||
}
|
||||
/* --------------------------------------------------------------- */
|
||||
Loading…
Reference in New Issue
Block a user