112 —
Dans l'équation de condition (1) ce terme est égal simplement à
. . . + d k . S~F . .. ,
le coefficient correspondant est donc immédiatement calculable.
L'obtention du tableau des coefficients CI est une phase essentielle du travail : ce ta
bleau est, compte-tenu de la colonne de termes tous connus (2ème nombre) et pour chaque ban
de, une matrice de : - (44 x 168) = 7.392 nombres.
3) - Normalisation et réduction.
Une des idées principales du programme consiste à faire apparaître, au niveau de la
jonction générale, une matrice normale déjà réduite de toutes les lignes relatives aux conditions
rhomboèdriques. Ceci implique que la construction des coefficients rectangles "CCI", qui
relient celles-ci (symbole C) aux conditions sur coordonnées (de symbole CI), soit aussitôt sui
vie de la réduction de l'effet de ces termes. En un mot, nous cherchons à obtenir, pour chaque
bande, un tableau des "CCI" réduits (y compris les termes de 2ème membre) après prise en
compte des lignes correspondant aux C. Cette double opération de construction et réduction des
"C C I" est facilitée du fait qu'elle est, tout comme la construction des C I eux-mêmes, une
énumération selon les nadirs; en outre, du fait que la matrice inverse du tableau des C est par
ticulièrement simple, les rhomboèdres étant "orthogonaux" entre eux.
Cette réduction préalable de la partie du système normal due à chaque bande est essen
tielle pour permettre l'implantation en mémoire, dans les limites de volume prescrites par le
matériel.
4) - Résolution du système normal.
Pour un système symétrique d'ordre 60, le problème de la résolution n'est pas consi
dérable avec les moyens modernes. Ont été utilisés avec succès un programme dû à Monsieur
Jacquard, qui met en oeuvre la méthode des gradients conjugués; également, un programme de
résolution classique par élimination, dû à Monsieur Lasvergeres.
5) - Programmation.
A l'exception des sous-programmes ci-dessus, la programmation est entièrement
nouvelle, réalisée par Monsieur Vosluisant en Fortran II, puis convertie en Fortran IV pour
passage en I B M 7040. Elle se subdivise dans les parties principales suivantes :
a) programme maître ou de coordination générale;
b) programme TRAIT, de traitement de chaque bande englobant la mise en place sur les
deux points principaux;
c) programme M A C I : construction et réduction préalable des tableaux de coefficients CI
et du tableau normal symétrique;
d) sortie des données de bande et des résultats intermédiaires de bande;
e) résolution du système symétrique de jonction;
f) rappel des données et résultats intermédiaires de bandes;
g) programme d'exploitation pour arriver aux corrections d a 2 : "DELDAL";
h) sortie des résultats finaux.
Environ 950 instructions FORTRAN, au total.
Pour un bloc de 3 bandes de 15 photographies, le passage demande environ 30 secondes
sur IBM 7094, 3 minutes sur IBM 7040.