These sound like just the solvers and preconditioners that I need. Does this also mean that flow is able to run in parallel?
Just for the record, I found out that Eigen features a number of iterative solvers though only CG and BiCGSTAB are supported. The others may be found here <http://eigen.tuxfamily.org/dox/unsupported/namespaceEigen.html>, together with much other functionality. On 20 May 2015 at 12:00, <opm-requ...@opm-project.org> wrote: > Send Opm mailing list submissions to > opm@opm-project.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.opm-project.org/mailman/listinfo/opm > or, via email, send a message with subject or body 'help' to > opm-requ...@opm-project.org > > You can reach the person managing the list at > opm-ow...@opm-project.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Opm digest..." > > > Today's Topics: > > 1. Iterative solvers for autodiff (Tobias Ritschel) > 2. Re: Iterative solvers for autodiff (Atgeirr Rasmussen) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 19 May 2015 22:09:15 +0200 > From: Tobias Ritschel <tobiasritsc...@gmail.com> > To: opm@opm-project.org > Subject: [OPM] Iterative solvers for autodiff > Message-ID: > <CAJcuCB3zD9wykjuWAcheh8xgu=-E8wz4= > sljybrfasbqe2w...@mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Hi again > > I want to solve the fully implicit residual equations using Newtons method > with an iterative method for solving the linear system. I am using the > AutoDiffBlock objects and thus the matrix and right-hand-side is in the V > and M formats, or Array and SparseMatrix in Eigen, respectively. > > As far as I can see, the only iterative solver for non-symmetric matrices > in Eigen is BiCGSTAB and it is a bit hard to work with. Besides, I would > like some freedom of choice with regard to algorithm. > > I can see that dune-istl has several solvers and eWoms as well (perhaps the > same?). Can you give some recommendation or advice as to what functions to > use? In addition, it would be preferable if the solver can run in parallel. > > Kind regards > Tobias Ritschel > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://www.opm-project.org/pipermail/opm/attachments/20150519/d9243fd7/attachment-0001.html > > > > ------------------------------ > > Message: 2 > Date: Tue, 19 May 2015 21:09:05 +0000 > From: Atgeirr Rasmussen <atgeirr.rasmus...@sintef.no> > To: OPM Mailing List <opm@opm-project.org> > Subject: Re: [OPM] Iterative solvers for autodiff > Message-ID: <1372bfd1-bfaa-45ab-878b-67798071b...@sintef.no> > Content-Type: text/plain; charset="us-ascii" > > > 19. mai 2015 kl. 22:09 skrev Tobias Ritschel <tobiasritsc...@gmail.com>: > > > Hi again > > > > I want to solve the fully implicit residual equations using Newtons > method with an iterative method for solving the linear system. I am using > the AutoDiffBlock objects and thus the matrix and right-hand-side is in the > V and M formats, or Array and SparseMatrix in Eigen, respectively. > > > > As far as I can see, the only iterative solver for non-symmetric > matrices in Eigen is BiCGSTAB and it is a bit hard to work with. Besides, I > would like some freedom of choice with regard to algorithm. > > > > I can see that dune-istl has several solvers and eWoms as well (perhaps > the same?). Can you give some recommendation or advice as to what functions > to use? In addition, it would be preferable if the solver can run in > parallel. > > Hi Tobias, > > The solvers of dune-istl are able to run in parallel, so I would recommend > that you try those. Which solver and preconditioner > to choose would depend on the problem, but we do have an implementation > using a CPR-type two-stage preconditioning for > fully implicit systems in the classes NewtonIterationBlackoilCPR (overall > solution) and CPRPreconditioner (preconditioner part) > that are used by the 'flow' fully implicit black-oil simulator. Those are > implemented with dune-istl, and have some choice of > options/solvers, so you may find them useful to look at. There is a class > (DuneMatrix) that allows you to construct a matrix that > can be used with the Dune solvers from an Eigen matrix, it could also be > useful to you. > > Good luck with your project! > > Atgeirr > > > > > > Kind regards > > Tobias Ritschel > > _______________________________________________ > > Opm mailing list > > Opm@opm-project.org > > http://www.opm-project.org/mailman/listinfo/opm > > > > > ------------------------------ > > Subject: Digest Footer > > _______________________________________________ > Opm mailing list > Opm@opm-project.org > http://www.opm-project.org/mailman/listinfo/opm > > > ------------------------------ > > End of Opm Digest, Vol 32, Issue 6 > ********************************** >
_______________________________________________ Opm mailing list Opm@opm-project.org http://www.opm-project.org/mailman/listinfo/opm