<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Mon, Oct 29, 2018 at 9:37 AM Fengwen Wang <<a href="mailto:fwan@mek.dtu.dk">fwan@mek.dtu.dk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div>


<div dir="ltr">
<div id="m_6998385602785042245x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif">
<p>Yes.  I would like to a rotational BCs(Twist bcs). </p></div></div></div></blockquote><div>You can do that with some work. I have implemented a Mobius strip (see DM_BOUNDARY_TWIST). However, </div><div>I have no idea how to put a function space down on a non-orientable manifold.</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><div dir="ltr"><div id="m_6998385602785042245x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif">
<p>Such a boundary is needed to reduce the computational domain.</p>
<p><br>
</p>
<p>Thanks lot</p>
<p><br>
</p>
<p>Best regards</p>
<p>Fengwen<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<div id="m_6998385602785042245x_Signature">
<div style="font-family:Tahoma;font-size:13px">
<div style="font-family:Tahoma;font-size:13px">
<div class="m_6998385602785042245x_BodyFragment"><font size="2">
<div class="m_6998385602785042245x_PlainText">____________________________________________________<br>
Senior Researcher<br>
Department of Mechanical Engineering, DTU<br>
Nils Koppels Allé       <br>
Building 404    <br>
2800  Kgs. Lyngby <br>
<a href="mailto:fwan@mek.dtu.dk" target="_blank">fwan@mek.dtu.dk</a></div>
</font></div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_6998385602785042245x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Lawrence Mitchell <<a href="mailto:wence@gmx.li" target="_blank">wence@gmx.li</a>><br>
<b>Sent:</b> Monday, October 29, 2018 1:17:34 PM<br>
<b>To:</b> Matt Knepley<br>
<b>Cc:</b> Fengwen Wang; PETSc; PETSc<br>
<b>Subject:</b> Re: [petsc-users] [petsc-maint] How to impose boundary conditions using DMDA</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt">
<div class="m_6998385602785042245PlainText"><br>
<br>
> On 29 Oct 2018, at 10:56, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br>
> <br>
> You can certainly map the right edge to the top edge (its topologically a square again), but that<br>
> mapping is not smooth, and I do not know how you would make a global basis for the approximation space.<br>
<br>
Yeah, I thought this was topologically just a cone, but that would have:<br>
<br>
u(x=1) = -v(y=1)<br>
<br>
and<br>
<br>
v(x=1) = u(y=1) (note matching sign).<br>
<br>
So if you wanted the v(x=1) = -u(y=1) there would be a twist somehow, and then a discontinuity in the mapping.<br>
<br>
Lawrence</div>
</span></font>
</div>

</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>