Sources Esope | Cast3M

orth10
C ORTH10    SOURCE    CHAT      05/01/13    02:06:25     5004
      SUBROUTINE ORTH10 (NUMOPT,IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR
     &                 ,EE,NN,IPYY)
************************************************************************
*
*                             O R T H 1 0
*                             -----------
*
* FONCTION:
* ---------
*
*     ORTHOGONALISER UN 'CHPOINT' "X" PAR RAPPORT A UNE SUITE DE
*     'CHPOINTS' "U(I)" ORTHOGONAUX ENTRE EUX, L'ORTHOGONALITE ETANT
*     DEFINIE AU MOYEN D'UNE 'RIGIDITE' "R" PAR:
*        XT.R.U(I) = 0
*     ("T" POUR "TRANSPOSE").
*
* MODE D'APPEL:
* -------------
*
*     CALL ORTH10 (NUMOPT,IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR
*    &           ,EE,NN,IPYY)
*
* PARAMETRES:   (E)=ENTREE   (S)=SORTIE
* -----------
*
*     NUMOPT  ENTIER    (E)  NUMERO D'OPTION.
*                            = 1  SI TOUS LES 'CHPOINTS' EN PRESENCE
*                            S'APPUIENT SUR LES MEMES POINTS, AVEC LES
*                            MEMES COMPOSANTES.
*     IPXX    ENTIER    (E)  POINTEUR SUR LE 'CHPOINT' A ORTHOGONALISER.
*     IPLIS1  ENTIER    (E)  POINTEUR SUR LE 'LISTCHPO' CONTENANT LES
*                            "U(I)".
*     IPLIS2  ENTIER    (E)  POINTEUR SUR LE 'LISTREEL' CONTENANT LES
*                            PRODUITS "U(I).R.U(I)".
*     IPLIS3  ENTIER    (E)  POINTEUR SUR LE 'LISTCHPO' CONTENANT LES
*                            PRODUITS "R.U(I)".
*                            = 0  SI CE 'LISTCHPO' N'A PAS ETE FOURNI.
*     LLIST   ENTIER    (E)  LONGUEUR DES LISTES DE POINTEURS "IPLIS1",
*                            "IPLIS2" ET "IPLIS3".
*     IPRR    ENTIER    (E)  POINTEUR SUR LA 'RIGIDITE' "R".
*     EE      REEL DP   (E)  PRECISION D'ORTHOGONALITE DEMANDEE.
*                            = 0  SI L'ON NE VEUT PAS DE TEST DE
*                            PRECISION.
*     NN      ENTIER    (E)  NOMBRE DE REPETITIONS MAXIMUM DE
*                            L'ORTHOGONALISATION.
*     IPYY    ENTIER    (S)  POINTEUR SUR LE 'CHPOINT' ORTHOGONALISE.
*
* LEXIQUE: (ORDRE ALPHABETIQUE)
* --------
*
*     COMBI2  SUBROUT.  SOUS-PROGRAMME DE COMBINAISON LINEAIRE "RAPIDE"
*                       (POUR DES 'CHPOINTS' SEMBLABLES).
*     COMB02  SUBROUT.  SOUS-PROGRAMME DE COMBINAISON LINEAIRE "GENERAL"
*                       (POUR DES 'CHPOINTS' DIFFERENTS LES UNS DES
*                       AUTRES).
*
* SOUS-PROGRAMMES APPELES:
* ------------------------
*
*     ORTHO1, ORTHO2.
*
* AUTEUR, DATE DE CREATION:
* -------------------------
*
*     PASCAL MANIGOT     18 AVRIL 1985
*
* LANGAGE:
* --------
*
*     FORTRAN77
*
************************************************************************
*
      IMPLICIT INTEGER(I-N)
      REAL*8 EE
*
      EXTERNAL COMBI2,COMB02
*
      IF (IPLIS3 .EQ. 0) THEN
*        RQ: LES PRODUITS "RR.U(I)" NE SONT PAS FAITS.
         IF (NUMOPT .EQ. 1) THEN
            CALL ORTHO2 (IPXX,IPLIS1,IPLIS2,LLIST,IPRR,EE,NN,COMBI2
     &                 ,IPYY)
         ELSE
            CALL ORTHO2 (IPXX,IPLIS1,IPLIS2,LLIST,IPRR,EE,NN,COMB02
     &                 ,IPYY)
         END IF
      ELSE
*        RQ: LES PRODUITS "RR.U(I)" SONT DEJA FAITS.
         IF (NUMOPT .EQ. 1) THEN
            CALL ORTHO1 (IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR,EE,NN
     &                 ,COMBI2,  IPYY)
         ELSE
            CALL ORTHO1 (IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR,EE,NN
     &                 ,COMB02,  IPYY)
         END IF
      END IF
*
      END