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

[Oliver Sander]

Am 27.05.2015 um 09:29 schrieb Markus Blatt:
> 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).

Why not?  Isn't ISTL all duck-typing?
Best,
Oliver



-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 473 bytes
Desc: OpenPGP digital signature
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150527/39353d01/attachment.sig>