Esope source files

   1 : C VRFMOD    SOURCE    CHAT      05/01/13    04:10:03     5004           
   2 :       SUBROUTINE VRFMOD (IPA,IPB,IPX,VALPP)                                     
   3 : ************************************************************************        
   4 : *                                                                               
   5 : *                             V R F M O D                                       
   6 : *                             -----------                                       
   7 : *                                                                               
   8 : * FONCTION:                                                                     
   9 : * ---------                                                                     
  10 : *                                                                               
  11 : *     PERMETTRE DE VERIFIER QU'UN MODE PROPRE (V,X) SOLUTION DE                 
  12 : *        |A|.(X) = V.|B|.(X)                                                    
  13 : *     A ETE CALCULE DE FACON EXACTE, EN FAISANT LES IMPRESSIONS DE              
  14 : *     |A|.(X) , V.|B|.(X)  ET  |A|.(X)-V.|B|.(X)                                
  15 : *                                                                               
  16 : * MODE D'APPEL:                                                                 
  17 : * -------------                                                                 
  18 : *                                                                               
  18 : *     CALL VRFMOD (IPA,IPB,IPX,VALPP)                                           
  19 : *                                                                               
  20 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE                                         
  21 : * -----------                                                                   
  22 : *                                                                               
  23 : *     IPA     ENTIER    (E)  POINTEUR SUR LA 'RIGIDITE' "A".                    
  24 : *     IPB     ENTIER    (E)  POINTEUR SUR LA 'RIGIDITE' "B".                    
  25 : *     IPX     ENTIER    (E)  POINTEUR SUR LE 'CHPOINT' "X".                     
  26 : *     VALPP   REEL DP   (E)  VALEUR PROPRE "V".                                 
  27 : *                                                                               
  28 : * LEXIQUE: (ORDRE ALPHABETIQUE)                                                 
  29 : * --------                                                                      
  30 : *                                                                               
  31 : *     IPAX    ENTIER    POINTEUR SUR LE 'CHPOINT' |A|.(X)                       
  32 : *     IPBX    ENTIER    POINTEUR SUR LE 'CHPOINT' |B|.(X)                       
  33 : *     IPVBX   ENTIER    POINTEUR SUR LE 'CHPOINT' V.|B|.(X)                     
  34 : *     IPZERO  ENTIER    POINTEUR SUR LE 'CHPOINT' |A|.(X) - V.|B|.(X)           
  35 : *                                                                               
  36 : * SOUS-PROGRAMMES APPELES:                                                      
  37 : * ------------------------                                                      
  38 : *                                                                               
  39 : *     ADCHPO, DTCHPO, ECCHPO, MUCHPO, MUCPRI.                                   
  40 : *                                                                               
  41 : * AUTEUR, DATE DE CREATION:                                                     
  42 : * -------------------------                                                     
  43 : *                                                                               
  44 : *     PASCAL MANIGOT     25 OCTOBRE 1984                                        
  45 : *                                                                               
  46 : * LANGAGE:                                                                      
  47 : * --------                                                                      
  48 : *                                                                               
  49 : *     FORTRAN77                                                                 
  50 : *                                                                               
  51 : ************************************************************************        
  52 : *                                                                               
  53 :       IMPLICIT INTEGER(I-N)
  54 : -INC CCOPTIO                                                                    
  55 : *                                                                               
  56 :       REAL*8 VALPP                                                              
  57 :       REAL*8 V                                                                  
  58 : *                                                                               
  59 : *     PRODUIT |A|.(X)                                                           
  59 :       CALL MUCPRI (IPX,IPA,  IPAX)                                              
  60 :       IF (IERR .NE. 0) RETURN                                                   
  61 : *                                                                               
  62 : *     IMPRESSION DE |A|.(X)                                                     
  63 :       WRITE (IOIMP,2010) '|A|.(X)            '                                  
  64 :  2010 FORMAT ('1 IMPRESSION DE ',A20,/)                                         
  64 :       CALL ECCHPO (IPAX,0)                                                        
  65 : *                                                                               
  66 : *     PRODUIT |B|.(X)                                                           
  66 :       CALL MUCPRI (IPX,IPB,  IPBX)                                              
  67 :       IF (IERR .NE. 0) RETURN                                                   
  68 : *                                                                               
  69 : *     PRODUIT V.|B|.(X)                                                         
  70 :       MULT = 1                                                                  
  71 :       V = VALPP                                                                 
  71 :       CALL MUCHPO (IPBX,V,  IPVBX,MULT)                                         
  72 :       IF (IERR .NE. 0) RETURN                                                   
  73 : *                                                                               
  74 : *     IMPRESSION DE V.|B|.(X)                                                   
  75 :       WRITE (IOIMP,2010) 'V.|B|.(X)          '                                  
  75 :       CALL ECCHPO (IPVBX,0)                                                       
  76 : *                                                                               
  77 : *     DIFFERENCE  |A|.(X) - V.|B|.(X)                                           
  78 :       MOINS = -1                                                                
  78 :       CALL ADCHPO (IPAX,IPVBX,  IPZERO,1.D0,-1.D0)
  79 :       IF (IERR .NE. 0) RETURN                                                   
  80 : *                                                                               
  81 : *     IMPRESSION DE |A|.(X) - V.|B|.(X)                                         
  82 :       WRITE (IOIMP,2010) '|A|.(X) - V.|B|.(X)'                                  
  82 :       CALL ECCHPO (IPZERO,0)                                                      
  83 : *                                                                               
  84 : *     DESTRUCTION DES OBJETS DE TRAVAIL:                                        
  84 :       CALL DTCHPO (IPZERO)                                                      
  84 :       CALL DTCHPO (IPAX)                                                        
  84 :       CALL DTCHPO (IPBX)                                                        
  84 :       CALL DTCHPO (IPVBX)                                                       
  85 : *                                                                               
  86 :       END                                                                       
  87 :  
  88 :