store coordinates at quadrature points or faster dg computation (40X faster for VortexInBox.dg) :-)

Common/gmshpy.i

@@ -5,6 +5,7 @@
 %include std_list.i
 %{
   #include "GModel.h"
+  #include "DefaultOptions.h"
   #include "fullMatrix.h"
   #include "function.h"
   #include "dofManager.h"
@@ -124,4 +125,5 @@ namespace std {
 %include "SPoint2.h"
 %include "GPoint.h"  
 %include "functionPython.h"
+%include "DefaultOptions.h"
 

Geo/GModel.h

@@ -404,7 +404,7 @@ class GModel
                    double angle);
   GEntity *extrude(GEntity *e, std::vector<double> p1, std::vector<double> p2);
   GEntity *addPipe(GEntity *e, std::vector<GEdge *>  edges);
-
+  
   void addRuledFaces(std::vector<std::vector<GEdge *> > edges);
   GFace *addFace(std::vector<GEdge *> edges, std::vector< std::vector<double > > points);
   GFace *addPlanarFace(std::vector<std::vector<GEdge *> > edges);