[dune-fem] Does anyone use fem::ISTLLinearOperator for rectangular matrices?

[Sacconi, Andrea]

Hi Andreas,

the problem is (or at least I think it is) related to the definition of the preconditioner.
When the matrix is square, everything goes fine, I have already solved Poisson and other stuff using ISTLLinearOperator class successfully.

When instead I need to assemble a rectangular matrix, since DomainType and RangeType are different, issues arise. In this case, preconditioning doesn't have much sense, and the compiler complains that methods are not compatible or some conversions can't be made. I attach a simple test case (with Makefile.am, main file, operator definition and output of the compilation).

Thanks for your help!
Andrea
__________________________________________________________

Andrea Sacconi
PhD student, Applied Mathematics
AMMP Section, Department of Mathematics, Imperial College London,
London SW7 2AZ, UK
a.sacconi11 at imperial.ac.uk

________________________________________
From: dune-fem-bounces+a.sacconi11=imperial.ac.uk at dune-project.org [dune-fem-bounces+a.sacconi11=imperial.ac.uk at dune-project.org] on behalf of Andreas Dedner [a.s.dedner at warwick.ac.uk]
Sent: 02 October 2014 13:27
To: dune-fem at dune-project.org
Subject: Re: [dune-fem] Does anyone use fem::ISTLLinearOperator for rectangular matrices?

Hi Andreas.
Whats the issue? I'm not doing it (but I thought I tried that in some
project).
What sort of error messages are you getting?
Andreas

On 02/10/14 11:35, Sacconi, Andrea wrote:
> Hi all,
>
> is there anyone in dune-fem community successfully using and assembling a rectangular operator (for example, divergence of the velocity in Stokes, or any operator with domain and range with different dimensions), with Dune::Fem::ISTLLinearOperator as linear operator counterpart (instead of Dune::Fem::SparseRowLinearOperator)?
>
> Cheers,
> Andrea
> __________________________________________________________
>
> Andrea Sacconi
> PhD student, Applied Mathematics
> AMMP Section, Department of Mathematics, Imperial College London,
> London SW7 2AZ, UK
> a.sacconi11 at imperial.ac.uk
>
> _______________________________________________
> dune-fem mailing list
> dune-fem at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-fem


_______________________________________________
dune-fem mailing list
dune-fem at dune-project.org
http://lists.dune-project.org/mailman/listinfo/dune-fem
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rectangularISTL.zip
Type: application/x-zip-compressed
Size: 18025 bytes
Desc: rectangularISTL.zip
URL: <https://lists.dune-project.org/pipermail/dune-fem/attachments/20141002/fa144519/attachment.bin>