qrrota2.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	qrrota2 (r22, r23, r32, r33, f11, f21, f12, f22, t11, t21, t12, t22, t13, t23, t14, t24, nel)

Function/Subroutine Documentation

qrrota2()

subroutine qrrota2	(		r22,
			r23,
			r32,
			r33,
			f11,
			f21,
			f12,
			f22,
			t11,
			t21,
			t12,
			t22,
			t13,
			t23,
			t14,
			t24,
		integer, intent(in)	nel )

Definition at line 28 of file qrrota2.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
C     REAL
      my_real
     .   r22(*), r23(*), r32(*), r33(*),
     .   f11(*), f21(*), f12(*), f22(*),
     .   t11(*), t21(*), t12(*), t22(*),
     .   t13(*), t23(*), t14(*), t24(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   y,z
C-----------------------------------------------
C     X=RX' <=> X'=t(R)X chgt de base.
      DO i=1,nel
        y=r22(i)*f11(i)+r32(i)*f21(i)
        z=r23(i)*f11(i)+r33(i)*f21(i)
        f11(i)=y
        f21(i)=z
      END DO
      DO i=1,nel
        y=r22(i)*f12(i)+r32(i)*f22(i)
        z=r23(i)*f12(i)+r33(i)*f22(i)
        f12(i)=y
        f22(i)=z
      END DO
      DO i=1,nel
        y=r22(i)*t11(i)+r32(i)*t21(i)
        z=r23(i)*t11(i)+r33(i)*t21(i)
        t11(i)=y
        t21(i)=z
      END DO
      DO i=1,nel
        y=r22(i)*t12(i)+r32(i)*t22(i)
        z=r23(i)*t12(i)+r33(i)*t22(i)
        t12(i)=y
        t22(i)=z
      END DO
      DO i=1,nel
        y=r22(i)*t13(i)+r32(i)*t23(i)
        z=r23(i)*t13(i)+r33(i)*t23(i)
        t13(i)=y
        t23(i)=z
      END DO
      DO i=1,nel
        y=r22(i)*t14(i)+r32(i)*t24(i)
        z=r23(i)*t14(i)+r33(i)*t24(i)
        t14(i)=y
        t24(i)=z
      END DO
c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
      RETURN