Télécharger la plaquette Cast3M


fr en

Télécharger gmomv.eso

Retour à la liste

Numérotation des lignes : 

gmomv

  1. C GMOMV     SOURCE    PV        20/09/26    21:17:08     10724          
  2.       SUBROUTINE GMOMV(IMVEC,TRANS,ALPHA,KMORS,KISA,X,BETA,Y)
  3.       IMPLICIT REAL*8 (A-H,O-Z)
  4.       IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM         : GMOMV
  7. C DESCRIPTION : Interface à DMOMV parallèle
  8. C
  9. C
  10. C
  11. C LANGAGE     : ESOPE
  12. C AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  13. C               mél : gounand@semt2.smts.cea.fr
  14. C***********************************************************************
  15. C APPELES          :
  16. C APPELES (E/S)    :
  17. C APPELES (BLAS)   :
  18. C APPELES (CALCUL) :
  19. C APPELE PAR       :
  20. C***********************************************************************
  21. C SYNTAXE GIBIANE    :
  22. C ENTREES            :
  23. C ENTREES/SORTIES    :
  24. C SORTIES            :
  25. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  26. C***********************************************************************
  27. C VERSION    : v1, 09/02/2006, version initiale
  28. C HISTORIQUE : v1, 09/02/2006, création
  29. C HISTORIQUE :     08/01/2008, ajout petit parallélisme
  30. C HISTORIQUE :
  31. C***********************************************************************
  32. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  33. C en cas de modification de ce sous-programme afin de faciliter
  34. C la maintenance !
  35. C***********************************************************************
  36.  
  37. -INC PPARAM
  38. -INC CCOPTIO
  40.       POINTEUR KMORS.PMORS
  41.       POINTEUR KISA.IZA
  42.       POINTEUR X.IZA,Y.IZA
  43. -INC CCASSIS
  44.       CHARACTER*1        TRANS
  45.       LOGICAL            LSAME
  46.       EXTERNAL           LSAME
  47.       EXTERNAL GMOMVI
  48.       COMMON/CGMOMV/xac,xbc,kmorsc,kisac,kxc,kyc,itran,nbthr,imvecc
  49. *
  50.       INTEGER IMPR,IRET
  51.       PARAMETER (NGRAND=100000)
  52.       LOGICAL LPARA
  53. *
  54. * Executable statements
  55. *
  56.       N =KMORS.IA(/1)-1
  57. *
  58.       LPARA=.TRUE.
  59.       if (.NOT.LSAME(TRANS,'N')) THEN
  60.          itran=1
  61.       else
  62.          itran=0
  63.       ENDIF
  64.       IF (ITRAN.EQ.1) LPARA=.FALSE.
  65. *      IF (N.LT.NGRAND) LPARA=.FALSE.
  66.       IF (IMVEC.EQ.0) LPARA=.FALSE.
  67. *      WRITE(IOIMP,*) 'LPARA=',LPARA
  68. *      IF (N.GT.NGRAND) THEN
  69.       IF (LPARA) THEN
  70. * Parallélisme recopié de chole3
  71.          xac=alpha
  72.          xbc=beta
  73.          kmorsc=kmors
  74.          kisac=kisa
  75.          kxc=x
  76.          kyc=y
  77.          imvecc=imvec
  78. *
  79.          nbthro=nbthrs
  80.          if (nbthro.gt.1) call threadii
  81.          nbthr=nbthro
  82.          if (nbthr.gt.1) then
  83.             do ithr=1,nbthr
  84.                call threadid(ithr,gmomvi)
  85.             enddo
  86.             do ithr=1,nbthr
  87.                call threadif(ithr)
  88.             enddo
  89.          else
  90.             do ithr=1,nbthr
  91.                call gmomvi(ithr)
  92.             enddo
  93.          endif
  94.          if (nbthr.gt.1) call threadis
  95.       ELSE
  96.          NNZ =KMORS.JA(/1)
  97.          CALL DMOMV(TRANS,N,ALPHA,
  98.      $        NNZ,KMORS.IA,KMORS.JA,KISA.A,
  99.      $        X.A,BETA,Y.A)
  100.       ENDIF
  101. *
  102. * Normal termination
  103. *
  104.       RETURN
  105. *
  106. * End of subroutine GMOMV
  107. *
  108.       END
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  

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