<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>