<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Feb 24, 2014 at 6:57 AM, De Groof, Vincent Frans Maria <span dir="ltr"><<a href="mailto:Vincent.De-Groof@uibk.ac.at" target="_blank">Vincent.De-Groof@uibk.ac.at</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Hi,
<div><br>
</div>
<div><br>
</div>
<div>I have an IS object, which lists per process certain DOFs in local numbering. I want to get the global numbering, so I apply ISLocalToGlobalMappingApplyIS. However, while I can view them in local numbering using the ISView(ISLocal, PETSC_VIEWER_STDOUT_WORLD),
 I cannot view them in global numbering using ISView(ISGlobal, VIEWER_STDOUT_WORLD). It complains that "Arguments must have same communicators". </div>
<div><br>
</div>
<div>Is this normal? Can someone explain what is happening? I assumed that the ISLocalToGlobalMappingApplyIS would just copy the communicator from the local IS.</div></div></div></blockquote><div><br></div><div>This is a bug. IS used to be purely serial, but then it evolved parallel abilities (it was bitten by a radioactive spider), so there</div>
<div>are some PETSC_COMM_SELFs hidden in the code. I will fix this, but you can just change the creation line in</div><div>src/vec/is/utils/isltog.c:ISLocalToGlobalMappingApplyIS() to use PetscObjectComm((PetscObject) is) instead.</div>
<div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div><div style="direction:ltr;font-size:10pt;font-family:Tahoma">
<div>Code snippet:</div>
<div><br>
</div>
<div>
<div><span style="white-space:pre-wrap"></span>IS<span style="white-space:pre-wrap">
</span>ISLocal, ISGlobal;</div>
<div><span style="white-space:pre-wrap"></span>ierr = ISCreateGeneral(PETSC_COMM_WORLD, nvtx, vtx_idx, PETSC_COPY_VALUES, &ISLocal); CHKERRQ(ierr);</div>
<div><span style="white-space:pre-wrap"></span>ierr = ISLocalToGlobalMappingApplyIS(ISg2lMap, ISLocal, &ISGlobal); CHKERRQ(ierr);</div>
<div><span style="white-space:pre-wrap"></span>ISView(ISLocal, PETSC_VIEWER_STDOUT_WORLD);</div>
<div><span style="white-space:pre-wrap"></span>ISView(ISGlobal, PETSC_VIEWER_STDOUT_WORLD);   <-- returns the error.</div>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>thanks,</div>
<div>Vincent</div>
</div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>