diff --git a/python/README.md b/python/README.md
index 182211d7a..a6eaed000 100644
--- a/python/README.md
+++ b/python/README.md
@@ -70,7 +70,7 @@ native python functions as callbacks. This is a simple example that loads
 two meshes and switches between the two when a key is pressed:
 
 ```python
-import igl
+import pyigl as igl
 
 V1 = igl.eigen.MatrixXd()
 F1 = igl.eigen.MatrixXi()
diff --git a/python/iglhelpers.py b/python/iglhelpers.py
index c5fcaeb26..8fa399718 100644
--- a/python/iglhelpers.py
+++ b/python/iglhelpers.py
@@ -1,6 +1,6 @@
 import numpy as np
 import scipy.sparse as sparse
-import igl
+import pyigl as igl
 
 def p2e(m):
     if isinstance(m, np.ndarray):
diff --git a/python/python_shared.cpp b/python/python_shared.cpp
index 7a154fd47..e541d5f67 100644
--- a/python/python_shared.cpp
+++ b/python/python_shared.cpp
@@ -14,9 +14,9 @@ extern void python_export_igl_viewer(py::module &);
 extern void python_export_igl_comiso(py::module &);
 #endif
 
-PYTHON_PLUGIN(igl) {
+PYTHON_PLUGIN(pyigl) {
     py::init_threading();
-    py::module m("igl", "Python wrappers for libigl");
+    py::module m("pyigl", "Python wrappers for libigl");
 
     python_export_vector(m);
     python_export_igl(m);
diff --git a/python/tcpviewer.py b/python/tcpviewer.py
index b89ec536c..f0743754a 100644
--- a/python/tcpviewer.py
+++ b/python/tcpviewer.py
@@ -1,6 +1,6 @@
 import socket
 import threading
-import igl
+import pyigl as igl
 import array
 import time
 
diff --git a/python/tutorial/101_FileIO.py b/python/tutorial/101_FileIO.py
index 17818444e..5553e51bd 100644
--- a/python/tutorial/101_FileIO.py
+++ b/python/tutorial/101_FileIO.py
@@ -3,7 +3,7 @@ from __future__ import print_function
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 # Load a mesh in OFF format
 V = igl.eigen.MatrixXd()
diff --git a/python/tutorial/102_DrawMesh.py b/python/tutorial/102_DrawMesh.py
index 1b1a0d1d1..15e40d86a 100755
--- a/python/tutorial/102_DrawMesh.py
+++ b/python/tutorial/102_DrawMesh.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 # Load a mesh in OFF format
 V = igl.eigen.MatrixXd()
diff --git a/python/tutorial/102_DrawMesh_TCP.py b/python/tutorial/102_DrawMesh_TCP.py
index e63987d54..60235e4ee 100644
--- a/python/tutorial/102_DrawMesh_TCP.py
+++ b/python/tutorial/102_DrawMesh_TCP.py
@@ -13,7 +13,7 @@ os.system("python ../tcpviewer.py&")
 # Wait for it to set up the socket
 time.sleep(1)
 
-import igl
+import pyigl as igl
 import tcpviewer
 
 # Read a mesh
diff --git a/python/tutorial/103_Events.py b/python/tutorial/103_Events.py
index a54c7671c..0b0462166 100755
--- a/python/tutorial/103_Events.py
+++ b/python/tutorial/103_Events.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V1 = igl.eigen.MatrixXd()
 F1 = igl.eigen.MatrixXi()
diff --git a/python/tutorial/104_Colors.py b/python/tutorial/104_Colors.py
index 26bd9507a..c850d2e78 100755
--- a/python/tutorial/104_Colors.py
+++ b/python/tutorial/104_Colors.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/105_Overlays.py b/python/tutorial/105_Overlays.py
index 680a52f7b..6bca2b39a 100755
--- a/python/tutorial/105_Overlays.py
+++ b/python/tutorial/105_Overlays.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/201_Normals.py b/python/tutorial/201_Normals.py
index 684ada726..b78a1577b 100755
--- a/python/tutorial/201_Normals.py
+++ b/python/tutorial/201_Normals.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/202_GaussianCurvature.py b/python/tutorial/202_GaussianCurvature.py
index 45f9a1206..189bd7db1 100644
--- a/python/tutorial/202_GaussianCurvature.py
+++ b/python/tutorial/202_GaussianCurvature.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 # Load mesh
 V = igl.eigen.MatrixXd()
diff --git a/python/tutorial/203_CurvatureDirections.py b/python/tutorial/203_CurvatureDirections.py
index bbe0386bf..27b98e9e8 100755
--- a/python/tutorial/203_CurvatureDirections.py
+++ b/python/tutorial/203_CurvatureDirections.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd();
 F = igl.eigen.MatrixXi();
diff --git a/python/tutorial/204_Gradient.py b/python/tutorial/204_Gradient.py
index 873f3462f..6e888352c 100755
--- a/python/tutorial/204_Gradient.py
+++ b/python/tutorial/204_Gradient.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/205_Laplacian.py b/python/tutorial/205_Laplacian.py
index 4a5e08080..2614ac26d 100755
--- a/python/tutorial/205_Laplacian.py
+++ b/python/tutorial/205_Laplacian.py
@@ -4,7 +4,7 @@ from __future__ import print_function
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 import math
 
 global V
diff --git a/python/tutorial/301_Slice.py b/python/tutorial/301_Slice.py
index 28df4011e..879b44b76 100644
--- a/python/tutorial/301_Slice.py
+++ b/python/tutorial/301_Slice.py
@@ -4,7 +4,7 @@ from __future__ import print_function
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/302_Sort.py b/python/tutorial/302_Sort.py
index 128e1672e..e0ba88c02 100644
--- a/python/tutorial/302_Sort.py
+++ b/python/tutorial/302_Sort.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/303_LaplaceEquation.py b/python/tutorial/303_LaplaceEquation.py
index d700009fa..1cb891e64 100644
--- a/python/tutorial/303_LaplaceEquation.py
+++ b/python/tutorial/303_LaplaceEquation.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/304_LinearEqualityConstraints.py b/python/tutorial/304_LinearEqualityConstraints.py
index bbff1d372..442c7f780 100644
--- a/python/tutorial/304_LinearEqualityConstraints.py
+++ b/python/tutorial/304_LinearEqualityConstraints.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/305_QuadraticProgramming.py b/python/tutorial/305_QuadraticProgramming.py
index 87a6453e2..3f9ff1039 100644
--- a/python/tutorial/305_QuadraticProgramming.py
+++ b/python/tutorial/305_QuadraticProgramming.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 b = igl.eigen.MatrixXi()
 B = igl.eigen.MatrixXd()
diff --git a/python/tutorial/306_EigenDecomposition.py b/python/tutorial/306_EigenDecomposition.py
index 3a7bbd565..b3bd3ab97 100644
--- a/python/tutorial/306_EigenDecomposition.py
+++ b/python/tutorial/306_EigenDecomposition.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 U = igl.eigen.MatrixXd()
diff --git a/python/tutorial/401_BiharmonicDeformation.py b/python/tutorial/401_BiharmonicDeformation.py
index 0b1cf75e7..68f6a0aac 100755
--- a/python/tutorial/401_BiharmonicDeformation.py
+++ b/python/tutorial/401_BiharmonicDeformation.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 global bc_frac, bc_dir,deformation_field, V, U, V_bc, U_bc, F, b
 bc_frac = 1.0
diff --git a/python/tutorial/402_PolyharmonicDeformation.py b/python/tutorial/402_PolyharmonicDeformation.py
index c991c91d4..1deacba13 100755
--- a/python/tutorial/402_PolyharmonicDeformation.py
+++ b/python/tutorial/402_PolyharmonicDeformation.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 global z_max, z_dir, k, resolve, V, U, Z, F, b, bc
 
diff --git a/python/tutorial/405_AsRigidAsPossible.py b/python/tutorial/405_AsRigidAsPossible.py
index 383ed6580..10cceb390 100755
--- a/python/tutorial/405_AsRigidAsPossible.py
+++ b/python/tutorial/405_AsRigidAsPossible.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 from math import sin,cos,pi
 
 sea_green = igl.eigen.MatrixXd([[70./255.,252./255.,167./255.]])
diff --git a/python/tutorial/501_HarmonicParam.py b/python/tutorial/501_HarmonicParam.py
index 0054df88d..50ab71673 100755
--- a/python/tutorial/501_HarmonicParam.py
+++ b/python/tutorial/501_HarmonicParam.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/502_LSCMParam.py b/python/tutorial/502_LSCMParam.py
index fca9ae35e..173c10b06 100755
--- a/python/tutorial/502_LSCMParam.py
+++ b/python/tutorial/502_LSCMParam.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/503_ARAPParam.py b/python/tutorial/503_ARAPParam.py
index 757a3816b..df93fe8c2 100755
--- a/python/tutorial/503_ARAPParam.py
+++ b/python/tutorial/503_ARAPParam.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 
 V = igl.eigen.MatrixXd()
 F = igl.eigen.MatrixXi()
diff --git a/python/tutorial/504_NRosyDesign.py b/python/tutorial/504_NRosyDesign.py
index b9fa50579..60993ee5c 100755
--- a/python/tutorial/504_NRosyDesign.py
+++ b/python/tutorial/504_NRosyDesign.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 from math import atan2,pi,cos,sin
 
 # Mesh
diff --git a/python/tutorial/505_MIQ.py b/python/tutorial/505_MIQ.py
index 7aa2db8d4..4b2c87a51 100755
--- a/python/tutorial/505_MIQ.py
+++ b/python/tutorial/505_MIQ.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 from math import pi
 
 V = igl.eigen.MatrixXd()
diff --git a/python/tutorial/507_PolyVectorField.py b/python/tutorial/507_PolyVectorField.py
index c4267de37..a7904c1ba 100755
--- a/python/tutorial/507_PolyVectorField.py
+++ b/python/tutorial/507_PolyVectorField.py
@@ -2,7 +2,7 @@
 import sys, os
 sys.path.insert(0, os.getcwd() + "/../")
 
-import igl
+import pyigl as igl
 import random
 from math import cos,sin,pi