OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
imp_setb.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!|| imp_setb ../engine/source/implicit/imp_setb.F
25!||--- called by ------------------------------------------------------
26!|| d_to_u ../engine/source/implicit/produt_v.F
27!|| imp_chkm ../engine/source/implicit/imp_solv.F
28!|| imp_dykv ../engine/source/implicit/imp_dyna.F
29!|| imp_dykv0 ../engine/source/implicit/imp_dyna.F
30!|| imp_solv ../engine/source/implicit/imp_solv.F
31!|| produt_u ../engine/source/implicit/produt_v.F
32!|| produt_u2 ../engine/source/implicit/produt_v.F
33!|| produt_uh ../engine/source/implicit/produt_v.F
34!|| produt_uh2 ../engine/source/implicit/produt_v.F
35!|| produt_uhp ../engine/source/implicit/produt_v.f
36!|| produt_uhp2 ../engine/source/implicit/produt_v.F
37!|| upd_rhs ../engine/source/implicit/upd_glob_k.F
38!||====================================================================
39 SUBROUTINE imp_setb(F ,M ,IDDL ,NDOF ,B )
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C C o m m o n B l o c k s
46C-----------------------------------------------
47#include "com04_c.inc"
48C-----------------------------------------------
49C D u m m y A r g u m e n t s
50C-----------------------------------------------
51 INTEGER IDDL(*),NDOF(*)
52C REAL
53 my_real
54 . f(3,*),m(3,*),b(*)
55C-----------------------------------------------
56C L o c a l V a r i a b l e s
57C-----------------------------------------------
58 INTEGER I,J,K,ND,ID,JR
59C------------------------------------------
60 DO i = 1,numnod
61 nd = iddl(i)
62 DO j =1,ndof(i)
63 id = nd + j
64 IF (j>3) THEN
65 jr=j-3
66 b(id) = m(jr,i)
67 ELSE
68 b(id) = f(j,i)
69 ENDIF
70 ENDDO
71 ENDDO
72C--------------------------------------------
73 RETURN
74 END
75!||====================================================================
76!|| imp_setbp ../engine/source/implicit/imp_setb.F
77!||--- called by ------------------------------------------------------
78!|| ext_rhs ../engine/source/implicit/upd_glob_k.F
79!|| get_fext ../engine/source/implicit/imp_solv.F
80!|| rer02 ../engine/source/implicit/upd_glob_k.F
81!|| upd_rhs ../engine/source/implicit/upd_glob_k.f
82!||====================================================================
83 SUBROUTINE imp_setbp(F ,M ,IDDL ,NDOF ,IKC ,
84 . INLOC ,B )
85C-----------------------------------------------
86C I m p l i c i t T y p e s
87C-----------------------------------------------
88#include "implicit_f.inc"
89C-----------------------------------------------
90C C o m m o n B l o c k s
91C-----------------------------------------------
92#include "com04_c.inc"
93C-----------------------------------------------
94C D u m m y A r g u m e n t s
95C-----------------------------------------------
96 INTEGER IDDL(*),NDOF(*),IKC(*),INLOC(*)
97C REAL
98 my_real
99 . f(3,*),m(3,*),b(*)
100C-----------------------------------------------
101C L o c a l V a r i a b l e s
102C-----------------------------------------------
103 INTEGER I,J,N,K,ND,ID,JR,NKC
104C------------------------------------------
105 nkc=0
106 DO n = 1,numnod
107 i=inloc(n)
108 DO j=1,ndof(i)
109 nd = iddl(i)+j
110 id = nd-nkc
111 IF (ikc(nd)<1) THEN
112 IF (j<=3) THEN
113 b(id)=b(id)+f(j,i)
114 ELSE
115 b(id)=b(id)+m(j-3,i)
116 ENDIF
117 ELSE
118 nkc=nkc+1
119 ENDIF
120 ENDDO
121 ENDDO
122C--------------------------------------------
123 RETURN
124 END
125!||====================================================================
126!|| imp_setba ../engine/source/implicit/imp_setb.F
127!||--- called by ------------------------------------------------------
128!|| ext_rhs ../engine/source/implicit/upd_glob_k.F
129!|| get_fext ../engine/source/implicit/imp_solv.F
130!|| imp_intfr ../engine/source/implicit/imp_solv.F
131!|| upd_rhs ../engine/source/implicit/upd_glob_k.F
132!||====================================================================
133 SUBROUTINE imp_setba(F ,M ,IDDL ,NDOF ,B ,
134 1 IFLAG )
135C-----------------------------------------------
136C I m p l i c i t T y p e s
137C-----------------------------------------------
138#include "implicit_f.inc"
139C-----------------------------------------------
140C C o m m o n B l o c k s
141C-----------------------------------------------
142#include "com04_c.inc"
143C-----------------------------------------------
144C D u m m y A r g u m e n t s
145C-----------------------------------------------
146 INTEGER IDDL(*),NDOF(*),IFLAG
147C REAL
148 my_real
149 . f(3,*),m(3,*),b(*)
150C-----------------------------------------------
151C L o c a l V a r i a b l e s
152C-----------------------------------------------
153 INTEGER I,J,K,ND,ID,JR
154C------------------------------------------
155 IF (iflag==1) THEN
156 DO i = 1,numnod
157 nd = iddl(i)
158 DO j =1,ndof(i)
159 id = nd + j
160 IF (j>3) THEN
161 jr=j-3
162 b(id) = b(id)+m(jr,i)
163 ELSE
164 b(id) = b(id)+f(j,i)
165 ENDIF
166 ENDDO
167 ENDDO
168 ELSE
169 DO i = 1,numnod
170 nd = iddl(i)
171 DO j =1,ndof(i)
172 id = nd + j
173 IF (j>3) THEN
174 jr=j-3
175 b(id) = b(id)-m(jr,i)
176 ELSE
177 b(id) = b(id)-f(j,i)
178 ENDIF
179 ENDDO
180 ENDDO
181 ENDIF
182C--------------------------------------------
183 RETURN
184 END
my_real
#define my_real
Definition cppsort.cpp:32
imp_setba
subroutine imp_setba(f, m, iddl, ndof, b, iflag)
Definition imp_setb.F:135
imp_setb
subroutine imp_setb(f, m, iddl, ndof, b)
Definition imp_setb.F:40
imp_setbp
subroutine imp_setbp(f, m, iddl, ndof, ikc, inloc, b)
Definition imp_setb.F:85
produt_uhp
subroutine produt_uhp(nddl0, nddl, iddl, ndof, ikc, dd, ddr, norm2, w_imp)
Definition produt_v.F:3359
produt_v
subroutine produt_v(nddl, x, y, r)
Definition produt_v.F:33
upd_rhs
subroutine upd_rhs(icodt, icodr, iskew, ibfv, xframe, rby, x, skew, lpby, npby, nrbyac, irbyac, nint2, iint2, ipari, intbuf_tab, ndof, iddl, ikc, nddl0, b, iupd, inloc, lj, a, ar, ac, acr, nt_rw, irflag, w_ddl, nddl, r02, idyna, v, vr, ms, in, irbe3, lrbe3, frbe3, weight, irbe2, lrbe2)
Definition upd_glob_k.F:664
upd_glob_k
subroutine upd_glob_k(icodt, icodr, iskew, ibfv, npc, tf, vel, xframe, rby, x, skew, lpby, npby, itab, weight, ms, in, nrbyac, irbyac, nsc, isij, nmc, imij, nss, iss, nint2, iint2, nsc2, isij2, nss2, iss2, ipari, intbuf_tab, nddl, nnz, iadk, jdik, diag_k, lt_k, ndof, iddl, ikc, ud, b, nkud, ikud, bkud, nmc2, imij2, nt_rw, rd, lj, irbe3, lrbe3, frbe3, iss3, irbe2, lrbe2, isb2, nsrb2)
Definition upd_glob_k.F:66