ajac3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ajac3_v (v1, v2, v3, aj1, aj2, aj3, ajv, nel)

Function/Subroutine Documentation

ajac3_v()

subroutine ajac3_v	(		v1,
			v2,
			v3,
			aj1,
			aj2,
			aj3,
			ajv,
		integer, intent(in)	nel )

Definition at line 28 of file ajac3.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-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   v1(mvsiz,3),v2(mvsiz,3),v3(mvsiz,3),aj1(mvsiz,3),aj2(mvsiz,3),aj3(mvsiz,3),ajv(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      DO i=1,nel
       aj1(i,1)=v2(i,2)*v3(i,3)-v2(i,3)*v3(i,2)
       aj1(i,2)=v2(i,3)*v3(i,1)-v2(i,1)*v3(i,3)
       aj1(i,3)=v2(i,1)*v3(i,2)-v2(i,2)*v3(i,1)
       aj2(i,1)=v3(i,2)*v1(i,3)-v3(i,3)*v1(i,2)
       aj2(i,2)=v3(i,3)*v1(i,1)-v3(i,1)*v1(i,3)
       aj2(i,3)=v3(i,1)*v1(i,2)-v3(i,2)*v1(i,1)
       aj3(i,1)=v1(i,2)*v2(i,3)-v1(i,3)*v2(i,2)
       aj3(i,2)=v1(i,3)*v2(i,1)-v1(i,1)*v2(i,3)
       aj3(i,3)=v1(i,1)*v2(i,2)-v1(i,2)*v2(i,1)
       ajv(i)=v1(i,1)*aj1(i,1)+v1(i,2)*aj1(i,2)+v1(i,3)*aj1(i,3)
      ENDDO
      RETURN