optimization for bilinear quads

libsrc/meshing/curvedelems.cpp

@@ -2254,6 +2254,18 @@ namespace netgen
             }
           break;
         }
+      case QUAD:
+        {
+          if (info.order >= 2) return false; // not yet supported
+          AutoDiff<2,T> lami[4] = { (1-x)*(1-y), x*(1-y), x*y, (1-x)*y };
+          for (int j = 0; j < 4; j++)
+            {
+              Point<3> p = mesh[el[j]];
+              for (int k = 0; k < DIM_SPACE; k++)
+                mapped_x[k] += p(k) * lami[j];
+            }
+          break;
+        }
       default:
         return false;
       }