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

[Christophe Geuzaine]

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