Thanks for your kind reply. That was exactly what I was searching for. Are there any plans to implement any functions for this kind of boundary-condition-definition in future versions? i.e. integrating user-defined (spatial) field-functions over elements and mapping of the equivalent values to their nodes? Even though I was able to do this with a little help of numpy, I assume that it would be faster to do this in gmsh "natively".

Hello all,

I am currently trying to impose boundary conditions from a field function on a mesh generated with the gmsh Python API. For calculating the equivalent nodal forces, it would be quite useful to be able to retrieve the global coordinates from a parametric point in an element.

An example: I want to calculate a given function f(x,y,z) at each integration point. To do this, I would need the x(u), y(v), and z(w) functions. I know that there is a method getValue, however it only returns the coordinate on a given entity (in my case a geometric line) and not that of the element on that entity (a 1D element).

Any help (or ideas for other approaches) would be very appreciated.

F.Frank