Fantastic! I did not try loading the 827 with ElmerGUI yet, but I will put in the patch so I don't notice a problem when I do try.
Edit: I did try, and it crashed. Looking at the code, it's not so simple to fix. For example there's a line
Code: Select all
int facenodes = degree * faceedges;
Which gives 8 for 820, but you can't get 9 from that. I can see how to add facemap9[] which would include the extra face point, but it takes a bit more TODO than I can figure out for the moment. I can always do 1st order, check my mesh with ElmerGUI, then go to second order and feed the mesh to ElmerSolver. For now, that's plenty good enough.