1 : C ORTH10 SOURCE CHAT 05/01/13 02:06:25 5004
2 : SUBROUTINE ORTH10 (NUMOPT,IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR
3 : & ,EE,NN,IPYY)
4 : ************************************************************************
5 : *
6 : * O R T H 1 0
7 : * -----------
8 : *
9 : * FONCTION:
10 : * ---------
11 : *
12 : * ORTHOGONALISER UN 'CHPOINT' "X" PAR RAPPORT A UNE SUITE DE
13 : * 'CHPOINTS' "U(I)" ORTHOGONAUX ENTRE EUX, L'ORTHOGONALITE ETANT
14 : * DEFINIE AU MOYEN D'UNE 'RIGIDITE' "R" PAR:
15 : * XT.R.U(I) = 0
16 : * ("T" POUR "TRANSPOSE").
17 : *
18 : * MODE D'APPEL:
19 : * -------------
20 : *
20 : * CALL ORTH10 (NUMOPT,IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR
21 : * & ,EE,NN,IPYY)
22 : *
23 : * PARAMETRES: (E)=ENTREE (S)=SORTIE
24 : * -----------
25 : *
26 : * NUMOPT ENTIER (E) NUMERO D'OPTION.
27 : * = 1 SI TOUS LES 'CHPOINTS' EN PRESENCE
28 : * S'APPUIENT SUR LES MEMES POINTS, AVEC LES
29 : * MEMES COMPOSANTES.
30 : * IPXX ENTIER (E) POINTEUR SUR LE 'CHPOINT' A ORTHOGONALISER.
31 : * IPLIS1 ENTIER (E) POINTEUR SUR LE 'LISTCHPO' CONTENANT LES
32 : * "U(I)".
33 : * IPLIS2 ENTIER (E) POINTEUR SUR LE 'LISTREEL' CONTENANT LES
34 : * PRODUITS "U(I).R.U(I)".
35 : * IPLIS3 ENTIER (E) POINTEUR SUR LE 'LISTCHPO' CONTENANT LES
36 : * PRODUITS "R.U(I)".
37 : * = 0 SI CE 'LISTCHPO' N'A PAS ETE FOURNI.
38 : * LLIST ENTIER (E) LONGUEUR DES LISTES DE POINTEURS "IPLIS1",
39 : * "IPLIS2" ET "IPLIS3".
40 : * IPRR ENTIER (E) POINTEUR SUR LA 'RIGIDITE' "R".
41 : * EE REEL DP (E) PRECISION D'ORTHOGONALITE DEMANDEE.
42 : * = 0 SI L'ON NE VEUT PAS DE TEST DE
43 : * PRECISION.
44 : * NN ENTIER (E) NOMBRE DE REPETITIONS MAXIMUM DE
45 : * L'ORTHOGONALISATION.
46 : * IPYY ENTIER (S) POINTEUR SUR LE 'CHPOINT' ORTHOGONALISE.
47 : *
48 : * LEXIQUE: (ORDRE ALPHABETIQUE)
49 : * --------
50 : *
51 : * COMBI2 SUBROUT. SOUS-PROGRAMME DE COMBINAISON LINEAIRE "RAPIDE"
52 : * (POUR DES 'CHPOINTS' SEMBLABLES).
53 : * COMB02 SUBROUT. SOUS-PROGRAMME DE COMBINAISON LINEAIRE "GENERAL"
54 : * (POUR DES 'CHPOINTS' DIFFERENTS LES UNS DES
55 : * AUTRES).
56 : *
57 : * SOUS-PROGRAMMES APPELES:
58 : * ------------------------
59 : *
60 : * ORTHO1, ORTHO2.
61 : *
62 : * AUTEUR, DATE DE CREATION:
63 : * -------------------------
64 : *
65 : * PASCAL MANIGOT 18 AVRIL 1985
66 : *
67 : * LANGAGE:
68 : * --------
69 : *
70 : * FORTRAN77
71 : *
72 : ************************************************************************
73 : *
74 : IMPLICIT INTEGER(I-N)
75 : REAL*8 EE
76 : *
77 : EXTERNAL COMBI2,COMB02
78 : *
79 : IF (IPLIS3 .EQ. 0) THEN
80 : * RQ: LES PRODUITS "RR.U(I)" NE SONT PAS FAITS.
81 : IF (NUMOPT .EQ. 1) THEN
81 : CALL ORTHO2 (IPXX,IPLIS1,IPLIS2,LLIST,IPRR,EE,NN,COMBI2
82 : & ,IPYY)
83 : ELSE
83 : CALL ORTHO2 (IPXX,IPLIS1,IPLIS2,LLIST,IPRR,EE,NN,COMB02
84 : & ,IPYY)
85 : END IF
86 : ELSE
87 : * RQ: LES PRODUITS "RR.U(I)" SONT DEJA FAITS.
88 : IF (NUMOPT .EQ. 1) THEN
88 : CALL ORTHO1 (IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR,EE,NN
89 : & ,COMBI2, IPYY)
90 : ELSE
90 : CALL ORTHO1 (IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR,EE,NN
91 : & ,COMB02, IPYY)
92 : END IF
93 : END IF
94 : *
95 : END
96 :
© Cast3M 2003 - All rights reserved.
Disclaimer