<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Jun 8, 2016 at 5:45 AM, Rongliang Chen <span dir="ltr"><<a href="mailto:rl.chen@siat.ac.cn" target="_blank">rl.chen@siat.ac.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div>Hi Matt,<br>
<br>
I tried your suggestions but I still can not get the refined
Vertex Sets. What I am doing looks like this:<br>
-----------------------------------------------------------------------<br>
if (!rank) {<br></div></div></blockquote><div><br></div><div>Put vertices in "Face Sets" instead</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="#FFFFFF" text="#000000"><div>
DMPlexSetLabelValue(dm, "Vertex Sets", xx, xx) </div></div></blockquote><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="#FFFFFF" text="#000000"><div> DMPlexSetLabelValue(dm, "Face Sets", xx, xx)<br></div></div></blockquote><div>// Now it has edges in it </div><div>DMPlexLabelComplete(dm, "Face Sets");</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="#FFFFFF" text="#000000"><div>
}<br>
<br>
DMPlexDistribute(dm, xx, xx)<br>
<br>
DMPlexGetStratumIS(dm, "Vertex Sets", xx, &vertexIS)<br>
ISGetLocalSize(vertexIS, &numVertexes1)<br>
<br>
DMRefine(dm, xx)<br>
<br>
DMPlexGetLabel(dm, "Face Sets", &facelabel)<br>
DMPlexLabelComplete(dm, facelabel)<br>
<br>
DMPlexGetStratumIS(dm, "Vertex Sets", xx, &vertexIS)<br>
ISGetLocalSize(vertexIS, &numVertexes2)<br>
-------------------------------------------------------------------<br>
<br>
I got "numVertexes1 == numVertexes2". But I expect to get
"numVertexes2 > numVertexes1". Did I miss something?<br></div></div></blockquote><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="#FFFFFF" text="#000000"><div>
Thanks,<br>
Rongliang<br>
<br>
<br>
<br>
On 06/08/2016 08:13 AM, Rongliang Chen wrote:<br>
</div>
<blockquote type="cite">
<div>Hi Matt,<br>
<br>
Thank you very much.<br>
<br>
Best,<br>
Rongliang<br>
<br>
<br>
On 06/08/2016 01:53 AM, Matthew Knepley wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>On Tue, Jun 7, 2016 at 5:18 PM, Rongliang Chen <span dir="ltr"><<a href="mailto:rl.chen@siat.ac.cn" target="_blank">rl.chen@siat.ac.cn</a>></span>
wrote:<br>
</div>
<div class="gmail_extra">
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div>Hi Matt,<br>
<br>
I know how to label vertex (ex_get_node_set) and
face (ex_get_side_set), but do not know how to label
the edge for a 3D mesh. Can you give me an example
for that?<br>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>You can label the faces, and then call
DMPlexLabelComplete().</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div> Thanks,<br>
Rongliang<br>
<pre cols="72">---------------------
Rongliang Chen, PhD
Associate Professor
Laboratory for Engineering and Scientific Computing
Shenzhen Institutes of Advanced Technology
Chinese Academy of Sciences
Address: 1068 Xueyuan Avenue, Shenzhen University Town, Shenzhen, Guangdong (518055), P. R. China
E-mail: <a href="mailto:rl.chen@siat.ac.cn" target="_blank">rl.chen@siat.ac.cn</a>
Phone: <a href="tel:%2B86-755-86392312" value="+8675586392312" target="_blank">+86-755-86392312</a></pre>
On 06/07/2016 10:06 PM, Matthew Knepley wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Tue, Jun 7, 2016 at
3:05 PM, Rongliang Chen <span dir="ltr"><<a href="mailto:rl.chen@siat.ac.cn" target="_blank">rl.chen@siat.ac.cn</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div>Hi Matt,<br>
<br>
Thanks for your reply. <br>
<br>
Yes, I only labels the vertex because I
am doing the finite element method where
I only need the vertex label. Do you
mean that I need to label the faces and
get the vertex label from the face
label?<br>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>You need to label vertices and edges.</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div> Thanks,<br>
Rongliang<br>
<br>
<br>
<br>
On 06/07/2016 09:51 PM, Matthew Knepley
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Tue, Jun
7, 2016 at 10:09 AM, Rongliang
Chen <span dir="ltr"><<a href="mailto:rl.chen@siat.ac.cn" target="_blank">rl.chen@siat.ac.cn</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear
All,<br>
<br>
I am using the DMRefine_Plex to
refine a mesh. In the coarse
mesh, it has a label "Vertex
Sets" to mark the boundary
nodes. I found that the
DMRefine_Plex did not refine
this label.<br>
<br>
For example, in the coarse mesh,
there are 100 nodes on the
boundary, that is the number of
nodes in the label is 100. In
the refined mesh (uniformly
refined), the number of nodes in
the label should be larger than
100, but I found that this
number is still 100, which means
the newly added nodes on the
boundary were not added to the
label.<br>
<br>
But I need a refined label to
set the boundary conditions.
Anyone can tell me how to get a
refined Vertex Sets label?<br>
</blockquote>
<div><br>
</div>
<div>All operations are
topological in Plex, not
geometric. When we refine, we
split edges (for instance) and
the new vertex inserted</div>
<div>inherits the labels from the
edge that was split. I am
guessing that you only labels
the vertices.</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> Best
regards,<br>
Rongliang<br>
<br>
</blockquote>
</div>
<br>
<br clear="all">
<span><font color="#888888"> <span><font color="#888888">
<div><br>
</div>
-- <br>
<div data-smartmail="gmail_signature">What
most experimenters take
for granted before they
begin their experiments is
infinitely more
interesting than any
results to which their
experiments lead.<br>
-- Norbert Wiener</div>
</font></span></font></span></div>
<span><font color="#888888"> </font></span></div>
<span><font color="#888888"> </font></span></blockquote>
<span><font color="#888888">
<br>
</font></span></div>
<span><font color="#888888">
</font></span></blockquote>
<span><font color="#888888"> </font></span></div>
<span><font color="#888888"> <br>
<br clear="all"><span class="HOEnZb"><font color="#888888">
<div><br>
</div>
-- <br>
<div data-smartmail="gmail_signature">What
most experimenters take for granted before
they begin their experiments is infinitely
more interesting than any results to which
their experiments lead.<br>
-- Norbert Wiener</div>
</font></span></font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
<br>
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div data-smartmail="gmail_signature">What most experimenters
take for granted before they begin their experiments is
infinitely more interesting than any results to which
their experiments lead.<br>
-- Norbert Wiener</div>
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
<br>
</font></span></blockquote>
<br>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>