[dune-fem] GeometryGrid and TwistUtility

[Martin Nolte]

Hi Andreas,

you are right, for AlbertaGrid LeafIntersection and LevelIntersection are of 
same type (actually, there is no level intersection unless the grid consists 
solely of the macro grid). Unfortunately, I never came around to implement the 
LevelIntersection.

Anyway, there are other grids, like SPGrid, whose leaf level is a "normal" level 
in the grid. In this case, LeafIntersection and LevelIntersection will always be 
of same type and we conclude that the implementation of the TwistUtility for 
GeometryGrid is plain wrong.

I hope to have fixed the bug in 68317f6a59382d0290eefd7fd1fa7249b911993d, but 
had not time to test it properly. Could you please have a look?

Best,

Martin

On 03/02/2015 10:03 PM, Andreas Dedner wrote:
> Hi.
> The problem seems to be that for Alberta LeafIntersection and LevelIntersection
> are the same type (is that right Martin?). Then the following two methods are
> the same:
> template< class HostGrid, class CoordFunction, class Allocator >
> struct TwistUtility< GeometryGrid< HostGrid, CoordFunction, Allocator > >
> {
>        static int twistInSelf ( const GridType &grid,
>                                 const LeafIntersection &intersection )
>        static int twistInSelf ( const GridType &grid,
>                                 const LevelIntersection &intersection )
> }
> and cause a problem.
>
> What happens if you use ALUGrid for surface and bulk?
>
> Andreas
>
>
> On 27/02/15 17:11, Agnese, Marco wrote:
>> Hi Duners,
>> in my code I have a GeometryGrid defined over an Alugrid and a lower dimension GeometryGrid over an Albertagrid. This works fine.
>>
>> Instead, if I use for both GeometryGrid-s an Alberta grid, it doesn't compile:  the error says that TwistUtility cannot be overloaded (see attachment).
>> What is the problem?
>>
>> Thank you very much,
>> cheers,
>> Marco.
>>
>>
>> _______________________________________________
>> dune-fem mailing list
>> dune-fem at dune-project.org
>> http://lists.dune-project.org/mailman/listinfo/dune-fem
>
>
>
> _______________________________________________
> dune-fem mailing list
> dune-fem at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-fem
>

-- 
Dr. Martin Nolte <nolte at mathematik.uni-freiburg.de>

Universität Freiburg                                   phone: +49-761-203-5630
Abteilung für angewandte Mathematik                    fax:   +49-761-203-5632
Hermann-Herder-Straße 10
79104 Freiburg, Germany