Consider using `basix::FiniteElement::push_forward` for evaluating basis functions.

[jorgensd]

Makes sense for non-standard elements (RT, N1curl etc), where expression would be very slow.
(https://github.com/FEniCS/basix/blob/eff3bb56cf80ba9507bc0495e371d5efb80924d1/cpp/basix/finite-element.h#L577)
Should be written in C++