OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
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.

37C-----------------------------------------------
38C I m p l i c i t T y p e s
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C C o m m o n B l o c k s
43C-----------------------------------------------
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 INTEGER, INTENT(IN) :: NEL
48C REAL
50 . aj1(*),aj2(*),aj3(*),
51 . aj4(*),aj5(*),aj6(*),
52 . aj7(*),aj8(*),aj9(*),
53 . det(*)
54C-----------------------------------------------
55C L o c a l V a r i a b l e s
56C-----------------------------------------------
57 INTEGER I
58C REAL
59C 12
61 . jac_59_68, jac_67_49, jac_48_57,
62 . aj11, aj22, aj33
63C-----------------------------------------------
64C JACOBIAN of [AJ]+[1]
65C
66 DO i=1,nel
67 aj11=aj1(i)+one
68 aj22=aj5(i)+one
69 aj33=aj9(i)+one
70 jac_59_68=aj22*aj33-aj6(i)*aj8(i)
71 jac_67_49=aj6(i)*aj7(i)-aj4(i)*aj33
72 jac_48_57=aj4(i)*aj8(i)-aj22*aj7(i)
73C
74 det(i)=aj11*jac_59_68+aj2(i)*jac_67_49+aj3(i)*jac_48_57
75 ENDDO
76C
77 RETURN
#define my_real
Definition cppsort.cpp:32