diff --git a/src/qiskit_symb/circuit/library/__init__.py b/src/qiskit_symb/circuit/library/__init__.py
index 6846512..b96f385 100644
--- a/src/qiskit_symb/circuit/library/__init__.py
+++ b/src/qiskit_symb/circuit/library/__init__.py
@@ -51,7 +51,8 @@
     'ctdg': CTdgGate,
     'u': UGate,
     'u3': UGate,
-    'cu': CUGate,
+    # https://github.com/SimoneGasperini/qiskit-symb/issues/4
+    # 'cu': CUGate,
     'cu3': CUGate,
     'rx': RXGate,
     'crx': CRXGate,
diff --git a/src/qiskit_symb/utils.py b/src/qiskit_symb/utils.py
index 24311ad..97de1e1 100644
--- a/src/qiskit_symb/utils.py
+++ b/src/qiskit_symb/utils.py
@@ -11,7 +11,7 @@ def get_init(name):
     """todo"""
     if name not in NAME_TO_INIT:
         raise NotImplementedError(
-            f'Gate "{name}" is not implemented in qiskit-symb')
+            f'Gate "{name}" is not implemented in qiskit-symb, use Qiskit transpiler!')
     return NAME_TO_INIT[name]
 
 
diff --git a/tests/__init__.py b/tests/__init__.py
new file mode 100644
index 0000000..e69de29
diff --git a/tests/test_ctrl_parametric_gates.py b/tests/test_ctrl_parametric_gates.py
index 61863c2..f068d14 100644
--- a/tests/test_ctrl_parametric_gates.py
+++ b/tests/test_ctrl_parametric_gates.py
@@ -5,7 +5,7 @@
 from qiskit.circuit import ParameterVector, Parameter
 from qiskit.quantum_info import Operator
 from qiskit.circuit.library import (
-    UGate, RXGate, RYGate, RZGate,
+    RXGate, RYGate, RZGate,
     PhaseGate, RGate,
     RXXGate, RYYGate,
     RZZGate, RZXGate
@@ -15,14 +15,14 @@
 
 val_range = {'min_value': -2*numpy.pi, 'max_value': 2*numpy.pi}
 
-
+# https://github.com/SimoneGasperini/qiskit-symb/issues/4
+"""
 @settings(deadline=None, max_examples=10)
 @given(theta=strategies.floats(**val_range),
        phi=strategies.floats(**val_range),
        lam=strategies.floats(**val_range),
        seed=strategies.integers(min_value=0))
 def test_cu(theta, phi, lam, seed):
-    """todo"""
     pars_vals = [theta, phi, lam]
     pars = ParameterVector(name='pars', length=len(pars_vals))
     circuit = get_random_controlled(base_gate=UGate(*pars), seed=seed)
@@ -34,6 +34,7 @@ def test_cu(theta, phi, lam, seed):
     arr2 = symb_Operator(
         circuit).subs({pars: pars_vals}).to_numpy()
     assert numpy.allclose(arr1, arr2)
+"""
 
 
 @settings(deadline=None, max_examples=10)