[dune-pdelab] Vector and Matrix backends for PETSc

Matthias Wohlmuth Matthias.Wohlmuth at ma.tum.de
Mon Sep 19 12:15:07 CEST 2011


  Hi,
the recently introduced petsc-backend added trouble vor codes using 
complex numbers because of

   typedef typename Backend::template Matrix<double> Matrix;
instead of
   typedef typename 
LFSU::Traits::FiniteElementType::Traits::LocalBasisType::Traits::RangeFieldType 
ElementType;
   typedef typename Backend::template Matrix<ElementType> Matrix;

Similar to several, other pieces of code in istl/common ....  it only 
works (and compiles!) with real numbers. Of course, this patch is not a 
systematic complex-numbers patch for the trunk, but perhaps this helps 
to be aware of this problem when adding new code.
The attached patch fixes the problems for the Accessor-Petsc class and 
similar places in the same directory.

Best,
Matthias





Am 12.09.2011 20:32, schrieb Steffen Müthing:
> If you want to look at the changes I did to the interface to shoehorn PETSc in there somehow,
> the three relevant commits are:
>
> r1607: Changes the way the local assembler accesses the matrix (the user perspective)
> r1613: Makes the ISTLMatrixBackend compatible with that new access pattern
> r1619: Not really important, but makes sure the default GridOperator actually does the right thing
>
> Steffen
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab
>


-- 
*****************************************************************
***************  Dipl.-Math. Matthias Wohlmuth
************
**********  Department of Mathematics (M2)
********  Technische Universität München (TU Munich)
******  Boltzmannstr. 3, 85748 Garching, Germany
****
**  http://www-m2.ma.tum.de/bin/view/Allgemeines/MatthiasWohlmuth
****  matthias.wohlmuth at ma.tum.de
******  Phone: +49 89 / 289 18441
********  Fax: +49 89 / 289 18435
**********  Office: 03.10.036
*****************************************************************

-------------- next part --------------
A non-text attachment was scrubbed...
Name: istlbackend_patches.diff
Type: text/x-patch
Size: 2408 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20110919/2894e7a9/attachment.bin>


More information about the dune-pdelab mailing list