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

Tobias Ritschel tobiasritschel at gmail.com
Wed May 27 12:51:16 CEST 2015


I see. I think I will copy the vector then :)

On 27 May 2015 at 11:59, Markus Blatt <markus at dr-blatt.de> wrote:

> 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
>
> _______________________________________________
> 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/20150527/6abf79fe/attachment.htm>


More information about the Dune mailing list