I coded loop subdivision and ran it on a cube with triangle faces. I noticed this:
I am using doubles to represent my vertex positions, I am not 100% sure these are caused by floating point imprecisions, but since the general shape looks ok I am not sure what else it could be.
Any suggestions on how to find the cause? My linear algebra library is glm and I am coding this in C++.