#1724 ISTL complains about non-matching right hand side together with UMFPack
Metadata
| Property | Value |
|---|---|
| Reported by | Christoph Grüninger (gruenich@iws.uni-stuttgart.de) |
| Reported at | Oct 20, 2015 15:53 |
| Type | Bug Report |
| Version | Git (pre3.0) |
| Operating System | Unspecified / All |
Description
I already have mentioned this problem in #1702 (closed). Now it occurred for one of our students and my office machine. The umfpacktest fails with Dune 3.0-git, but works fine with 2.4. I could not figure out any relevant changes since the 2.4 release.
dune-istl/build-cmake/dune/istl/test> make umfpacktest && ./umfpacktest
[100%] Built target umfpacktest
testing for N=100 BS=1
UMFPACK V5.7.1 (Oct 10, 2014): OK
[UMFPack Decomposition]
Wallclock Time taken: 0.0205744 (CPU Time: 0.0205744)
Flops taken: 2.35676e+07
Peak Memory Usage: 5.14298e+06 bytes
Condition number estimate: 4.6416
Numbers of non-zeroes in decomposition: L: 206332 U: 206332
UMFPACK V5.7.1 (Oct 10, 2014): OK
[UMFPack Solve]
Wallclock Time: 0.00166421 (CPU Time: 0.00166421)
Flops Taken: 2.22706e+06
Iterative Refinement steps taken: 1
Error Estimate: 1.32357e-16 resp. 0
Dune reported error: Dune::ISTLError [apply:/temp/gruenich/dune/complete/dune-istl/dune/istl/umfpack.hh:318]: Size of right-hand-side vector b does not match the number of matrix rows!