[Dune] Initialization of BlockVector with C array or Eigen array

Markus Blatt markus at dr-blatt.de
Wed May 27 09:29:35 CEST 2015


On Wed, May 27, 2015 at 09:05:48AM +0200, Oliver Sander wrote:
> Am 26.05.2015 um 21:38 schrieb Tobias Ritschel:
> > Hi
> > 
> > I would like to use dune-istl to solve linear systems in a fully implicit
> > two-phase solver and I use Eigen to generate the Jacobi matrix and
> > right-hand-side. Is there a way to initialize a
> > BlockVector<FieldVector<double,1>> with a Eigen array or a C array without
> > copying?
> 
> IMO the most elegant solution would be to write a new istl vector class that behaves
> like BlockVector<FieldVector> but uses Eigen under the hood. 

Please be aware that using this appoach some of the more powerful
preconditioners might not be usable (e.g. AMG).

Markus

-- 
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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 836 bytes
Desc: Digital signature
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150527/df6fe128/attachment.sig>


More information about the Dune mailing list