<div dir="ltr"><div><div><div><div>Dear all,<br><br></div>my last mail was too big for some mailservers, so I'm sending attachment only containing relevant information.<br><br><br></div>Thanks!<br><br><br></div>Best regards,<br><br></div>Jelena Grižić<br></div><div class="gmail_extra"><br><div class="gmail_quote">2015-11-13 13:09 GMT+01:00 Jelena Grižić <span dir="ltr"><<a href="mailto:jelena.grizic@gmail.com" target="_blank">jelena.grizic@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi Andreas,<br><br></div>file with all requested information is in the attachment. I'm using module dune-alugrid.<br></div><div>I can also share dgf file on Google Disc or Dropbox, if necessarily.<br><br><br></div><div>Thank you!<br><br></div><div>Best regards,<br><br></div><div>Jelena <br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">2015-11-13 12:00 GMT+01:00  <span dir="ltr"><<a href="mailto:dune-request@dune-project.org" target="_blank">dune-request@dune-project.org</a>></span>:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">Send Dune mailing list submissions to<br>
        <a href="mailto:dune@dune-project.org" target="_blank">dune@dune-project.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://lists.dune-project.org/mailman/listinfo/dune" rel="noreferrer" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:dune-request@dune-project.org" target="_blank">dune-request@dune-project.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:dune-owner@dune-project.org" target="_blank">dune-owner@dune-project.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Dune digest..."<br>
<br></div></div>Today's Topics:<br>
<br>
   1. dune-alugrid: Assertion `DetDf > 0' failed (Jelena Gri?i?)<br>
   2. Re: dune-alugrid: Assertion `DetDf > 0' failed (Andreas Dedner)<br>
<br><br>---------- Proslijeđena poruka ----------<br>From: "Jelena Grižić" <<a href="mailto:jelena.grizic@gmail.com" target="_blank">jelena.grizic@gmail.com</a>><br>To: <a href="mailto:dune@dune-project.org" target="_blank">dune@dune-project.org</a><br>Cc: <br>Date: Thu, 12 Nov 2015 18:11:38 +0100<br>Subject: [Dune] dune-alugrid: Assertion `DetDf > 0' failed<br><div dir="ltr"><div><div><div><div><div><div>Dear all,<br><br><br></div>while constructing mesh from dgf file, I receive message :<br><br>dune-alugrid/dune/alugrid/3d/mappings_imp.cc:187: Dune::alu3d_ctype Dune::TrilinearMapping::det(const coord_t&): Assertion `DetDf > 0' failed.<br><br></div>Could you explain what that exactly means and what could be problem with my mesh?<br><br></div>Unfortunately, dgf file is to big (56,7 MB) to be included in mail.<br><br><br></div>Thank you.<br><br></div>Best regards,<br><br></div>Jelena Grižić<br></div>
<br><br>---------- Proslijeđena poruka ----------<br>From: Andreas Dedner <<a href="mailto:a.s.dedner@warwick.ac.uk" target="_blank">a.s.dedner@warwick.ac.uk</a>><br>To: <<a href="mailto:dune@dune-project.org" target="_blank">dune@dune-project.org</a>><br>Cc: <br>Date: Thu, 12 Nov 2015 18:44:49 +0000<br>Subject: Re: [Dune] dune-alugrid: Assertion `DetDf > 0' failed<br>Hi.<br>
The determinant in the alugrid TrilinearMapping class where the<br>
assertion is triggered, is used to compute for example the integration<br>
element or the volume of the cells in the grid. Alugrid assumes a<br>
certain orientation of the grid elements (determinant > 0). I would have<br>
thought that would be checked during the grid generation process and the<br>
elements local vertex numbering changed if necessary but you might have<br>
found a bad case here. I would need (at least):<br>
0) which alugrid type are you using, e.g., Cube<3> or Cube<2,3> or...<br>
1) the value of detDf (is it zero or negative)<br>
2) the call stack for the error<br>
3) The class also contains a variable 'a' which describes the mapping<br>
   and can be used to figure out which element in your dgf file is<br>
   causing the problem<br>
4) If that doesn't work then the dgf file...<br>
Andreas<br>
<br>
<br>
On 12/11/15 17:11, Jelena Grižić wrote:<br>
> Dear all,<br>
><br>
><br>
> while constructing mesh from dgf file, I receive message :<br>
><br>
> dune-alugrid/dune/alugrid/3d/mappings_imp.cc:187: Dune::alu3d_ctype<br>
> Dune::TrilinearMapping::det(const coord_t&): Assertion `DetDf > 0' failed.<br>
><br>
> Could you explain what that exactly means and what could be problem with<br>
> my mesh?<br>
><br>
> Unfortunately, dgf file is to big (56,7 MB) to be included in mail.<br>
><br>
><br>
> Thank you.<br>
><br>
> Best regards,<br>
><br>
> Jelena Grižić<br>
><br>
><br>
> _______________________________________________<br>
> Dune mailing list<br>
> <a href="mailto:Dune@dune-project.org" target="_blank">Dune@dune-project.org</a><br>
> <a href="http://lists.dune-project.org/mailman/listinfo/dune" rel="noreferrer" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
><br>
<br>
<br>
<br>_______________________________________________<br>
Dune mailing list<br>
<a href="mailto:Dune@dune-project.org" target="_blank">Dune@dune-project.org</a><br>
<a href="http://lists.dune-project.org/mailman/listinfo/dune" rel="noreferrer" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>