#1727 kamgtest fails after throwing an exception
Metadata
Property | Value |
---|---|
Reported by | Elias Pipping (elias.pipping@fu-berlin.de) |
Reported at | Oct 25, 2015 12:09 |
Type | Bug Report |
Version | Git (pre3.0) |
Operating System | Unspecified / All |
Last edited by | Elias Pipping (elias.pipping@fu-berlin.de) |
Last edited at | Oct 25, 2015 13:45 |
Closed by | Elias Pipping (elias.pipping@fu-berlin.de) |
Closed at | Oct 25, 2015 13:45 |
Closed in version | Unknown |
Resolution | Duplicate |
Comment | This looks like a duplicate of bug #1724 (closed). Sorry about that. |
Description
The test ./dune/istl/paamg/test/kamgtest fails for me.
I've modified dune-istl slightly to make things a tad easier to debug, namely in two ways: (1) make the test print the exception that is thrown (a size mismatch) and (2) make umfpack print those numbers.
Here's the resulting output:
% ./dune/istl/paamg/test/kamgtest N=100 coarsenTarget=1200 maxlevel=10 Level 0 has 10000 unknowns, 10000 unknowns per proc (procs=1) aggregating finished. Level 1 has 2500 unknowns, 2500 unknowns per proc (procs=1) aggregating finished. Level 2 has 600 unknowns, 600 unknowns per proc (procs=1) operator complexity: 1.30319 Using a direct coarse solver (UMFPACK) Building hierarchy of 3 levels (inclusive coarse solver) took 0.0168803 seconds. Building hierarchy took 0.0169085 seconds === CGSolver Iter Defect Rate 0 70.7254 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 1 1.99283 0.0281771 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 2 0.666525 0.334461 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 3 0.146794 0.220237 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 4 0.0441693 0.300894 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 5 0.0138672 0.313957 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 6 0.00395074 0.284897 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 7 0.00142139 0.359778 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 8 0.000299216 0.21051 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 9 0.000119257 0.398566 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 umfpackMatrix_.N() = 600, b.size() = 600 umfpackMatrix_.M() = 600, x.size() = 600 10 2.57486e-05 0.215908 === rate=0.227048, T=0.0347132, TIT=0.00347132, IT=10 AMG solving took 0.0347341 seconds AMG building took 4.8709 iterations AMG building together with solving took 0.0516426 N=100 coarsenTarget=1200 maxlevel=10 Level 0 has 10000 unknowns, 10000 unknowns per proc (procs=1) aggregating finished. Level 1 has 2500 unknowns, 2500 unknowns per proc (procs=1) aggregating finished. Level 2 has 600 unknowns, 600 unknowns per proc (procs=1) operator complexity: 1.30319 Using a direct coarse solver (UMFPACK) Building hierarchy of 3 levels (inclusive coarse solver) took 0.0159227 seconds. Building hierarchy took 0.0159388 seconds === CGSolver Iter Defect Rate 0 99.8816 umfpackMatrix_.N() = 1200, b.size() = 600 umfpackMatrix_.M() = 1200, x.size() = 600 Dune::ISTLError [apply:[..]/dune/istl/umfpack.hh:323]: Size of right-hand-side vector b does not match the number of matrix rows! %