OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
mintp_rt.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_rt ../common_source/eos/mintp_rt.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!|| m26th ../engine/source/materials/mat/mat026/m26th.F
29!|| sesa20 ../engine/source/materials/mat/mat026/sesa20.F
30!|| sesame ../common_source/eos/sesame.F
31!||--- calls -----------------------------------------------------
32!|| minter1d_rat ../common_source/eos/minter1d_rat.F
33!||====================================================================
34 SUBROUTINE mintp_rt(XX,NX,YY,NY,ZZ,X,Y,Z,DZDX)
35C=======================================================================
36C
37C | | | |
38C --+---------+--------*-------+----------+-- YY(IY+2)
39C INPUT : | | Z3 | |
40C XX(NX) | | | |
41C NX | | | |
42C YY(NY) | | | |
43C ZZ(NX,NY) | Z(IX,IY+1) Z(IX+1,IY+1) |
44C X --+---------+--------*-------+----------+-- YY(IY+1)
45C Y | | Z2 | |
46C | | | |
47C OUTPUT: | | | |
48C Z | | Z(X,Y) | |
49C | | * - - - | - - - - -| - - - Y
50C | | . | |
51C | | . | |
52C | Z(IX,IY) . Z(IX+1,IY) |
53C --+---------+--------*-------+----------+-- YY(IY)
54C | | Z1 | |
55C | | . | |
56C | | . | |
57C | | Z0 | |
58C --+---------+--------*-------+----------+-- YY(IY-1)
59C | | . | |
60C XX(IX-1) XX(IX) . XX(IX+1) XX(IX+2)
61C X
62C-----------------------------------------------
63C I m p l i c i t T y p e s
64C-----------------------------------------------
65#include "implicit_f.inc"
66C-----------------------------------------------
67C D u m m y A r g u m e n t s
68C-----------------------------------------------
69 INTEGER NX, NY
70C REAL
71 my_real
72 . x, y, z, dzdx
73C REAL
74 my_real
75 . xx(nx), yy(ny), zz(nx,ny)
76C-----------------------------------------------
77C L o c a l V a r i a b l e s
78C-----------------------------------------------
79 INTEGER IX, IY, IXM1, IXP1, IXP2, IYM1, IYP1, IYP2
80C REAL
81 my_real
82 . xx0, xx1, xx2, xx3, zz00, zz10, zz20, zz30, zz01, zz11, zz21,
83 . zz31, zz02, zz12, zz22, zz32, zz03, zz13, zz23, zz33, z0,
84 . dzdx0, z1, dzdx1, z2, dzdx2, z3, dzdx3, dzdy, d2zdxdy
85C-----------------------------------------------
86C E x t e r n a l F u n c t i o n s
87C-----------------------------------------------
88 INTEGER MINDEX
89C-----------------------------------------------
90 ix = mindex(xx,1,nx,x)
91 iy = mindex(yy,1,ny,y)
92C
93 ix = max(1,min(nx-1,ix))
94 iy = max(1,min(ny-1,iy))
95C
96 ix = max(1,min(nx-1,ix))
97 iy = max(1,min(ny-1,iy))
98C
99 ixm1 = max(1,ix-1)
100 ixp1 = ix+1
101 ixp2 = min(nx,ix+2)
102C
103 iym1 = max(1,iy-1)
104 iyp1 = iy+1
105 iyp2 = min(ny,iy+2)
106C
107 xx0 = xx(ixm1)
108 xx1 = xx(ix)
109 xx2 = xx(ixp1)
110 xx3 = xx(ixp2)
111C
112 zz00 = zz(ixm1,iym1)
113 zz10 = zz(ix ,iym1)
114 zz20 = zz(ixp1,iym1)
115 zz30 = zz(ixp2,iym1)
116C
117 zz01 = zz(ixm1,iy)
118 zz11 = zz(ix ,iy)
119 zz21 = zz(ixp1,iy)
120 zz31 = zz(ixp2,iy)
121C
122 zz02 = zz(ixm1,iyp1)
123 zz12 = zz(ix ,iyp1)
124 zz22 = zz(ixp1,iyp1)
125 zz32 = zz(ixp2,iyp1)
126C
127 zz03 = zz(ixm1,iyp2)
128 zz13 = zz(ix ,iyp2)
129 zz23 = zz(ixp1,iyp2)
130 zz33 = zz(ixp2,iyp2)
131C
132 CALL minter1d_rat(
133 + xx0 ,xx1 ,xx2 ,xx3 ,zz00 ,zz10 ,zz20 ,zz30 ,x ,z0 ,dzdx0,ix,nx)
134 CALL minter1d_rat(
135 + xx0 ,xx1 ,xx2 ,xx3 ,zz01 ,zz11 ,zz21 ,zz31 ,x ,z1 ,dzdx1,ix,nx)
136 CALL minter1d_rat(
137 + xx0 ,xx1 ,xx2 ,xx3 ,zz02 ,zz12 ,zz22 ,zz32 ,x ,z2 ,dzdx2,ix,nx)
138 CALL minter1d_rat(
139 + xx0 ,xx1 ,xx2 ,xx3 ,zz03 ,zz13 ,zz23 ,zz33 ,x ,z3 ,dzdx3,ix,nx)
140C
141 CALL minter1d_rat(
142 + yy(iym1) ,yy(iy) ,yy(iyp1) ,yy(iyp2),
143 + z0 ,z1 ,z2 ,z3 ,
144 + y ,z ,dzdy ,iy ,ny)
145C
146 CALL minter1d_rat(
147 + yy(iym1) ,yy(iy) ,yy(iyp1) ,yy(iyp2),
148 + dzdx0 ,dzdx1 ,dzdx2 ,dzdx3 ,
149 + y ,dzdx ,d2zdxdy ,iy ,ny)
150C
151 RETURN
152 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_rt
subroutine mintp_rt(xx, nx, yy, ny, zz, x, y, z, dzdx)
Definition mintp_rt.F:35