#include <stdio.h>
#include <stdlib.h>
#define STRINGIFY(x) #x
#define TOSTRING(x) STRINGIFY(x)
#define EXCHECK(funcall) \
do { \
int error = (funcall); \
printf("after %s, error = %d\n", TOSTRING(funcall), error); \
if (error != EX_NOERR) { \
fprintf(stderr, "Error calling %s\n", TOSTRING(funcall)); \
ex_close(exoid); \
exit(-1); \
} \
} while (0)
int main(int argc, char **argv)
{
int i, j, k, m;
int node_list[100], elem_list[100], side_list[100];
int num_glo_vars, num_nod_vars, num_ele_vars, num_nset_vars, num_sset_vars;
int *truth_tab, *nset_tab, *sset_tab;
int whole_time_step, num_time_steps;
int prop_array[2];
float *glob_var_vals, *nodal_var_vals, *elem_var_vals, *nset_var_vals, *sset_var_vals;
float time_value;
float dist_fact[100];
char *var_names[7];
int CPU_word_size = 0;
int IO_word_size = 4;
&CPU_word_size,
&IO_word_size);
printf("after ex_create for test.exo, exoid = %d\n", exoid);
printf(" cpu word size: %d io word size: %d\n", CPU_word_size, IO_word_size);
int num_dim = 3;
int num_nodes = 33;
int num_elem = 7;
int num_elem_blk = 7;
int num_node_sets = 2;
int num_side_sets = 5;
num_node_sets, num_side_sets));
int num_qa_rec = 2;
char *qa_record[2][4];
qa_record[0][0] = "TESTWT";
qa_record[0][1] = "testwt";
qa_record[0][2] = "07/07/93";
qa_record[0][3] = "15:41:33";
qa_record[1][0] = "";
qa_record[1][1] = " ";
qa_record[1][2] = "";
qa_record[1][3] = " ";
int num_info = 3;
char *info[3];
info[0] = "This is the first information record.";
info[1] = "";
info[2] = " ";
float x[100], y[100], z[100];
x[0] = 0.0;
y[0] = 0.0;
z[0] = 0.0;
x[1] = 1.0;
y[1] = 0.0;
z[1] = 0.0;
x[2] = 1.0;
y[2] = 1.0;
z[2] = 0.0;
x[3] = 0.0;
y[3] = 1.0;
z[3] = 0.0;
x[4] = 1.0;
y[4] = 0.0;
z[4] = 0.0;
x[5] = 2.0;
y[5] = 0.0;
z[5] = 0.0;
x[6] = 2.0;
y[6] = 1.0;
z[6] = 0.0;
x[7] = 1.0;
y[7] = 1.0;
z[7] = 0.0;
x[8] = 0.0;
y[8] = 0.0;
z[8] = 0.0;
x[9] = 10.0;
y[9] = 0.0;
z[9] = 0.0;
x[10] = 10.0;
y[10] = 0.0;
z[10] = -10.0;
x[11] = 1.0;
y[11] = 0.0;
z[11] = -10.0;
x[12] = 1.0;
y[12] = 10.0;
z[12] = 0.0;
x[13] = 10.0;
y[13] = 10.0;
z[13] = 0.0;
x[14] = 10.0;
y[14] = 10.0;
z[14] = -10.0;
x[15] = 1.0;
y[15] = 10.0;
z[15] = -10.0;
x[16] = 0.0;
y[16] = 0.0;
z[16] = 0.0;
x[17] = 1.0;
y[17] = 0.0;
z[17] = 5.0;
x[18] = 10.0;
y[18] = 0.0;
z[18] = 2.0;
x[19] = 7.0;
y[19] = 5.0;
z[19] = 3.0;
x[20] = 3.0;
y[20] = 0.0;
z[20] = 6.0;
x[21] = 6.0;
y[21] = 0.0;
z[21] = 0.0;
x[22] = 0.0;
y[22] = 0.0;
z[22] = 0.0;
x[23] = 3.0;
y[23] = 2.0;
z[23] = 6.0;
x[24] = 6.0;
y[24] = 2.0;
z[24] = 2.0;
x[25] = 0.0;
y[25] = 2.0;
z[25] = 0.0;
x[26] = 2.7;
y[26] = 1.7;
z[26] = 2.7;
x[27] = 6.0;
y[27] = 1.7;
z[27] = 3.3;
x[28] = 5.7;
y[28] = 1.7;
z[28] = 1.7;
x[29] = 3.7;
y[29] = 0.0;
z[29] = 2.3;
x[30] = 0.0;
y[30] = 0.0;
z[30] = 0.0;
x[31] = 10.0;
y[31] = 0.0;
z[31] = 0.0;
x[32] = 10.0;
y[32] = 10.0;
z[32] = 10.0;
char *coord_names[] = {"xcoor", "ycoor", "zcoor"};
int *elem_map = (int *)calloc(num_elem, sizeof(int));
for (i = 1; i <= num_elem; i++) {
elem_map[i - 1] = i;
}
free(elem_map);
int num_elem_in_block[] = {1, 1, 1, 1, 1, 1, 1};
int num_nodes_per_elem[] = {4,
4,
8,
4,
6,
8,
3};
int ebids[] = {10, 11, 12, 13, 14, 15, 16};
int nattr[] = {1, 1, 1, 1, 1, 1, 1};
char *eb_type[] = {"quad", "quad", "hex", "tetra", "wedge", "tetra", "tri"};
nattr, 0));
char *prop_names[] = {"MATL", "DENSITY"};
int *connect = (int *)calloc(8, sizeof(int));
connect[0] = 1;
connect[1] = 2;
connect[2] = 3;
connect[3] = 4;
connect[0] = 5;
connect[1] = 6;
connect[2] = 7;
connect[3] = 8;
connect[0] = 9;
connect[1] = 10;
connect[2] = 11;
connect[3] = 12;
connect[4] = 13;
connect[5] = 14;
connect[6] = 15;
connect[7] = 16;
connect[0] = 17;
connect[1] = 18;
connect[2] = 19;
connect[3] = 20;
connect[0] = 21;
connect[1] = 22;
connect[2] = 23;
connect[3] = 24;
connect[4] = 25;
connect[5] = 26;
connect[0] = 17;
connect[1] = 18;
connect[2] = 19;
connect[3] = 20;
connect[4] = 27;
connect[5] = 28;
connect[6] = 30;
connect[7] = 29;
connect[0] = 31;
connect[1] = 32;
connect[2] = 33;
free(connect);
float attrib[] = {3.14159};
attrib[0] = 6.14159;
int nsids[] = {20, 21};
int num_nodes_per_set[] = {5, 3};
{
int num_df_per_set[] = {5, 3};
}
node_list[0] = 10;
node_list[1] = 11;
node_list[2] = 12;
node_list[3] = 13;
node_list[4] = 14;
dist_fact[0] = 1.0;
dist_fact[1] = 2.0;
dist_fact[2] = 3.0;
dist_fact[3] = 4.0;
dist_fact[4] = 5.0;
node_list[0] = 20;
node_list[1] = 21;
node_list[2] = 22;
dist_fact[0] = 1.1;
dist_fact[1] = 2.1;
dist_fact[2] = 3.1;
prop_array[0] = 1000;
prop_array[1] = 2000;
int ssids[] = {30, 31, 32, 33, 34};
int num_elem_per_set[] = {2, 2, 7, 8, 10};
{
int num_df_per_set[] = {4, 4, 0, 0, 0};
}
elem_list[0] = 2;
elem_list[1] = 2;
side_list[0] = 4;
side_list[1] = 2;
dist_fact[0] = 30.0;
dist_fact[1] = 30.1;
dist_fact[2] = 30.2;
dist_fact[3] = 30.3;
elem_list[0] = 1;
elem_list[1] = 2;
side_list[0] = 2;
side_list[1] = 3;
dist_fact[0] = 31.0;
dist_fact[1] = 31.1;
dist_fact[2] = 31.2;
dist_fact[3] = 31.3;
elem_list[0] = 3;
elem_list[1] = 3;
elem_list[2] = 3;
elem_list[3] = 3;
elem_list[4] = 3;
elem_list[5] = 3;
elem_list[6] = 3;
side_list[0] = 5;
side_list[1] = 3;
side_list[2] = 3;
side_list[3] = 2;
side_list[4] = 4;
side_list[5] = 1;
side_list[6] = 6;
elem_list[0] = 4;
elem_list[1] = 4;
elem_list[2] = 4;
elem_list[3] = 4;
elem_list[4] = 6;
elem_list[5] = 6;
elem_list[6] = 6;
elem_list[7] = 6;
side_list[0] = 1;
side_list[1] = 2;
side_list[2] = 3;
side_list[3] = 4;
side_list[4] = 1;
side_list[5] = 2;
side_list[6] = 3;
side_list[7] = 4;
elem_list[0] = 5;
elem_list[1] = 5;
elem_list[2] = 5;
elem_list[3] = 5;
elem_list[4] = 5;
elem_list[5] = 7;
elem_list[6] = 7;
elem_list[7] = 7;
elem_list[8] = 7;
elem_list[9] = 7;
side_list[0] = 1;
side_list[1] = 2;
side_list[2] = 3;
side_list[3] = 4;
side_list[4] = 5;
side_list[5] = 1;
side_list[6] = 2;
side_list[7] = 3;
side_list[8] = 4;
side_list[9] = 5;
num_glo_vars = 1;
num_nod_vars = 2;
num_ele_vars = 3;
num_nset_vars = 4;
num_sset_vars = 7;
truth_tab = (int *)calloc((num_elem_blk * num_ele_vars), sizeof(int));
nset_tab = (int *)calloc((num_node_sets * num_nset_vars), sizeof(int));
sset_tab = (int *)calloc((num_side_sets * num_sset_vars), sizeof(int));
k = 0;
for (i = 0; i < num_elem_blk; i++) {
for (j = 0; j < num_ele_vars; j++) {
truth_tab[k++] = 1;
}
}
k = 0;
for (i = 0; i < num_node_sets; i++) {
for (j = 0; j < num_nset_vars; j++) {
if (k % 2 == 0) {
nset_tab[k++] = 1;
}
else {
nset_tab[k++] = 0;
}
}
}
k = 0;
for (i = 0; i < num_side_sets; i++) {
for (j = 0; j < num_sset_vars; j++) {
if (k % 2 == 0) {
sset_tab[k++] = 0;
}
else {
sset_tab[k++] = 1;
}
}
}
num_nset_vars, nset_tab, num_sset_vars, sset_tab));
free(truth_tab);
free(nset_tab);
free(sset_tab);
var_names[0] = "glo_vars";
var_names[0] = "node_variable_a_very_long_name_0";
var_names[1] = "nod_var1";
var_names[0] = "ele_var0";
var_names[1] = "ele_var1";
var_names[2] = "ele_var2";
var_names[0] = "nset_var0";
var_names[1] = "nset_var1";
var_names[2] = "nset_var2";
var_names[3] = "nset_var3";
var_names[0] = "sset_var0";
var_names[1] = "sset_var1";
var_names[2] = "sset_var2";
var_names[3] = "sset_var3";
var_names[4] = "sset_var4";
var_names[5] = "sset_var5";
var_names[6] = "sset_var6";
whole_time_step = 1;
num_time_steps = 10;
glob_var_vals = (float *)calloc(num_glo_vars, CPU_word_size);
nodal_var_vals = (float *)calloc(num_nodes, CPU_word_size);
elem_var_vals = (float *)calloc(4, CPU_word_size);
nset_var_vals = (float *)calloc(5, CPU_word_size);
sset_var_vals = (float *)calloc(10, CPU_word_size);
for (i = 0; i < num_time_steps; i++) {
time_value = (float)(i + 1) / 100.;
for (j = 0; j < num_glo_vars; j++) {
glob_var_vals[j] = (float)(j + 2) * time_value;
}
for (k = 1; k <= num_nod_vars; k++) {
for (j = 0; j < num_nodes; j++) {
nodal_var_vals[j] = (float)k + ((float)(j + 1) * time_value);
}
}
for (k = 1; k <= num_ele_vars; k++) {
for (j = 0; j < num_elem_blk; j++) {
for (m = 0; m < num_elem_in_block[j]; m++) {
elem_var_vals[m] = (float)(k + 1) + (float)(j + 2) + ((float)(m + 1) * time_value);
}
elem_var_vals));
}
}
int kk = 0;
for (j = 0; j < num_node_sets; j++) {
for (k = 0; k < num_nset_vars; k++) {
if (kk++ % 2 == 0) {
for (m = 0; m < num_nodes_per_set[j]; m++) {
nset_var_vals[m] = (float)(k + 1) + (float)(j + 2) + ((float)(m + 1) * time_value);
}
num_nodes_per_set[j], nset_var_vals));
}
}
}
kk = 0;
for (j = 0; j < num_side_sets; j++) {
for (k = 0; k < num_sset_vars; k++) {
if (kk++ % 2 != 0) {
for (m = 0; m < num_elem_per_set[j]; m++) {
sset_var_vals[m] = (float)(k + 1) + (float)(j + 2) + ((float)(m + 1) * time_value);
}
num_elem_per_set[j], sset_var_vals));
}
}
}
whole_time_step++;
}
free(glob_var_vals);
free(nodal_var_vals);
free(elem_var_vals);
free(nset_var_vals);
free(sset_var_vals);
return 0;
}
#define EXCHECK(funcall)
Definition ex_copy.c:15
@ EX_SIDE_SET
Definition exodusII.h:278
@ EX_NODE_SET
Definition exodusII.h:270
@ EX_NODAL
Definition exodusII.h:268
@ EX_GLOBAL
Definition exodusII.h:285
@ EX_ELEM_MAP
Definition exodusII.h:280
@ EX_ELEM_BLOCK
Definition exodusII.h:275
@ EX_ABORT
Definition exodusII.h:392
@ EX_VERBOSE
Definition exodusII.h:390
#define EX_CLOBBER
Definition exodusII.h:101
int ex_put_conn(int exoid, ex_entity_type blk_type, ex_entity_id blk_id, const void_int *node_conn, const void_int *elem_edge_conn, const void_int *elem_face_conn)
Definition ex_put_conn.c:44
int ex_put_prop_names(int exoid, ex_entity_type obj_type, int num_props, char **prop_names)
Definition ex_put_prop_names.c:97
int ex_put_coord_names(int exoid, char *const coord_names[])
Definition ex_put_coord_names.c:46
int ex_put_coord(int exoid, const void *x_coor, const void *y_coor, const void *z_coor)
Definition ex_put_coord.c:87
int ex_put_prop(int exoid, ex_entity_type obj_type, ex_entity_id obj_id, const char *prop_name, ex_entity_id value)
Definition ex_put_prop.c:77
int ex_put_concat_elem_block(int exoid, const void_int *elem_blk_id, char *const elem_type[], const void_int *num_elem_this_blk, const void_int *num_nodes_per_elem, const void_int *num_attr_this_blk, int define_maps)
Definition ex_put_concat_elem_block.c:36
int ex_put_attr(int exoid, ex_entity_type blk_type, ex_entity_id blk_id, const void *attrib)
Definition ex_put_attr.c:37
int ex_put_prop_array(int exoid, ex_entity_type obj_type, const char *prop_name, const void_int *values)
Definition ex_put_prop_array.c:71
int ex_put_init(int exoid, const char *title, int64_t num_dim, int64_t num_nodes, int64_t num_elem, int64_t num_elem_blk, int64_t num_node_sets, int64_t num_side_sets)
Definition ex_put_init.c:53
int ex_put_set_dist_fact(int exoid, ex_entity_type set_type, ex_entity_id set_id, const void *set_dist_fact)
Definition ex_put_set_dist_fact.c:37
int ex_put_set(int exoid, ex_entity_type set_type, ex_entity_id set_id, const void_int *set_entry_list, const void_int *set_extra_list)
Definition ex_put_set.c:41
int ex_put_concat_sets(int exoid, ex_entity_type set_type, const struct ex_set_specs *set_specs)
Definition ex_put_concat_sets.c:37
int ex_put_id_map(int exoid, ex_entity_type map_type, const void_int *map)
Definition ex_put_id_map.c:37
int ex_put_time(int exoid, int time_step, const void *time_value)
Definition ex_put_time.c:51
int ex_put_variable_names(int exoid, ex_entity_type obj_type, int num_vars, char *const var_names[])
Definition ex_put_variable_names.c:125
int ex_put_all_var_param(int exoid, int num_g, int num_n, int num_e, const int *elem_var_tab, int num_m, const int *nset_var_tab, int num_s, const int *sset_var_tab)
Definition ex_put_all_var_param.c:45
int ex_put_var(int exoid, int time_step, ex_entity_type var_type, int var_index, ex_entity_id obj_id, int64_t num_entries_this_obj, const void *var_vals)
Definition ex_put_var.c:77
#define ex_create(path, mode, comp_ws, io_ws)
Definition exodusII.h:591
int ex_put_qa(int exoid, int num_qa_records, char *qa_record[][4])
Definition ex_put_qa.c:63
int ex_opts(int options)
Definition ex_opts.c:56
int ex_put_info(int exoid, int num_info, char *const info[])
Definition ex_put_info.c:73
int ex_close(int exoid)
Definition ex_close.c:47
int ex_update(int exoid)
Definition ex_update.c:30
Definition exodusII.h:537
void_int * sets_dist_index
Definition exodusII.h:542
void_int * sets_ids
Definition exodusII.h:538
void_int * num_dist_per_set
Definition exodusII.h:540
void * sets_dist_fact
Definition exodusII.h:545
void_int * sets_extra_list
Definition exodusII.h:544
void_int * num_entries_per_set
Definition exodusII.h:539
void_int * sets_entry_list
Definition exodusII.h:543
void_int * sets_entry_index
Definition exodusII.h:541