[dune-pdelab] Fwd: Re: PGFS with 1 component

[Oliver Sander]

Hi Martin,
please send such mails to the dune-pdelab list.  I personally have very little knowledge
of the pdelab internals -- I only know how to write helpful bug reports.
Cheers,
Oliver


-------- Original-Nachricht --------
Betreff: Re: [dune-pdelab] PGFS with 1 component
Datum: Mon, 17 Feb 2014 16:55:49 +0100
Von: Martin Rückl <mrueckl at physik.hu-berlin.de>
An: Oliver Sander <sander at igpm.rwth-aachen.de>

Hi Oliver,
I have hijacked one of your pdelab course examples (attached). Might be
I'm doing something wrong with the PGFS setup, but at least it compiles
fine for one and two component PGFS. However with 1 component execution
results in

mrueckl at dong30:/home/mrueckl/dune-git/dune-pdelab-howto/build-cmake/src/course-examples$
example05 1 2 3 4
example05:
/home/mrueckl/dune-git/dune-pdelab/dune/pdelab/backend/istl/vectorhelpers.hh:27:
typename Block::field_type&
Dune::PDELab::istl::access_vector_element(Dune::PDELab::istl::tags::field_vector_1,
Block&, const CI&, int) [with CI = Dune::PDELab::MultiIndex<long
unsigned int, 2ul>; Block = Dune::FieldVector<double, 1>; typename
Block::field_type = double]: Assertion `i == -1' failed.
...
/
/This also means, its not an issue of my

/little bit of extra functionality around the ISTLBlockVectorConainer/

;-)

Martin



On 17/02/14 15:39, Oliver Sander wrote:
> Am 17.02.2014 15:27, schrieb 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
>>
> Hi Martin,
> I think the most helpful thing would be if you could provide a small test program
> that triggers the bug, and which only relies on pdelab functionality.
> Cheers,
> Oliver
>




-------------- next part --------------
A non-text attachment was scrubbed...
Name: example05_QkQk.hh
Type: text/x-c++hdr
Size: 3110 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20140217/6c187453/attachment.hh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 534 bytes
Desc: OpenPGP digital signature
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20140217/6c187453/attachment.sig>