[Dune-devel] dune-istl umfpack support

[Christoph Grüninger]

Hi Dominic, hi Markus,
I gave the UMFPACK branch a try. Thanks for the work. I wonder whether
UMFPACK is faster compared to superLU? Did you any test?

On my openSuse UMFPACK was not found, because I needed additional libraries
  -lcholmod -lcolamd -lamd -lsuitesparseconfig -llapack
which I hacked into to m4 test and dune/istl/test/Makefile.am. Maybe
this is a bug in the openSuse package as libumfpack should be linked to
these libraries. Had anyone similar problems?

Find attached a patch to get rid of two compiler warnings and gitignore
the test.

With Clang 3.3 it worked for me, GCC 4.7 and 4.8 segfaulted. I attached
the back trace:

> #0  0x00007ffff6483a3e in umfpack_di_free_symbolic () from /usr/lib64/libumfpack-5.6.1.so
> #1  0x00000000004053c0 in Dune::UMFPackMethodChooser<double>::free_symbolic<void**> () at ../../../dune/istl/umfpack.hh:73
> #2  0x0000000000403553 in Dune::UMFPack<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > > >::free (
>     this=0x7fffffffd4c0) at ../../../dune/istl/umfpack.hh:401
> #3  0x0000000000404880 in Dune::UMFPack<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > > >::setMatrix (this=0x7fffffffd4c0, _mat=...) at ../../../dune/istl/umfpack.hh:360
> #4  0x000000000040394d in Dune::UMFPack<Dune::BCRSMatrix<Dune::FieldMatrix<double, 1, 1>, std::allocator<Dune::FieldMatrix<double, 1, 1> > > >::UMFPack
>     (this=0x7fffffffd4c0, mat_=..., file=0x40af50 "umfpack_decomp", verbose=0) at ../../../dune/istl/umfpack.hh:256
> #5  0x0000000000401fce in main (argc=1, argv=0x7fffffffd7f8) at umfpacktest.cc:70

Ciao
Christoph

-- 
Une science n'était vraiment développée que quand elle
pouvait utiliser les mathématiques.    (Paul Lafargue)
*********************************************
CMWR 2014: 10th - 13th June 2014 in Stuttgart
         Please visit www.cmwr14.de
*********************************************
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 0001-umfpack-Get-rid-of-warnings.patch
Type: text/x-patch
Size: 1505 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune-devel/attachments/20131017/2e8dcd9d/attachment.bin>