TetaTricks/chap/Fortran.tex

366 lines
11 KiB
TeX
Raw Normal View History

2020-12-12 04:13:54 +01:00
\chapter{Fortran} \label{chap:Fortran} \index{Fortran}
2023-03-04 12:48:05 +01:00
\begin{quote}
2020-12-12 04:13:54 +01:00
Fortran (\textsc{FORmula TRANslator}) est, avec le Cobol\index{Cobol},
un des premiers langages évolués. Il a été conçu après la seconde
guerre mondiale, en 1954, par
John Warner Backus, ingénieur en radiophonie chez IBM\index{IBM},
2023-03-04 12:48:05 +01:00
\end{quote}
2020-12-12 04:13:54 +01:00
2023-03-21 00:16:54 +01:00
Pourquoi j'ai redécouvert un langage de ma jeunesse, et pourquoi
vous devriez vous aussi y jeter un regard attentif. Il y a plein
de bonnes choses, la septième va vous tenir compagnie.
2021-12-23 20:06:52 +01:00
% ========================================
2023-03-21 15:03:09 +01:00
%
% un petit historique ?
%
\section{On the big irons}
TODO\index{XXX}
2023-03-04 12:48:05 +01:00
\section{FORTRAN77}
2020-12-12 04:13:54 +01:00
Dernière version à exiger le format '80 colonnes' des cartes perforées,
cette version est aussi (aux yeux de vieux dinos) la plus pure,
à cause d'instructions comme le \textsl{computed goto} ou le
2022-12-16 20:41:32 +01:00
mythique \textsl{comefrom}\footnote{Keyword: Datamatic}.
Nom officiel : \texttt{ANSI X3.9-1978}.
2020-12-12 04:13:54 +01:00
2023-03-04 12:48:05 +01:00
L'incarnation \texttt{g77}\index{g77} n'existe plus, meis peut être
à-priori émulable en jouant avec les options de \texttt{gfortran}
et un brin de shell.
2022-01-01 00:04:23 +01:00
% ============****************===============
% https://linuxfr.org/users/vmagnin/journaux/taptempo-fortran
2023-06-29 11:23:11 +02:00
% https://en.wikipedia.org/wiki/Fortran_95_language_features
2021-12-23 20:06:52 +01:00
\section{Gfortan} \index{gfortran}
2020-12-12 04:13:54 +01:00
Puis ce langage du passé a subi une évolution lente vers quelque
2021-12-23 20:06:52 +01:00
chose de bien plus moderne : le Fortran 90.
2022-12-16 20:41:32 +01:00
Les dinos vont devoir s'adapter.
Et songer à utiliser un \texttt{gfortran} moderne,
2022-02-16 00:53:49 +01:00
un élément de la Gnu Compiler Collection\index{gcc}.
2021-12-23 20:06:52 +01:00
2022-02-16 00:53:49 +01:00
Comme pour tous les langages de programmation,
il est obligatoire d'avoir l'exemple canonique, donc le voilà.
2021-12-23 20:06:52 +01:00
Il faut l'enregister dans un fichier nommé \texttt{hello.f90}
2022-01-01 00:04:23 +01:00
pour que le compilateur puisse admettre que c'est écrit
en \textsl{free form format}.
2021-12-23 20:06:52 +01:00
2023-03-04 12:48:05 +01:00
\begin{lstlisting}
2022-01-01 00:04:23 +01:00
$ cat hello.f90
2021-12-23 20:06:52 +01:00
program hello
implicit none
print *, "hello world"
end
2022-01-01 00:04:23 +01:00
$ gfortran -Wall hello.f90 -o hello && ./hello
hello world
$
2023-03-04 12:48:05 +01:00
\end{lstlisting}
2021-12-23 20:06:52 +01:00
2022-01-19 23:31:51 +01:00
% --------------------------------------------------------
2022-04-03 11:56:49 +02:00
%
% https://fortranwiki.org/fortran/show/Command-line+arguments
%
2022-01-19 23:31:51 +01:00
\section{Args \& Env}
2023-03-21 00:16:54 +01:00
Un fort beau programme que ce \textsc{helloworld}, mais il est un peu
fermé au monde extérieur.
2022-01-19 23:31:51 +01:00
Quand on lance un programme, il peut recevoir des instructions
par (mais pas que) deux canaux : les arguments de la ligne de commande
et les variables d'environnement.
2023-03-04 12:48:05 +01:00
\begin{lstlisting}
2022-01-19 23:31:51 +01:00
program rum_me
implicit none
integer :: nbarg, foo
character(len=32) :: arg
nbarg = IARGC()
if (nbarg .GT. 0) then
do foo=1, nbarg
call GETARG(foo, arg)
print *, foo," ", arg
enddo
endif
end
2023-03-04 12:48:05 +01:00
\end{lstlisting}
2022-01-19 23:31:51 +01:00
2022-02-16 00:53:49 +01:00
Si un des arguments doir être vu comme une valeur numérique,
2023-03-04 12:48:05 +01:00
il faut convertir la chaine de caractères avant usage.
Exemple pour un nombre flottant~:
2022-02-16 00:53:49 +01:00
2023-03-04 12:48:05 +01:00
\begin{lstlisting}
2022-12-16 20:41:32 +01:00
character(len=89) :: arg
real :: cx
call getarg(2, string)
read (string, *) cx
2023-03-04 12:48:05 +01:00
\end{lstlisting}
2022-02-16 00:53:49 +01:00
2022-01-01 00:04:23 +01:00
% --------------------------------------------------------
2021-12-23 20:06:52 +01:00
2022-04-03 11:56:49 +02:00
\section{Tableaux}
2021-12-23 20:06:52 +01:00
2022-01-19 23:31:51 +01:00
Rank, Size, \textsc{allocatable}, toussa\dots
2022-01-01 00:04:23 +01:00
2022-04-03 11:56:49 +02:00
Amies du C, soyez d'entrée prévenues, en Fortran, l'indice
2022-12-16 20:41:32 +01:00
par défaut du premier élément d'un tableau est \textbf{1},
2022-04-03 11:56:49 +02:00
mais cette valeur peut être modifiée à la déclaration
2023-03-04 12:48:05 +01:00
du tableau. Ce qui permet des \textsl{of-by-anynumber}
du plus bel effet.
% --------------------------------------------------------
% new 29 decembre 2022
\section{Nombres complexes}
Partie imaginaire : \texttt{ipart = aimag(cmplx)}.
% --------------------------------------------------------
% nouveau 4 mars 2023
\section {Les structures}\index{type}
Une « structure », c'est comme une petite boite dans laquelle on peut
mettre des variables de différentes natures.
\begin{lstlisting}
type t_pixrgb
integer :: r, g, b
integer :: alpha = 0
end type
\end{lstlisting}
Nous verrons plus loin l'importance de cette notion.
2022-04-03 11:56:49 +02:00
2022-01-01 00:04:23 +01:00
% --------------------------------------------------------
2022-04-03 11:56:49 +02:00
\section{Fichiers}
2022-02-16 00:53:49 +01:00
Pour qui vient du C et de la philosophie Unix, la gestion
des fichiers en Fortran est assez déroutante%
\footnote{Par contre, pour un dino des \textsl{big irons}\dots}.
2022-12-31 15:39:58 +01:00
Il semble qu'il y ait deux concepts orthogonaux :
\begin{itemize}
\item Direct access \textit{vs.} Sequential access
\item Formtated \textit{vs.} Unformated
\end{itemize}
2022-02-16 00:53:49 +01:00
2022-04-01 14:23:45 +02:00
\textit{Open, Read, Write, Inquire, toussa\dots}
2022-02-16 00:53:49 +01:00
2022-12-16 20:41:32 +01:00
\subsection{Sequentiel texte}
XXX\index{XXX}
2023-03-21 00:16:54 +01:00
Pour faire simple, on accède à chaque unité de donnée l'une après l'autre,
comme quand on déroule une bande 6250bpi.
2022-12-16 20:41:32 +01:00
\subsection{Indexé} \index{index}
XXX\index{XXX}
2023-03-21 00:16:54 +01:00
C'est pas pareil.
2022-02-16 00:53:49 +01:00
% --------------------------------------------------------
2022-04-03 11:56:49 +02:00
\section{Modules}
2022-01-01 00:04:23 +01:00
Une notion essentielle, à mon avis, pour écrire des programmes
2022-01-19 23:31:51 +01:00
au-delà de 104 lignes. Oui, j'avoue, j'ai bien changé depuis
2022-02-16 00:53:49 +01:00
l'époque où je confectionnait des procédures \textsl{monobloc}
de 1337 lignes, avec une bonne centaine de \texttt{GOTO} dedans.
2021-12-23 20:06:52 +01:00
2022-12-16 20:41:32 +01:00
Il serait bon de voir un exemple du monde réel. Ou juste un
exemple basique. Mais avec des explications. Parce que ce n'est pas
vraiment évident. Mais c'est un concept balaize.
2023-03-04 12:48:05 +01:00
Nous allons donc créer un module qui ne fait quasiment rien, ce qui
limitera le nombre de pannes possibles dans la suite des opérations.
\begin{lstlisting}
module dummy
implicit none
contains
subroutine print_dummy
print *, 'this is the dummy subroutine'
end subroutine
end module
\end{lstlisting}
Rien de bien compliqué~:
Un module peut être vu comme une boite qui contient (\texttt{contains})
des procédures (\texttt{function} ou \texttt{subroutine},
et qui sera utilisé avec ce genre de petit programme~:
\begin{lstlisting}
program t
use dummy
implicit none
print *, '=== programme de test ==='
call print_dummy
end program
\end{lstlisting}
2022-04-01 14:23:45 +02:00
% --------------------------------------------------------
2022-04-03 11:56:49 +02:00
\section{Random et Aléa}
2022-04-01 14:23:45 +02:00
La dernière ligne va vous étonner.
2023-03-04 12:48:05 +01:00
\begin{lstlisting}
2022-04-01 14:23:45 +02:00
integer, dimension(3) :: tarray
integer :: t3
real :: dummy
call itime(tarray)
t3 = 3600*tarray(1) + 60*tarray(2) + tarray(3)
write(0, '(A,3I3,A,I6)') "sranding: ", tarray, " --> ", t3
call srand(t3)
! after initializing the random generator engine,
! you MUST use it for initializing the initializer
dummy = rand()
2023-03-04 12:48:05 +01:00
\end{lstlisting}
2022-01-01 00:04:23 +01:00
2022-12-16 20:41:32 +01:00
% --------------------------------------------------------
% new Thu 24 Nov 2022 02:27:05 AM CET
\section{Options de gfortran} \index{gfortran}
Comme d'habitude avec gcc\footnote{The Gnu Compiler Collection}, les
options sont innombrables. J'ai essayé d'en dégager les plus importantes,
mais ce ne sont probablement pas les seules.
2022-01-01 00:04:23 +01:00
% --------------------------------------------------------
2022-04-03 11:56:49 +02:00
\section{Images}
2022-01-01 00:04:23 +01:00
2022-06-09 21:08:46 +02:00
La bibliothèque \textbf{Plplot} semble être un bon choix pour commencer,
mais il existe bien d'autres options.
2023-03-21 00:16:54 +01:00
Par exemple \textbf{g2} (vue en page \pageref{g2}) semble être utilisable
2022-06-24 03:19:32 +02:00
avec les vieux F77, et possiblement avec un Fortran moderne, il faut juste
un peu de \textsc{rtfm}.
2023-03-21 15:03:09 +01:00
Nous allons donc commencer par PlPlot, pour envisager ensuite les alternatives.
2022-06-24 03:19:32 +02:00
% --------------------------------------------------------
2023-03-04 12:48:05 +01:00
%
% $ dillo /usr/share/doc/plplot-doc/html/index.html &
%
2022-06-09 21:08:46 +02:00
\subsection{Plplot} \index{plplot} \label{ploplt}
2022-06-24 03:19:32 +02:00
Commençons directement par un EMC\footnote{Exemple Minimal Compilable}
qui ne va pas faire grand chose, à part initialiser la machinerie
interne de Plplot qui est assez avancée.
Il y a beaucoup de paramètres réglables, mais aussi des fonctions
qui gèrent pour vous les valeurs « par défaut ».
Voyez cette exemple comme un \textsl{boilerplate} simple
ou un \textsc{HelloWorld}~:
2023-03-04 12:48:05 +01:00
\lstinputlisting{code/fortran/plplotting.f90}
2022-06-24 03:19:32 +02:00
2023-03-07 20:49:36 +01:00
Il ne se passe pas grand chose, sauf qu'on a une proposition de
choisir le type de sortie.
2022-06-24 03:19:32 +02:00
Dans la version que j'ai (XXX\index{XXX}), on a le choix entre
2023-03-07 20:49:36 +01:00
X-Window, PostScript mono ou couleur, Xfig, PNG, SVG, et bien
d'autres dont certains, pour moi, assez ésotériques.
2023-03-21 00:16:54 +01:00
D'autre part, avec cette méthode, il semble impossible de
2023-03-07 20:49:36 +01:00
préciser la taille de l'image.
2023-03-21 00:16:54 +01:00
Mais la démonstration est faite.
Il faut quand même regarder de plus près la subroutine
\texttt{plenv()} (\textsl{Set up standard window and draw box})
qui initialise simplement un contexte graphique 2D somme toute
assez classique. Elle prend six paramètres. Les quatre premiers
sont les coordonnées flottantes de deux coins du papier,
le cinquième controle l'échelle des axes, et le dernier
controle le tracé des boites autour du graphique.
Nous avons les bases, commençons donc à dessiner.
Premièrement, nous allons chosir de démarrer automatiquement
dans une fenètre X11,
en gardant le reste de la configuration avec ses valeurs par défaut.
Pour cela, il suffit de rajouter \texttt{call plsdev('xwin')} juste avant
2023-03-21 15:03:09 +01:00
l'appel à \texttt{plinit},
2023-03-21 00:16:54 +01:00
Ensuite, pour mettre en évidence notre art,
nous allons choisir une couleur pour l'encre (qui est d'un
2023-03-07 20:49:36 +01:00
rouge du meilleur effet par défaut en mode "xwin"),
puis écrire quelques légendes canoniques, donc inutiles, .
\begin{verbatim}
call plcol0 (15) ! pure white
call pllab ("Fuzzfactor", "Yoyodines", "Some nice plots")
\end{verbatim}
Il est maintenant temps de poser un premier trait de crayon numérique
sur notre feuille blanche numérique\footnote{Actuellement noire, mais
nous trouverons bien comment changer ça}.
La forme de ce tracé sera donnée par une suite de coordonnées
flottantes \textsl{x/y}
2023-03-21 00:16:54 +01:00
stockées dans deux tableaux parallèles que nous allons
immédiatement garnir de données pertinentes~:
2023-03-07 20:49:36 +01:00
\begin{verbatim}
integer, parameter :: lg = 500
real :: x(lg), y(lg)
real :: k
integer :: i
k = 0.1
do i = 1, lg
x(i) = k * sin(real(i)/3.0)
2023-03-21 00:16:54 +01:00
y(i) = k * cos(real(i)/5.0)
2023-03-07 20:49:36 +01:00
k = k * 1.006021
enddo
\end{verbatim}
2023-03-21 00:16:54 +01:00
Une fois en possession de ces données conceptuelles, que l'on peut
aisément voire comme un cheminement, nous allons demander
au crayon de se promener sur notre feuille virtuelle en suivant
2023-03-21 15:03:09 +01:00
ce chemin, et enfin de s'arreter de dessiner~:
2023-03-21 00:16:54 +01:00
\begin{verbatim}
call plcol0 (12)
call plline (x, y)
2023-03-21 15:03:09 +01:00
call plend ()
2023-03-21 00:16:54 +01:00
\end{verbatim}
Au passage, vous pouvez remarquer la magie du Fortran moderne~:
il n'est plus nécessaire de signifier la dimension des tableaux
lors de l'appel de la procédure.
2022-06-09 21:08:46 +02:00
2023-03-21 15:03:09 +01:00
Dans ce mode de fonctionnement interactif,
il vous faudra utiliser la touche
\textsl{<enter>} pour sortir\footnote{Une de mes prochaines tâche sera
de regarder comment quitter avec le '\textbf{q}' canonique}.
Nous verrons pas la suite comment écrire des fichiers 'image'
vectoriels ou bitmap en mode \textsl{batch}.
Enfin une ouverture vers les Gifs animées\index{gif89a}.
2022-06-09 21:08:46 +02:00
% --------------------------------------------------------
2023-03-07 20:49:36 +01:00
\section{Questions en attente}
2022-06-09 21:08:46 +02:00
2022-06-24 03:19:32 +02:00
\index{XXX}
2022-06-09 21:08:46 +02:00
2022-06-24 03:19:32 +02:00
\begin{itemize}
\item Comment déclarer une variable statique dans une procédure ?
\item Pourquoi ne pas parler de DISLIN ?
\item Est-il possible de causer à \texttt{libsndfile} ?
2022-12-16 20:41:32 +01:00
\item Comment caler une chaine à gauche avec un \textsc{format} ?
2023-02-09 23:41:16 +01:00
\item Is there a \texttt{-fortran} option for making comments with
2023-03-07 20:49:36 +01:00
\texttt{boxes} + \texttt{figlet}\index{figlet} ?
2022-06-24 03:19:32 +02:00
\end{itemize}
2022-01-01 00:04:23 +01:00
2023-02-09 23:41:16 +01:00
% --------------------------------------------------------
2023-02-11 12:19:00 +01:00
\section{Réponses}
2023-02-09 23:41:16 +01:00
2023-02-11 12:19:00 +01:00
\textsl{<Pas de réponse>}
2023-02-09 23:41:16 +01:00
2023-03-07 20:49:36 +01:00
Voilà, c'est tout pour le moment\dots
2023-02-09 23:41:16 +01:00
2022-01-01 00:04:23 +01:00
% ========================================