[Dune] Dune Digest, Vol 125, Issue 7

[Andreas Dedner]

Hi.
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
diff.
Andreas

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