[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