Télécharger la plaquette Cast3M


fr en

Télécharger kres15.eso

Retour à la liste

Numérotation des lignes : 

kres15

  1. C KRES15    SOURCE    GOUNAND   25/04/30    21:15:12     12258          
  2.       SUBROUTINE KRES15(KMORS,KIZA)
  3.       IMPLICIT REAL*8 (A-H,O-Z)
  4.       IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM         : KRES15
  7. C DESCRIPTION : - Ordonnancement des colonnes d'une matrice Morse
  8. C
  9. C
  10. C LANGAGE     : ESOPE
  11. C AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  12. C               mél : gounand@semt2.smts.cea.fr
  13. C***********************************************************************
  14. C VERSION    : v1, 10/04/2025, version initiale
  15. C HISTORIQUE : v1, 10/04/2025, création
  16. C HISTORIQUE :
  17. C HISTORIQUE :
  18. C***********************************************************************
  19.       REAL*8 XKT,PREC
  20. -INC SMRIGID
  21. -INC SMMATRI
  22.  
  23. -INC PPARAM
  24. -INC CCOPTIO
  25. -INC SMLENTI
  26.       POINTEUR IWORK.MLENTI
  27.       SEGMENT PMORS
  28.       INTEGER IA (NTT+1)
  29.       INTEGER JA (NJA)
  30.       ENDSEGMENT
  31.       SEGMENT IZA
  32.       REAL*8  A(NBVA)
  33.       ENDSEGMENT
  34.       POINTEUR KMORU.PMORS,KMORUT.PMORS,KMORL.PMORS,KMORS.PMORS
  35.       POINTEUR KIZAU.IZA,KIZAUT.IZA,KIZAL.IZA,KIZA.IZA
  36. *
  37. *      WRITE(IOIMP,*) 'COUCOU KRES15'
  38. *
  39. * Et on ordonne les colonnes si besoin
  40. *  Les colonnes sont-elles ordonnées ?
  41. * Réponse : pas toujours !
  42.       PMORS=KMORS
  43.       IZA=KIZA
  44.       SEGACT PMORS
  45.       SEGACT IZA
  46.       NTT=IA(/1)-1
  47.       DO ITT=1,NTT
  48.          IAD=IA(ITT)
  49.          IAF=IA(ITT+1)-1
  50.          JINI=0
  51.          DO IJA=IAD,IAF
  52.             JCOU=JA(IJA)
  53.             IF (JCOU.LT.JINI) GOTO 30
  54.             JINI=JCOU
  55.          ENDDO
  56.       ENDDO
  57. *      WRITE(IOIMP,*) 'Les colonnes sont ordonnees'
  58.       GOTO 40
  59. *
  60.  30   CONTINUE
  61. *      WRITE(IOIMP,*) 'Les colonnes ne sont pas ordonnees'
  62.       SEGACT PMORS*MOD
  63.       SEGACT IZA*MOD
  64.       NTT=IA(/1)-1
  65.       NJA=JA(/1)
  66.       JG=MAX(NTT+1,2*NJA)
  67.       SEGINI IWORK
  68.       CALL CSORT(IA(/1)-1,A,JA,IA,
  69.      $     IWORK.LECT,.TRUE.)
  70.       SEGSUP IWORK
  71. C      WRITE(IOIMP,*) 'Matrice Morse  ordonnée'
  72. C      CALL ECMORS(PMORS,IZA,4)
  73.  40   CONTINUE
  74. C
  75.       RETURN
  76.       END
  77.  
  78.  

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