<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Hi Lawrence and Matt,<br>
<br>
Many thanks for your help.<br>
<br>
Best regards,<br>
Rongliang<br>
<br>
On 06/08/2016 04:36 PM, Matthew Knepley wrote:<br>
</div>
<blockquote
cite="mid:CAMYG4GmYwn09QgtZXAAZcpU0wMZaCfNpVVz5+gdoRCth5jZ-ZA@mail.gmail.com"
type="cite">
<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
moz-do-not-send="true" 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
moz-do-not-send="true"
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 moz-do-not-send="true" href="mailto:rl.chen@siat.ac.cn" target="_blank">rl.chen@siat.ac.cn</a>
Phone: <a moz-do-not-send="true" 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
moz-do-not-send="true"
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
moz-do-not-send="true" 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>
</blockquote>
<br>
</body>
</html>