[Dune] Projection direction in dune-grid-glue with psurface

[Miguel de Benito Delgado]

Hi,

On 14 May, 2013, at 16:05, Oliver Sander <sander at igpm.rwth-aachen.de> wrote:

> Patches for the tests and for the code are welcome.

I had to change ReferenceElements to GenericReferenceElements, as in the attachment. Otherwise things seem to work just fine, though I haven't really tested much yet.

> Consequently, these two vectors fields are not exactly the same.
> If you want to obtain b) through the interface you have to use
> 
> intersection->geometryInOutside().global(pos) - intersection->geometry().global(pos)


I thought I had done this, but instead I have

    auto ref = GenericReferenceElements<ctype, dim>::general (is->inside()->type());
    auto pos = is->geometryInInside().local (ref.position (i, dim));
    auto nor = is->geometryOutside().global (pos) - is->geometry().global (pos);

I'm not sure this is the best way to accomplish what I need, though.

Best,
--
Miguel de Benito.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fix.diff
Type: application/octet-stream
Size: 550 bytes
Desc: not available
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20130516/60c49382/attachment.obj>