
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">

<p>Hi Andreas,</p>

<p>I will try to move to the newer version of dune with the dune-alugrid module.</p>

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

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

<p>Thanks again,</p>

<p><br>

</p>

<p>Marco<br>

</p>

<br>

<br>

<div style="color: rgb(0, 0, 0);">

<div>

<hr tabindex="-1" style="display:inline-block; width:98%">

<div id="x_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>Da:</b> Andreas Dedner <a.s.dedner@warwick.ac.uk><br>

<b>Inviato:</b> lunedė 23 novembre 2015 20.24<br>

<b>A:</b> dune@dune-project.org<br>

<b>Oggetto:</b> Re: [Dune] weird intersection centers</font>

<div> </div>

</div>

</div>

<font size="2"><span style="font-size:10pt;">

<div class="PlainText">Hi Marco.<br>

For your email I can't decide what the mistake could be. A picture plus<br>

some details on the input would be needed. But actually it would be even<br>

better if you would test this again with dune 2.3 and<br>

dune-alugrid: You could try<br>

git clone -b releases/2.3.1<br>

<a id="LPlnk904987" href="https://gitlab.dune-project.org/extensions/dune-alugrid.git">https://gitlab.dune-project.org/extensions/dune-alugrid.git</a><br>

and then add dune-alugrid to the Suggestions of your module. And also<br>

remove the with-alugrid from your config file.<br>

Best<br>

Andreas<br>

<br>

<br>

On 23/11/15 19:12, Marco Cisternino wrote:<br>

> Hi duners,<br>

> I'm using dune 2.3 and alugrid 1.52 and I see a weird behaviour when I ask dune for intersection centers.<br>

> The element I'm considering should have interfaces 4 and 5 centers with the same z-coordinate and different x and y coordinates.<br>

> But it is not.<br>

> This statement<br>

> liit->geometry().center();<br>

> where liit is a iterator over element intersections<br>

> gives me the same x coordinate for the two faces. And this coordinate is equal to the element center one and it shouldn't.<br>

> Plotting the element using the vtk writer gives me the correct element.<br>

> I'm wondering if the source of error could be the call of center() method of the geometry in that way.<br>

> I want to remark that faces 2 and 3 are in an analogous situation but their centers are correct.<br>

> Faces 0 and 1 are correct too, but in a different configuration.<br>

> I can provide you with a plot of the element if it could help.<br>

> <br>

> Thanks for your attention.<br>

> Best regards,<br>

> <br>

> Marco Cisternino<br>

> <br>

> _______________________________________________<br>

> Dune mailing list<br>

> Dune@dune-project.org<br>

> <a href="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</a><br>

> <br>

<br>

<br>

</div>

</span></font></div>

</div>

</body>

</html>