Re: [OPM] Iterative solvers for autodiff

2015-05-21 Thread Atgeirr Rasmussen
21. mai 2015 kl. 10:25 skrev Tobias Ritschel tobiasritsc...@gmail.com: Thanks for the great feedback. I will proceed with the DUNE solvers and I might compare to the Eigen solvers for a simple setup. I also considered PETSc as I use TAO (part of PETSc) to carry out optimization, though I

Re: [OPM] Iterative solvers for autodiff

2015-05-20 Thread Arne Morten Kvarving
On 20/05/15 14:02, Markus Blatt wrote: Like I said I am biased as a DUNE developer (parallel iterative solvers). But if you do not want to use DUNE, then I would rather recommend PETSc, Trilinos or you will find others that might be more suitable from Jack Dongarra's list

Re: [OPM] Iterative solvers for autodiff

2015-05-20 Thread Atgeirr Rasmussen
20. mai 2015 kl. 12:06 skrev Tobias Ritschel tobiasritsc...@gmail.com: These sound like just the solvers and preconditioners that I need. Does this also mean that flow is able to run in parallel? No, it is not. But a variant (flow_cp) is, however it is still in heavy development and

Re: [OPM] Iterative solvers for autodiff

2015-05-20 Thread Markus Blatt
On Wed, May 20, 2015 at 12:06:41PM +0200, Tobias Ritschel wrote: 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

Re: [OPM] Iterative solvers for autodiff

2015-05-20 Thread Tobias Ritschel
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

[OPM] Iterative solvers for autodiff

2015-05-19 Thread Tobias Ritschel
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