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