Télécharger orth10.eso

Retour à la liste

Numérotation des lignes :

  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. *
  21. * CALL ORTH10 (NUMOPT,IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR
  22. * & ,EE,NN,IPYY)
  23. *
  24. * PARAMETRES: (E)=ENTREE (S)=SORTIE
  25. * -----------
  26. *
  27. * NUMOPT ENTIER (E) NUMERO D'OPTION.
  28. * = 1 SI TOUS LES 'CHPOINTS' EN PRESENCE
  29. * S'APPUIENT SUR LES MEMES POINTS, AVEC LES
  30. * MEMES COMPOSANTES.
  31. * IPXX ENTIER (E) POINTEUR SUR LE 'CHPOINT' A ORTHOGONALISER.
  32. * IPLIS1 ENTIER (E) POINTEUR SUR LE 'LISTCHPO' CONTENANT LES
  33. * "U(I)".
  34. * IPLIS2 ENTIER (E) POINTEUR SUR LE 'LISTREEL' CONTENANT LES
  35. * PRODUITS "U(I).R.U(I)".
  36. * IPLIS3 ENTIER (E) POINTEUR SUR LE 'LISTCHPO' CONTENANT LES
  37. * PRODUITS "R.U(I)".
  38. * = 0 SI CE 'LISTCHPO' N'A PAS ETE FOURNI.
  39. * LLIST ENTIER (E) LONGUEUR DES LISTES DE POINTEURS "IPLIS1",
  40. * "IPLIS2" ET "IPLIS3".
  41. * IPRR ENTIER (E) POINTEUR SUR LA 'RIGIDITE' "R".
  42. * EE REEL DP (E) PRECISION D'ORTHOGONALITE DEMANDEE.
  43. * = 0 SI L'ON NE VEUT PAS DE TEST DE
  44. * PRECISION.
  45. * NN ENTIER (E) NOMBRE DE REPETITIONS MAXIMUM DE
  46. * L'ORTHOGONALISATION.
  47. * IPYY ENTIER (S) POINTEUR SUR LE 'CHPOINT' ORTHOGONALISE.
  48. *
  49. * LEXIQUE: (ORDRE ALPHABETIQUE)
  50. * --------
  51. *
  52. * COMBI2 SUBROUT. SOUS-PROGRAMME DE COMBINAISON LINEAIRE "RAPIDE"
  53. * (POUR DES 'CHPOINTS' SEMBLABLES).
  54. * COMB02 SUBROUT. SOUS-PROGRAMME DE COMBINAISON LINEAIRE "GENERAL"
  55. * (POUR DES 'CHPOINTS' DIFFERENTS LES UNS DES
  56. * AUTRES).
  57. *
  58. * SOUS-PROGRAMMES APPELES:
  59. * ------------------------
  60. *
  61. * ORTHO1, ORTHO2.
  62. *
  63. * AUTEUR, DATE DE CREATION:
  64. * -------------------------
  65. *
  66. * PASCAL MANIGOT 18 AVRIL 1985
  67. *
  68. * LANGAGE:
  69. * --------
  70. *
  71. * FORTRAN77
  72. *
  73. ************************************************************************
  74. *
  75. IMPLICIT INTEGER(I-N)
  76. REAL*8 EE
  77. *
  78. EXTERNAL COMBI2,COMB02
  79. *
  80. IF (IPLIS3 .EQ. 0) THEN
  81. * RQ: LES PRODUITS "RR.U(I)" NE SONT PAS FAITS.
  82. IF (NUMOPT .EQ. 1) THEN
  83. CALL ORTHO2 (IPXX,IPLIS1,IPLIS2,LLIST,IPRR,EE,NN,COMBI2
  84. & ,IPYY)
  85. ELSE
  86. CALL ORTHO2 (IPXX,IPLIS1,IPLIS2,LLIST,IPRR,EE,NN,COMB02
  87. & ,IPYY)
  88. END IF
  89. ELSE
  90. * RQ: LES PRODUITS "RR.U(I)" SONT DEJA FAITS.
  91. IF (NUMOPT .EQ. 1) THEN
  92. CALL ORTHO1 (IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR,EE,NN
  93. & ,COMBI2, IPYY)
  94. ELSE
  95. CALL ORTHO1 (IPXX,IPLIS1,IPLIS2,IPLIS3,LLIST,IPRR,EE,NN
  96. & ,COMB02, IPYY)
  97. END IF
  98. END IF
  99. *
  100. END
  101.  
  102.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales