[Dune] Runtime error with parallel linear solver

[Tobias Ritschel]

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.

On 11 June 2015 at 11:36, Markus Blatt <markus at dr-blatt.de> wrote:

> On Thu, Jun 11, 2015 at 10:38:10AM +0200, Tobias Ritschel wrote:
> > Ah, this is clearly my error. I am using the constructor
> >
> > OwnerOverlapCopyCommunication
> > <
> http://www.dune-project.org/doc/doxygen/dune-istl-html/class_dune_1_1_owner_overlap_copy_communication.html#a93f446939bd4d6623bdf2e544f183f8c
> >
> > (MPI_Comm
> > comm_, SolverCategory::Category
> > <
> http://www.dune-project.org/doc/doxygen/dune-istl-html/struct_dune_1_1_solver_category.html#ae061380ac961490be6ee353cf0dc1733
> >
> > cat_=SolverCategory::overlapping,
> > bool freecomm_=false)
> >
> > and used the created object without an index set. Is the information
> > automatically set up when constructing an IndexInfoFromGrid object, or
> > should they be set with AddLocalIndex and AddRemoteIndex? And is this the
> > most convenient way to set up the indices?
>
> Don't use IndexInfoGrid. I assume that is a left over from ancient
> dune-disc times. You can access the underlying index set by
> OwnerOverlapCopyCommunication::indexSet() and use that to add the
> indices with the correct tag. Later on get the remote indices with
> OwnerOverlapCopyCommunication::remoteIndices() and call rebuild on it,
>
> You will find more information on the index sets and how they are used
> in dune-istl in the "Description of the parallel communication
> interface" on https://dune-project.org/doc/index.html. You need one
> index per degree of freedom and use the attribute owner to mark a
> partitioning. In the easiest case you need to mark connected degrees
> of freedom that are not owner on the process as copy.
>
> If you are using dune-cornerpoint with one unknown per cell, then you
> can use ParallelISTLInformation::copyValuesTo from opm-core to do the job.
>
> Otherwise there is some manual work to do. E.g. compute the owner
> region using grid communication.
>
> Markus
>
> >
> > On 11 June 2015 at 09:59, Markus Blatt <markus at dr-blatt.de> wrote:
> >
> > > Hi,
> > >
> > > On Wed, Jun 10, 2015 at 09:00:36PM +0200, Tobias Ritschel wrote:
> > > >
> > > > terminate called after throwing an instance of
> > > > 'Dune::InterfaceBuilder::RemotexIndicesStateError'
> > > >
> > > > As far as I can see, it is called when I call the member function
> apply
> > > on
> > > > the RestartedGMRES solver object.
> > > >
> > >
> > > This only happens if the information in the RemoteIndices is not in
> > > sync with the underlying IndexSets. This can be case if these were
> > > changed, but there was not call to RemoteIndices::rebuild().
> > > Thi is also reflected in the error message: "RemoteIndices is not in
> > > sync with the index set. Call RemoteIndices::rebuild first!"
> > >
> > > So the question is: How do you set up the index sets and the remote
> > > information that you are using in OwnerOverlapCopyCommunication?
> > >
> > > Cheers,
> > >
> > > Markus
> > >
> > >
> > > Dr. Markus Blatt - HPC-Simulation-Software & Services
> > > http://www.dr-blatt.de
> > > Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany
> > > Tel.: +49 (0) 160 97590858
> > >
> > > _______________________________________________
> > > Dune mailing list
> > > Dune at dune-project.org
> > > http://lists.dune-project.org/mailman/listinfo/dune
> > >
> > >
>
> --
> Do you need more support with DUNE or HPC in general?
>
> Dr. Markus Blatt - HPC-Simulation-Software & Services
> http://www.dr-blatt.de
> Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany
> Tel.: +49 (0) 160 97590858
>
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150611/2ec059ea/attachment.htm>