<div dir="ltr">I believe this is a change in the API between Xdmf2 and 3.<div><br></div><div>In Xdmf2 the VTK level reader would break up a single regular grid automatically for you. In Xdmf3 I didn't implement that. </div><div><br></div><div>Workarounds include modifying your xdmf file to manually partition the underlying hdf5 grid, or reverting to xdmf2.</div><div><br></div><div>Others have come across this already before. This and the new boost dependency are the reasons why I have not yet deprecated xdmf2 and made xdmf3 the default in VTK and ParaView. I'ld love to fix this issue, but have not yet found the funding to devote the time for it.</div><div><br></div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature">David E DeMarle<br>Kitware, Inc.<br>R&D Engineer<br>21 Corporate Drive<br>Clifton Park, NY 12065-8662<br>Phone: 518-881-4909</div></div>
<br><div class="gmail_quote">On Thu, Jun 18, 2015 at 3:15 PM, Nico Schlömer <span dir="ltr"><<a href="mailto:nico.schloemer@gmail.com" target="_blank">nico.schloemer@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I've created a minimal working example that shows what's going wrong; find it attached.<div><br></div><div>I'd be glad about any hints here.</div><div><br></div><div>Cheers,</div><div>Nico<div><div class="h5"><br><br><div class="gmail_quote"><div dir="ltr">On Thu, Jun 18, 2015 at 8:37 PM Nico Schlömer <<a href="mailto:nico.schloemer@gmail.com" target="_blank">nico.schloemer@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">It's quite the same: All nodes are on all processes. (Or so it seems: `m->GetNumberOfPoints()` returns the number of all points on all nodes.)</div><div dir="ltr"><br><div><br></div><div>--Nico</div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Jun 17, 2015 at 6:00 PM Berk Geveci <<a href="mailto:berk.geveci@kitware.com" target="_blank">berk.geveci@kitware.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Nico,<div><br></div><div>What happens if you call</div><div><br></div><div>xreader->UpdateInformation();</div><div><br></div><div>before you call SetUpdateExtent() ?</div><div><br></div><div>Best,</div><div>-berk</div></div><div class="gmail_extra"><br><div class="gmail_quote"></div></div><div class="gmail_extra"><div class="gmail_quote">On Tue, Jun 16, 2015 at 9:06 AM, Nico Schlömer <span dir="ltr"><<a href="mailto:nico.schloemer@gmail.com" target="_blank">nico.schloemer@gmail.com</a>></span> wrote:<br></div></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">From [1], I find that parallel reading of XDMF files is as simple as<br>
```<br>
  vtkMPIController *contr = vtkMPIController::New();<br>
  contr->Initialize();<br>
  int proc = contr->GetLocalProcessId();<br>
  int numprocs = contr->GetNumberOfProcesses();<br>
  contr->Barrier();<br>
  vtkXdmf3Reader *xreader = vtkXdmf3Reader::New();<br>
  xreader->SetFileName("data/pacman.xmf");<br>
  contr->Barrier();<br>
  xreader->SetUpdateExtent(proc, numprocs, 0);<br>
  xreader->Update();<br>
```<br>
How can I check if this really worked?<br>
After retrieving the output<br>
```<br>
  vtkDataObject *rOutput = xreader->GetOutputDataObject(0);<br>
  TEUCHOS_ASSERT(rOutput != nullptr);<br>
  vtkUnstructuredGrid * m = dynamic_cast<vtkUnstructuredGrid*>(rOutput);<br>
  TEUCHOS_ASSERT(m != nullptr);<br>
<br>
  std::cout << m->GetNumberOfCells() << " " << m->GetNumberOfPoints()<br>
<< std::endl;<br>
```<br>
it seems that the all procs have all points and cells.<br>
<br>
Is there a way to only retrieve the locally owned cells/points?<br>
<br>
Cheers,<br>
Nico<br></blockquote></div></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">
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/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" rel="noreferrer" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=vtkusers" rel="noreferrer" target="_blank">http://markmail.org/search/?q=vtkusers</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/vtkusers" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/vtkusers</a><br>
</blockquote></div><br></div>
</blockquote></div></blockquote></div></div></div></div></div>
<br>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/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" rel="noreferrer" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=vtkusers" rel="noreferrer" target="_blank">http://markmail.org/search/?q=vtkusers</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/vtkusers" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/vtkusers</a><br>
<br></blockquote></div><br></div>