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
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
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
#define my_real
Definition cppsort.cpp:32
subroutine imp_setba(f, m, iddl, ndof, b, iflag)
Definition imp_setb.F:135
subroutine imp_setb(f, m, iddl, ndof, b)
Definition imp_setb.F:40
subroutine imp_setbp(f, m, iddl, ndof, ikc, inloc, b)
Definition imp_setb.F:85
subroutine imp_solv(output, timers, python, icode, iskew, iskwn, ipart, ixtg, ixs, ixq, ixc, ixt, ixp, ixr, ixtg1, itab, itabm1, npc, ibcl, ibfv, sensor_tab, nnlink, lnlink, iparg, igrv, ipari, intbuf_tab, nprw, iconx, npby, lpby, lrivet, nstrf, ljoint, icodt, icodr, isky, adsky, iads_f, ilink, llink, weight, itask, ibvel, lbvel, fbvel, x, d, v, vr, dr, thke, damp, ms, in, pm, skews, geo, eani, bufmat, bufgeo, bufsf, tf, forc, vel, fsav, agrv, fr_wave, parts0, elbuf, rby, rivet, fr_elem, iad_elem, wa, a, ar, stifn, stifr, partsav, fsky, fskyi, iframe, xframe, w16, iactiv, fskym, igeo, ipm, wfext, nodft, nodlt, nint7, num_imp, ns_imp, ne_imp, ind_imp, it, rwbuf, lprw, fr_wall, nbintc, intlist, fopt, rwsav, fsavd, graphe, fac_k, ipiv_k, nkcond, nsensor, monvol, igrsurf, fr_mv, volmon, dirul, nodglob, mumps_par, cddlp, isendto, irecvfrom, newfront, imsch, i2msch, isizxv, ilenxv, islen7, irlen7, islen11, irlen11, islen17, irlen17, irlen7t, islen7t, kinet, num_imp1, temp, dt2prev, waint, lgrav, sh4tree, sh3tree, irlen20, islen20, irlen20t, islen20t, irlen20e, islen20e, irbe3, lrbe3, frbe3, fr_i2m, iad_i2m, fr_rbe3m, iad_rbe3m, frwl6, irbe2, lrbe2, intbuf_tab_c, ikine, diag_sms, icfield, lcfield, cfield, count_remslv, count_remslve, elbuf_tab, elbuf_imp, xdp, weight_md, stack, dimfb, fbsav6, stabsen, tabsensor, drape_sh4n, drape_sh3n, h3d_data, multi_fvm, igrbric, igrsh4n, igrsh3n, igrbeam, forneqs, maxdgap, nddl0, nnzk0, it_t, impbuf_tab, cptreac, fthreac, nodreac, drapeg, interfaces, th_surf, dpl0cld, vel0cld, snpc, stf, glob_therm, wfext_md)
Definition imp_solv.F:175
subroutine imp_intfr(num_imp, ns_imp, ne_imp, ipari, intbuf_tab, npby, lpby, itab, nrbyac, irbyac, nint2, iint2, iddl, ikc, ndof, inloc, nsrem, nsl, nbintc, intlist, x, ibfv, dirul, skew, xframe, iskew, icodt, de, d_imp, lb, ifdis, nddl, dr_imp, iddli, irbe3, lrbe3, frbe3, irbe2, lrbe2, dd, ddr, a, ar, ac, acr, ms, v, nddl0, r02, rby, icodr, nt_rw, w_ddl, weight, irflag)
Definition imp_solv.F:7331