<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body text="#000000" bgcolor="#FFFFFF">
    <div class="moz-cite-prefix">A final suggestion I could make:<br>
      start a new project on the dune user wiki.<br>
      Move the geometrygrid files from dune-grid/dune/grid/geometrygrid
      <br>
      to newproject/dune/grid/geometrygrid<br>
      Now you should be able to test everything including
      dune-geometrygrid as<br>
      before and everything can be easily copied back to dune-grid when
      completed.<br>
      Andreas<br>
      <br>
      On 23.01.2014 14:04, Benedikt Oswald wrote:<br>
    </div>
    <blockquote cite="mid:E50ABFD0-0404-4F32-B0A5-A47656095D2F@lspr.ch"
      type="cite">
      <pre wrap="">hello DUNE

so, reflecting the comments the approach to get curvilinear meshes
should be along those lines:

1) extend the gush reader to read tetrahedral grids up to 3rd order
    gmsh at present works up to 3rd order, higher order mesh generation
    crashes

2) use the grid factory facility to create the grid in dune

3) implement the curvilinear geometry grid in

        a) multilineargeometry.hh
        b) geometry-grid


4) since we would like to integrate this into Dune as soon as it is mature,
    hearing all the suggestions, where should we put the work ?

    keep it on our machines ? keep it in a project on the Dune user wiki ?

 Greetings, Benedikt

   





On Jan 23, 2014, at 11:48 AM, Oliver Sander <a class="moz-txt-link-rfc2396E" href="mailto:sander@igpm.rwth-aachen.de"><sander@igpm.rwth-aachen.de></a> wrote:

</pre>
      <blockquote type="cite">
        <pre wrap="">Am 23.01.2014 10:25, schrieb Andreas Dedner:
</pre>
        <blockquote type="cite">
          <pre wrap="">What I would suggest would be to setup a project on the dune user wiki and
</pre>
        </blockquote>
        <pre wrap="">
A good idea, but while doing that I suggest to keep in mind that in the long
run it may be preferable to merge your code into GeometryGrid proper.
You are implementing a more general version of GeometryGrid; it is not obvious
why it should be a completely separate implementation (except to get you started).

</pre>
        <blockquote type="cite">
          <pre wrap="">giving it public access - that hopefully will give you some help from the
experienced dune developers if you get stuck.
I'm terrible with naming so will not make a suggestion but
CurvilinearGeometryGrid and simply CurvilinearGrid comes to mind.
</pre>
        </blockquote>
        <pre wrap="">
-1, I think this is too unspecific.  If you intend to implement polynomial
geometries how about PolynomialGrid?  Or QuadraticGrid or SecondOrderGrid
if you only want to do quadratic geometries?

Best,
Oliver

</pre>
        <blockquote type="cite">
          <pre wrap="">Best
Andreas

On 23.01.2014 10:13, Benedikt Oswald wrote:
</pre>
          <blockquote type="cite">
            <pre wrap="">Hello Andreas, thanks.

Then, I shall do so!

However, not being an expert in Dune internals programming,
where would I start ?

would the extension of dune / grid / geometry grid / coordfunction be
the right place ?

And, more technically speaking, should copy the geometry grid
directory into my
own source tree, implement the things there and, once it works, submit
back to Dune?

What about naming the new implementation ?

Greetings, Benedikt




--------------------------------------------------------------------------------------------------------------------------------------------------------------
Dr. sc. techn. Benedikt Oswald - first engineer - LSPR AG - phone -
+41 43 366 90 74
Technoparkstrasse 1, CH-8005 Zürich, <a class="moz-txt-link-abbreviated" href="mailto:benedikt.oswald@lspr.ch">benedikt.oswald@lspr.ch</a>
<a class="moz-txt-link-rfc2396E" href="mailto:benedikt.oswald@lspr.ch"><mailto:benedikt.oswald@lspr.ch></a>
--------------------------------------------------------------------------------------------------------------------------------------------------------------

On Jan 23, 2014, at 10:04 AM, Andreas Dedner <<a class="moz-txt-link-abbreviated" href="mailto:a.s.dedner@warwick.ac.uk">a.s.dedner@warwick.ac.uk</a>
<a class="moz-txt-link-rfc2396E" href="mailto:a.s.dedner@warwick.ac.uk"><mailto:a.s.dedner@warwick.ac.uk></a>> wrote:

</pre>
            <blockquote type="cite">
              <pre wrap="">Hi Benedikt.
You would need to extend GeometryGrid (something many of us would
really appreciate). Basically GeometryGrid replaces the geometries of
the host grid
using the GenericGeometries which are only "linear". This part of
GeometryGrid
would need to be replaced with something higher order, e.g., some
curvilinear approximation.
Best
Andreas

On 23.01.2014 08:32, Benedikt Oswald wrote:
</pre>
              <blockquote type="cite">
                <pre wrap="">Dear Dune

as it happens we are forced to consider drastic steps to increase
the efficiency of our electromagnetic solvers,
especially in the frequency domain.

Recently, Peter Bastian suggested we use curvilinear elements in
order to model our curved geometry and domain
boundaries. In fact, we believe that this will help us enormously
since the element count will be much reduced.

Peter also suggested we use the geometry grid on top of ALUGrid in
order to achieve this.

Now, my concrete questions is: given an existing curvilinear mesh
(e.g. from gush), how would we do this ?

I read that geometry grid manages the vertices of the underlying
grid, but, a curvilinear tetrahedral element has
vertices not only in the corners but also on the edges, faces and
inside the volume.
Perhaps my question is trivial, nevertheless, I'd appreciate your
suggestions.

I should say that gmsh is capable of generating quite nice
curvilinear tetrahedral meshes, at least up to 3rd order.

Greetings and thanks in advance, Benedikt



--------------------------------------------------------------------------------------------------------------------------------------------------------------
Dr. sc. techn. Benedikt Oswald - first engineer - LSPR AG - phone -
+41 43 366 90 74
Technoparkstrasse 1, CH-8005 Zürich, <a class="moz-txt-link-abbreviated" href="mailto:benedikt.oswald@lspr.ch">benedikt.oswald@lspr.ch</a>
<a class="moz-txt-link-rfc2396E" href="mailto:benedikt.oswald@lspr.ch"><mailto:benedikt.oswald@lspr.ch></a>
--------------------------------------------------------------------------------------------------------------------------------------------------------------



_______________________________________________
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>
              <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-rfc2396E" href="mailto:Dune@dune-project.org"><mailto: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>
            <pre wrap="">
</pre>
          </blockquote>
          <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>
        <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>
      <pre wrap="">
</pre>
      <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>