<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">I have done this a few tme and it did
seem to work (but possibly by chance).<br>
But if you keep the tests within the new project it will
definitely work.<br>
<br>
<br>
On 23.01.2014 15:32, Markus Blatt wrote:<br>
</div>
<blockquote cite="mid:20140123153232.GI4407@dr-blatt.de" type="cite">
<pre wrap="">On Thu, Jan 23, 2014 at 01:40:49PM +0000, Andreas Dedner wrote:
</pre>
<blockquote type="cite">
<pre wrap="">A final suggestion I could make:
start a new project on the dune user wiki.
Move the geometrygrid files from dune-grid/dune/grid/geometrygrid
to newproject/dune/grid/geometrygrid
Now you should be able to test everything including dune-geometrygrid as
before and everything can be easily copied back to dune-grid when completed.
Andreas
</pre>
</blockquote>
<pre wrap="">
Actually I am not sure that this would work. Pretty much depends on
how the include directories take precedence.
The include path when compiling dune-grid-howto is
-I/include/path/dune-grid-howto
-I/include/path//dune-common
-I/include/path//dune-common
-I/include/path//dune-geometry
-I/include/path/dune-geometry
-I/include/path/dune-grid
-I/include/path/dune-grid
-I/include/path/dune-istl
-I/include/path/dune-istl
-I/include/path/dune-grid-howto
So this will probably not work if you another module that depends on
your new one.
</pre>
<blockquote type="cite">
<pre wrap="">On 23.01.2014 14:04, Benedikt Oswald wrote:
</pre>
<blockquote 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>
</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="">
_______________________________________________
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>
</blockquote>
<pre wrap="">
</pre>
<blockquote type="cite">
<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>