Sources Esope | Cast3M

lspfrm
C LSPFRM    SOURCE    CHAT      05/01/13    01:25:38     5004
      SUBROUTINE LSPFRM(IWRK,KERRE,VECT,ICAS)
C-----------------------------------------------------------------------
C  ROUTINE DE REMISE A JOUR DU VECTEUR DEFINISSANT LE REPERE LOCAL
C                     POUR LE LINESPRING
C                HYPOTHESE DES PETITES ROTATIONS
C-----------------------------------------------------------------------
C  ENTREE
C          IWRK    POINTEUR SUR SEGMENT DE TRAVAIL
C          VECT    LE VECTEUR A REACTUALISER
C          ICAS    INDICE DU POINT A TRAITER
C  SORTIE
C          KERRE   INDICE D'ERREUR  ( 0 SI TOUT EST OK )
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
      SEGMENT/IWRK/(XDDL(LRE)*D,XE(3,NBNN)*D,WORK(LW)*D)
      DIMENSION OMEG(3),VECT(*)
C
      KERRE=0
C
      IF(ICAS.EQ.1) THEN
      OMEG(1)=0.5D0*(XDDL(4)+XDDL(22))
      OMEG(2)=0.5D0*(XDDL(5)+XDDL(23))
      OMEG(3)=0.5D0*(XDDL(6)+XDDL(24))
      ELSE IF(ICAS.EQ.2) THEN
      OMEG(1)=0.5D0*(XDDL(10)+XDDL(16))
      OMEG(2)=0.5D0*(XDDL(11)+XDDL(17))
      OMEG(3)=0.5D0*(XDDL(12)+XDDL(18))
      ENDIF
      CALL ROTVEC(VECT,OMEG,IDIM,KERRE)
      RETURN
      END