OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
mintp_re.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| mintp_re ../common_source/eos/mintp_re.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_eos_sesame ../starter/source/materials/eos/hm_read_eos_sesame.F
27!|| hm_read_mat26 ../starter/source/materials/mat/mat026/hm_read_mat26.F
28!|| sesa20 ../engine/source/materials/mat/mat026/sesa20.F
29!|| sesame ../common_source/eos/sesame.F
30!||--- calls -----------------------------------------------------
31!|| minter1d_rat ../common_source/eos/minter1d_rat.F
32!||====================================================================
33 SUBROUTINE mintp_re(XX,NX,YY,NY,ZZ,X,Y,Z,DYDZ)
34C=======================================================================
35C
36C | | | |
37C --+---------+--------*-------+----------+-- YY(IY+2)
38C INPUT : | | Z3 | |
39C XX(NX) | | | |
40C NX | | | |
41C YY(NY) | | | |
42C ZZ(NX,NY) | Z(IX,IY+1) Z(IX+1,IY+1) |
43C X --+---------+--------*-------+----------+-- YY(IY+1)
44C Z | | Z2 | |
45C | | | |
46C OUTPUT: | | | |
47C Y | | Z(X,Y) | |
48C | | * - - - | - - - - -| - - - Y
49C | | . | |
50C | | . | |
51C | Z(IX,IY) . Z(IX+1,IY) |
52C --+---------+--------*-------+----------+-- YY(IY)
53C | | Z1 | |
54C | | . | |
55C | | . | |
56C | | . | |
57C | | Z0 | |
58C --+---------+--------*-------+----------+-- YY(IY-1)
59C | | . | |
60C | | . | |
61C XX(IX-1) XX(IX) . XX(IX+1) XX(IX+2)
62C X
63C
64C=======================================================================
65C-----------------------------------------------
66C I m p l i c i t T y p e s
67C-----------------------------------------------
68#include "implicit_f.inc"
69C-----------------------------------------------
70C D u m m y A r g u m e n t s
71C-----------------------------------------------
72 INTEGER NX, NY
73 my_real
74 . x, y, z, dydz
75 my_real
76 . xx(nx), yy(ny), zz(nx,ny)
77C-----------------------------------------------
78C L o c a l V a r i a b l e s
79C-----------------------------------------------
80 INTEGER IX, IY, IXM1, IXP1, IXP2, IYM1, IYP1, IYP2
81C REAL
82 my_real
83 . xx0, xx1, xx2, xx3, zz00, zz10, zz20, zz30, zz01, zz11, zz21,
84 . zz31, zz02, zz12, zz22, zz32, zz03, zz13, zz23, zz33, z0,
85 . dzdx0, z1, dzdx1, z2, dzdx2, z3, dzdx3, dzdx, d2zdzdx
86C-----------------------------------------------
87C E x t e r n a l F u n c t i o n s
88C-----------------------------------------------
89 INTEGER MINDEX
90C-----------------------------------------------
91 ix = mindex(xx,1,nx,x)
92 iy = mindex(zz(ix,1),nx,ny,z)
93C
94 ix = max(1,min(nx-1,ix))
95 iy = max(1,min(ny-1,iy))
96C
97 ixm1 = max(1,ix-1)
98 ixp1 = ix+1
99 ixp2 = min(nx,ix+2)
100C
101 iym1 = max(1,iy-1)
102 iyp1 = iy+1
103 iyp2 = min(ny,iy+2)
104C
105 xx0 = xx(ixm1)
106 xx1 = xx(ix)
107 xx2 = xx(ixp1)
108 xx3 = xx(ixp2)
109C
110 zz00 = zz(ixm1,iym1)
111 zz10 = zz(ix ,iym1)
112 zz20 = zz(ixp1,iym1)
113 zz30 = zz(ixp2,iym1)
114C
115 zz01 = zz(ixm1,iy)
116 zz11 = zz(ix ,iy)
117 zz21 = zz(ixp1,iy)
118 zz31 = zz(ixp2,iy)
119C
120 zz02 = zz(ixm1,iyp1)
121 zz12 = zz(ix ,iyp1)
122 zz22 = zz(ixp1,iyp1)
123 zz32 = zz(ixp2,iyp1)
124C
125 zz03 = zz(ixm1,iyp2)
126 zz13 = zz(ix ,iyp2)
127 zz23 = zz(ixp1,iyp2)
128 zz33 = zz(ixp2,iyp2)
129C
130 CALL minter1d_rat(
131 + xx0 ,xx1 ,xx2 ,xx3 ,zz00 ,zz10 ,zz20 ,zz30 ,x ,z0 ,dzdx0,ix,nx)
132 CALL minter1d_rat(
133 + xx0 ,xx1 ,xx2 ,xx3 ,zz01 ,zz11 ,zz21 ,zz31 ,x ,z1 ,dzdx1,ix,nx)
134 CALL minter1d_rat(
135 + xx0 ,xx1 ,xx2 ,xx3 ,zz02 ,zz12 ,zz22 ,zz32 ,x ,z2 ,dzdx2,ix,nx)
136 CALL minter1d_rat(
137 + xx0 ,xx1 ,xx2 ,xx3 ,zz03 ,zz13 ,zz23 ,zz33 ,x ,z3 ,dzdx3,ix,nx)
138C
139 CALL minter1d_rat(
140 + z0 ,z1 ,z2 ,z3 ,
141 + yy(iym1) ,yy(iy) ,yy(iyp1) ,yy(iyp2),
142 + z ,y ,dydz ,iy ,ny)
143 CALL minter1d_rat(
144 + z0 ,z1 ,z2 ,z3 ,
145 + dzdx0 ,dzdx1 ,dzdx2 ,dzdx3 ,
146 + z ,dzdx ,d2zdzdx ,iy ,ny)
147C
148 RETURN
149 END
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
minter1d_rat
subroutine minter1d_rat(x0, x1, x2, x3, y0, y1, y2, y3, x, y, yp, i, n)
Definition minter1d_rat.F:33
mintp_re
subroutine mintp_re(xx, nx, yy, ny, zz, x, y, z, dydz)
Definition mintp_re.F:34