[dune-pdelab] PGFS with 1 component

Christian Engwer christian.engwer at uni-muenster.de
Tue Feb 18 07:58:05 CET 2014


Hi Pavel,

>your code works, but if you change Dune::PDELab::ISTLVectorBackend<>
>(default no blocking) to
>Dune::PDELab::ISTLVectorBackend<Dune::PDELab::ISTLParameters::static_blocking>
>(as in Martin's code), it breaks.
>
>Does it make a sense use static_blocking with blocksize 1?

In generic cover it does. Consider a general reaction diffusion system... Why should you're code fail for one component?!

Or your stokes code... Why should it fall in 1d?.. At least for debugging it is still interesting!

Christian




More information about the dune-pdelab mailing list