[dune-pdelab] PGFS with 1 component

[Martin Rückl]

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 :-(

HtH!
Martin Rueckl




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20140217/f72631d2/attachment.htm>