FrontISTR 5.2.0
Large-scale structural analysis program with finit element method
Loading...
Searching...
No Matches
hecmw_debug_write_dist.c
Go to the documentation of this file.
1/*****************************************************************************
2 * Copyright (c) 2019 FrontISTR Commons
3 * This software is released under the MIT License, see LICENSE.txt
4 *****************************************************************************/
5
6#include <stdio.h>
7#include "hecmw_struct.h"
8#include "hecmw_util.h"
10
11/*============================================================================*/
12/* */
13/* node information */
14/* */
15/*============================================================================*/
16/*----------------------------------------------------------------------------*/
17/* number of nodes < n_node > */
18/*----------------------------------------------------------------------------*/
19extern void HECMW_dbg_n_node_(struct hecmwST_local_mesh *local_mesh, char *file,
20 int line) {
21 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
22 local_mesh->n_node);
23}
24
25/*----------------------------------------------------------------------------*/
26/* number of nodes including extra overlap after refinement < n_node_gross > */
27/*----------------------------------------------------------------------------*/
28extern void HECMW_dbg_n_node_gross_(struct hecmwST_local_mesh *local_mesh,
29 char *file, int line) {
30 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node_gross = %d", file, line,
31 local_mesh->n_node_gross);
32}
33
34/*----------------------------------------------------------------------------*/
35/* number of nodes without extra slave nodes < nn_middle > */
36/*----------------------------------------------------------------------------*/
37extern void HECMW_dbg_nn_middle_(struct hecmwST_local_mesh *local_mesh,
38 char *file, int line) {
39 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: nn_middle = %d", file, line,
40 local_mesh->nn_middle);
41}
42
43/*----------------------------------------------------------------------------*/
44/* number of internal nodes < nn_internal > */
45/*----------------------------------------------------------------------------*/
46extern void HECMW_dbg_nn_internal_(struct hecmwST_local_mesh *local_mesh,
47 char *file, int line) {
48 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: nn_internal = %d", file, line,
49 local_mesh->nn_internal);
50}
51
52/*----------------------------------------------------------------------------*/
53/* maximal number of DOF < n_dof > */
54/*----------------------------------------------------------------------------*/
55extern void HECMW_dbg_n_dof_(struct hecmwST_local_mesh *local_mesh, char *file,
56 int line) {
57 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_dof = %d", file, line,
58 local_mesh->n_dof);
59}
60
61/*----------------------------------------------------------------------------*/
62/* number of DOF groups < n_dof_grp > */
63/*----------------------------------------------------------------------------*/
64extern void HECMW_dbg_n_dof_grp_(struct hecmwST_local_mesh *local_mesh,
65 char *file, int line) {
66 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_dof_grp = %d", file, line,
67 local_mesh->n_dof_grp);
68}
69
70/*----------------------------------------------------------------------------*/
71/* nodal coordinates < node > */
72/*----------------------------------------------------------------------------*/
73extern void HECMW_dbg_node_(struct hecmwST_local_mesh *local_mesh, char *file,
74 int line) {
75 int i;
76 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
77 local_mesh->n_node);
78
79 if (local_mesh->node) {
80 for (i = 0; i < local_mesh->n_node; i++) {
83 "%s:%d: i = %d, node[3*i] = %E, node[3*i+1] = %E, node[3*i+2] = %E",
84 file, line, i, local_mesh->node[3 * i], local_mesh->node[3 * i + 1],
85 local_mesh->node[3 * i + 2]);
86 }
87
88 } else {
89 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node = NULL");
90 }
91}
92
93/*----------------------------------------------------------------------------*/
94/* global node id < global_node_ID > */
95/*----------------------------------------------------------------------------*/
96extern void HECMW_dbg_global_node_id_(struct hecmwST_local_mesh *local_mesh,
97 char *file, int line) {
98 int i;
99 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
100 local_mesh->n_node);
101
102 if (local_mesh->global_node_ID) {
103 for (i = 0; i < local_mesh->n_node; i++) {
104 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, global_node_ID[i] = %d", file,
105 line, i, local_mesh->global_node_ID[i]);
106 }
107
108 } else {
109 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: global_node_ID = NULL", file, line);
110 }
111}
112
113/*----------------------------------------------------------------------------*/
114/* local node id & belonging domain of node < node_ID > */
115/*----------------------------------------------------------------------------*/
116extern void HECMW_dbg_node_id_(struct hecmwST_local_mesh *local_mesh,
117 char *file, int line) {
118 int i;
119 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
120 local_mesh->n_node);
121
122 if (local_mesh->node_ID) {
123 for (i = 0; i < local_mesh->n_node; i++) {
125 "%s:%d: i = %d, node_ID[2*i] = %d, node_ID[2*i+1] = %d", file,
126 line, i, local_mesh->node_ID[2 * i],
127 local_mesh->node_ID[2 * i + 1]);
128 }
129
130 } else {
131 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_ID = NULL", file, line);
132 }
133}
134
135/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
136/* local node id < node_ID[2*i] > */
137/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
138extern void HECMW_dbg_node_id_lid_(struct hecmwST_local_mesh *local_mesh,
139 char *file, int line) {
140 int i;
141 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
142 local_mesh->n_node);
143
144 if (local_mesh->node_ID) {
145 for (i = 0; i < local_mesh->n_node; i++) {
146 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, node_ID[2*i] = %d", file, line,
147 i, local_mesh->node_ID[2 * i]);
148 }
149
150 } else {
151 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_ID = NULL", file, line);
152 }
153}
154
155/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
156/* belonging domain of node < node_ID[2*i+1] > */
157/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
158extern void HECMW_dbg_node_id_domain_(struct hecmwST_local_mesh *local_mesh,
159 char *file, int line) {
160 int i;
161 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
162 local_mesh->n_node);
163
164 if (local_mesh->node_ID) {
165 for (i = 0; i < local_mesh->n_node; i++) {
166 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, node_ID[2*i+1] = %d", file,
167 line, i, local_mesh->node_ID[2 * i + 1]);
168 }
169
170 } else {
171 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_ID = NULL", file, line);
172 }
173}
174
175/*----------------------------------------------------------------------------*/
176/* node DOF group < node_dof_index, node_dof_item > */
177/*----------------------------------------------------------------------------*/
178extern void HECMW_dbg_node_dof_item_(struct hecmwST_local_mesh *local_mesh,
179 char *file, int line) {
180 int i;
181 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_dof_grp = %d", file, line,
182 local_mesh->n_dof_grp);
183
184 if (local_mesh->node_dof_index && local_mesh->node_dof_item) {
185 for (i = 0; i < local_mesh->n_dof_grp; i++) {
187 "%s:%d: i = %d, node_dof_index[i] = %d, node_dof_index[i+1] = "
188 "%d, node_dof_item[i] = %d",
189 file, line, i, local_mesh->node_dof_index[i],
190 local_mesh->node_dof_index[i + 1],
191 local_mesh->node_dof_item[i]);
192 }
193
194 } else {
195 if (local_mesh->node_dof_index == NULL) {
196 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_dof_index = NULL", file, line);
197 }
198
199 if (local_mesh->node_dof_item == NULL) {
200 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_dof_item = NULL", file, line);
201 }
202 }
203}
204
205/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
206/* node DOF group < node_dof_index > */
207/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
208extern void HECMW_dbg_node_dof_index_(struct hecmwST_local_mesh *local_mesh,
209 char *file, int line) {
210 int i;
211 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_dof_grp = %d", file, line,
212 local_mesh->n_dof_grp);
213
214 if (local_mesh->node_dof_index) {
215 for (i = 0; i < local_mesh->n_dof_grp; i++) {
216 HECMW_log(
218 "%s:%d: i = %d, node_dof_index[i] = %d, node_dof_index[i+1] = %d",
219 file, line, i, local_mesh->node_dof_index[i],
220 local_mesh->node_dof_index[i + 1]);
221 }
222
223 } else {
224 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_dof_index = NULL", file, line);
225 }
226}
227
228/*----------------------------------------------------------------------------*/
229/* initial condition of node < node_init_val_index, node_init_val_item > */
230/*----------------------------------------------------------------------------*/
232 char *file, int line) {
233 int i, j;
234 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
235 local_mesh->n_node);
236
237 if (local_mesh->node_init_val_index && local_mesh->node_init_val_item) {
238 for (i = 0; i < local_mesh->n_node; i++) {
240 "%s:%d: i = %d, node_init_val_index[i] = %d, "
241 "node_init_val_index[i+1] = %d",
242 file, line, i, local_mesh->node_init_val_index[i],
243 local_mesh->node_init_val_index[i + 1]);
244
245 for (j = local_mesh->node_init_val_index[i];
246 j < local_mesh->node_init_val_index[i + 1]; j++) {
247 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: j = %d, node_init_val_item[j] = %E",
248 file, line, j, local_mesh->node_init_val_item[j]);
249 }
250 }
251
252 } else {
253 if (local_mesh->node_init_val_index == NULL) {
254 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_init_val_index = NULL", file,
255 line);
256 }
257
258 if (local_mesh->node_init_val_item == NULL) {
259 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_init_val_item = NULL", file,
260 line);
261 }
262 }
263}
264
265/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
266/* initial condition of node < node_init_val_index > */
267/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
269 struct hecmwST_local_mesh *local_mesh, char *file, int line) {
270 int i;
271 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_node = %d", file, line,
272 local_mesh->n_node);
273
274 if (local_mesh->node_init_val_index) {
275 for (i = 0; i < local_mesh->n_node; i++) {
277 "%s:%d: i = %d, node_init_val_index[i] = %d, "
278 "node_init_val_index[i+1] = %d",
279 file, line, i, local_mesh->node_init_val_index[i],
280 local_mesh->node_init_val_index[i + 1]);
281 }
282
283 } else {
284 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_init_val_index = NULL", file, line);
285 }
286}
287
288/*============================================================================*/
289/* */
290/* element information */
291/* */
292/*============================================================================*/
293/*----------------------------------------------------------------------------*/
294/* number of elements < n_elem > */
295/*----------------------------------------------------------------------------*/
296extern void HECMW_dbg_n_elem_(struct hecmwST_local_mesh *local_mesh, char *file,
297 int line) {
298 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
299 local_mesh->n_elem);
300}
301
302/*----------------------------------------------------------------------------*/
303/* number of internal elements < ne_internal > */
304/*----------------------------------------------------------------------------*/
305extern void HECMW_dbg_ne_internal_(struct hecmwST_local_mesh *local_mesh,
306 char *file, int line) {
307 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: ne_internal = %d", file, line,
308 local_mesh->ne_internal);
309}
310
311/*----------------------------------------------------------------------------*/
312/* number of finite element types < n_elem_type > */
313/*----------------------------------------------------------------------------*/
314extern void HECMW_dbg_n_elem_type_(struct hecmwST_local_mesh *local_mesh,
315 char *file, int line) {
316 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem_type = %d", file, line,
317 local_mesh->n_elem_type);
318}
319
320/*----------------------------------------------------------------------------*/
321/* finite element type < elem_type_index, elem_type_item > */
322/*----------------------------------------------------------------------------*/
323extern void HECMW_dbg_elem_type_item_(struct hecmwST_local_mesh *local_mesh,
324 char *file, int line) {
325 int i;
326 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem_type = %d", file, line,
327 local_mesh->n_elem_type);
328
329 if (local_mesh->elem_type_index && local_mesh->elem_type_item) {
330 for (i = 0; i < local_mesh->n_elem_type; i++) {
332 "%s:%d: i = %d, elem_type_index[i] = %d, elem_type_index[i+1] "
333 "= %d: elem_type_item[i] = %d",
334 file, line, i, local_mesh->elem_type_index[i],
335 local_mesh->elem_type_index[i + 1],
336 local_mesh->elem_type_item[i]);
337 }
338
339 } else {
340 if (local_mesh->elem_type_index == NULL) {
341 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_type_index = NULL", file, line);
342 }
343
344 if (local_mesh->elem_type_item == NULL) {
345 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_type_item = NULL", file, line);
346 }
347 }
348}
349
350/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
351/* finite element type < elem_type_index > */
352/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
353extern void HECMW_dbg_elem_type_index_(struct hecmwST_local_mesh *local_mesh,
354 char *file, int line) {
355 int i;
356 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem_type = %d", file, line,
357 local_mesh->n_elem_type);
358
359 if (local_mesh->elem_type_index) {
360 for (i = 0; i < local_mesh->n_elem_type; i++) {
361 HECMW_log(
363 "%s:%d: i = %d, elem_type_index[i] = %d, elem_type_index[i+1] = %d",
364 file, line, i, local_mesh->elem_type_index[i],
365 local_mesh->elem_type_index[i + 1]);
366 }
367
368 } else {
369 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_type_index = NULL", file, line);
370 }
371}
372
373/*----------------------------------------------------------------------------*/
374/* finite element type < elem_type > */
375/*----------------------------------------------------------------------------*/
376extern void HECMW_dbg_elem_type_(struct hecmwST_local_mesh *local_mesh,
377 char *file, int line) {
378 int i;
379 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
380 local_mesh->n_elem);
381
382 if (local_mesh->elem_type) {
383 for (i = 0; i < local_mesh->n_elem; i++) {
384 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, elem_type[i] = %d", file, line,
385 i, local_mesh->elem_type[i]);
386 }
387
388 } else {
389 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_type = NULL", file, line);
390 }
391}
392
393/*----------------------------------------------------------------------------*/
394/* component node of element < elem_node_index, elem_node_item > */
395/*----------------------------------------------------------------------------*/
396extern void HECMW_dbg_elem_node_item_(struct hecmwST_local_mesh *local_mesh,
397 char *file, int line) {
398 int i, j;
399 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
400 local_mesh->n_elem);
401
402 if (local_mesh->elem_node_index && local_mesh->elem_node_item) {
403 for (i = 0; i < local_mesh->n_elem; i++) {
404 HECMW_log(
406 "%s:%d: i = %d, elem_node_index[i] = %d, elem_node_index[i+1] = %d",
407 file, line, i, local_mesh->elem_node_index[i],
408 local_mesh->elem_node_index[i + 1]);
409
410 for (j = local_mesh->elem_node_index[i];
411 j < local_mesh->elem_node_index[i + 1]; j++) {
412 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: j = %d, elem_node_item[j] = %d",
413 file, line, j, local_mesh->elem_node_item[j]);
414 }
415 }
416
417 } else {
418 if (local_mesh->elem_node_index == NULL) {
419 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_node_index = NULL", file, line);
420 }
421
422 if (local_mesh->elem_node_item == NULL) {
423 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_node_item = NULL", file, line);
424 }
425 }
426}
427
428/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
429/* component node of element < elem_node_index > */
430/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
431extern void HECMW_dbg_elem_node_index_(struct hecmwST_local_mesh *local_mesh,
432 char *file, int line) {
433 int i;
434 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
435 local_mesh->n_elem);
436
437 if (local_mesh->elem_node_index) {
438 for (i = 0; i < local_mesh->n_elem; i++) {
439 HECMW_log(
441 "%s:%d: i = %d, elem_node_index[i] = %d, elem_node_index[i+1] = %d",
442 file, line, i, local_mesh->elem_node_index[i],
443 local_mesh->elem_node_index[i + 1]);
444 }
445
446 } else {
447 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_node_index = NULL", file, line);
448 }
449}
450
451/*----------------------------------------------------------------------------*/
452/* local element id & belonging domain of element < elem_ID > */
453/*----------------------------------------------------------------------------*/
454extern void HECMW_dbg_elem_id_(struct hecmwST_local_mesh *local_mesh,
455 char *file, int line) {
456 int i;
457 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
458 local_mesh->n_elem);
459
460 if (local_mesh->elem_ID) {
461 for (i = 0; i < local_mesh->n_elem; i++) {
463 "%s:%d: i = %d, elem_ID[2*i] = %d, elem_ID[2*i+1] = %d", file,
464 line, i, local_mesh->elem_ID[2 * i],
465 local_mesh->elem_ID[2 * i + 1]);
466 }
467
468 } else {
469 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_ID = NULL", file, line);
470 }
471}
472
473/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
474/* local element id < elem_ID[2*i] > */
475/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
476extern void HECMW_dbg_elem_id_lid_(struct hecmwST_local_mesh *local_mesh,
477 char *file, int line) {
478 int i;
479 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
480 local_mesh->n_elem);
481
482 if (local_mesh->elem_ID) {
483 for (i = 0; i < local_mesh->n_elem; i++) {
484 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, elem_ID[2*i] = %d", file, line,
485 i, local_mesh->elem_ID[2 * i]);
486 }
487
488 } else {
489 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_ID = NULL", file, line);
490 }
491}
492
493/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
494/* belonging domain of element < elem_ID[2*i+1] > */
495/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
496extern void HECMW_dbg_elem_id_domain_(struct hecmwST_local_mesh *local_mesh,
497 char *file, int line) {
498 int i;
499 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
500 local_mesh->n_elem);
501
502 if (local_mesh->elem_ID) {
503 for (i = 0; i < local_mesh->n_elem; i++) {
504 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, elem_ID[2*i+1] = %d", file,
505 line, i, local_mesh->elem_ID[2 * i + 1]);
506 }
507
508 } else {
509 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_ID = NULL", file, line);
510 }
511}
512
513/*----------------------------------------------------------------------------*/
514/* global element id < global_element_ID > */
515/*----------------------------------------------------------------------------*/
516extern void HECMW_dbg_global_elem_id_(struct hecmwST_local_mesh *local_mesh,
517 char *file, int line) {
518 int i;
519 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
520 local_mesh->n_elem);
521
522 if (local_mesh->global_elem_ID) {
523 for (i = 0; i < local_mesh->n_elem; i++) {
524 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, global_elem_ID[i] = %d", file,
525 line, i, local_mesh->global_elem_ID[i]);
526 }
527
528 } else {
529 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: global_elem_ID = NULL", file, line);
530 }
531}
532
533/*----------------------------------------------------------------------------*/
534/* list of internal element < elem_internal_list > */
535/*----------------------------------------------------------------------------*/
537 char *file, int line) {
538 int i;
539 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: ne_internal = %d", file, line,
540 local_mesh->ne_internal);
541
542 if (local_mesh->elem_internal_list) {
543 for (i = 0; i < local_mesh->ne_internal; i++) {
544 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, elem_internal_list[i] = %d",
545 file, line, i, local_mesh->elem_internal_list[i]);
546 }
547
548 } else {
549 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_internal_list = NULL", file, line);
550 }
551}
552
553/*----------------------------------------------------------------------------*/
554/* section id < section_ID > */
555/*----------------------------------------------------------------------------*/
556extern void HECMW_dbg_section_id_(struct hecmwST_local_mesh *local_mesh,
557 char *file, int line) {
558 int i;
559 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
560 local_mesh->n_elem);
561
562 if (local_mesh->section_ID) {
563 for (i = 0; i < local_mesh->n_elem; i++) {
564 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, section_ID[i] = %d", file,
565 line, i, local_mesh->section_ID[i]);
566 }
567
568 } else {
569 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section_ID = NULL", file, line);
570 }
571}
572
573/*----------------------------------------------------------------------------*/
574/* number of material id < n_elem_mat_ID > */
575/*----------------------------------------------------------------------------*/
576extern void HECMW_dbg_n_elem_mat_id_(struct hecmwST_local_mesh *local_mesh,
577 char *file, int line) {
578 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem_mat_ID = %d", file, line,
579 local_mesh->n_elem_mat_ID);
580}
581
582/*----------------------------------------------------------------------------*/
583/* material id < elem_mat_ID_index, elem_mat_ID_item > */
584/*----------------------------------------------------------------------------*/
585extern void HECMW_dbg_elem_mat_id_item_(struct hecmwST_local_mesh *local_mesh,
586 char *file, int line) {
587 int i, j;
588 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
589 local_mesh->n_elem);
590
591 if (local_mesh->elem_mat_ID_index && local_mesh->elem_mat_ID_item) {
592 for (i = 0; i < local_mesh->n_elem; i++) {
594 "%s:%d: i = %d, elem_mat_ID_index[i] = %d, "
595 "elem_mat_ID_index[i+1] = %d",
596 file, line, i, local_mesh->elem_mat_ID_index[i],
597 local_mesh->elem_mat_ID_index[i + 1]);
598
599 for (j = local_mesh->elem_mat_ID_index[i];
600 j < local_mesh->elem_mat_ID_index[i + 1]; j++) {
601 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: j = %d, elem_mat_ID_item[j] = %d",
602 file, line, j, local_mesh->elem_mat_ID_item[j]);
603 }
604 }
605
606 } else {
607 if (local_mesh->elem_mat_ID_index == NULL) {
608 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_mat_ID_index = NULL", file, line);
609 }
610
611 if (local_mesh->elem_mat_ID_item == NULL) {
612 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_mat_ID_item = NULL", file, line);
613 }
614 }
615}
616
617/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
618/* material id < elem_mat_ID_index > */
619/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
621 char *file, int line) {
622 int i;
623 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_elem = %d", file, line,
624 local_mesh->n_elem);
625
626 if (local_mesh->elem_mat_ID_index) {
627 for (i = 0; i < local_mesh->n_elem; i++) {
629 "%s:%d: i = %d, elem_mat_ID_index[i] = %d, "
630 "elem_mat_ID_index[i+1] = %d",
631 file, line, i, local_mesh->elem_mat_ID_index[i],
632 local_mesh->elem_mat_ID_index[i + 1]);
633 }
634
635 } else {
636 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_mat_ID_index = NULL", file, line);
637 }
638}
639
640/*============================================================================*/
641/* */
642/* parallel & communication table information */
643/* */
644/*============================================================================*/
645/*----------------------------------------------------------------------------*/
646/* number of neighboring domains < n_neighbor_pe > */
647/*----------------------------------------------------------------------------*/
648extern void HECMW_dbg_n_neighbor_pe_(struct hecmwST_local_mesh *local_mesh,
649 char *file, int line) {
650 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
651 local_mesh->n_neighbor_pe);
652}
653
654/*----------------------------------------------------------------------------*/
655/* neighbor domain < neighbor_pe > */
656/*----------------------------------------------------------------------------*/
657extern void HECMW_dbg_neighbor_pe_(struct hecmwST_local_mesh *local_mesh,
658 char *file, int line) {
659 int i;
660 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
661 local_mesh->n_neighbor_pe);
662
663 if (local_mesh->neighbor_pe) {
664 for (i = 0; i < local_mesh->n_neighbor_pe; i++) {
665 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: i = %d, neighbor_pe[i] = %d", file,
666 line, i, local_mesh->neighbor_pe[i]);
667 }
668
669 } else {
670 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: neighbor_pe = NULL", file, line);
671 }
672}
673
674/*----------------------------------------------------------------------------*/
675/* import information < import_index, import_item > */
676/*----------------------------------------------------------------------------*/
677extern void HECMW_dbg_import_item_(struct hecmwST_local_mesh *local_mesh,
678 char *file, int line) {
679 int i, j;
680 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
681 local_mesh->n_neighbor_pe);
682
683 if (local_mesh->import_index && local_mesh->import_item) {
684 for (i = 0; i < local_mesh->n_neighbor_pe; i++) {
686 "%s:%d: i = %d, neighbor_pe[i] = %d, import_index[i] = %d, "
687 "import_index[i+1] = %d",
688 file, line, i, local_mesh->neighbor_pe[i],
689 local_mesh->import_index[i], local_mesh->import_index[i + 1]);
690
691 for (j = local_mesh->import_index[i]; j < local_mesh->import_index[i + 1];
692 j++) {
693 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: j = %d, import_item[j] = %d", file,
694 line, j, local_mesh->import_item[j]);
695 }
696 }
697
698 } else {
699 if (local_mesh->import_index == NULL) {
700 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: import_index = NULL", file, line);
701 }
702
703 if (local_mesh->import_item == NULL) {
704 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: import_item = NULL", file, line);
705 }
706 }
707}
708
709/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
710/* import information < import_index > */
711/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
712extern void HECMW_dbg_import_index_(struct hecmwST_local_mesh *local_mesh,
713 char *file, int line) {
714 int i;
715 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
716 local_mesh->n_neighbor_pe);
717
718 if (local_mesh->import_index) {
719 for (i = 0; i < local_mesh->n_neighbor_pe; i++) {
721 "%s:%d: i = %d, neighbor_pe[i] = %d, import_index[i] = %d, "
722 "import_index[i+1] = %d",
723 file, line, i, local_mesh->neighbor_pe[i],
724 local_mesh->import_index[i], local_mesh->import_index[i + 1]);
725 }
726
727 } else {
728 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: import_index = NULL", file, line);
729 }
730}
731
732/*----------------------------------------------------------------------------*/
733/* export information < export_index, export_item > */
734/*----------------------------------------------------------------------------*/
735extern void HECMW_dbg_export_item_(struct hecmwST_local_mesh *local_mesh,
736 char *file, int line) {
737 int i, j;
738 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
739 local_mesh->n_neighbor_pe);
740
741 if (local_mesh->export_index && local_mesh->export_item) {
742 for (i = 0; i < local_mesh->n_neighbor_pe; i++) {
744 "%s:%d: i = %d, neighbor_pe[i] = %d, export_index[i] = %d, "
745 "export_index[i+1] = %d",
746 file, line, i, local_mesh->neighbor_pe[i],
747 local_mesh->export_index[i], local_mesh->export_index[i + 1]);
748
749 for (j = local_mesh->export_index[i]; j < local_mesh->export_index[i + 1];
750 j++) {
751 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: j = %d, export_item[j] = %d", file,
752 line, j, local_mesh->export_item[j]);
753 }
754 }
755
756 } else {
757 if (local_mesh->export_index == NULL) {
758 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: export_index = NULL", file, line);
759 }
760
761 if (local_mesh->export_item == NULL) {
762 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: export_item = NULL", file, line);
763 }
764 }
765}
766
767/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
768/* export information < export_index > */
769/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
770extern void HECMW_dbg_export_index_(struct hecmwST_local_mesh *local_mesh,
771 char *file, int line) {
772 int i;
773 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
774 local_mesh->n_neighbor_pe);
775
776 if (local_mesh->export_index) {
777 for (i = 0; i < local_mesh->n_neighbor_pe; i++) {
779 "%s:%d: i = %d, neighbor_pe[i] = %d, export_index[i] = %d, "
780 "export_index[i+1] = %d",
781 file, line, i, local_mesh->neighbor_pe[i],
782 local_mesh->export_index[i], local_mesh->export_index[i + 1]);
783 }
784
785 } else {
786 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: export_index = NULL", file, line);
787 }
788}
789
790/*----------------------------------------------------------------------------*/
791/* shared information < shared_index, export_item > */
792/*----------------------------------------------------------------------------*/
793extern void HECMW_dbg_shared_item_(struct hecmwST_local_mesh *local_mesh,
794 char *file, int line) {
795 int i, j;
796 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
797 local_mesh->n_neighbor_pe);
798
799 if (local_mesh->shared_index && local_mesh->shared_item) {
800 for (i = 0; i < local_mesh->n_neighbor_pe; i++) {
802 "%s:%d: i = %d, neighbor_pe[i] = %d, shared_index[i] = %d, "
803 "shared_index[i+1] = %d",
804 file, line, i, local_mesh->neighbor_pe[i],
805 local_mesh->shared_index[i], local_mesh->shared_index[i + 1]);
806
807 for (j = local_mesh->shared_index[i]; j < local_mesh->shared_index[i + 1];
808 j++) {
809 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: j = %d, shared_item[j] = %d", file,
810 line, j, local_mesh->shared_item[j]);
811 }
812 }
813
814 } else {
815 if (local_mesh->shared_index == NULL) {
816 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: shared_index = NULL", file, line);
817 }
818
819 if (local_mesh->shared_item == NULL) {
820 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: shared_item = NULL", file, line);
821 }
822 }
823}
824
825/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
826/* shared information < shared_index > */
827/*- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
828extern void HECMW_dbg_shared_index_(struct hecmwST_local_mesh *local_mesh,
829 char *file, int line) {
830 int i;
831 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: n_neighbor_pe = %d", file, line,
832 local_mesh->n_neighbor_pe);
833
834 if (local_mesh->shared_index) {
835 for (i = 0; i < local_mesh->n_neighbor_pe; i++) {
837 "%s:%d: i = %d, neighbor_pe[i] = %d, shared_index[i] = %d, "
838 "shared_index[i+1] = %d",
839 file, line, i, local_mesh->neighbor_pe[i],
840 local_mesh->shared_index[i], local_mesh->shared_index[i + 1]);
841 }
842
843 } else {
844 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: shared_index = NULL", file, line);
845 }
846}
847
848/*============================================================================*/
849/* */
850/* section information */
851/* */
852/*============================================================================*/
853/*----------------------------------------------------------------------------*/
854/* number of sections < section->n_sect > */
855/*----------------------------------------------------------------------------*/
856extern void HECMW_dbg_n_sect_(struct hecmwST_local_mesh *local_mesh, char *file,
857 int line) {
858 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->n_sect = %d", file, line,
859 local_mesh->section->n_sect);
860}
861
862/*----------------------------------------------------------------------------*/
863/* section information */
864/* < section->n_sect, section->sect_type, section->sect_opt, */
865/* section->sect_mat_ID_index, section->sect_mat_ID_item, */
866/* section->sect_I_index, section->sect_I_item, */
867/* section->sect_R_index, section->sect_R_item > */
868/*----------------------------------------------------------------------------*/
869extern void HECMW_dbg_section_(struct hecmwST_local_mesh *local_mesh,
870 char *file, int line) {
871 int i, j;
872 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->n_sect = %d", file, line,
873 local_mesh->section->n_sect);
874
875 if (local_mesh->section->sect_type && local_mesh->section->sect_opt) {
876 for (i = 0; i < local_mesh->section->n_sect; i++) {
878 "%s:%d: i = %d, section->sect_type[i] = %d: "
879 "section->sect_opt[i] = %d",
880 file, line, i, local_mesh->section->sect_type[i],
881 local_mesh->section->sect_opt[i]);
882 }
883
884 } else {
885 if (local_mesh->section->sect_type == NULL) {
886 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_type = NULL", file,
887 line);
888 }
889
890 if (local_mesh->section->sect_opt == NULL) {
891 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_opt = NULL", file, line);
892 }
893 }
894
895 if (local_mesh->section->sect_mat_ID_index &&
896 local_mesh->section->sect_mat_ID_item) {
897 for (i = 0; i < local_mesh->section->n_sect; i++) {
899 "%s:%d: i = %d, section->sect_mat_ID_index[i] = %d, "
900 "section->sect_mat_ID_index[i+1] = %d",
901 file, line, i, local_mesh->section->sect_mat_ID_index[i],
902 local_mesh->section->sect_mat_ID_index[i + 1]);
903
904 for (j = local_mesh->section->sect_mat_ID_index[i];
905 j < local_mesh->section->sect_mat_ID_index[i + 1]; j++) {
907 "%s:%d: j = %d, section->sect_mat_ID_item[j] = %d", file,
908 line, j, local_mesh->section->sect_mat_ID_item[j]);
909 }
910 }
911
912 } else {
913 if (local_mesh->section->sect_mat_ID_index == NULL) {
914 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_mat_ID_index = NULL",
915 file, line);
916 }
917
918 if (local_mesh->section->sect_mat_ID_item == NULL) {
919 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_mat_ID_item = NULL",
920 file, line);
921 }
922 }
923
924 if (local_mesh->section->sect_I_index && local_mesh->section->sect_I_item) {
925 for (i = 0; i < local_mesh->section->n_sect; i++) {
927 "%s:%d: i = %d, section->sect_I_index[i] = %d, "
928 "section->sect_I_index[i+1] = %d",
929 file, line, i, local_mesh->section->sect_I_index[i],
930 local_mesh->section->sect_I_index[i + 1]);
931
932 for (j = local_mesh->section->sect_I_index[i];
933 j < local_mesh->section->sect_I_index[i + 1]; j++) {
935 "%s:%d: j = %d, section->sect_I_item[j] = %d", file, line, j,
936 local_mesh->section->sect_I_item[j]);
937 }
938 }
939
940 } else {
941 if (local_mesh->section->sect_I_index == NULL) {
942 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_I_index = NULL", file,
943 line);
944 }
945
946 if (local_mesh->section->sect_I_item == NULL) {
947 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_I_item = NULL", file,
948 line);
949 }
950 }
951
952 if (local_mesh->section->sect_R_index && local_mesh->section->sect_R_item) {
953 for (i = 0; i < local_mesh->section->n_sect; i++) {
955 "%s:%d: i = %d, section->sect_R_index[i] = %d, "
956 "section->sect_R_index[i+1] = %d",
957 file, line, i, local_mesh->section->sect_R_index[i],
958 local_mesh->section->sect_R_index[i + 1]);
959
960 for (j = local_mesh->section->sect_R_index[i];
961 j < local_mesh->section->sect_R_index[i + 1]; j++) {
963 "%s:%d: j = %d, section->sect_R_item[j] = %E", file, line, j,
964 local_mesh->section->sect_R_item[j]);
965 }
966 }
967
968 } else {
969 if (local_mesh->section->sect_R_index == NULL) {
970 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_R_index = NULL", file,
971 line);
972 }
973
974 if (local_mesh->section->sect_R_item == NULL) {
975 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: section->sect_R_item = NULL", file,
976 line);
977 }
978 }
979}
980
981/*============================================================================*/
982/* */
983/* material information */
984/* */
985/*============================================================================*/
986/*----------------------------------------------------------------------------*/
987/* number of materials < material->n_mat > */
988/*----------------------------------------------------------------------------*/
989extern void HECMW_dbg_n_mat_(struct hecmwST_local_mesh *local_mesh, char *file,
990 int line) {
991 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->n_mat = %d", file, line,
992 local_mesh->material->n_mat);
993}
994
995/*----------------------------------------------------------------------------*/
996/* material information */
997/* < material->n_mat, material->n_mat_item, material->n_mat_subitem, */
998/* material->n_mat_table, material->mat_name, material->mat_item_index, */
999/* material->mat_subitem_index, material->mat_table_index, */
1000/* material->mat_val, material->mat_temp > */
1001/*----------------------------------------------------------------------------*/
1002extern void HECMW_dbg_material_(struct hecmwST_local_mesh *local_mesh,
1003 char *file, int line) {
1004 int i, j;
1005 /* material->n_mat */
1006 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->n_mat = %d", file, line,
1007 local_mesh->material->n_mat);
1008
1009 /* material->mat_item_index */
1010 if (local_mesh->material->mat_item_index) {
1011 for (i = 0; i < local_mesh->material->n_mat; i++) {
1013 "%s:%d: i = %d, material->mat_item_index[i] = %d, "
1014 "material->mat_item_index[i+1] = %d",
1015 file, line, i, local_mesh->material->mat_item_index[i],
1016 local_mesh->material->mat_item_index[i + 1]);
1017 }
1018
1019 } else {
1020 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->mat_item_index = NULL", file,
1021 line);
1022 }
1023
1024 /* material->mat_name */
1025 if (local_mesh->material->mat_name) {
1026 for (i = 0; i < local_mesh->material->n_mat; i++) {
1028 "%s:%d: i = %d, material->mat_name[i] = \"%s\"", file, line, i,
1029 local_mesh->material->mat_name[i]);
1030 }
1031
1032 } else {
1033 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->mat_name = NULL", file, line);
1034 }
1035
1036 /* material->n_mat_item */
1037 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->n_mat_item = %d", file, line,
1038 local_mesh->material->n_mat_item);
1039
1040 /* material->mat_subitem_index */
1041 if (local_mesh->material->mat_subitem_index) {
1042 for (i = 0; i < local_mesh->material->n_mat_item; i++) {
1044 "%s:%d: i = %d, material->mat_subitem_index[i] = %d, "
1045 "material->mat_subitem_index[i+1] = %d",
1046 file, line, i, local_mesh->material->mat_subitem_index[i],
1047 local_mesh->material->mat_subitem_index[i + 1]);
1048 }
1049
1050 } else {
1051 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->mat_subitem_index = NULL",
1052 file, line);
1053 }
1054
1055 /* material->n_mat_subitem */
1056 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->n_mat_subitem = %d", file, line,
1057 local_mesh->material->n_mat_subitem);
1058
1059 /* material->mat_table_index, material->mat_val, material->mat_temp */
1060 if (local_mesh->material->mat_table_index && local_mesh->material->mat_val &&
1061 local_mesh->material->mat_temp) {
1062 for (i = 0; i < local_mesh->material->n_mat_subitem; i++) {
1064 "%s:%d: i = %d, material->mat_table_index[i] = %d, "
1065 "material->mat_table_index[i+1] = %d",
1066 file, line, i, local_mesh->material->mat_table_index[i],
1067 local_mesh->material->mat_table_index[i + 1]);
1068
1069 for (j = local_mesh->material->mat_table_index[i];
1070 j < local_mesh->material->mat_table_index[i + 1]; j++) {
1072 "%s:%d: j = %d, material->mat_val[j] = %E: "
1073 "material->mat_temp[j] = %E",
1074 file, line, j, local_mesh->material->mat_val[j],
1075 local_mesh->material->mat_temp[j]);
1076 }
1077 }
1078
1079 } else {
1080 if (local_mesh->material->mat_table_index == NULL) {
1081 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->mat_table_index = NULL",
1082 file, line);
1083 }
1084
1085 if (local_mesh->material->mat_val == NULL) {
1086 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->mat_val = NULL", file, line);
1087 }
1088
1089 if (local_mesh->material->mat_temp == NULL) {
1090 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->mat_temp = NULL", file,
1091 line);
1092 }
1093 }
1094
1095 /* material->n_mat_table */
1096 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: material->n_mat_table = %d", file, line,
1097 local_mesh->material->n_mat_table);
1098}
1099
1100/*============================================================================*/
1101/* */
1102/* MPC group information */
1103/* */
1104/*============================================================================*/
1105/*----------------------------------------------------------------------------*/
1106/* number of MPC groups < mpc->n_mpc > */
1107/*----------------------------------------------------------------------------*/
1108extern void HECMW_dbg_n_mpc_(struct hecmwST_local_mesh *local_mesh, char *file,
1109 int line) {
1110 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: mpc->n_mpc = %d", file, line,
1111 local_mesh->mpc->n_mpc);
1112}
1113
1114/*----------------------------------------------------------------------------*/
1115/* MPC group information */
1116/* < mpc->n_mpc, mpc->mpc_index, */
1117/* mpc->mpc_item, mpc->mpc_dof, mpc->mpc_val > */
1118/*----------------------------------------------------------------------------*/
1119extern void HECMW_dbg_mpc_(struct hecmwST_local_mesh *local_mesh, char *file,
1120 int line) {
1121 int i, j;
1122 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: mpc->n_mpc = %d", file, line,
1123 local_mesh->mpc->n_mpc);
1124
1125 if (local_mesh->mpc->mpc_index && local_mesh->mpc->mpc_item &&
1126 local_mesh->mpc->mpc_dof && local_mesh->mpc->mpc_val) {
1127 for (i = 0; i < local_mesh->mpc->n_mpc; i++) {
1128 HECMW_log(
1130 "%s:%d: i = %d, mpc->mpc_index[i] = %d, mpc->mpc_index[i+1] = %d",
1131 file, line, i, local_mesh->mpc->mpc_index[i],
1132 local_mesh->mpc->mpc_index[i + 1]);
1133
1134 for (j = local_mesh->mpc->mpc_index[i];
1135 j < local_mesh->mpc->mpc_index[i + 1]; j++) {
1137 "%s:%d: j = %d, mpc->mpc_item[j] = %d: mpc->mpc_dof[j] = %d: "
1138 "mpc->mpc_val[j] = %E",
1139 file, line, j, local_mesh->mpc->mpc_item[j],
1140 local_mesh->mpc->mpc_dof[j], local_mesh->mpc->mpc_val[j]);
1141 }
1142 }
1143
1144 } else {
1145 if (local_mesh->mpc->mpc_index == NULL) {
1146 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: mpc->mpc_index = NULL", file, line);
1147 }
1148
1149 if (local_mesh->mpc->mpc_item == NULL) {
1150 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: mpc->mpc_item = NULL", file, line);
1151 }
1152
1153 if (local_mesh->mpc->mpc_dof == NULL) {
1154 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: mpc->mpc_dof = NULL", file, line);
1155 }
1156
1157 if (local_mesh->mpc->mpc_val == NULL) {
1158 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: mpc->mpc_val = NULL", file, line);
1159 }
1160 }
1161}
1162
1163/*============================================================================*/
1164/* */
1165/* amplitude information */
1166/* */
1167/*============================================================================*/
1168/*----------------------------------------------------------------------------*/
1169/* number of amplitude groups < amp->n_amp > */
1170/*----------------------------------------------------------------------------*/
1171extern void HECMW_dbg_n_amp_(struct hecmwST_local_mesh *local_mesh, char *file,
1172 int line) {
1173 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->n_amp = %d", file, line,
1174 local_mesh->amp->n_amp);
1175}
1176
1177/*----------------------------------------------------------------------------*/
1178/* amplitude group information */
1179/* < amp->n_amp, amp->amp_index, */
1180/* amp->amp_type_definition, amp->amp_type_time, amp->amp_type_value, */
1181/* amp->amp_val, amp->amp_table > */
1182/*----------------------------------------------------------------------------*/
1183extern void HECMW_dbg_amp_(struct hecmwST_local_mesh *local_mesh, char *file,
1184 int line) {
1185 int i, j;
1186 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->n_amp = %d", file, line,
1187 local_mesh->amp->n_amp);
1188
1189 if (local_mesh->amp->amp_index && local_mesh->amp->amp_type_definition &&
1190 local_mesh->amp->amp_type_time && local_mesh->amp->amp_type_value &&
1191 local_mesh->amp->amp_val && local_mesh->amp->amp_table) {
1192 for (i = 0; i < local_mesh->amp->n_amp; i++) {
1193 HECMW_log(
1195 "%s:%d: i = %d, amp->amp_index[i] = %d, amp->amp_index[i+1] = %d",
1196 file, line, i, local_mesh->amp->amp_index[i],
1197 local_mesh->amp->amp_index[i + 1]);
1199 "%s:%d: i = %d, amp->amp_type_definition[i] = %d: "
1200 "amp->amp_type_time[i] = %d: amp->amp_type_value[i] = %d",
1201 file, line, i, local_mesh->amp->amp_type_definition[i],
1202 local_mesh->amp->amp_type_time[i],
1203 local_mesh->amp->amp_type_value[i]);
1204
1205 for (j = local_mesh->amp->amp_index[i];
1206 j < local_mesh->amp->amp_index[i + 1]; j++) {
1208 "%s:%d: j = %d, amp->amp_val[j] = %E, amp_table[j] = %E",
1209 file, line, j, local_mesh->amp->amp_val[j],
1210 local_mesh->amp->amp_table[j]);
1211 }
1212 }
1213
1214 } else {
1215 if (local_mesh->amp->amp_index == NULL) {
1216 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->amp_index = NULL", file, line);
1217 }
1218
1219 if (local_mesh->amp->amp_type_definition == NULL) {
1220 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->amp_type_definition = NULL", file,
1221 line);
1222 }
1223
1224 if (local_mesh->amp->amp_type_time == NULL) {
1225 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->amp_type_time = NULL", file,
1226 line);
1227 }
1228
1229 if (local_mesh->amp->amp_type_value == NULL) {
1230 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->amp_type_value = NULL", file,
1231 line);
1232 }
1233
1234 if (local_mesh->amp->amp_val == NULL) {
1235 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->amp_val = NULL", file, line);
1236 }
1237
1238 if (local_mesh->amp->amp_table == NULL) {
1239 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: amp->amp_table = NULL", file, line);
1240 }
1241 }
1242}
1243
1244/*============================================================================*/
1245/* */
1246/* node group information */
1247/* */
1248/*============================================================================*/
1249/*----------------------------------------------------------------------------*/
1250/* number of node groups < node_group->n_grp > */
1251/*----------------------------------------------------------------------------*/
1252extern void HECMW_dbg_n_node_grp_(struct hecmwST_local_mesh *local_mesh,
1253 char *file, int line) {
1254 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_group->n_grp = %d", file, line,
1255 local_mesh->node_group->n_grp);
1256}
1257
1258/*----------------------------------------------------------------------------*/
1259/* node group information */
1260/* < node_group->n_grp, node_group->grp_index, */
1261/* node_group->grp_name, node_group->grp_item > */
1262/*----------------------------------------------------------------------------*/
1263extern void HECMW_dbg_node_group_(struct hecmwST_local_mesh *local_mesh,
1264 char *file, int line) {
1265 int i, j;
1266 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_group->n_grp = %d", file, line,
1267 local_mesh->node_group->n_grp);
1268
1269 if (local_mesh->node_group->grp_index && local_mesh->node_group->grp_name &&
1270 local_mesh->node_group->grp_item) {
1271 for (i = 0; i < local_mesh->node_group->n_grp; i++) {
1273 "%s:%d: i = %d, node_group->grp_index[i] = %d, "
1274 "node_group->grp_index[i+1] = %d",
1275 file, line, i, local_mesh->node_group->grp_index[i],
1276 local_mesh->node_group->grp_index[i + 1]);
1278 "%s:%d: i = %d, node_group->grp_name[i] = \"%s\"", file, line,
1279 i, local_mesh->node_group->grp_name[i]);
1280
1281 for (j = local_mesh->node_group->grp_index[i];
1282 j < local_mesh->node_group->grp_index[i + 1]; j++) {
1284 "%s:%d: j = %d, node_group->grp_item[j] = %d", file, line, j,
1285 local_mesh->node_group->grp_item[j]);
1286 }
1287 }
1288
1289 } else {
1290 if (local_mesh->node_group->grp_index == NULL) {
1291 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_group->grp_index == NULL", file,
1292 line);
1293 }
1294
1295 if (local_mesh->node_group->grp_name == NULL) {
1296 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_group->grp_name == NULL", file,
1297 line);
1298 }
1299
1300 if (local_mesh->node_group->grp_item == NULL) {
1301 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: node_group->grp_item == NULL", file,
1302 line);
1303 }
1304 }
1305}
1306
1307/*============================================================================*/
1308/* */
1309/* element group information */
1310/* */
1311/*============================================================================*/
1312/*----------------------------------------------------------------------------*/
1313/* number of element groups < elem_group->n_grp > */
1314/*----------------------------------------------------------------------------*/
1315extern void HECMW_dbg_n_elem_grp_(struct hecmwST_local_mesh *local_mesh,
1316 char *file, int line) {
1317 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_group->n_grp = %d", file, line,
1318 local_mesh->elem_group->n_grp);
1319}
1320
1321/*----------------------------------------------------------------------------*/
1322/* element group information */
1323/* < elem_group->n_grp, elem_group->grp_index, */
1324/* elem_group->grp_name, elem_group->grp_item > */
1325/*----------------------------------------------------------------------------*/
1326extern void HECMW_dbg_elem_group_(struct hecmwST_local_mesh *local_mesh,
1327 char *file, int line) {
1328 int i, j;
1329 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_group->n_grp = %d", file, line,
1330 local_mesh->elem_group->n_grp);
1331
1332 if (local_mesh->elem_group->grp_index && local_mesh->elem_group->grp_name &&
1333 local_mesh->elem_group->grp_item) {
1334 for (i = 0; i < local_mesh->elem_group->n_grp; i++) {
1336 "%s:%d: i = %d, elem_group->grp_index[i] = %d, "
1337 "elem_group->grp_index[i+1] = %d",
1338 file, line, i, local_mesh->elem_group->grp_index[i],
1339 local_mesh->elem_group->grp_index[i + 1]);
1341 "%s:%d: i = %d, elem_group->grp_name[i] = \"%s\"", file, line,
1342 i, local_mesh->elem_group->grp_name[i]);
1343
1344 for (j = local_mesh->elem_group->grp_index[i];
1345 j < local_mesh->elem_group->grp_index[i + 1]; j++) {
1347 "%s:%d: j = %d, elem_group->grp_item[j] = %d", file, line, j,
1348 local_mesh->elem_group->grp_item[j]);
1349 }
1350 }
1351
1352 } else {
1353 if (local_mesh->elem_group->grp_index == NULL) {
1354 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_group->grp_index = NULL", file,
1355 line);
1356 }
1357
1358 if (local_mesh->elem_group->grp_name == NULL) {
1359 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_group->grp_name = NULL", file,
1360 line);
1361 }
1362
1363 if (local_mesh->elem_group->grp_item == NULL) {
1364 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: elem_group->grp_item = NULL", file,
1365 line);
1366 }
1367 }
1368}
1369
1370/*============================================================================*/
1371/* */
1372/* surface group information */
1373/* */
1374/*============================================================================*/
1375/*----------------------------------------------------------------------------*/
1376/* number of surface groups < surf_group->n_grp > */
1377/*----------------------------------------------------------------------------*/
1378extern void HECMW_dbg_n_surf_grp_(struct hecmwST_local_mesh *local_mesh,
1379 char *file, int line) {
1380 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: surf_group->n_grp = %d", file, line,
1381 local_mesh->surf_group->n_grp);
1382}
1383
1384/*----------------------------------------------------------------------------*/
1385/* surface group information */
1386/* < surf_group->n_grp, surf_group->grp_index, */
1387/* surf_group->grp_name, surf_group->grp_item > */
1388/*----------------------------------------------------------------------------*/
1389extern void HECMW_dbg_surf_group_(struct hecmwST_local_mesh *local_mesh,
1390 char *file, int line) {
1391 int i, j;
1392 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: surf_group->n_grp = %d", file, line,
1393 local_mesh->surf_group->n_grp);
1394
1395 if (local_mesh->surf_group->grp_index && local_mesh->surf_group->grp_name &&
1396 local_mesh->surf_group->grp_item) {
1397 for (i = 0; i < local_mesh->surf_group->n_grp; i++) {
1399 "%s:%d: i = %d, surf_group->grp_index[i] = %d, "
1400 "surf_group->grp_index[i+1] = %d",
1401 file, line, i, local_mesh->surf_group->grp_index[i],
1402 local_mesh->surf_group->grp_index[i + 1]);
1404 "%s:%d: i = %d, surf_group->grp_name[i] = \"%s\"", file, line,
1405 i, local_mesh->surf_group->grp_name[i]);
1406
1407 for (j = local_mesh->surf_group->grp_index[i];
1408 j < local_mesh->surf_group->grp_index[i + 1]; j++) {
1410 "%s:%d: j = %d, surf_group->grp_item[2*i] = %d, "
1411 "surf_group->grp_item[2*i+1] = %d",
1412 file, line, j, local_mesh->surf_group->grp_item[2 * j],
1413 local_mesh->surf_group->grp_item[2 * j + 1]);
1414 }
1415 }
1416
1417 } else {
1418 if (local_mesh->surf_group->grp_index == NULL) {
1419 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: surf_group->grp_index = NULL", file,
1420 line);
1421 }
1422
1423 if (local_mesh->surf_group->grp_name == NULL) {
1424 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: surf_group->grp_name = NULL", file,
1425 line);
1426 }
1427
1428 if (local_mesh->surf_group->grp_item == NULL) {
1429 HECMW_log(HECMW_LOG_DEBUG, "%s:%d: surf_group->grp_item = NULL", file,
1430 line);
1431 }
1432 }
1433}
void HECMW_dbg_section_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_neighbor_pe_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_surf_grp_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_import_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_init_val_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_node_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_export_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_amp_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_global_node_id_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_node_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_export_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_group_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_elem_mat_id_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_type_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_mat_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_elem_type_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_node_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_elem_grp_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_id_domain_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_surf_group_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_internal_list_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_mpc_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_init_val_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_id_lid_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_id_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_dof_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_id_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_type_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_id_domain_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_node_gross_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_dof_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_mat_id_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_material_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_type_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_elem_mat_id_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_shared_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_global_elem_id_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_sect_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_dof_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_dof_grp_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_node_grp_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_neighbor_pe_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_nn_internal_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_id_lid_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_section_id_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_shared_index_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_import_item_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_amp_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_ne_internal_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_nn_middle_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_n_elem_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_node_group_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
void HECMW_dbg_mpc_(struct hecmwST_local_mesh *local_mesh, char *file, int line)
#define NULL
int HECMW_log(int loglv, const char *fmt,...)
Definition: hecmw_log.c:260
#define HECMW_LOG_DEBUG
Definition: hecmw_log.h:21
int * amp_type_definition
Definition: hecmw_struct.h:61
double * amp_table
Definition: hecmw_struct.h:72
struct hecmwST_section * section
Definition: hecmw_struct.h:244
struct hecmwST_amplitude * amp
Definition: hecmw_struct.h:247
struct hecmwST_material * material
Definition: hecmw_struct.h:245
struct hecmwST_mpc * mpc
Definition: hecmw_struct.h:246
struct hecmwST_node_grp * node_group
Definition: hecmw_struct.h:248
double * node_init_val_item
Definition: hecmw_struct.h:180
struct hecmwST_surf_grp * surf_group
Definition: hecmw_struct.h:250
struct hecmwST_elem_grp * elem_group
Definition: hecmw_struct.h:249
int * mat_subitem_index
Definition: hecmw_struct.h:42
double * mat_val
Definition: hecmw_struct.h:44
double * mat_temp
Definition: hecmw_struct.h:45
int * mpc_dof
Definition: hecmw_struct.h:52
double * mpc_val
Definition: hecmw_struct.h:53
int * mpc_index
Definition: hecmw_struct.h:50
int * mpc_item
Definition: hecmw_struct.h:51
double * sect_R_item
Definition: hecmw_struct.h:32
int * sect_mat_ID_index
Definition: hecmw_struct.h:27
int * sect_mat_ID_item
Definition: hecmw_struct.h:28