
<html>

  <head>

    <meta content="text/html; charset=ISO-8859-1"

      http-equiv="Content-Type">

  </head>

  <body bgcolor="#FFFFFF" text="#000000">

    Hi PDELab team,<br>

    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.: <br>

    <br>

    <small><small><font face="monospace"><big><big><big>    typedef

                Dune::PDELab::PowerGridFunctionSpace<GFS,1,VBE,Dune::PDELab::EntityBlockedOrderingTag>

                PGFS;</big></big></big></font></small></small><br>

    <br>

    while the above is still possible, i get an assertion error when

    interpolating an initial solution into the dof vector:<br>

    <br>

    <big><font face="monospace">   

        Dune::PDELab::interpolate(initialSolution,pgfs,u);</font></big><br>

    <br>

    -> assert(i==-1) in dune/pdelab/backend/istl/vectorhelpers.hh

    line #27 fails<br>

    <br>

    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.<br>

    <br>

    Since i added a little bit of extra functionality around the

    ISTLBlockVectorConainer, I’m not 100% sure whether it is an PDELab

    issue :-(<br>

    <br>

    HtH!<br>

    Martin Rueckl<br>

    <br>

    <br>

    <br>

    <br>

  </body>

</html>