Esope source files

   1 : C DBST      SOURCE    CHAT      05/01/12    22:37:28     5004           
   2 :       SUBROUTINE DBST(BGENE,DDHOOK,XDDL,LRE,NSTRS,XSTRS)                        
   3 : C======================================================================         
   4 : C                                                                               
   5 : C  CALCULE LE PRODUIT DDHOOK*BGENE=XSTRS                                        
   6 : C                                                                               
   7 : C  ENTREE                                                                       
   8 : C    BGENE(NSTRS,LRE)   =MATRICE B                                              
   9 : C    DDHOOK(NSTRS,NSTRS)=MATRICE DE HOOKE                                       
  10 : C    XDDL(LRE)          =D D L DE LA FORMULATION                                
  11 : C    LRE                =NOMBRE DE D D L                                        
  12 : C    NSTRS              =NOMBRE DE CONTRAINTES                                  
  13 : C  SORTIE                                                                       
  14 : C    XSTRS(NSTRS)       =CONTRAINTES                                            
  15 : C                                                                               
  16 : C   CODE EBERSOLT JUILLET 84                                                    
  17 : C=======================================================================        
  18 :       IMPLICIT INTEGER(I-N)
  19 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  20 :       DIMENSION BGENE(NSTRS,1),DDHOOK(NSTRS,1),XSTRS(1),XDDL(1)                 
  21 :       DIMENSION XWRK(10)                                                        
  22 :       DO 100 IA=1,NSTRS                                                         
  23 :       CC=0.D0                                                                   
  24 :       DO 200 IB=1,LRE                                                           
  25 :       CC=CC+XDDL(IB)*BGENE(IA,IB)                                               
  26 :   200 CONTINUE                                                                  
  27 :       XWRK(IA)=CC                                                               
  28 :   100 CONTINUE                                                                  
  29 :       DO 300 IA=1,NSTRS                                                         
  30 :       CC=0.D0                                                                   
  31 :       DO 400 IB=1,NSTRS                                                         
  32 :       CC=CC+DDHOOK(IA,IB)*XWRK(IB)                                              
  33 :   400 CONTINUE                                                                  
  34 :       XSTRS(IA)=CC                                                              
  35 :   300 CONTINUE                                                                  
  36 :       RETURN                                                                    
  37 :       END                                                                       
  38 :  
  39 :