Esope source files

   1 : C PMIXT     SOURCE    PV        07/11/23    21:18:35     5978           
   2 : C   CALCULE LE PRODUIT MIXTE DE IDIM VECTEURS                                   
   3 : C                                                                               
   4 :       SUBROUTINE PMIXT                                                          
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT real*8 (a-h,o-z)
   7 : -INC CCOPTIO                                                                    
   8 : -INC SMCOORD                                                                    
   9 :       REAL*8 PM                                                                 
   9 :       CALL LIROBJ('POINT ',IP1,1,IRETOU)                                        
   9 :       CALL LIROBJ('POINT ',IP2,1,IRETOU)                                        
  10 :       IF (IERR.NE.0) RETURN                                                     
  11 :       SEGACT MCOORD                                                             
  12 :       IREF1=(IDIM+1)*(IP1-1)                                                    
  13 :       IREF2=(IDIM+1)*(IP2-1)                                                    
  14 :       X1=XCOOR(IREF1+1)                                                         
  15 :       X2=XCOOR(IREF2+1)                                                         
  16 :       Y1=XCOOR(IREF1+2)                                                         
  17 :       Y2=XCOOR(IREF2+2)                                                         
  18 :       IF (IDIM.EQ.3) GOTO 10                                                    
  19 :       PM=X1*Y2-Y1*X2                                                            
  20 :       GOTO 20                                                                   
  20 :   10  CALL LIROBJ('POINT ',IP3,1,IRETOU)                                        
  21 :       IF (IERR.NE.0) RETURN                                                     
  22 :       IREF3=(IDIM+1)*(IP3-1)                                                    
  23 :       Z1=XCOOR(IREF1+3)                                                         
  24 :       Z2=XCOOR(IREF2+3)                                                         
  25 :       X3=XCOOR(IREF3+1)                                                         
  26 :       Y3=XCOOR(IREF3+2)                                                         
  27 :       Z3=XCOOR(IREF3+3)                                                         
  28 :       PM=X1*(Y2*Z3-Y3*Z2)+Y1*(Z2*X3-Z3*X2)+Z1*(X2*Y3-X3*Y2)                     
  29 :   20  CONTINUE                                                                  
  29 :       CALL ECRREE(PM)                                                           
  30 :       RETURN                                                                    
  31 :       END                                                                       
  32 :  
  33 :