OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
c_h3d_update_solid_tensor.cpp
Go to the documentation of this file.
1//Copyright> OpenRadioss
2//Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3//Copyright>
4//Copyright> This program is free software: you can redistribute it and/or modify
5//Copyright> it under the terms of the GNU Affero General Public License as published by
6//Copyright> the Free Software Foundation, either version 3 of the License, or
7//Copyright> (at your option) any later version.
8//Copyright>
9//Copyright> This program is distributed in the hope that it will be useful,
10//Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11//Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12//Copyright> GNU Affero General Public License for more details.
13//Copyright>
14//Copyright> You should have received a copy of the GNU Affero General Public License
15//Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16//Copyright>
17//Copyright>
18//Copyright> Commercial Alternative: Altair Radioss Software
19//Copyright>
20//Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21//Copyright> software under a commercial license. Contact Altair to discuss further if the
22//Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23//
24#include <stdio.h>
25#include <string.h>
26#include <stdlib.h>
27#include <math.h>
28#include <fcntl.h>
29
30
31
32#ifdef _WIN32
33/* Windows includes */
34#include <windows.h>
35#include <process.h>
36#include <io.h>
37#include <sys\types.h>
38#include <sys/stat.h>
39
40
41#elif 1
42
43/* Linux includes */
44#include <sys/resource.h>
45#include <sys/types.h>
46#include <time.h>
47#include <sys/stat.h>
48#include <unistd.h>
49#include <dlfcn.h>
50#define _FCALL
51#include <math.h>
52#include <stdbool.h>
53
54#endif
55
56#include "h3dpublic_defs.h"
57#include "h3dpublic_export.h"
58
59#define _FCALL
60
61#include "h3d_values.h"
62
63#ifdef MYREAL4
64#define my_real float
65#endif
66
67#ifdef MYREAL8
68#define my_real double
69#endif
70
71extern "C"
72/*=================================================================*/
73{
74
75/*=================================================================*/
76/* C_H3D_UPDATE_SOLID_TENSOR */
77/*=================================================================*/
78
79void c_h3d_update_solid_tensor_(my_real *TT,int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS,
80 my_real *FUNC ,int *ID_ELEM,int *CPT_DATATYPE,
81 int *ITY_ELEM, int *NUMELQ , int *NUMELT , int *NUMELP , int *NUMELR, int *IS_WRITTEN,
82 my_real *FUNC_CORNER, int *IS_CORNER_DATA,int *ISOLNOD,int *MAXNOD)
83{
84 int i,j,element_nb_nodes;
85 int offset;
86 unsigned int NNOD = *MAXNOD;
87 H3D_ID elem_id;
88 H3D_ID comp_id;
89//
90 // initialize
91
92 try {
93 // create Subcase (Loadcase)
94 unsigned int max_sims = 10;
95 unsigned int sub_count = 1;
96 float elem_result[6] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f };
97 float* elem_result1 = nullptr;
98 elem_result1 = new float[6 * NNOD];
99
100 for( i = 0; i < 6*NNOD; i++ )
101 {
102 elem_result1[i] = 0.0f;
103 }
104
105 // create Result Data sets
106 unsigned int num_corners = 0;
107 unsigned int num_modes = 0;
108 bool complex = false;
109 float value[1] ;
110
111 sim_idx = *IH3D;
112
113 if(*IS_CORNER_DATA == 1)
114 {
115 if(*NUMELS != 0)
116 {
117 rc = Hyper3DDatasetBegin(h3d_file, *NUMELS, sim_idx, subcase_id, H3D_DS_ELEM,
118 H3D_DS_TENSOR3D, NNOD , num_modes, *CPT_DATATYPE,
119 0, solid_poolname_id, complex);
120 if( !rc ) throw rc;
121
122 for( i = 0; i < *NUMELS; i++ )
123 {
124 if( ITY_ELEM[i] == 1 && IS_WRITTEN[i] == 1)
125 {
126 elem_id = ID_ELEM[i];
127 element_nb_nodes = ISOLNOD[i];
128// 4 nodes & 10 nodes
129 if( ISOLNOD[i] == 4)
130 {
131 for( j = 0; j < 4; j++ )
132 {
133 elem_result1[j*6] = FUNC_CORNER[6*(NNOD*(i)+j)];
134 elem_result1[j*6+1] = FUNC_CORNER[6*(NNOD*(i)+j)+1];
135 elem_result1[j*6+2] = FUNC_CORNER[6*(NNOD*(i)+j)+2];
136 elem_result1[j*6+3] = FUNC_CORNER[6*(NNOD*(i)+j)+3];
137 elem_result1[j*6+4] = FUNC_CORNER[6*(NNOD*(i)+j)+4];
138 elem_result1[j*6+5] = FUNC_CORNER[6*(NNOD*(i)+j)+5];
139 }
140 }
141// 6 nodes
142 if( ISOLNOD[i] == 6)
143 {
144 for( j = 0; j < 6; j++ )
145 {
146 elem_result1[j*6] = FUNC_CORNER[6*(NNOD*(i)+j)];
147 elem_result1[j*6+1] = FUNC_CORNER[6*(NNOD*(i)+j)+1];
148 elem_result1[j*6+2] = FUNC_CORNER[6*(NNOD*(i)+j)+2];
149 elem_result1[j*6+3] = FUNC_CORNER[6*(NNOD*(i)+j)+3];
150 elem_result1[j*6+4] = FUNC_CORNER[6*(NNOD*(i)+j)+4];
151 elem_result1[j*6+5] = FUNC_CORNER[6*(NNOD*(i)+j)+5];
152 }
153 }
154// 8 nodes
155 if( ISOLNOD[i] == 8)
156 {
157 for( j = 0; j < 8; j++ )
158 {
159 elem_result1[j*6] = FUNC_CORNER[6*(NNOD*(i)+j)];
160 elem_result1[j*6+1] = FUNC_CORNER[6*(NNOD*(i)+j)+1];
161 elem_result1[j*6+2] = FUNC_CORNER[6*(NNOD*(i)+j)+2];
162 elem_result1[j*6+3] = FUNC_CORNER[6*(NNOD*(i)+j)+3];
163 elem_result1[j*6+4] = FUNC_CORNER[6*(NNOD*(i)+j)+4];
164 elem_result1[j*6+5] = FUNC_CORNER[6*(NNOD*(i)+j)+5];
165 }
166 }
167
168// 10 nodes
169 if( ISOLNOD[i] == 10)
170 {
171 for( j = 0; j < 10; j++ )
172 {
173 elem_result1[j*6] = FUNC_CORNER[6*(NNOD*(i)+j)];
174 elem_result1[j*6+1] = FUNC_CORNER[6*(NNOD*(i)+j)+1];
175 elem_result1[j*6+2] = FUNC_CORNER[6*(NNOD*(i)+j)+2];
176 elem_result1[j*6+3] = FUNC_CORNER[6*(NNOD*(i)+j)+3];
177 elem_result1[j*6+4] = FUNC_CORNER[6*(NNOD*(i)+j)+4];
178 elem_result1[j*6+5] = FUNC_CORNER[6*(NNOD*(i)+j)+5];
179 }
180 }
181// 16 nodes : not available for the moment
182/*
183 if( ISOLNOD[i] == 16)
184 {
185 for( j = 0; j < 16; j++ )
186 {
187 elem_result1[j*6] = FUNC_CORNER[6*(20*(i)+j)];
188 elem_result1[j*6+1] = FUNC_CORNER[6*(20*(i)+j)+1];
189 elem_result1[j*6+2] = FUNC_CORNER[6*(20*(i)+j)+2];
190 elem_result1[j*6+3] = FUNC_CORNER[6*(20*(i)+j)+3];
191 elem_result1[j*6+4] = FUNC_CORNER[6*(20*(i)+j)+4];
192 elem_result1[j*6+5] = FUNC_CORNER[6*(20*(i)+j)+5];
193 }
194
195 }
196*/
197// 20 nodes
198/*
199
200 if( ISOLNOD[i] == 20)
201 {
202 for( j = 0; j < 20; j++ )
203 {
204 elem_result1[j*6] = FUNC_CORNER[6*(20*(i)+j)];
205 elem_result1[j*6+1] = FUNC_CORNER[6*(20*(i)+j)+1];
206 elem_result1[j*6+2] = FUNC_CORNER[6*(20*(i)+j)+2];
207 elem_result1[j*6+3] = FUNC_CORNER[6*(20*(i)+j)+3];
208 elem_result1[j*6+4] = FUNC_CORNER[6*(20*(i)+j)+4];
209 elem_result1[j*6+5] = FUNC_CORNER[6*(20*(i)+j)+5];
210 }
211 }
212*/
213 rc = Hyper3DDatasetWrite(h3d_file, elem_id, &elem_result1[0]);
214 }
215 }
216 rc = Hyper3DDatasetEnd(h3d_file);
217 if( !rc ) throw rc;
218 }
219 }
220
221 if(*NUMELS != 0)
222 {
223 rc = Hyper3DDatasetBegin(h3d_file, *NUMELS, sim_idx, subcase_id, H3D_DS_ELEM,
224 H3D_DS_TENSOR3D, num_corners, num_modes, *CPT_DATATYPE,
225 0, solid_poolname_id, complex);
226 if( !rc ) throw rc;
227
228 for( i = 0; i < *NUMELS; i++ )
229 {
230 if( ITY_ELEM[i] == 1 && IS_WRITTEN[i] == 1)
231 {
232 elem_id = ID_ELEM[i];
233 elem_result[0] = FUNC[6*i];
234 elem_result[1] = FUNC[6*i+1];
235 elem_result[2] = FUNC[6*i+2];
236 elem_result[3] = FUNC[6*i+3];
237 elem_result[4] = FUNC[6*i+4];
238 elem_result[5] = FUNC[6*i+5];
239 rc = Hyper3DDatasetWrite(h3d_file, elem_id, &elem_result[0]);
240 IS_WRITTEN[i] = 0;
241 }
242 }
243 rc = Hyper3DDatasetEnd(h3d_file);
244 if( !rc ) throw rc;
245 }
246
247 // When done, deallocate the memory
248 delete[] elem_result1;
249 elem_result1 = nullptr; // Optional: set pointer to nullptr after deletion
250
251 } // end of try
252
253 catch(...) {
254 Hyper3DExportClearError(h3d_file);
255 }
256
257}
258
259void _FCALL C_H3D_UPDATE_SOLID_TENSOR(my_real *TT,int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS,
260 my_real *FUNC ,int *ID_ELEM,int *CPT_DATATYPE,
261 int *ITY_ELEM, int *NUMELQ , int *NUMELT , int *NUMELP , int *NUMELR, int *IS_WRITTEN,
262 my_real *FUNC_CORNER, int *IS_CORNER_DATA,int *ISOLNOD,int *MAXNOD)
263{c_h3d_update_solid_tensor_ (TT,IH3D,ITAB,NUMNOD,IXS,NIXS,NUMELS,IPARTS,FUNC,ID_ELEM,CPT_DATATYPE,ITY_ELEM,
264 NUMELQ,NUMELT,NUMELP,NUMELR,IS_WRITTEN,FUNC_CORNER,IS_CORNER_DATA,ISOLNOD, MAXNOD);}
265
266void c_h3d_update_solid_tensor__ (my_real *TT,int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS,
267 my_real *FUNC ,int *ID_ELEM,int *CPT_DATATYPE,
268 int *ITY_ELEM, int *NUMELQ , int *NUMELT , int *NUMELP , int *NUMELR, int *IS_WRITTEN,
269 my_real *FUNC_CORNER, int *IS_CORNER_DATA,int *ISOLNOD,int *MAXNOD)
270{c_h3d_update_solid_tensor_ (TT,IH3D,ITAB,NUMNOD,IXS,NIXS,NUMELS,IPARTS,FUNC,ID_ELEM,CPT_DATATYPE,ITY_ELEM,
271 NUMELQ,NUMELT,NUMELP,NUMELR,IS_WRITTEN,FUNC_CORNER,IS_CORNER_DATA,ISOLNOD,MAXNOD);}
272
273void c_h3d_update_solid_tensor (my_real *TT,int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS,
274 my_real *FUNC ,int *ID_ELEM,int *CPT_DATATYPE,
275 int *ITY_ELEM, int *NUMELQ , int *NUMELT , int *NUMELP , int *NUMELR, int *IS_WRITTEN,
276 my_real *FUNC_CORNER, int *IS_CORNER_DATA,int *ISOLNOD,int *MAXNOD)
277{c_h3d_update_solid_tensor_ (TT,IH3D,ITAB,NUMNOD,IXS,NIXS,NUMELS,IPARTS,FUNC,ID_ELEM,CPT_DATATYPE,ITY_ELEM,
278 NUMELQ,NUMELT,NUMELP,NUMELR,IS_WRITTEN,FUNC_CORNER,IS_CORNER_DATA,ISOLNOD,MAXNOD);}
279
280
281}
h3d_file
H3DFileInfo * h3d_file
Definition c_h3d_open_file.cpp:69
rc
bool rc
Definition c_h3d_open_file.cpp:70
comp_id
H3D_ID comp_id
Definition c_h3d_open_file.cpp:92
solid_poolname_id
H3D_ID solid_poolname_id
Definition c_h3d_open_file.cpp:118
sim_idx
H3D_SIM_IDX sim_idx
Definition c_h3d_open_file.cpp:95
subcase_id
H3D_ID subcase_id
Definition c_h3d_open_file.cpp:71
C_H3D_UPDATE_SOLID_TENSOR
void _FCALL C_H3D_UPDATE_SOLID_TENSOR(my_real *TT, int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS, my_real *FUNC, int *ID_ELEM, int *CPT_DATATYPE, int *ITY_ELEM, int *NUMELQ, int *NUMELT, int *NUMELP, int *NUMELR, int *IS_WRITTEN, my_real *FUNC_CORNER, int *IS_CORNER_DATA, int *ISOLNOD, int *MAXNOD)
Definition c_h3d_update_solid_tensor.cpp:259
c_h3d_update_solid_tensor
void c_h3d_update_solid_tensor(my_real *TT, int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS, my_real *FUNC, int *ID_ELEM, int *CPT_DATATYPE, int *ITY_ELEM, int *NUMELQ, int *NUMELT, int *NUMELP, int *NUMELR, int *IS_WRITTEN, my_real *FUNC_CORNER, int *IS_CORNER_DATA, int *ISOLNOD, int *MAXNOD)
Definition c_h3d_update_solid_tensor.cpp:273
c_h3d_update_solid_tensor__
void c_h3d_update_solid_tensor__(my_real *TT, int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS, my_real *FUNC, int *ID_ELEM, int *CPT_DATATYPE, int *ITY_ELEM, int *NUMELQ, int *NUMELT, int *NUMELP, int *NUMELR, int *IS_WRITTEN, my_real *FUNC_CORNER, int *IS_CORNER_DATA, int *ISOLNOD, int *MAXNOD)
Definition c_h3d_update_solid_tensor.cpp:266
c_h3d_update_solid_tensor_
void c_h3d_update_solid_tensor_(my_real *TT, int *IH3D, int *ITAB, int *NUMNOD, int *IXS, int *NIXS, int *NUMELS, int *IPARTS, my_real *FUNC, int *ID_ELEM, int *CPT_DATATYPE, int *ITY_ELEM, int *NUMELQ, int *NUMELT, int *NUMELP, int *NUMELR, int *IS_WRITTEN, my_real *FUNC_CORNER, int *IS_CORNER_DATA, int *ISOLNOD, int *MAXNOD)
Definition c_h3d_update_solid_tensor.cpp:79
FUNC
#define FUNC
Definition clamov.c:8
my_real
#define my_real
Definition cppsort.cpp:32
Hyper3DDatasetBegin
bool Hyper3DDatasetBegin(H3DFileInfo *h3d_file, unsigned int count, H3D_SIM_IDX idx, H3D_ID subcase_id, H3D_DS_TYPE type, H3D_DS_FORMAT format, unsigned int num_corners, unsigned int num_modes, H3D_ID dt_id, int layer_idx, H3D_ID data_poolname_id, bool complex)
Definition h3d_dl.c:1352
Hyper3DExportClearError
bool Hyper3DExportClearError(H3DFileInfo *h3d_file)
Definition h3d_dl.c:939
Hyper3DDatasetEnd
bool Hyper3DDatasetEnd(H3DFileInfo *h3d_file)
Definition h3d_dl.c:1400
Hyper3DDatasetWrite
bool Hyper3DDatasetWrite(H3DFileInfo *h3d_file, H3D_ID id, const float *data)
Definition h3d_dl.c:1372
_FCALL
#define _FCALL
Definition memory_use_c.c:25