Esope source files

   1 : C DIMNO1    SOURCE    CHAT      09/10/09    21:17:15     6519           
   2 :       SUBROUTINE DIMNO1 (IPRIGI,NOYAU)                                          
   3 : ************************************************************************        
   4 : *                                                                               
   5 : *                             D I M N O 1                                       
   6 : *                             -----------                                       
   7 : *                                                                               
   8 : * FONCTION:                                                                     
   9 : * ---------                                                                     
  10 : *                                                                               
  11 : *     DIMENSION DU NOYAU D'UN OBJET 'RIGIDITE'.                                 
  12 : *                                                                               
  13 : * MODE D'APPEL:                                                                 
  14 : * -------------                                                                 
  15 : *                                                                               
  15 : *     CALL DIMNO1 (IPRIGI,NOYAU)                                                
  16 : *                                                                               
  17 : * PARAMETRES:  (E)=ENTREE   (S)=SORTIE                                          
  18 : * -----------                                                                   
  19 : *                                                                               
  20 : *     IPRIGI  ENTIER    (E)  POINTEUR D'UNE 'RIGIDITE'.                         
  21 : *     NOYAU   ENTIER    (S)  DIMENSION DU NOYAU DE LA 'RIGIDITE'.               
  22 : *                                                                               
  23 : * REMARQUES:                                                                    
  24 : * ----------                                                                    
  25 : *                                                                               
  26 : *     LA VALEUR "NOYAU" RETOURNEE PEUT ETRE SUPERIEURE A LA REALITE:            
  27 : *     VOIR LA DEFINITION DE "NENS", CALCULE DANS LE SOUS-PROGRAMME              
  28 : *     TRAITANT DE LA DECOMPOSITION EN "L.D.LT".                                 
  29 : *     (LA DETECTION D'UN PIVOT NUL -- AU LIEU DE PIVOT MAXIMUM NUL --           
  30 : *     EST ABUSIVEMENT ASSIMILEE A UNE SINGULARITE DE MATRICE)                   
  31 : *                                                                               
  32 : * AUTEUR, DATE DE CREATION:                                                     
  33 : * -------------------------                                                     
  34 : *                                                                               
  35 : *     PASCAL MANIGOT     10 SEPTEMBRE 1985                                      
  36 : *                                                                               
  37 : * LANGAGE:                                                                      
  38 : * --------                                                                      
  39 : *                                                                               
  40 : *     ESOPE + FORTRAN77                                                         
  41 : *                                                                               
  42 : ************************************************************************        
  43 : *                                                                               
  44 :       IMPLICIT INTEGER(I-N)
  45 : -INC CCOPTIO                                                                    
  46 : -INC SMMATRI                                                                    
  47 : -INC SMRIGID                                                                    
  48 : *
  49 :       integer insym
  50 :       insym = 0  
  51 :       MRIGID = IPRIGI                                                           
  52 :       SEGACT,MRIGID  
  53 :       NRG = IRIGEL(/1)
  54 :       NBR = IRIGEL(/2)
  55 :       IPMATR = ICHOLE 
  56 :       IF(NORINC.GT.0  .AND. NORIND.GT.0) THEN
  57 :          INSYM = 1
  58 :       ENDIF
  59 :       IF (NRG.GE.7) THEN
  60 :          DO  9 IN = 1,NBR
  61 :             IANTI=IRIGEL(7,IN)
  62 :             IF(IANTI.GT.0) THEN
  63 :                INSYM = 1
  64 :             ENDIF
  65 :   9      CONTINUE
  66 :       ENDIF                                                          
  67 :       SEGDES,MRIGID                                                       
  68 : *                                                                               
  69 :       IF (IPMATR .EQ. 0) THEN 
  70 :         IF (INSYM .EQ. 0) THEN   
  70 :          CALL TRIANG (IPRIGI,1D-18)
  71 :         ELSE
  71 :          CALL ldmt1(IPRIGI,1d-18)
  72 :         ENDIF
  73 :          IF (IERR .NE. 0) RETURN                                                
  74 :          MRIGID = IPRIGI                                                        
  75 :          SEGACT,MRIGID                                                          
  76 :          IPMATR = ICHOLE                                                        
  77 :          SEGDES,MRIGID                                                    
  78 :       END IF                                                                    
  79 : *                                                                               
  80 :       MMATRI = IPMATR                                                           
  81 :       SEGACT,MMATRI                                                             
  82 :       NOYAU = NENS                                                              
  83 :       SEGDES,MMATRI                                                       
  84 : *                                                                               
  85 :       END                                                                       
  86 :  
  87 :  
  88 :  
  89 :  
  90 :  
  91 :  
  92 :  
  93 :  
  94 :  
  95 :  
  96 :  
  97 :