Télécharger la plaquette Cast3M


fr en

Télécharger psolve.eso

Retour à la liste

Numérotation des lignes : 

psolve

  1. C PSOLVE    SOURCE    PV        20/09/26    21:19:38     10724          
  2.       SUBROUTINE PSOLVE(TRANS,KPREC,KMORS,KISA,INVDIA,ILUM,ILUI,Z,R)
  3.       IMPLICIT REAL*8 (A-H,O-Z)
  4.       IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM         : PSOLVE
  7. C DESCRIPTION : Preconditionner solve : résoud un système avec
  8. C               un préconditionneur
  9. C               Fait l'interface pour les différents types de
  10. C               préconditionneurs
  11. C
  12. C LANGAGE     : ESOPE
  13. C AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  14. C               mél : gounand@semt2.smts.cea.fr
  15. C***********************************************************************
  16. C APPELES          :
  17. C APPELES (E/S)    :
  18. C APPELES (BLAS)   :
  19. C APPELES (CALCUL) :
  20. C APPELE PAR       :
  21. C***********************************************************************
  22. C SYNTAXE GIBIANE    :
  23. C ENTREES            :
  24. C ENTREES/SORTIES    :
  25. C SORTIES            :
  26. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  27. C***********************************************************************
  28. C VERSION    : v1, 09/02/2006, version initiale
  29. C HISTORIQUE : v1, 09/02/2006, création
  30. C HISTORIQUE :
  31. C HISTORIQUE :
  32. C***********************************************************************
  33. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  34. C en cas de modification de ce sous-programme afin de faciliter
  35. C la maintenance !
  36. C***********************************************************************
  37.  
  38. -INC PPARAM
  39. -INC CCOPTIO
  41.       POINTEUR KMORS.PMORS,KISA.IZA
  42.       POINTEUR INVDIA.IZA,ILUM.PMORS,ILUI.IZA
  43.       POINTEUR Z.IZA,R.IZA
  44. *
  45.       CHARACTER*1        TRANS
  46.       INTEGER IMPR,IRET
  47. *
  48. * Executable statements
  49. *
  50.       IF (KPREC.EQ.0) THEN
  51.          CALL GCOPY(R,Z)
  52.       ELSEIF (KPREC.EQ.1) THEN
  53.          N=INVDIA.A(/1)
  54.          CALL PSDIAG(N,INVDIA.A,Z.A,R.A)
  55.       ELSEIF (KPREC.EQ.2) THEN
  56.          N=KMORS.IA(/1)-1
  57.          NNZ=KMORS.JA(/1)
  58.          CALL PSDILU(N,NNZ,KMORS.IA,KMORS.JA,KISA.A,
  59.      $        INVDIA.A,
  60.      $        Z.A,R.A)
  61.       ELSEIF (KPREC.GE.3.AND.KPREC.LE.10) THEN
  62.          N=KMORS.IA(/1)-1
  63.          NNZ=KMORS.JA(/1)
  64. *         N=ILUM.IA(/1)-1
  65. *         NNZ=ILUM.JA(/1)
  66.          CALL PSILU0(TRANS,N,NNZ,ILUM.IA,ILUM.JA,ILUI.A,
  67.      $        Z.A,R.A)
  68.       ELSE
  69.          WRITE(IOIMP,*) 'KPREC=',KPREC,' not allowed'
  70.          GOTO 9999
  71.       ENDIF
  72. *
  73. * Normal termination
  74. *
  75.       RETURN
  76. *
  77. * Format handling
  78. *
  79. *
  80. * Error handling
  81. *
  82.  9999 CONTINUE
  83.       CALL ERREUR(5)
  84.       WRITE(IOIMP,*) 'An error was detected in subroutine psolve'
  85.       RETURN
  86. *
  87. * End of subroutine PSOLVE
  88. *
  89.       END
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  

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