diff --git a/validation/scripts/test_anaklasis.py b/validation/scripts/test_anaklasis.py
index 0103571..f8d719a 100644
--- a/validation/scripts/test_anaklasis.py
+++ b/validation/scripts/test_anaklasis.py
@@ -2,7 +2,6 @@
 
 import numpy as np
 from test_discovery import get_test_data
-from anaklasis import ref
 
 
 tests = list(get_test_data())
@@ -56,6 +55,8 @@ def anaklasis_layer_matrix(slabs):
 
 def kernel_test(slabs, data):
     # test unsmeared reflectivity calculation
+    from anaklasis import ref
+
     q = data[:, 0]
     R = data[:, 1]
 
@@ -79,6 +80,8 @@ def kernel_test(slabs, data):
 
 def resolution_test(slabs, data):
     # test unsmeared reflectivity calculation
+    from anaklasis import ref
+
     q = data[:, 0]
     R = data[:, 1]
     # dq is a standard deviation but with anaklasis dQ represents the FWHM
diff --git a/validation/scripts/test_bornagain.py b/validation/scripts/test_bornagain.py
index 030878d..8cb7f58 100644
--- a/validation/scripts/test_bornagain.py
+++ b/validation/scripts/test_bornagain.py
@@ -2,14 +2,14 @@
 import numpy as np
 from test_discovery import get_test_data
 
-import bornagain as ba
-from bornagain import angstrom
-
 
 def get_sample(slabs):
     """
     Defines sample and returns it. Note that SLD-based materials are used.
     """
+    import bornagain as ba
+    from bornagain import angstrom
+
     # creating materials
     multi_layer = ba.MultiLayer()
 
@@ -43,6 +43,8 @@ def get_simulation(qzs):
     with a qz-defined beam
     """
     # bornagain requires Qz in nm
+    import bornagain as ba
+
     scan = ba.QzScan(qzs * 10.0)
     simulation = ba.SpecularSimulation()
     simulation.setScan(scan)
@@ -55,6 +57,8 @@ def get_simulation_smeared(qzs, dqzs):
     with a qz-defined beam
     """
     # 3.5 sigma to sync with refnx
+    import bornagain as ba
+
     n_sig = 3.5
     n_samples = 21
     distr = ba.RangedDistributionGaussian(n_samples, n_sig)