Hi everyone,
I have a similar problem. With the Basix documentation, I can calculate the points and weights for the master element (for example, unit triangle). But what should I do to apply this to all triangles? I want to find the Gaussian points and weights of any triangle.
With this code, I get the correct point and basis function values for the reference triangle, but what should I do to generalize this to any triangle?
import basic
import numpy as np
from basix import CellType, ElementFamily, LagrangeVariant
wht = basix.create_element(basix.ElementFamily.N1E, CellType.triangle, 1)
print(wht.dim)
points = np.array([[0.0, 0.0], [0.1, 0.1], [0.2, 0.3], [0.3, 0.6], [0.4, 1.0]])
tab = wht.tabulate(0, points)
print(tab)
print(tab.shape)
Thanks