<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Hi,<br>
<br>
I am encountering an issue with the GMSH reader which seems to ignore vertex sets. I am attaching a small c example, as well as two meshes, both contain cell, face, and vertex sets, but somehow the vertex sets are not read by the GMSH reader.
<br>
I am particularly confused since the plexgmsh.c seems to have code to handle vertex sets.
<br>
<br>
Regards,<br>
Blaise<br>
<br>
<br>
</div>
</span></font></div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText"><br>
<br>
For reference, I am attaching the .geo file for the GMSH file and a dump of the exodus version:<br>
bblaptop:Tests (dmplex)$ ncdump ../TestMeshes/SquareFaceSet.gen<br>
netcdf SquareFaceSet {<br>
dimensions:<br>
len_name = 256 ;<br>
time_step = UNLIMITED ; // (0 currently)<br>
num_dim = 2 ;<br>
num_nodes = 9 ;<br>
num_elem = 8 ;<br>
num_el_blk = 1 ;<br>
num_node_sets = 2 ;<br>
num_side_sets = 3 ;<br>
num_el_in_blk1 = 8 ;<br>
num_nod_per_el1 = 3 ;<br>
num_side_ss1 = 2 ;<br>
num_df_ss1 = 4 ;<br>
num_side_ss2 = 2 ;<br>
num_df_ss2 = 4 ;<br>
num_side_ss3 = 4 ;<br>
num_df_ss3 = 8 ;<br>
num_qa_rec = 1 ;<br>
four = 4 ;<br>
len_string = 33 ;<br>
num_nod_ns1 = 3 ;<br>
num_nod_ns2 = 3 ;<br>
variables:<br>
double time_whole(time_step) ;<br>
int eb_status(num_el_blk) ;<br>
int eb_prop1(num_el_blk) ;<br>
eb_prop1:name = "ID" ;<br>
int ns_status(num_node_sets) ;<br>
int ns_prop1(num_node_sets) ;<br>
ns_prop1:name = "ID" ;<br>
int ss_status(num_side_sets) ;<br>
int ss_prop1(num_side_sets) ;<br>
ss_prop1:name = "ID" ;<br>
double coordx(num_nodes) ;<br>
double coordy(num_nodes) ;<br>
char eb_names(num_el_blk, len_name) ;<br>
char ns_names(num_node_sets, len_name) ;<br>
char ss_names(num_side_sets, len_name) ;<br>
char coor_names(num_dim, len_name) ;<br>
int connect1(num_el_in_blk1, num_nod_per_el1) ;<br>
connect1:elem_type = "TRI3" ;<br>
int elem_ss1(num_side_ss1) ;<br>
int side_ss1(num_side_ss1) ;<br>
double dist_fact_ss1(num_df_ss1) ;<br>
int elem_ss2(num_side_ss2) ;<br>
int side_ss2(num_side_ss2) ;<br>
double dist_fact_ss2(num_df_ss2) ;<br>
int elem_ss3(num_side_ss3) ;<br>
int side_ss3(num_side_ss3) ;<br>
double dist_fact_ss3(num_df_ss3) ;<br>
char qa_records(num_qa_rec, four, len_string) ;<br>
int elem_map(num_elem) ;<br>
int elem_num_map(num_elem) ;<br>
int node_num_map(num_nodes) ;<br>
int node_ns1(num_nod_ns1) ;<br>
double dist_fact_ns1(num_nod_ns1) ;<br>
int node_ns2(num_nod_ns2) ;<br>
double dist_fact_ns2(num_nod_ns2) ;<br>
<br>
// global attributes:<br>
:api_version = 8.03f ;<br>
:version = 8.03f ;<br>
:floating_point_word_size = 8 ;<br>
:file_size = 1 ;<br>
:maximum_name_length = 32 ;<br>
:int64_status = 0 ;<br>
:title = "cubit(SquareFaceSet.gen): 07/22/2022: 16:01:07" ;<br>
data:<br>
<br>
eb_status = 1 ;<br>
<br>
eb_prop1 = 1 ;<br>
<br>
ns_status = 1, 1 ;<br>
<br>
ns_prop1 = 10, 20 ;<br>
<br>
ss_status = 1, 1, 1 ;<br>
<br>
ss_prop1 = 20, 21, 22 ;<br>
<br>
coordx = 0, 0.5, -0.166666666666667, 0, 0.5, -0.5, -0.5, -0.5, 0.5 ;<br>
<br>
coordy = -0.5, 0, 0, 0.5, 0.5, 0, 0.5, -0.5, -0.5 ;<br>
<br>
eb_names =<br>
"" ;<br>
<br>
ns_names =<br>
"",<br>
"" ;<br>
<br>
ss_names =<br>
"",<br>
"",<br>
"" ;<br>
<br>
coor_names =<br>
"x",<br>
"y" ;<br>
<br>
connect1 =<br>
1, 2, 3,<br>
2, 4, 3,<br>
5, 4, 2,<br>
6, 3, 7,<br>
1, 3, 8,<br>
9, 2, 1,<br>
3, 6, 8,<br>
3, 4, 7 ;<br>
<br>
elem_ss1 = 5, 6 ;<br>
<br>
side_ss1 = 3, 3 ;<br>
<br>
dist_fact_ss1 = 1, 1, 1, 1 ;<br>
<br>
elem_ss2 = 3, 8 ;<br>
<br>
side_ss2 = 1, 2 ;<br>
<br>
dist_fact_ss2 = 1, 1, 1, 1 ;<br>
<br>
elem_ss3 = 6, 3, 4, 7 ;<br>
<br>
side_ss3 = 1, 3, 3, 2 ;<br>
<br>
dist_fact_ss3 = 1, 1, 1, 1, 1, 1, 1, 1 ;<br>
<br>
qa_records =<br>
"CUBIT",<br>
"2021.11",<br>
"07/22/2022",<br>
"16:01:07" ;<br>
<br>
elem_map = 1, 2, 3, 4, 5, 6, 7, 8 ;<br>
<br>
elem_num_map = 1, 2, 3, 4, 5, 6, 7, 8 ;<br>
<br>
node_num_map = 8, 3, 9, 5, 2, 7, 4, 6, 1 ;<br>
<br>
node_ns1 = 5, 7, 4 ;<br>
<br>
dist_fact_ns1 = 1, 1, 1 ;<br>
<br>
node_ns2 = 8, 9, 1 ;<br>
<br>
dist_fact_ns2 = 1, 1, 1 ;<br>
}<br>
<br>
— <br>
Canada Research Chair in Mathematical and Computational Aspects of Solid Mechanics (Tier 1)<br>
Professor, Department of Mathematics & Statistics<br>
Hamilton Hall room 409A, McMaster University<br>
1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada <br>
<a href="https://www.math.mcmaster.ca/bourdin">https://www.math.mcmaster.ca/bourdin</a> | +1 (905) 525 9140 ext. 27243<br>
<br>
</div>
</span></font></div>
</body>
</html>