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

[Buesing, Henrik]

> > [Buesing, Henrik] I get the following output, which is not really
> helpful...
> >
> > Catchpoint 1 (exception thrown), 0x0000003b292bcca0 in __cxa_throw ()
> > from /usr/lib64/libstdc++.so.6 Missing separate debuginfos, use:
> > debuginfo-install atlas-3.8.4-2.el6.x86_64 blas-3.2.1-4.el6.x86_64
> > glibc-2.12-1.149.el6_6.7.x86_64 libgcc-4.4.7-11.el6.x86_64
> > libgfortran-4.4.7-11.el6.x86_64 libstdc++-4.4.7-11.el6.x86_64
> 
> Look at the backtrace.

[Buesing, Henrik] Here it is!

@Markus: Does that say sth. to you? 

(gdb) backtrace
#0  0x0000003b292bcca0 in __cxa_throw () from /usr/lib64/libstdc++.so.6
#1  0x0000000000487d60 in Dune::Amg::AMG<Dune::MatrixAdapter<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::SeqSSOR<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, 1>, Dune::Amg::SequentialInformation, std::allocator<Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > > >::mgc(Dune::Amg::AMG<Dune::MatrixAdapter<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::SeqSSOR<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, 1>, Dune::Amg::SequentialInformation, std::allocator<Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > > >::LevelContext&) ()
#2  0x0000000000485074 in Dune::Amg::AMG<Dune::MatrixAdapter<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::SeqSSOR<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >, 1>, Dune::Amg::SequentialInformation, std::allocator<Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > > >::apply(Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >&, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > const&) ()
#3  0x0000000000461ae2 in Dune::GeneralizedPCGSolver<Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > > >::apply(Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >&, Dune::BlockVector<Dune::FieldVector<double, 1>, std::allocator<Dune::FieldVector<double, 1> > >&, Dune::InverseOperatorResult&) ()
#4  0x000000000045dd3e in main ()


> 
> 
> >
> >>> Henrik
> >>>
> >>>
> >>> %%MatrixMarket
> >>> %%MatrixMarket
> >>> Level 0 has 6.3058e+06 unknowns, 6.3058e+06 unknowns per proc
> >>> (procs=1) aggregating finished.
> >>> Level 1 has 2.45706e+06 unknowns, 2.45706e+06 unknowns per proc
> >>> (procs=1) aggregating finished.
> >>> Level 2 has 904947 unknowns, 904947 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Level 3 has 500740 unknowns, 500740 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Level 4 has 179584 unknowns, 179584 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Level 5 has 85507 unknowns, 85507 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Level 6 has 39132 unknowns, 39132 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Level 7 has 21563 unknowns, 21563 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Level 8 has 13195 unknowns, 13195 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Level 9 has 10225 unknowns, 10225 unknowns per proc (procs=1)
> >>> aggregating finished.
> >>> Stopped coarsening because of rate breakdown 10225/9256=1.10469<1.2
> >>> operator complexity: 1.62054 Building hierarchy of 10 levels
> >>> (inclusive coarse solver) took 223.879 seconds.
> >>> Building hierarchy took 223.88 seconds === GeneralizedPCGSolver
> >>>  Iter          Defect            Rate
> >>>     0        0.0181497
> >>> terminate called after throwing an instance of 'Dune::MathError'
> >>>
> >>>
> >>>
> >>> _______________________________________________
> >>> Dune mailing list
> >>> Dune at dune-project.org
> >>> http://lists.dune-project.org/mailman/listinfo/dune
> >>
> >> _______________________________________________
> >> Dune mailing list
> >> Dune at dune-project.org
> >> http://lists.dune-project.org/mailman/listinfo/dune
> > _______________________________________________
> > Dune mailing list
> > Dune at dune-project.org
> > http://lists.dune-project.org/mailman/listinfo/dune
> 
> 
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune