<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Anton,<br>
I think that what you want can be changed in your code. Moab does not need to do anything for it.<br>
<br>
Let's take a simple case;<br>
You can have a 1d mesh, from node 0 to node N-1, along a line, with x varying from 0 to N-1, for each node<br>
You can repeat that , so you will have some nodes from -N to -1, and from N to 2N+1 (periodicity "N")<br>
<br>
<br>
In your code, you just need to assemble you matrices carefully; <br>
you will have only N degrees of freedom<br>
<br>
nodes at positions -N, 0, and N will be mapped (assembled) to dof 0<br>
at dof 1, you will have nodes at position -N+1, 1, N+1, and so on (kN+1)<br>
<br>
Your loop over nodes maybe start at position x=-1 (for one layer of ghosts), and end at N+1 , for one layer at the other end<br>
<br>
So , from your description, at dof 0, will correspond to both input and output flow.<br>
<br>
Of course, you will have to generalize for 3d, the same idea.<br>
<br>
Would that work for you?<br>
<br>
Iulian<br>
<br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF386061"><font color="#000000" face="Tahoma" size="2"><b>From:</b> kanaev@ibrae.ac.ru [kanaev@ibrae.ac.ru]<br>
<b>Sent:</b> Wednesday, April 09, 2014 4:18 PM<br>
<b>To:</b> Grindeanu, Iulian R.<br>
<b>Subject:</b> RE: PBC<br>
</font><br>
</div>
<div></div>
<div>
<p style="font-family:'Lucida Grande',Verdana,Arial,Helvetica,sans-serif"> For example, i've got the tube with inlet surface and outlet surface. 2d meshing of inlet and outlet surfaces is equal node to node. Having PBCs set on inlet and outlet surfaces we can
model an infinite flow by getting the data from outlet to the inlet of the tube. The most obvious way of implementing that is to make inlet and outlet surfaces shared, like a surface between two neighbouring processors and make a ghost layer of cells "before"
inlet surface (doubling layer of cells containing outlet surface) and a ghost layer of cells "after" outlet surface (doubling layer of cells containing inlet surface).
</p>
<p style="font-family:'Lucida Grande',Verdana,Arial,Helvetica,sans-serif">I've looked briefly through resolve_shared_entities function and saw that shared entities are resolved by finding nodes with the same GLOBAL_ID tag on the different processors. I've tried
to do the trick, did a simple brick geometry and manually changed nodes ids on the outlet to make them equal to node ids on the inlet, called resolve shared entities while loading mesh in parallel but it didn't resolve that vertices as shared))
</p>
<p style="font-family:'Lucida Grande',Verdana,Arial,Helvetica,sans-serif"> </p>
<p style="font-family:'Lucida Grande',Verdana,Arial,Helvetica,sans-serif">Can you think out some way of implementing that with what moab provides?
</p>
<br>
On Mon, 7 Apr 2014 03:21:24 +0000, "Grindeanu, Iulian R." wrote:<br>
<blockquote style="padding-left:5px; border-left-color:#1010ff; border-left-width:2px; border-left-style:solid; margin-left:5px; width:100%">
P {margin-top:0;margin-bottom:0;}
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi Anton,<br>
I was on vacation last week; I will look at your problems; I do not understand what you need yet, maybe after I read more carefully<br>
Iulian<br>
<div style="font-family:'Times New Roman'; color:#000000; font-size:16px">
<hr>
<div id="divRpF595248" style="direction:ltr"><font color="#000000" face="Tahoma" size="2"><strong>From:</strong> kanaev@ibrae.ac.ru [kanaev@ibrae.ac.ru]<br>
<strong>Sent:</strong> Wednesday, April 02, 2014 3:23 AM<br>
<strong>To:</strong> Grindeanu, Iulian R.<br>
<strong>Subject:</strong> PBC<br>
</font><br>
</div>
<div>
<p>Hello Iulian! </p>
<p>I've sent a mail about PBC to moab-dev and i saw that in <span style="background-color:#ebebeb; font-family:'Lucida Grande',Verdana,Arial,Helvetica,sans-serif; font-size:11px; font-weight:bold">
moab-dev Digest, Vol 65, Issue 32</span> </p>
<p>Did it go there automatically or someone answered? </p>
<p>I can see it like this: </p>
<p><span style="font-family:monospace; white-space:pre">Today's Topics:</span> </p>
<p><span style="font-family:monospace; white-space:pre">1. PBC (</span><a href="mailto:kanaev@ibrae.ac.ru" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">kanaev@ibrae.ac.ru</a><span style="font-family:monospace; white-space:pre">)
2. commit/MOAB: danwu: Minor changes to WriteNC.cpp (code formatting, remove some compiler warnings). Added a new test for writing variable gw (gauss weights). Running this test shows "Writer with name NC for file testgw.nc using extension nc was unsuccessful.
Using default writer WriteHDF5Parallel for file testgw.nc". Will update NC writer later to make it work on variable gw. (</span><a href="mailto:commits-noreply@bitbucket.org" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">commits-noreply@bitbucket.org</a><span style="font-family:monospace; white-space:pre">)
---------------------------------------------------------------------- Message: 1 Date: Tue, 01 Apr 2014 15:34:28 -0400 From: <</span><a href="mailto:kanaev@ibrae.ac.ru" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">kanaev@ibrae.ac.ru</a><span style="font-family:monospace; white-space:pre">>
To: MOAB dev <</span><a href="mailto:MOAB-dev@mcs.anl.gov" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">MOAB-dev@mcs.anl.gov</a><span style="font-family:monospace; white-space:pre">> Subject: [MOAB-dev] PBC Message-ID:
<</span><a href="mailto:bc4cbd0d4c5454e36c766953a8694287@mailx.ibrae.ac.ru" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">bc4cbd0d4c5454e36c766953a8694287@mailx.ibrae.ac.ru</a><span style="font-family:monospace; white-space:pre">>
Content-Type: text/plain; charset="utf-8" Hello MOAB-dev, What can be the way of implementing periodical boundary conditions using MOAB? As i can see, i need to make a connectivity table for opposing entities of e.g. inlet/outlet (which generally are on different
partitions) Thanks Anton -------------- next part -------------- An HTML attachment was scrubbed... URL: <</span><a href="http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20140401/ca7d50dc/attachment-0001.html" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20140401/ca7d50dc/attachment-0001.html</a><span style="font-family:monospace; white-space:pre">>
------------------------------ Message: 2 Date: Tue, 01 Apr 2014 21:11:38 -0000 From:
</span><a href="mailto:commits-noreply@bitbucket.org" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">commits-noreply@bitbucket.org</a><span style="font-family:monospace; white-space:pre"> To:
</span><a href="mailto:moab-dev@mcs.anl.gov" style="color:#0000cc; outline:none; font-family:monospace; white-space:pre" target="_blank">moab-dev@mcs.anl.gov</a><span style="font-family:monospace; white-space:pre"> Subject: [MOAB-dev] commit/MOAB: danwu: Minor
changes to WriteNC.cpp</span> </p>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</body>
</html>