<div dir="ltr">Thanks, I will have a look at the documentation. I solve for both pressure and saturation in each cell, so I I'll create the index set.</div><div class="gmail_extra"><br><div class="gmail_quote">On 11 June 2015 at 11:36, Markus Blatt <span dir="ltr"><<a href="mailto:markus@dr-blatt.de" target="_blank">markus@dr-blatt.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On Thu, Jun 11, 2015 at 10:38:10AM +0200, Tobias Ritschel wrote:<br>
> Ah, this is clearly my error. I am using the constructor<br>
><br>
> OwnerOverlapCopyCommunication<br>
</span>> <<a href="http://www.dune-project.org/doc/doxygen/dune-istl-html/class_dune_1_1_owner_overlap_copy_communication.html#a93f446939bd4d6623bdf2e544f183f8c" target="_blank">http://www.dune-project.org/doc/doxygen/dune-istl-html/class_dune_1_1_owner_overlap_copy_communication.html#a93f446939bd4d6623bdf2e544f183f8c</a>><br>
> (MPI_Comm<br>
> comm_, SolverCategory::Category<br>
> <<a href="http://www.dune-project.org/doc/doxygen/dune-istl-html/struct_dune_1_1_solver_category.html#ae061380ac961490be6ee353cf0dc1733" target="_blank">http://www.dune-project.org/doc/doxygen/dune-istl-html/struct_dune_1_1_solver_category.html#ae061380ac961490be6ee353cf0dc1733</a>><br>
<span class="">> cat_=SolverCategory::overlapping,<br>
> bool freecomm_=false)<br>
><br>
> and used the created object without an index set. Is the information<br>
> automatically set up when constructing an IndexInfoFromGrid object, or<br>
> should they be set with AddLocalIndex and AddRemoteIndex? And is this the<br>
> most convenient way to set up the indices?<br>
<br>
</span>Don't use IndexInfoGrid. I assume that is a left over from ancient<br>
dune-disc times. You can access the underlying index set by<br>
OwnerOverlapCopyCommunication::indexSet() and use that to add the<br>
indices with the correct tag. Later on get the remote indices with<br>
OwnerOverlapCopyCommunication::remoteIndices() and call rebuild on it,<br>
<br>
You will find more information on the index sets and how they are used<br>
in dune-istl in the "Description of the parallel communication<br>
interface" on <a href="https://dune-project.org/doc/index.html" target="_blank">https://dune-project.org/doc/index.html</a>. You need one<br>
index per degree of freedom and use the attribute owner to mark a<br>
partitioning. In the easiest case you need to mark connected degrees<br>
of freedom that are not owner on the process as copy.<br>
<br>
If you are using dune-cornerpoint with one unknown per cell, then you<br>
can use ParallelISTLInformation::copyValuesTo from opm-core to do the job.<br>
<br>
Otherwise there is some manual work to do. E.g. compute the owner<br>
region using grid communication.<br>
<span class="HOEnZb"><font color="#888888"><br>
Markus<br>
</font></span><span class="im HOEnZb"><br>
><br>
> On 11 June 2015 at 09:59, Markus Blatt <<a href="mailto:markus@dr-blatt.de">markus@dr-blatt.de</a>> wrote:<br>
><br>
> > Hi,<br>
> ><br>
> > On Wed, Jun 10, 2015 at 09:00:36PM +0200, Tobias Ritschel wrote:<br>
> > ><br>
> > > terminate called after throwing an instance of<br>
> > > 'Dune::InterfaceBuilder::RemotexIndicesStateError'<br>
> > ><br>
> > > As far as I can see, it is called when I call the member function apply<br>
> > on<br>
> > > the RestartedGMRES solver object.<br>
> > ><br>
> ><br>
> > This only happens if the information in the RemoteIndices is not in<br>
> > sync with the underlying IndexSets. This can be case if these were<br>
> > changed, but there was not call to RemoteIndices::rebuild().<br>
> > Thi is also reflected in the error message: "RemoteIndices is not in<br>
> > sync with the index set. Call RemoteIndices::rebuild first!"<br>
> ><br>
> > So the question is: How do you set up the index sets and the remote<br>
> > information that you are using in OwnerOverlapCopyCommunication?<br>
> ><br>
> > Cheers,<br>
> ><br>
> > Markus<br>
> ><br>
> ><br>
</span><div class="HOEnZb"><div class="h5">> > Dr. Markus Blatt - HPC-Simulation-Software & Services<br>
> > <a href="http://www.dr-blatt.de" target="_blank">http://www.dr-blatt.de</a><br>
> > Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany<br>
> > Tel.: <a href="tel:%2B49%20%280%29%20160%2097590858" value="+4916097590858">+49 (0) 160 97590858</a><br>
> ><br>
> > _______________________________________________<br>
> > Dune mailing list<br>
> > <a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
> > <a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
> ><br>
> ><br>
<br>
--<br>
Do you need more support with DUNE or HPC in general?<br>
<br>
Dr. Markus Blatt - HPC-Simulation-Software & Services <a href="http://www.dr-blatt.de" target="_blank">http://www.dr-blatt.de</a><br>
Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany<br>
Tel.: <a href="tel:%2B49%20%280%29%20160%2097590858" value="+4916097590858">+49 (0) 160 97590858</a><br>
</div></div><br>_______________________________________________<br>
Dune mailing list<br>
<a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
<a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
<br></blockquote></div><br></div>