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 ?

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

Dimitri