Hi, I was wondering if there is a reason to use a left-handed coordinate system for the 2D-elements in the function `SF::get_transformed_pts()`. For the 3D-elements, the coordinates are just copied so they are still given in a conventional right-handed coordinate system which leads to inconsistencies between the element formulations.
Would be a pretty easy fix by just exchanging `y = cross(x, z)` with `y = cross(z, x)` but I didn't want to do that without knowing whether there is a reason to do it the way it's implemented.