diff --git a/shenfun/legendre/bases.py b/shenfun/legendre/bases.py
index b2e4b7f1..cc57218c 100644
--- a/shenfun/legendre/bases.py
+++ b/shenfun/legendre/bases.py
@@ -50,7 +50,6 @@
 from shenfun.utilities import n
 from shenfun.jacobi import JacobiBase
 from .lobatto import legendre_lobatto_nodes_and_weights
-from . import fastgl
 
 
 bases = ['Orthogonal',
@@ -120,6 +119,7 @@ def family():
         return 'legendre'
 
     def points_and_weights(self, N=None, map_true_domain=False, weighted=True, **kw):
+        from . import fastgl
         if N is None:
             N = self.shape(False)
         if self.quad == "LG":
diff --git a/shenfun/legendre/dlt.py b/shenfun/legendre/dlt.py
index 4fbcb79a..896c8d90 100644
--- a/shenfun/legendre/dlt.py
+++ b/shenfun/legendre/dlt.py
@@ -10,7 +10,6 @@
 from shenfun.optimization import cython
 from shenfun.spectralbase import islicedict, slicedict
 from shenfun.forms.arguments import FunctionSpace
-from . import fastgl
 
 __all__ = ['DLT', 'leg2cheb', 'cheb2leg', 'Leg2chebHaleTownsend',
            'Leg2Cheb', 'Cheb2Leg', 'FMMLeg2Cheb', 'FMMCheb2Leg']
@@ -100,6 +99,7 @@ class DLT:
     """
     def __init__(self, input_array, s=None, axes=(-1,), threads=1, kind='backward',
                  flags=(FFTW_MEASURE, FFTW_PRESERVE_INPUT), output_array=None):
+        from . import fastgl
         if isinstance(axes, tuple):
             assert len(axes) == 1
             axis = axes[-1]