[dune-pdelab] PGFS with 1 component

Steffen Müthing steffen.muething at iwr.uni-heidelberg.de
Mon Feb 17 23:59:29 CET 2014


Hi Martin,

Am 17.02.2014 um 06:27 schrieb Martin Rückl <mrueckl at physik.hu-berlin.de>:

> Hi PDELab team,
> i am using the Dune::PDELab::PowerGridFunctionSpace in one of my generic classes. At the moment I’m porting lots of my code to the new Dune/Pdelab version. I think in the old version of PDElab (like one year ago, before the changes with the dof handling) it was possible to create a PowerGridFunctionSpace with only a single component, i.e.: 
> 
>     typedef Dune::PDELab::PowerGridFunctionSpace<GFS,1,VBE,Dune::PDELab::EntityBlockedOrderingTag> PGFS;
> 
> while the above is still possible, i get an assertion error when interpolating an initial solution into the dof vector:
> 
>     Dune::PDELab::interpolate(initialSolution,pgfs,u);
> 
> -> assert(i==-1) in dune/pdelab/backend/istl/vectorhelpers.hh line #27 fails
> 
> i had a short glimpse into the callstack, and  I think the cache type had size 1 hence i resulted in 0. If you like, i can provide a full call stack(which looks REALLY ugly :-() and/or fill a bug report.
> 
> Since i added a little bit of extra functionality around the ISTLBlockVectorConainer, I’m not 100% sure whether it is an PDELab issue :-(

that’s weird, I just clobbered together a small example (see attached file), and that works without a problem.
Is your initialSolution a PowerGridFunction?

Steffen

-------------- next part --------------
A non-text attachment was scrubbed...
Name: testgridfunctionspace.cc
Type: application/octet-stream
Size: 2759 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20140217/ac67019a/attachment.obj>
-------------- next part --------------

> 
> HtH!
> Martin Rueckl
> 
> 
> 
> 
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 495 bytes
Desc: Message signed with OpenPGP using GPGMail
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20140217/ac67019a/attachment.sig>


More information about the dune-pdelab mailing list