Hello to all experts,
I am wondering if it is possible to define a new finite element type, where the interpolation points (as well as the base) depend not only on the current triangle, but also on the neighbours. This is the case if we consider for example the gravity center of each triangle, join them to find the intersection with an edge, as on the figure
Dear Frédéric,
Thank you for your answer. Upon your advice, I have found a bypass, by using the fact that my space is a subspace of a known fespace on the split mesh Ths=splitmesh6PowellSabin(Th0).
I have performed a “space reduction” similar to what is done with on(tgv=-1): I replaced some well chosen lines in my matrix by constraints defining the subspace.