diff --git a/lbmpy/methods/__init__.py b/lbmpy/methods/__init__.py
index 9fa48fdd027109c3667308021fbe296cbac3c797..3be0dbc3e7bf8970d6d4cb494e7d5ca0443cbf5e 100644
--- a/lbmpy/methods/__init__.py
+++ b/lbmpy/methods/__init__.py
@@ -5,7 +5,7 @@ from lbmpy.methods.creationfunctions import (
     create_centered_cumulant_model, create_with_default_polynomial_cumulants,
     create_with_polynomial_cumulants, create_with_monomial_cumulants)
 
-from lbmpy.methods.default_moment_sets import mrt_orthogonal_modes_literature
+from lbmpy.methods.default_moment_sets import mrt_orthogonal_modes_literature, cascaded_moment_sets_literature
 
 from lbmpy.methods.abstractlbmethod import AbstractLbMethod, RelaxationInfo
 from lbmpy.methods.conservedquantitycomputation import AbstractConservedQuantityComputation
@@ -17,6 +17,6 @@ __all__ = ['RelaxationInfo', 'AbstractLbMethod',
            'create_srt', 'create_trt', 'create_trt_with_magic_number', 'create_trt_kbc',
            'create_mrt_orthogonal', 'create_mrt_raw', 'create_central_moment',
            'create_with_continuous_maxwellian_eq_moments', 'create_with_discrete_maxwellian_eq_moments',
-           'mrt_orthogonal_modes_literature', 'create_centered_cumulant_model',
-           'create_with_default_polynomial_cumulants', 'create_with_polynomial_cumulants',
-           'create_with_monomial_cumulants']
+           'mrt_orthogonal_modes_literature', 'cascaded_moment_sets_literature',
+           'create_centered_cumulant_model', 'create_with_default_polynomial_cumulants',
+           'create_with_polynomial_cumulants', 'create_with_monomial_cumulants']