diff --git a/CMakeLists.txt b/CMakeLists.txt
index e9078711ecc7634139e86ed9a62150dcb8ed69a0..0cdf4d7352f2773c3c4f99fecbe7d635159e05ad 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -72,7 +72,7 @@ set(GMSH_API
   ${CMAKE_CURRENT_BINARY_DIR}/Common/GmshVersion.h
   Common/Gmsh.h Common/Context.h Common/GmshDefines.h Common/GmshMessage.h
     Common/VertexArray.h Common/Octree.h Common/OctreeInternals.h
-  Numeric/Numeric.h Numeric/Gauss.h Numeric/polynomialBasis.h
+  Numeric/Numeric.h Numeric/GaussIntegration.h Numeric/polynomialBasis.h
     Numeric/JacobianBasis.h Numeric/fullMatrix.h
     Numeric/simpleFunction.h Numeric/cartesian.h
   Geo/GModel.h Geo/GEntity.h Geo/GPoint.h Geo/GVertex.h Geo/GEdge.h 
diff --git a/Geo/MElement.h b/Geo/MElement.h
index f4711dc1ba87d6d6e855bdd9e085083f82ebd106..63b58ccc74c91ffb8dd86192e48df852b63a36ca 100644
--- a/Geo/MElement.h
+++ b/Geo/MElement.h
@@ -16,9 +16,7 @@
 #include "MFace.h"
 #include "polynomialBasis.h"
 #include "JacobianBasis.h"
-#include "Gauss.h"
-
-//class GFace; ??
+#include "GaussIntegration.h"
 
 // A mesh element.
 class MElement
diff --git a/Mesh/Levy3D.cpp b/Mesh/Levy3D.cpp
index 61d146b059fccc37c213d5c1af3431cd9484992d..44ab810f3d77bc2ed867c8cfae5caafd3dfddfec 100755
--- a/Mesh/Levy3D.cpp
+++ b/Mesh/Levy3D.cpp
@@ -12,7 +12,7 @@
 
 #include "Levy3D.h"
 #include "polynomialBasis.h"
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GModel.h"
 #include "MElement.h"
 #include "MElementOctree.h"
diff --git a/Numeric/CMakeLists.txt b/Numeric/CMakeLists.txt
index 88cd8089b565425a8b8bea09eb559b87637f9a87..388f7259c01074e44a21f5923a6c7a702a9e5070 100644
--- a/Numeric/CMakeLists.txt
+++ b/Numeric/CMakeLists.txt
@@ -8,14 +8,14 @@ set(SRC
     fullMatrix.cpp
     polynomialBasis.cpp
     JacobianBasis.cpp
-  GaussQuadratureLin.cpp
+  GaussIntegration.cpp
+    GaussQuadratureLin.cpp
     GaussQuadratureTri.cpp
     GaussQuadratureQuad.cpp
     GaussQuadratureTet.cpp
     GaussQuadratureHex.cpp
     GaussQuadraturePri.cpp
     GaussLegendreSimplex.cpp
-    Gauss.cpp
   robustPredicates.cpp
   mathEvaluator.cpp
   Iso.cpp
diff --git a/Numeric/Gauss.cpp b/Numeric/GaussIntegration.cpp
similarity index 98%
rename from Numeric/Gauss.cpp
rename to Numeric/GaussIntegration.cpp
index ca0c88b2e914f7375f94903212698885ae31dad9..9e1f471012dc4ff8b9f20438726d4e4c19ce4caa 100644
--- a/Numeric/Gauss.cpp
+++ b/Numeric/GaussIntegration.cpp
@@ -1,4 +1,4 @@
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GmshDefines.h"
 
 static void pts2fullMatrix(int npts, IntPt *pts, fullMatrix<double> &pMatrix, 
diff --git a/Numeric/Gauss.h b/Numeric/GaussIntegration.h
similarity index 95%
rename from Numeric/Gauss.h
rename to Numeric/GaussIntegration.h
index bc2b35595933ed6c4a3ca193f0d58ba6dc778fdd..099b6b22af8d8734f7721f970a736f0e5c505d07 100644
--- a/Numeric/Gauss.h
+++ b/Numeric/GaussIntegration.h
@@ -3,8 +3,8 @@
 // See the LICENSE.txt file for license information. Please report all
 // bugs and problems to <gmsh@geuz.org>.
 
-#ifndef _GAUSS_H_
-#define _GAUSS_H_
+#ifndef _GAUSS_INTEGRATION_H_
+#define _GAUSS_INTEGRATION_H_
 
 #include "fullMatrix.h"
 
diff --git a/Numeric/GaussLegendreSimplex.cpp b/Numeric/GaussLegendreSimplex.cpp
index aeed112f422f93a09bb3aa5699901c95f5478911..8e24a370e200d50c5956949670afb6bdc684464c 100644
--- a/Numeric/GaussLegendreSimplex.cpp
+++ b/Numeric/GaussLegendreSimplex.cpp
@@ -4,7 +4,7 @@
 // bugs and problems to <gmsh@geuz.org>.
 
 #include <math.h>
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GaussLegendre1D.h"
 
 static void brickToTet(double xi, double eta, double zeta,
diff --git a/Numeric/GaussQuadratureHex.cpp b/Numeric/GaussQuadratureHex.cpp
index 808aae08b0949af58d6b9ba576e57af901f05410..43fee450ec649467116c95601e82b4c4dd9e72ac 100644
--- a/Numeric/GaussQuadratureHex.cpp
+++ b/Numeric/GaussQuadratureHex.cpp
@@ -3,7 +3,7 @@
 // See the LICENSE.txt file for license information. Please report all
 // bugs and problems to <gmsh@geuz.org>.
 
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GaussLegendre1D.h"
 
 const double a1  = 0.40824826;
diff --git a/Numeric/GaussQuadratureLin.cpp b/Numeric/GaussQuadratureLin.cpp
index 3562a03b0b99516a4de2bbe83d3d108e3676df71..e747b4cced54e6f7e69a27622c87fbac9e79663a 100644
--- a/Numeric/GaussQuadratureLin.cpp
+++ b/Numeric/GaussQuadratureLin.cpp
@@ -3,7 +3,7 @@
 // See the LICENSE.txt file for license information. Please report all
 // bugs and problems to <gmsh@geuz.org>.
 
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GaussLegendre1D.h"
 
 IntPt * GQL[20] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
diff --git a/Numeric/GaussQuadraturePri.cpp b/Numeric/GaussQuadraturePri.cpp
index 1a328b750c181ef4c4ae9f86d1b2009fcc4655ff..c7bddb484caae1f2bcf19d4272e1a694d6426835 100644
--- a/Numeric/GaussQuadraturePri.cpp
+++ b/Numeric/GaussQuadraturePri.cpp
@@ -3,7 +3,7 @@
 // See the LICENSE.txt file for license information. Please report all
 // bugs and problems to <gmsh@geuz.org>.
 
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GaussLegendre1D.h"
 #include "stdio.h"
 
diff --git a/Numeric/GaussQuadratureQuad.cpp b/Numeric/GaussQuadratureQuad.cpp
index ddfde8c4f7854640869148d59c5772ac0d571519..6e4987d65c0f33da226ddcaa0c7434ba72fb4d18 100644
--- a/Numeric/GaussQuadratureQuad.cpp
+++ b/Numeric/GaussQuadratureQuad.cpp
@@ -3,7 +3,7 @@
 // See the LICENSE.txt file for license information. Please report all
 // bugs and problems to <gmsh@geuz.org>.
 
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GaussLegendre1D.h"
 
 IntPt GQQ1[1] = {
diff --git a/Numeric/GaussQuadratureTet.cpp b/Numeric/GaussQuadratureTet.cpp
index c30f455cfd8bfe32ffed8b802775ede50776a897..ceb35dc72d11c40b447ed8a90dea7f2160c38725 100644
--- a/Numeric/GaussQuadratureTet.cpp
+++ b/Numeric/GaussQuadratureTet.cpp
@@ -3,7 +3,7 @@
 // See the LICENSE.txt file for license information. Please report all
 // bugs and problems to <gmsh@geuz.org>.
 
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GaussLegendre1D.h"
 
 /* constants for 4-point rule */
diff --git a/Numeric/GaussQuadratureTri.cpp b/Numeric/GaussQuadratureTri.cpp
index 7409586c3f36a4e8d270b0abf39faa6c60c01260..7f3a7feacb853cabd05337e41b53000f2b9ab363 100644
--- a/Numeric/GaussQuadratureTri.cpp
+++ b/Numeric/GaussQuadratureTri.cpp
@@ -3,7 +3,7 @@
 // See the LICENSE.txt file for license information. Please report all
 // bugs and problems to <gmsh@geuz.org>.
 
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "GaussLegendre1D.h"
 
 IntPt GQT1[1] = { 
diff --git a/Numeric/polynomialBasis.cpp b/Numeric/polynomialBasis.cpp
index 18b2caf491d553d135483d7afdc16629baa8c3fc..aa84c6c4e4fa0840436b4b0f4dd8528f1888566a 100644
--- a/Numeric/polynomialBasis.cpp
+++ b/Numeric/polynomialBasis.cpp
@@ -12,7 +12,7 @@
 #include "GmshDefines.h"
 #include "GmshMessage.h"
 #include "polynomialBasis.h"
-#include "Gauss.h"
+#include "GaussIntegration.h"
 
 static void printClosure(polynomialBasis::clCont &fullClosure, std::vector<int> &closureRef,
                          polynomialBasis::clCont &closures)
diff --git a/Solver/quadratureRules.h b/Solver/quadratureRules.h
index d84d4b1dc8b43e3ad26c98d8ac8560d80db557fc..bf5fe5b399541dfb6e5eee1316b70c554b940e17 100644
--- a/Solver/quadratureRules.h
+++ b/Solver/quadratureRules.h
@@ -14,7 +14,7 @@
 #ifndef _QUADRATURERULES_H_
 #define _QUADRATURERULES_H_
 
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "MElement.h"
 
 class QuadratureBase
diff --git a/contrib/DiscreteIntegration/Integration3D.cpp b/contrib/DiscreteIntegration/Integration3D.cpp
index a9f8ed9c36855ab477e3874636c596a4ddb3466b..8219446b601dfe4e7f140807e64c0d1ef445542e 100644
--- a/contrib/DiscreteIntegration/Integration3D.cpp
+++ b/contrib/DiscreteIntegration/Integration3D.cpp
@@ -5,7 +5,7 @@
 
 #include "Integration3D.h"
 #include "recurCut.h"
-#include "Gauss.h"
+#include "GaussIntegration.h"
 #include "polynomialBasis.h"
 #include "GmshDefines.h" 
 
diff --git a/gmshpy/gmshNumeric.i b/gmshpy/gmshNumeric.i
index 77c064d0d287fad0d818f630cdcd39e55ba3f540..6d0a1653e63591412c8e6d0a4f75d5230664f139 100644
--- a/gmshpy/gmshNumeric.i
+++ b/gmshpy/gmshNumeric.i
@@ -8,14 +8,14 @@
   #include "GmshConfig.h"
 
   #include "DivideAndConquer.h"
-  #include "Gauss.h"
+  #include "GaussIntegration.h"
   #include "JacobianBasis.h"
   #include "fullMatrix.h"
   #include "polynomialBasis.h"
 %}
 
 %include "DivideAndConquer.h"
-%include "Gauss.h"
+%include "GaussIntegration.h"
 %include "JacobianBasis.h"  
 %include "fullMatrix.h"
 %template(fullMatrixDouble) fullMatrix<double>;