On 28/10/14 13:19, Atgeirr Rasmussen wrote:
hi,
my .02c;
However, I am afraid that it
may end up as "Yet Another Sparse Matrix" and only add to any complexity or
confusion there already is. Abandoning the Eigen
SparseMatrix for opm-autodiff will require a substantial effort, that I do not
think we'll make (at this point at least), so I think that will remain.
this. no matter which is chosen, conversions will be the consequence and
these are nontrivial for large systems. having the flexibility to use
the native type of whatever underlying solver library used should be the
goal. that is, a thin wrapper on top of native matrix types. istl got
this right with its LinearOperators (sadly, but obviously, only as long
as no algebraic preconditioners are used).
arnem
_______________________________________________
Opm mailing list
[email protected]
http://www.opm-project.org/mailman/listinfo/opm
