jacob_j33.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	jacob_j33 (det, aj1, aj2, aj3, aj4, aj5, aj6, aj7, aj8, aj9, nel)

Function/Subroutine Documentation

jacob_j33()

subroutine jacob_j33	(		det,
			aj1,
			aj2,
			aj3,
			aj4,
			aj5,
			aj6,
			aj7,
			aj8,
			aj9,
		integer, intent(in)	nel )

Definition at line 33 of file jacob_j33.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.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
C     REAL
      my_real
     .   aj1(*),aj2(*),aj3(*),
     .   aj4(*),aj5(*),aj6(*),
     .   aj7(*),aj8(*),aj9(*),
     .   det(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
C                                                                     12
      my_real
     .   jac_59_68, jac_67_49, jac_48_57,
     .   aj11, aj22, aj33
C-----------------------------------------------
C     JACOBIAN of [AJ]+[1]
C
      DO i=1,nel
       aj11=aj1(i)+one
       aj22=aj5(i)+one
       aj33=aj9(i)+one
       jac_59_68=aj22*aj33-aj6(i)*aj8(i)
       jac_67_49=aj6(i)*aj7(i)-aj4(i)*aj33
       jac_48_57=aj4(i)*aj8(i)-aj22*aj7(i)
C
       det(i)=aj11*jac_59_68+aj2(i)*jac_67_49+aj3(i)*jac_48_57
      ENDDO
C
      RETURN