Télécharger ordm01.eso

Retour à la liste

Numérotation des lignes :

  1. C ORDM01 SOURCE JC220346 14/12/16 21:15:03 8324
  2. ************************************************************************
  3. * NOM : ORDM01
  4. * DESCRIPTION : Tri par fusion (Merge Sort) pour une liste de reels
  5. ************************************************************************
  6. * APPELE PAR : pilot.eso
  7. ************************************************************************
  8. * ENTREES :: A = liste de reels a trier
  9. * N = longueur de la liste A
  10. * T = tableau de travail de reels de longueur (N+1)/2
  11. * CROISS = booleen indiquant l'ordre du tri
  12. * SORTIES :: A = liste triee
  13. ************************************************************************
  14. RECURSIVE SUBROUTINE ORDM01(A,N,T,CROISS)
  15. *
  16. INTEGER I,J,K,N,NA,NB
  17. REAL*8 A,T,V
  18. DIMENSION A(N),T((N+1)/2)
  19. LOGICAL CROISS,DECROI
  20. *
  21. * On arrete la recursion pour N=1 ou N=2
  22. IF (N.LT.2) RETURN
  23. DECROI = .NOT.CROISS
  24. IF (N.EQ.2) THEN
  25. IF ((CROISS.AND.A(1).GT.A(2)).OR.
  26. & (DECROI.AND.A(1).LT.A(2))) THEN
  27. V = A(1)
  28. A(1) = A(2)
  29. A(2) = V
  30. ENDIF
  31. RETURN
  32. ENDIF
  33. *
  34. * Sinon, on subdivise en deux sous-listes plus petites que l'on
  35. * va ordonner separement
  36. NA=(N+1)/2
  37. NB=N-NA
  38. CALL ORDM01(A,NA,T,CROISS)
  39. CALL ORDM01(A(NA+1),NB,T,CROISS)
  40. *
  41. * On melange les deux petites sous-listes en une grande liste
  42. * ordonnee (avec un peu de chance, il n'y a rien a faire...)
  43. IF ((CROISS.AND.A(NA).GT.A(NA+1)).OR.
  44. & (DECROI.AND.A(NA).LT.A(NA+1))) THEN
  45. T(1:NA)=A(1:NA)
  46. I = 1
  47. J = 1
  48. K = 1
  49. DO WHILE (I.LE.NA.AND.J.LE.NB)
  50. V = A(NA+J)
  51. IF ((CROISS.AND.T(I).LE.V).OR.
  52. & (DECROI.AND.T(I).GE.V)) THEN
  53. A(K) = T(I)
  54. I = I+1
  55. ELSE
  56. A(K) = V
  57. J = J+1
  58. ENDIF
  59. K = K + 1
  60. ENDDO
  61. *
  62. DO WHILE (I.LE.NA)
  63. A(K) = T(I)
  64. I = I + 1
  65. K = K + 1
  66. ENDDO
  67. ENDIF
  68. *
  69. RETURN
  70. *
  71. END
  72. *
  73.  

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