[Dune] Read in matrix and rhs, solve linear system

[Buesing, Henrik]

> -----Ursprüngliche Nachricht-----
> Von: dune-bounces+hbuesing=eonerc.rwth-aachen.de at dune-project.org
> [mailto:dune-bounces+hbuesing=eonerc.rwth-aachen.de at dune-project.org] Im
> Auftrag von Markus Blatt
> Gesendet: 01 September 2015 14:08
> An: dune at dune-project.org
> Betreff: Re: [Dune] Read in matrix and rhs, solve linear system
> 
> Hi,
> 
> On Tue, Sep 01, 2015 at 08:27:16AM +0000, Buesing, Henrik wrote:
> > [Buesing, Henrik] I now get the following output (see below). The
> solver does not crash anymore, but gets NaNs instead...
> 
> 
> That is weired and I have no idea why this happens. Can you try another
> solver (e.g. BiCGSTABSolver instead of GeneralizedPCGSolver)?
> 
[Buesing, Henrik] Yes! I have substituted Dune::GeneralizedPCGSolver with Dune::BiCGSTABSolver in matrixmarket-amg.cc. Still NaNs (for output see below. Note, this is the output for the smaller matrix with the same behavior.)

> Is it possible to publish your test matrix and right hand side for
> testing?
> 
[Buesing, Henrik] Yes. I have uploaded a smaller matrix and rhs, which has the same behavior, but allows for easier testing. 
The following link is valid for one week:
https://gigamove.rz.rwth-aachen.de/download/id/WbADRKiCzjuCp8


Output of matrixmarket-amg

matrixmarket-amg           
%%MatrixMarket                                                                                
%%MatrixMarket                                                                                
Level 0 has 196608 unknowns, 196608 unknowns per proc (procs=1)                               
aggregating finished.                                                                         
Level 1 has 72552 unknowns, 72552 unknowns per proc (procs=1)
aggregating finished.
Level 2 has 23068 unknowns, 23068 unknowns per proc (procs=1)
aggregating finished.
Level 3 has 11228 unknowns, 11228 unknowns per proc (procs=1)
aggregating finished.
Level 4 has 4197 unknowns, 4197 unknowns per proc (procs=1)
aggregating finished.
Level 5 has 2039 unknowns, 2039 unknowns per proc (procs=1)
aggregating finished.
Level 6 has 895 unknowns, 895 unknowns per proc (procs=1)
operator complexity: 1.53953
Using a direct coarse solver (UMFPACK)
Building hierarchy of 7 levels (inclusive coarse solver) took 8.79153 seconds.
Building hierarchy took 8.79157 seconds
=== BiCGSTABSolver
 Iter          Defect            Rate
    0        0.0390357
  0.5              nan              nan
    1              nan              nan
  1.5              nan              nan
    2              nan              nan
  2.5              nan              nan
    3              nan              nan 

-- 
Dipl.-Math. Henrik Büsing
Institute for Applied Geophysics and Geothermal Energy
E.ON Energy Research Center
RWTH Aachen University
------------------------------------------------------
Mathieustr. 10            |    Tel +49 (0)241 80 49907
52074 Aachen, Germany     |    Fax +49 (0)241 80 49889
------------------------------------------------------
http://www.eonerc.rwth-aachen.de/GGE
hbuesing at eonerc.rwth-aachen.de
------------------------------------------------------