Télécharger pmixt.eso

Retour à la liste

Numérotation des lignes :

  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
  10. CALL LIROBJ('POINT ',IP1,1,IRETOU)
  11. CALL LIROBJ('POINT ',IP2,1,IRETOU)
  12. IF (IERR.NE.0) RETURN
  13. SEGACT MCOORD
  14. IREF1=(IDIM+1)*(IP1-1)
  15. IREF2=(IDIM+1)*(IP2-1)
  16. X1=XCOOR(IREF1+1)
  17. X2=XCOOR(IREF2+1)
  18. Y1=XCOOR(IREF1+2)
  19. Y2=XCOOR(IREF2+2)
  20. IF (IDIM.EQ.3) GOTO 10
  21. PM=X1*Y2-Y1*X2
  22. GOTO 20
  23. 10 CALL LIROBJ('POINT ',IP3,1,IRETOU)
  24. IF (IERR.NE.0) RETURN
  25. IREF3=(IDIM+1)*(IP3-1)
  26. Z1=XCOOR(IREF1+3)
  27. Z2=XCOOR(IREF2+3)
  28. X3=XCOOR(IREF3+1)
  29. Y3=XCOOR(IREF3+2)
  30. Z3=XCOOR(IREF3+3)
  31. PM=X1*(Y2*Z3-Y3*Z2)+Y1*(Z2*X3-Z3*X2)+Z1*(X2*Y3-X3*Y2)
  32. 20 CONTINUE
  33. CALL ECRREE(PM)
  34. RETURN
  35. END
  36.  
  37.  
  38.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales