[Dune] weird intersection centers

[Marco Cisternino]

Hi Andreas,

I will try to move to the newer version of dune with the dune-alugrid module.

Do you think my code, working in the version I told, will still work?

I use some experimental stuffs, like the method hostEntity() of the entity in geometric grid...

Thanks again,


Marco


________________________________
Da: Andreas Dedner <a.s.dedner at warwick.ac.uk>
Inviato: lunedì 23 novembre 2015 20.24
A: dune at dune-project.org
Oggetto: Re: [Dune] weird intersection centers

Hi Marco.
For your email I can't decide what the mistake could be. A picture plus
some details on the input would be needed. But actually it would be even
better if you would test this again with dune 2.3 and
dune-alugrid: You could try
git clone -b releases/2.3.1
https://gitlab.dune-project.org/extensions/dune-alugrid.git
and then add dune-alugrid to the Suggestions of your module. And also
remove the with-alugrid from your config file.
Best
Andreas


On 23/11/15 19:12, Marco Cisternino wrote:
> Hi duners,
> I'm using dune 2.3 and alugrid 1.52 and I see a weird behaviour when I ask dune for intersection centers.
> The element I'm considering should have interfaces 4 and 5 centers with the same z-coordinate and different x and y coordinates.
> But it is not.
> This statement
> liit->geometry().center();
> where liit is a iterator over element intersections
> gives me the same x coordinate for the two faces. And this coordinate is equal to the element center one and it shouldn't.
> Plotting the element using the vtk writer gives me the correct element.
> I'm wondering if the source of error could be the call of center() method of the geometry in that way.
> I want to remark that faces 2 and 3 are in an analogous situation but their centers are correct.
> Faces 0 and 1 are correct too, but in a different configuration.
> I can provide you with a plot of the element if it could help.
>
> Thanks for your attention.
> Best regards,
>
> Marco Cisternino
>
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune
>


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20151126/273df82c/attachment.htm>