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

[Tobias Malkmus]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hi Andrea, Andreas

Since the example of Andrea does not compile, there are some files
missing, it wrote a small example code for the case domainSpace !=
rangeSpace with ISTLLinearOeprator, main file is attached.

The problem is that the preconditioners (istl and fem) are hard linked
into the matrix class. Most of the preconditioners expect domainSpace
== rangeSpace.

In a private branch i have extracted them into a separate class some
time ago, i will see whether i can merge this changes into the master.

Andrea, this means you will have to change to the master branch of
dune-fem or will have to wait for the next release.

I think its time to review the matrix implementations and bindings.

Best Tobias

On 10/02/2014 02:43 PM, Sacconi, Andrea wrote:
> 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
> 
> 
> 
> _______________________________________________ dune-fem mailing
> list dune-fem at dune-project.org 
> http://lists.dune-project.org/mailman/listinfo/dune-fem
> 


- -- 
Tobias Malkmus                 <tomalk at mathematik.uni-freiburg.de>

Mathematisches Institut               Tel: +49 761 203 5627
Abt. für Angewandte Mathematik        Universität Freiburg
Hermann-Herder-Str. 10
79104 Freiburg

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v2

iQEcBAEBAgAGBQJULWxVAAoJEBMY7DHvcWNgPBgH/RCy6iLjPnl7BbXsawvSZWtL
FUI9zi9nufkkJyewgwGEOO5t5VVM5sLP6Vq4Y7Q4CZpYRK9111L5JfnitKAzjee1
5Vvhqv63wDJBcLCDk+Q/MlT8E4OgLE7tx6aPxZH9dRSl3QAemHp7rTA/c3yWNaqv
lCB4EMkMNB06g1MDfTuO2W6iWxkcBcDm5+RCJDg6BYJ34OqBIBTQ94RxszJejNu9
sjZ7qsHzIEPlfRirUCO0GJNrNX3+O8P6CNTfxVWXWTwGDqQ3I5so1ykmBmf3SCc1
UNS7CRK1XJZD6cH9oUC5KM9WmAkgb+bLdFmd25j3HUAHMd7qzyO0Oqdku3jrcRo=
=TsFV
-----END PGP SIGNATURE-----
-------------- next part --------------
A non-text attachment was scrubbed...
Name: main.cc
Type: text/x-c
Size: 2772 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune-fem/attachments/20141002/f353b8de/attachment.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: main.cc.sig
Type: application/pgp-signature
Size: 287 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune-fem/attachments/20141002/f353b8de/attachment.sig>