[Getdp] [Gmsh] Desirable Feature for getdp: higher order, at least quadratic shape functions also in 3D

Christophe Geuzaine cgeuzaine at ulg.ac.be
Fri Jan 13 13:57:37 CET 2012


On 13 Jan 2012, at 12:09, moritz braun wrote:

> Dear Christophe
> 
> I am only looking at higher interpolation order for the fields
> and for it seems that getdp cannot  handle higher then linear for the
> tetrehedral elements
> i.e. 3D
> 

It can actually handle second order elements (and third order for a subset of element types). Here's a simple example.

Hope this helps,

Christophe

-------------- next part --------------
A non-text attachment was scrubbed...
Name: example.pro
Type: application/octet-stream
Size: 2378 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20120113/71d55e79/attachment.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: example.geo
Type: application/octet-stream
Size: 398 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20120113/71d55e79/attachment.geo>
-------------- next part --------------




> Thanks for your info
> 
> I will have a look at the code and see whether I understand it
> I would also be nice to find a list of all allowed combinations of element
> geometries and orders of interpolations
> 
> regards
> 
> Moritz
> 
> On Fri, Jan 13, 2012 at 1:02 PM, Christophe Geuzaine
> <cgeuzaine at ulg.ac.be> wrote:
>> 
>> On 13 Jan 2012, at 11:31, moritz braun wrote:
>> 
>>> Dear Colleagues
>>> 
>>> In order to make use of getdp for
>>> work leadting to potentially publishable atomic / molecular physics results
>>> I would appreciate being able to use higher order shape functions
>>> in getdp. Since I am a reasonably good programmer I would
>>> also be glad to help, but a look at the source code for getdp showed
>>> me that it is not easy to see where the shape functiojns  are actually defined.
>>> 
>> 
>> Hi Moritz,
>> 
>> In GetDP the shape functions used for the geometry and the interpolation of the fields are defined separately. At the moment, you can use high-order basis functions for the interpolation of the fields, but only linear ones for the geometry.
>> 
>> Do you need high order shape functions for the geometry (i.e., curved elements) ? Or only for the interpolation of the fields ?
>> 
>> (In GetDP the basis functions are hardcoded in the getdp/Legacy/BF_*.cpp files)
>> 
>> Hope this helps,
>> 
>> Christophe
>> 
>>> regards
>>> 
>>> Moritz
>>> 
>>> 
>>> 
>>> --
>>> Prof M Braun         Tel.:27-12-4298006/8027
>>> Physics Department  Fax.: 27-12-4293643
>>> University of South Africa (UNISA)
>>> moritz.braun at gmail.com
>>> P.O. Box 392
>>> 0003
>>> UNISA
>>>  South Africa
>>> 
>>> _______________________________________________
>>> gmsh mailing list
>>> gmsh at geuz.org
>>> http://www.geuz.org/mailman/listinfo/gmsh
>> 
>> --
>> Prof. Christophe Geuzaine
>> University of Liege, Electrical Engineering and Computer Science
>> http://www.montefiore.ulg.ac.be/~geuzaine
>> 
>> 
>> 
> 
> 
> 
> -- 
> Prof M Braun         Tel.:27-12-4298006/8027
> Physics Department  Fax.: 27-12-4293643
> University of South Africa (UNISA)
> moritz.braun at gmail.com
> P.O. Box 392
> 0003
> UNISA
>  South Africa

-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 
http://www.montefiore.ulg.ac.be/~geuzaine