<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<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">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>