[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 Oliver Sander
> Gesendet: 28 August 2015 16:22
> An: dune at dune-project.org
> Betreff: Re: [Dune] Read in matrix and rhs, solve linear system
> 
> 
> >
> > @Markus: Does that say sth. to you?
> 
> If you build your program with -g then the backtrace will contain line
> numbers.
> That will be more informative.

[Buesing, Henrik] I tried to add 

CFLAGS='-g'

to my DUNE options file, but that seemed to have no effect.


> 
> >
> > (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<do
> uble, 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::Field!
>  Ma!
> >  trix<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<do
> uble, 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::allocato!
>  r<!
> >  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
> > _______________________________________________
> > 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