i12rot3.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    i12rot3    ../engine/source/interfaces/inter3d/i12rot3.F
!||--- called by ------------------------------------------------------
!||    intti12a   ../engine/source/interfaces/interf/intti12.F
!||====================================================================
      SUBROUTINE i12rot3(XS,XM,X,SKEW,ISKWN,ILEV,IREF,ICEN,NSN,NMN,
     &                   TETS,TETM)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ILEV,IREF,ICEN,ISKWN(LISKN,*),NSN,NMN
      my_real
     .   xs(3,*),xm(3,*),x(3,*),skew(lskew,*),tets(*),tetm(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,ORIG
      my_real
     . xx(3),xo(3),xr(3),rr,tet,p(9)
C
      IF(ilev==1)THEN
C        ORIG=ISKWN(1,IREF+1)
        IF(icen>0)THEN
          xo(1)=x(1,icen)
          xo(2)=x(2,icen)
          xo(3)=x(3,icen)
        ELSE
          xo(1)=zero
          xo(2)=zero
          xo(3)=zero
        ENDIF
        DO i=1,9
          p(i)=skew(i,iref+1)
        ENDDO
        DO i=1,nsn
          xx(1)=xs(1,i)-xo(1)
          xx(2)=xs(2,i)-xo(2)
          xx(3)=xs(3,i)-xo(3)
          xr(1)=p(1)*xx(1)+p(2)*xx(2)+p(3)*xx(3)
          xr(2)=p(4)*xx(1)+p(5)*xx(2)+p(6)*xx(3)
          xr(3)=p(7)*xx(1)+p(8)*xx(2)+p(9)*xx(3)
          rr=sqrt(xr(2)**2+xr(3)**2)
          IF(xr(3)/=zero.AND.xr(2)/=zero) THEN
            tet=atan2(xr(3),xr(2))
          ELSE
            tet=zero
          ENDIF
          xs(1,i)=xr(1)
          xs(2,i)=rr
          xs(3,i)=tet
          tets(i)=tet
        ENDDO
        DO i=1,nmn
          xx(1)=xm(1,i)-xo(1)
          xx(2)=xm(2,i)-xo(2)
          xx(3)=xm(3,i)-xo(3)
          xr(1)=p(1)*xx(1)+p(2)*xx(2)+p(3)*xx(3)
          xr(2)=p(4)*xx(1)+p(5)*xx(2)+p(6)*xx(3)
          xr(3)=p(7)*xx(1)+p(8)*xx(2)+p(9)*xx(3)
          rr=sqrt(xr(2)**2+xr(3)**2)
          IF(xr(3)/=zero.AND.xr(2)/=zero) THEN
            tet=atan2(xr(3),xr(2))
          ELSE
            tet=zero
          ENDIF
          xm(1,i)=xr(1)
          xm(2,i)=rr
          xm(3,i)=tet
          tetm(i)=tet
        ENDDO
      ENDIF
 
      RETURN
      END