<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><blockquote type="cite" class=""><div dir="ltr" class=""><div class="gmail_quote"><div class=""> I think the real question is, how do layers interact.</div></div></div></blockquote><div class="">The main source of interaction between cells is the energy between directly connected nearest-neighbor atom pairs, so crudely stiffness in some form. The layers are set up such each atom in one layer is linked to a unique atom in the layer above or below (basically one layer sits directly on top of the other, except for a possible xy shift, communication-wise it is identical to in-plane connections though). </div><div class=""><br class=""></div><div class="">The wrinkle with this however, is that I need to extend this beyond nearest-neighbor where I can set an interaction radius, where any cell connected to the current cell (not necessarily directly, can be through multiple jumps) within that radius also interacts with my current cell.</div><div class=""><br class=""></div><div class="">It seems a combination of DMNetwork + DMPlex seems the most fitting here then, since it seems that DMNetwork is built on plex.</div><div class=""><br class=""></div><div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div>Best regards,<br class=""><br class="">Jacob Faibussowitsch<br class="">(Jacob Fai - booss - oh - vitch)<br class="">Cell: (312) 694-3391</div></div>
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">On Jan 28, 2020, at 3:07 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div dir="ltr" class="">On Tue, Jan 28, 2020 at 2:32 PM Jacob Faibussowitsch <<a href="mailto:jacob.fai@gmail.com" class="">jacob.fai@gmail.com</a>> wrote:<br class=""></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word" class="">Hello All,<div class=""><br class=""></div><div class="">I am currently working on rewriting an existing code into PETSc. This particular code performs tight binding and works on a tiled hexagonal carbon lattice as a mesh. I am wondering what the best DM format is to represent this? Some potentially relevant info:</div><div class=""><br class=""></div><div class="">1. General setup looks like this (<a href="https://www.iconfinder.com/icons/3417099/carbon_structure_geometric_pattern_graphene_technology_hexagonal_lattice_molecular_structure_icon" target="_blank" class="">https://www.iconfinder.com/icons/3417099/carbon_structure_geometric_pattern_graphene_technology_hexagonal_lattice_molecular_structure_icon</a>), initially in 2D, however it should be extensible to 3D in the Z direction by layering these configurations on top of one another. They are not necessarily layered nicely, i.e. having one layer directly above the other, but should be able to handle offsets in the X or Y directions.</div><div class="">2. A “Cell” is a pair of atoms, so in the above link 2 connected circles represents 1 cell.</div><div class="">3. Connectivity between cells is always structured.</div><div class="">4. The mesh should only be interpolated (in the DMPlex sense) once, i.e. I need the 2 atom “cell" to store some geometric information about its internal bond.</div><div class="">5. Computational stencil should always be pointing “outwards”, i.e. it should not loop back on itself.</div><div class="">6. The computational kernel is a band-diagonal Hamiltonian matrix generated from the tight binding wave function (function of vectors between atoms, some known constants, and stencil size)</div><div class=""><br class=""></div><div class="">Granted some of these are not necessarily things that PETSc should facilitate on its own (most likely point #4 and #5) but there may be DM implementations that make writing this by hand easier than others. From looking into this some potential contenders I identified were DMNetwork, (a modified) DMDA, DMStag, DMPlex, or DMMOAB, but I welcome any and all input here!</div></div></blockquote><div class=""><br class=""></div><div class="">I don't think you can do this with DMDA/DMStag since they are based on tensor-product meshes, and I don't</div><div class="">see how to phrase this in that way. You can do this in DMPlex (and probably MOAB, I don't know) but since you</div><div class="">only really use 1D information, DMNetwork might be easier. I think the real question is, how do layers interact. One</div><div class="">layer is probably best as a Network. Multiple layers might still be better, or might need something more general.</div><div class=""><br class=""></div><div class=""> Thanks,</div><div class=""><br class=""></div><div class=""> Matt</div><div class=""> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word" class=""><div class=""><div class="">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; word-wrap: break-word;" class=""><div class="">Best regards,<br class=""><br class="">Jacob Faibussowitsch<br class="">(Jacob Fai - booss - oh - vitch)<br class="">Cell: (312) 694-3391</div></div>
</div>
<br class=""></div></div></blockquote></div><br clear="all" class=""><div class=""><br class=""></div>-- <br class=""><div dir="ltr" class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">-- Norbert Wiener</div><div class=""><br class=""></div><div class=""><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank" class="">https://www.cse.buffalo.edu/~knepley/</a><br class=""></div></div></div></div></div></div></div></div>
</div></blockquote></div><br class=""></body></html>