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
Yes, i is a little tricky,
But before, it is for what,
They are lot way to solve this problem,
For exemple store the point A,B,C in finite element P0 or in P0edgedc ,
So you can make a visio begin of january after the 8 /01 morning.
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.
Nevertheless I would be interested in learning to define fespaces, and understanding the meaning of all the parameters involved. Probably other users would be interested also (for example What does the array of ``data[ ]'' represent during the process of adding a new finite element? ). I don’t think there is a tuto on that.
