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

Markus Blatt markus at dr-blatt.de
Wed May 27 11:59:44 CEST 2015


On Tue, May 26, 2015 at 09:38:22PM +0200, Tobias Ritschel wrote:

> 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?

Just a comment: There is not so much overhead when copying a
vector. It does not need much time. In any case dune-istl will create
quite a few vectors depending on the solver/preconditioner
(e.g. BiCGStab creates 5, AMG approximately another 5).

Saving just one for the cost of spending valuable thesis time might
not be a good investment.

And in OPM there is already quite some copying ;)

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/25206630/attachment.sig>


More information about the Dune mailing list