<div><div dir="auto">Andy,</div><div dir="auto"><br></div><div dir="auto">My own code. I’ve always had those Piece lines in my src from day one, and always used it via mpi, so don’t need added. All procs call all lines in my file writer. There is no if-master-section.  Literally just a vtk version number change that causes the problem. Can’t get the data files or the cut down source until tomorrow. </div><div dir="auto"><br></div><div dir="auto">I shared the other two links as those in addition to this post are the only places I can find that reference such a major change to the library that is undocumented and causes such a significant change to user code. Are there any further references on this change that I can read about to try to resolve this?</div><div dir="auto"><br></div><div dir="auto">I find it strange that this has not come up more. It would mean all users of 7.1 without the additional vtkmpicontroller lines would fail when calling a parallel writer from a parallel mpi-user-code....</div><div dir="auto"><br></div><div dir="auto">Cheers,</div><div dir="auto">Andy</div><br><div class="gmail_quote"><div>On Mon, 5 Feb 2018 at 18:07, Andy Bauer <<a href="mailto:andy.bauer@kitware.com">andy.bauer@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><div><div><div>Are you using libMesh or another simulation code? You'll definitely need the following to tell each process what piece it has:<br>pwriter->SetNumberOfPieces(nprocs);                                                                   <br>pwriter->SetStartPiece(rank);                                                                                                     <br>pwriter->SetEndPiece(rank);                                                                                                       <br>pwriter->Write();                <br><br></div>Can you share the portion of your code that uses vtkXMLPUnstructuredGridWriter? Also, looking at the data files from 6.3 and 7.1.1 may help.<br><br></div>Thanks,<br></div>Andy                          <br><div><div><div><div><div><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 5, 2018 at 12:08 PM, Andrew Parker <span><<a href="mailto:andy.john.parker@googlemail.com" target="_blank">andy.john.parker@googlemail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div><div><div><div>Andy,<br><br></div>Just moved from vtk 6.3.0 to vtk 7.1.1.  Cannot for the life of me get your modifications to work from your vtk_test.C.  I've included the vtkmpiController stuff after I init mpi myself.  I found the migration problem mentioned here:<br><br><a href="https://github.com/libMesh/libmesh/issues/1179" target="_blank">https://github.com/libMesh/libmesh/issues/1179</a><br><br></div>and here:<br><br><a href="https://gitlab.kitware.com/vtk/vtk/issues/16924" target="_blank">https://gitlab.kitware.com/vtk/vtk/issues/16924</a><br><br></div>as well as this post.  My parallel unstructured grid now only contains one source within the piece section of the pvtu, as opposed to n-processor sources (all processors have real data to write and this is not a case of empty meshes on these processors).  I can re-link and recompile against 6.3 and the problem goes away and I get the correct number of sources written to the pvtu file as before.  Can you help?  I make extensive use of pvtu files from my mpi-based code and now all grid files (pvtu and pvtp) have stopped working just by moving to 7.1.1 (which has been built with mpi-enabled).<br><br></div>Any thoughts?   <br><br></div><div>Thanks,<br></div>Andy<br><br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_5133892900030491468h5">On 6 December 2016 at 16:12, Andy Bauer <span><<a href="mailto:andy.bauer@kitware.com" target="_blank">andy.bauer@kitware.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_5133892900030491468h5"><div><div><div><div><div>Hi John,<br><br></div>The issue is indeed due to changes that I mentioned before (i.e. only data sets that have data get written out and the .pvtu file is adjusted accordingly). With this change though interprocess communication is required which is done the the vtkMultiProcessController::GetGlobalController() object which essentially acts as MPI_COMM_WORLD within VTK. You need to create a vtkMPIController which derives from vtkMultiProcessController. See the attached files for how that is done. This should also be backwards compatible but if it isn't, please let us know.<br><br></div>Best,<br></div>Andy<br><br></div>ps. Ignore the ADIOS stuff in there (it's commented out). I was trying to do two things at once and that's why the ADIOS stuff is in the CMakeLists.txt.<br></div><div class="m_5133892900030491468m_2367226779657850128HOEnZb"><div class="m_5133892900030491468m_2367226779657850128h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Dec 5, 2016 at 6:06 PM, John Peterson <span><<a href="mailto:jwpeterson@gmail.com" target="_blank">jwpeterson@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><br><div class="gmail_extra"><br><div class="gmail_quote"><span>On Mon, Dec 5, 2016 at 12:45 PM, Andy Bauer <span><<a href="mailto:andy.bauer@kitware.com" target="_blank">andy.bauer@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><div>Beyond this I can't think of any change off the top of my head that would cause your issue. I looked at the repo but can't figure out the VTK configurations from that. If you can share a simple test case which demonstrates the issue I may be able to easily diagnose it.<br></div></div></div></blockquote><div><br></div></span><div>Here's a standalone test code that reproduces the issue for me:</div><div><br></div><div><a href="https://gist.github.com/jwpeterson/92f4b8422d6fbb1056e848c9b14ee1d7" target="_blank">https://gist.github.com/jwpeterson/92f4b8422d6fbb1056e848c9b14ee1d7</a></div><div><br></div><div>Run on two procs with: mpiexec -np 2 and I get the following test.pvtu file:<br></div><div><br></div><div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><?xml version="1.0"?><span><br><VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor"><br>  <PUnstructuredGrid GhostLevel="1"><br>    <PPointData><br></span><span>      <PDataArray type="Float64" Name="u"/><br>    </PPointData><br>    <PCellData><br></span>      <PDataArray type="Int32" Name="elem_id"/><span><br>      <PDataArray type="Int32" Name="subdomain_id"/><br>      <PDataArray type="Int32" Name="processor_id"/><br>    </PCellData><br>    <PPoints><br></span>      <PDataArray type="Float64" Name="Points" NumberOfComponents="3"/><br>    </PPoints><br>    <Piece Source="test_0.vtu"/><br>  </PUnstructuredGrid><br></VTKFile></blockquote></div><div><br></div><div><br></div><div>As you can see, only one "Piece Source" is listed, but the output of the program is actually all 3 files (test.pvtu, test_0.vtu, test_1.vtu) and both the _0 and _1 files have nodes in them.  Furthermore, if I simply add a second "Piece Source" line corresponding to "test_1.vtu", the whole thing opens up just fine in Paraview.  I am a novice VTK programmer so it's very possible I'm doing something wrong, but this code definitely worked in VTK 6.x, and 7.0.</div><div><br></div><div>Thanks for your help,</div></div><div class="m_5133892900030491468m_2367226779657850128m_-8802502105114214918m_-3042913779814129432gmail_signature">John</div>
</div></div>
</blockquote></div><br></div>
</div></div><br></div></div><span>_______________________________________________<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></span></blockquote></div><br></div>
</blockquote></div><br></div>
</blockquote></div></div>