Esope source files

   1 : C MATVE1    SOURCE    CHAT      05/01/13    01:35:41     5004           
   2 :       SUBROUTINE MATVE1(A,V,NL,NC,X,ICLE)                                       
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       DIMENSION A(NL,*),V(*),X(*)                                               
   6 : *                                                                               
   7 : * MULTIPLICATION DE:   X(NC)=V(NL)*A(NL,NC) SI ICLE=1                           
   8 : *                      X(NL)=A(NL,NC)*V(NC) SI ICLE=2                           
   9 : *  X CONTIENT LE RESULTAT                                                       
  10 : *  LE PREMIER INDICE EST LE NB DE LIGNE                                         
  11 :       IF(ICLE.EQ.1) THEN                                                        
  12 :         DO 1 I=1,NC                                                             
  13 :           XX= 0.D0                                                              
  14 :           DO 2 J=1,NL                                                           
  15 :           XX = A(J,I)*V(J)  +XX                                                 
  16 :     2     CONTINUE                                                              
  17 :           X(I)=XX                                                               
  18 :     1   CONTINUE                                                                
  19 :       ELSEIF(ICLE.EQ.2) THEN                                                    
  20 :         DO 3 I=1,NL                                                             
  21 :           XX= 0.D0                                                              
  22 :           DO 4 J=1,NC                                                           
  23 :           XX = A(I,J)*V(J) + XX                                                 
  24 :     4     CONTINUE                                                              
  25 :           X(I)=XX                                                               
  26 :     3     CONTINUE                                                              
  27 :       ELSE                                                                      
  27 :           CALL  ERREUR(5)                                                       
  28 :       ENDIF                                                                     
  29 :       RETURN                                                                    
  30 :       END                                                                       
  31 :