[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