<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">I just realized that the TrilinearMap
class does not only store the a[.][.] (which describes the actual
mapping) but<br>
also the actual coordinates of the cube in
p0,p1,p2,p3,p4,p5,p6,p7.<br>
Having those of the offending cube would actually be even more
interesting. So if you could add those to your output<br>
that would simplify figuring out what the problem is.<br>
Andreas<br>
<br>
On 17/11/15 16:02, Andreas Dedner wrote:<br>
</div>
<blockquote cite="mid:564B4F87.4020800@warwick.ac.uk" type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<div class="moz-cite-prefix">Hi.<br>
To get that output of detDf and a[i][j] you must have added some
lines to the code - could you send me the corresponding<br>
diff.<br>
Andreas<br>
<br>
On 17/11/15 15:44, Jelena Grižić wrote:<br>
</div>
<blockquote
cite="mid:CABBkLKgMSw=kU5fBQ86ZknEBDcT0Hw4CV=ie-aym_q_CJJcSWw@mail.gmail.com"
type="cite">
<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 moz-do-not-send="true"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:jelena.grizic@gmail.com"
target="_blank">jelena.grizic@gmail.com</a>><br>
To: <a moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:a.s.dedner@warwick.ac.uk"
target="_blank">a.s.dedner@warwick.ac.uk</a>><br>
To: <<a moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:Dune@dune-project.org"
target="_blank">Dune@dune-project.org</a><br>
> <a moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:Dune@dune-project.org"
target="_blank">Dune@dune-project.org</a><br>
<a moz-do-not-send="true"
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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Dune mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:Dune@dune-project.org">Dune@dune-project.org</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</a>
</pre>
</blockquote>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Dune mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Dune@dune-project.org">Dune@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</a>
</pre>
</blockquote>
<br>
</body>
</html>