Hi Dimitri,

On 04/24/2016 08:48 PM, Dimitri Balasoiu wrote:

Hello everyone,

I have a question about sfepy's definition of regions.

I want to define a region Omega_1, wich is a subregion of Omega. So I type :

regions = { 'Omega' : 'all', 'Omega_1' : ('vertices of group 1', 'cell', Omega),}

, where group 1 is a group of vertices.

Question : wich field element P1 will be created over the region Omega_1 ?

More precisely, I have a triangle (A,B,C), with node A in group 1, but the other ones outside of group 1.

Will the finite element wich is worth 1 at A, and 0 at B and C be in my field element basis ?

If the field was defined in Omega_1 (assuming it is the default, volume, field and not a surface field), it will use only the entire cells in Omega_1 - so no, cells that only touch the region in a vertex or edge/face are not in the field.

But you can define a field over the whole domain Omega, and integrate over Omega_1. What do you try to achieve?

In wich file are the basis functions for the fe fields defined ?

See sfepy/discrete/fem/poly_spaces.py.

r.