Télécharger la plaquette Cast3M


fr en

Télécharger hbmordo.eso

Retour à la liste

Numérotation des lignes : 

hbmordo

  1. C HBMORDO   SOURCE    OF166741  26/05/11    21:15:13     12538          
  2.  
  3. *=======================================================================
  4. *  tri des valeurs propres par partie imaginaire croisante en module
  5. *=======================================================================
  6.  
  7.       SUBROUTINE HBMORDO(NM,NDDL,OMEG,EXPIM,INDEIG)
  8.  
  9.       IMPLICIT INTEGER(I-N)
  10.       IMPLICIT REAL*8(A-H,O-Z)
  11.  
  12.       INTEGER NM,NDDL
  13. C-      DIMENSION INDEIG(2*NM)
  14.       INTEGER INDEIG(*)
  15.       REAL*8 OMEG
  16. C-      DIMENSION*8 EXPIM(NM)
  17.       REAL*8 EXPIM(*)
  18.  
  19.       INTEGER I
  20. *      INTEGER MINT
  21. *      REAL*8 EXPIW(NM) <- Alloc. dynamique ! MLREEL ?
  22.  
  23. c      DO I=1,2*NDDL
  24. c        INDEIG(I)=0
  25. c      ENDDO
  26.       DO I=1,NM
  27.         INDEIG(I)=I
  28.       ENDDO
  29.  
  30. *c  Copier EXPIM dans le tableau de travail EXPIW
  31. *      DO I = 1,NM
  32. *        EXPIW(I) = ABS(EXPIM(I))
  33. *      ENDDO
  34. *c  Pour l'instant, on prend simplement les 2*NDDL VP avec les parties
  35. *c  imaginaires les plus petites en module
  36. *      DO I=1,2*NDDL
  37. *        MINT = MINLOC(EXPIW,DIM=1)
  38. *        INDEIG(I) = MINT
  39. *        EXPIW(MINT) = 1.E6
  40. *      ENDDO
  41.  
  42. c   bp: on remplace ci-dessus par le classement de la liste
  43. c       on recuperera seulement les 1eres valeurs
  44.       CALL ORDO13(EXPIM,NM,.TRUE.,INDEIG)
  45.  
  46.       RETURN
  47.       END
  48.  
  49.  
  50.  

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