<div dir="ltr">Hi Jeremy,<div><br></div><div style>To the best of my knowledge, partition interfaces in VTK need to abut.</div><div style><br></div><div style>Although I do understand the reasoning of what you are trying to achieve, from a finite difference code perspective, leaving a gap is problematic b/c it creates an implicit cell that could span two or more partitions/zones. For most post-processing/visualization operations, the algorithms need to loop over cells, e.g., consider contouring or interpolating the data over a line. Not having an explicit definition of all the cells in a partition makes this more complicated -- it would require building that information when reading in the partitioned grid, e.g., by growing the grid to the right and communicating the data (assuming we run in parallel).</div>
<div style><br></div><div style>More importantly, VTK structured grids can also store values at cell-centers. The output scheme that you describe would only work for node-centered data.</div><div style><br></div><div style>
So, to best support node-centered and cell-centered data, partition interfaces need to abut.</div><div style><br></div><div style>Hope this helps.</div><div style><br></div><div style>Best,</div><div style>George</div><div style>
<br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jun 27, 2013 at 4:04 PM, Kozdon, Jeremy (CIV) <span dir="ltr"><<a href="mailto:jekozdon@nps.edu" target="_blank">jekozdon@nps.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
I am trying to get vtk to properly work with a parallel finite difference code. What I would like to do is have each processor write the data they "own" to a separate file. For instance<br>
processor 0 would write local extends "0 9 0 10 0 0" to file file0.vts<br>
processor 1 would write local extends "10 19 0 10 0 0" to file file1.vts<br>
processor 2 would write local extends "20 29 0 10 0 0" to file file2.vts<br>
...<br>
I thought that when I wrote the pvts file to combine the data this would be ghostlevel="0" since there is no overlap of the data, but I get errors on a lack of extents. For instance in the above it would be for lack of extents "9 10 0 10 0 0" and "19 20 0 10 0 0", i.e., where the parallel boundary occurs.<br>
<br>
I am wondering if there is a better way to do what I am trying to accomplish / whether it is possible to do such formatting with vtk. I would prefer not having processors write duplicated data if possible.<br>
<br>
Thanks!<br>
Jeremy<br>
______________________________<u></u>_________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/<u></u>opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_<u></u>FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/<u></u>listinfo/vtkusers</a><br>
</blockquote></div><br></div>