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

[Oliver Sander]

>
> @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.
Best,
Oliver

>
> (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::FieldMa!
>  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<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
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune