Esope source files

   1 : C POLYN2    SOURCE    PV        07/11/23    21:18:38     5978           
   2 :       FUNCTION POLYN2 (COEFP,LLIG,LCOL,V,U)                                     
   3 : ************************************************************************        
   4 : *                                                                               
   5 : *                             P O L Y N 2                                       
   6 : *                             -----------                                       
   7 : *                                                                               
   8 : * FONCTION:                                                                     
   9 : * ---------                                                                     
  10 : *                                                                               
  11 : *     DONNER LA VALEUR D'UN POLYNOME A 2 VARIABLES "U ET "V", TEL QUE:          
  12 : *                                                                               
  13 : *                            | P00 P01 P02 ... |   | 1  |                       
  14 : *               2  3         | P10 P11 P12 ... |   | U  |                       
  15 : *         (1 V V  V  ... ) X | P20 P21 P22 ... | X | U2 |                       
  16 : *                            | P30 P31 P32 ... |   | .. |                       
  17 : *                            |       ...       |                                
  18 : *                                                                               
  19 : *     POUR UN COUPLE (U,V) DONNE.                                               
  20 : *                                                                               
  21 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE   (+ = CONTENU DANS UN COMMUN)          
  22 : * -----------                                                                   
  23 : *                                                                               
  24 : *     COEFP   (E)  COEFFICIENTS DU POLYNOME (MATRICE RECTANGULAIRE).            
  25 : *     LLIG    (E)  NOMBRE DE LIGNES DE LA MATRICE.                              
  26 : *     LCOL    (E)  NOMBRE DE COLONNES DE LA MATRICE.                            
  27 : *     V,U     (E)  COUPLE EN LEQUEL ON CHERCHE LA VALEUR DU POLYNOME.           
  28 : *     POLYN2  (S)  VALEUR DU POLYNOME.                                          
  29 : *                  = 0  SI LE TABLEAU "COEFP" EST VIDE.                         
  30 : *                                                                               
  31 :       IMPLICIT INTEGER(I-N)
  32 :       IMPLICIT real*8 (a-h,o-z)
  33 :       INTEGER LLIG,LCOL                                                         
  34 :       REAL*8 COEFP(LLIG,LCOL),V,U                                               
  35 :       REAL*8 POLYN2                                                             
  36 : *                                                                               
  37 : * VARIABLES:                                                                    
  38 : * ----------                                                                    
  39 : *                                                                               
  40 :       REAL*8 ZZ                                                                 
  41 : *                                                                               
  42 : * FONCTIONS:                                                                    
  43 : * ----------                                                                    
  44 : *                                                                               
  45 :       REAL*8 POLYNO                                                             
  46 : *                                                                               
  47 : * AUTEUR, DATE DE CREATION:                                                     
  48 : * -------------------------                                                     
  49 : *                                                                               
  50 : *     PASCAL MANIGOT     24 FEVRIER 1987                                        
  51 : *                                                                               
  52 : * LANGAGE:                                                                      
  53 : * --------                                                                      
  54 : *                                                                               
  55 : *     FORTRAN77 + EXTENSION: DECLARATION "REAL*8".                              
  56 : *                                                                               
  57 : ************************************************************************        
  58 : *                                                                               
  59 : *     ON FAIT D'ABORD LA MULTIPLICATION A GAUCHE, POUR EVITER LES GRANDS        
  60 : *     DEPLACEMENTS DE POINTEUR.                                                 
  61 : *                                                                               
  62 :       ZZ = 0.D0                                                                 
  63 :       DO 100 IB=LCOL,1,-1                                                       
  64 :          ZZ = ZZ * U + POLYNO (COEFP(1,IB),LLIG,1,V)                            
  65 :   100    CONTINUE                                                               
  66 : *     END DO                                                                    
  67 :       POLYN2 = ZZ                                                               
  68 : *                                                                               
  69 :       END                                                                       
  70 :  
  71 :