<div dir="ltr"><div><div><div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">a simpler, relatively straighforward approach could be<br>
to read your polygons into gmsh and then have gmsh<br>
triangulate them, into triangles.</blockquote><div> </div>@Benedikt - Yes, that is precisely what I thought I should do but did not know you could do the triangulation in Gmsh. Any pointers? I still have to use quadrature on polygons in my code and for that I am going to need the triangulation anyway.<br></div>@Andreas, Christian: I have been having trouble building opm-cornerpoint module on my laptop. Is it best to seek help on their mailing list to get it working with Dune?<br><br></div>Thanks<br></div>Ganesh<br><div><div><div><div><br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Feb 26, 2015 at 2:31 PM, Christian Engwer <span dir="ltr"><<a href="mailto:christian.engwer@uni-muenster.de" target="_blank">christian.engwer@uni-muenster.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span>On Thu, Feb 26, 2015 at 03:05:13PM +0100, Benedikt Oswald wrote:<br>
> aha. But still, what about triangulating into triangles and assigning a<br>
> certain tag to a group of elements<br>
> that form the original polygon, then evaluating the algorthm based on<br>
> triangles with that same tag ? Benedikt<br>
<br>
</span>It will make handling of DOFs nearly impossible. Often you have things<br>
like one DOF per intersection/face (in the case of conforming polygons<br>
every intersection corresponfs to a face), but now you have polygon<br>
internal faces and resolving this is all but simple.<br>
<span><font color="#888888"><br>
Christian<br>
</font></span><div><div><br>
> Am 26/02/15 um 15:01 schrieb Christian Engwer:<br>
> > Hi Benedikt,<br>
> ><br>
> > he says he wants to use the virtual element method and there the main<br>
> > idea is to have a method which wokrs on arbitrary polygons.<br>
> ><br>
> > Christian<br>
> ><br>
> > On Thu, Feb 26, 2015 at 02:56:01PM +0100, Benedikt Oswald wrote:<br>
> >> Hello Ganesh<br>
> >><br>
> >> a simpler, relatively straighforward approach could be<br>
> >> to read your polygons into gmsh and then have gmsh<br>
> >> triangulate them, into triangles. The same approach<br>
> >> would also work in 3D. Then, there is not need to<br>
> >> write a new grid. Greetings, Benedikt<br>
> >><br>
> >> Am 26/02/15 um 14:49 schrieb Christian Engwer:<br>
> >>> Hi,<br>
> >>><br>
> >>> currently no grid can handle such elements. Once upon a time there was<br>
> >>> a simple implementation, which was able to handle arbitrary polygonal<br>
> >>> elements in 2D with out refinement. But this doesn't work with the<br>
> >>> current Dune anymore. On the other hand I believe it wouldn't be too<br>
> >>> difficult to implement such a mesh, as we have introducent several<br>
> >>> changes in the interface to allow for non-standard meshes (this was<br>
> >>> all with the cornerpoint mesh in mind).<br>
> >>><br>
> >>> Ciao<br>
> >>> Christian<br>
> >>><br>
> >>> On Thu, Feb 26, 2015 at 12:02:02PM +0000, Andreas Dedner wrote:<br>
> >>>> Hi.<br>
> >>>> The main problem is not the grid type or the factory but that you<br>
> >>>> need an actual<br>
> >>>> grid implementation that can handle these general meshes. ALU and<br>
> >>>> SPGrid etc.<br>
> >>>> can only handle cubes and simplex elements and that is the only<br>
> >>>> geometry types<br>
> >>>> they can understand in their GridFactory.<br>
> >>>><br>
> >>>> So for something more general you need to provide a new grid<br>
> >>>> implementation that<br>
> >>>> allows for more general grid types. The geometry type of these<br>
> >>>> elements would be<br>
> >>>> "none" and probably you would also be implementing only part of the<br>
> >>>> geometry interface.<br>
> >>>> Something along those lines was done in the CornerPointGrid implementation.<br>
> >>>> <a href="https://github.com/OPM/dune-cornerpoint" target="_blank">https://github.com/OPM/dune-cornerpoint</a><br>
> >>>><br>
> >>>> Best<br>
> >>>> Andreas<br>
> >>>><br>
> >>>><br>
> >>>><br>
> >>>> On 26/02/15 11:38, Ganesh Diwan wrote:<br>
> >>>>> Dear list,<br>
> >>>>><br>
> >>>>> I am working on implementing virtual element computations inside<br>
> >>>>> Dune. This will involve implementing the algorithm on arbitrary<br>
> >>>>> shaped elements. Is there a way to read arbitrary polygonal (2d)<br>
> >>>>> or polyhedral (3d) meshes in Dune? Currently implemented Grid<br>
> >>>>> managers (ALUGrid, SPGrid etc.) identify either simplices or cube<br>
> >>>>> types.<br>
> >>>>><br>
> >>>>> Can the GridFactory class be used for this purpose? Grid-howto<br>
> >>>>> document says when building a new gridtype Dune may have a problem<br>
> >>>>> for unidentified element type. Say, I have an element with five<br>
> >>>>> edges (a pentagon). In this case, is Dune likely to throw an<br>
> >>>>> exception for unidentified element type?<br>
> >>>>><br>
> >>>>> I wanted to get the Grid guru's thoughts before investing time in this.<br>
> >>>>><br>
> >>>>> Thank you!<br>
> >>>>> Ganesh<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" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><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" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
> >> --<br>
> >> ------------------------------------------------------------------------------------------------<br>
> >> Dr. sc. techn. Benedikt Oswald - first engineer - LSPR AG - phone - <a href="tel:%2B41%2043%20366%2090%2074" value="+41433669074" target="_blank">+41 43 366 90 74</a><br>
> >> Technoparkstrasse 1, CH-8005 Zürich, <a href="mailto:benedikt.oswald@lspr.ch" target="_blank">benedikt.oswald@lspr.ch</a> - labor vincit omnia improbus<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" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
> >><br>
><br>
> --<br>
> ------------------------------------------------------------------------------------------------<br>
> Dr. sc. techn. Benedikt Oswald - first engineer - LSPR AG - phone - <a href="tel:%2B41%2043%20366%2090%2074" value="+41433669074" target="_blank">+41 43 366 90 74</a><br>
> Technoparkstrasse 1, CH-8005 Zürich, <a href="mailto:benedikt.oswald@lspr.ch" target="_blank">benedikt.oswald@lspr.ch</a> - labor vincit omnia improbus<br>
> ------------------------------------------------------------------------------------------------<br>
><br>
><br>
<br>
</div></div><span>--<br>
Prof. Dr. Christian Engwer<br>
Institut für Numerische und Angewandte Mathematik<br>
Fachbereich Mathematik und Informatik der Universität Münster<br>
Einsteinstrasse 62<br>
48149 Münster<br>
<br>
E-Mail  <a href="mailto:christian.engwer@uni-muenster.de" target="_blank">christian.engwer@uni-muenster.de</a><br>
Telefon <a href="tel:%2B49%20251%2083-35067" value="+492518335067" target="_blank">+49 251 83-35067</a><br>
FAX             <a href="tel:%2B49%20251%2083-32729" value="+492518332729" target="_blank">+49 251 83-32729</a><br>
<br>
</span><div><div>_______________________________________________<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" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
</div></div></blockquote></div><br></div></div></div></div></div></div>