diff --git a/autoemulate/emulators/__init__.py b/autoemulate/emulators/__init__.py
index 868a4f20..01f29542 100644
--- a/autoemulate/emulators/__init__.py
+++ b/autoemulate/emulators/__init__.py
@@ -1,16 +1,16 @@
 from .base import Emulator
 from .gaussian_process import GaussianProcess
 from .gaussian_process_sk import GaussianProcessSk
-from .neural_network import NeuralNetwork
+from .neural_net_sk import NeuralNetSk
 from .random_forest import RandomForest
 from .radial_basis import RadialBasis
-from .neural_net_pt import SkorchMLPRegressor
+from .neural_net_torch import NeuralNetTorch
 
 MODEL_REGISTRY = {
     # "GaussianProcess": GaussianProcess,
     "GaussianProcessSk": GaussianProcessSk,
-    "NeuralNetwork": NeuralNetwork,
+    "NeuralNetSk": NeuralNetSk,
     "RandomForest": RandomForest,
     "RadialBasis": RadialBasis,
-    "SkorchMLPRegressor": SkorchMLPRegressor,
+    # "NeuralNetTorch": NeuralNetTorch,
 }
diff --git a/autoemulate/emulators/gaussian_process.py b/autoemulate/emulators/gaussian_process.py
index 5271d3f6..25ddec3c 100644
--- a/autoemulate/emulators/gaussian_process.py
+++ b/autoemulate/emulators/gaussian_process.py
@@ -33,7 +33,7 @@ def fit(self, X, y):
         X, y = check_X_y(X, y, multi_output=False, y_numeric=True)
         self.n_features_in_ = X.shape[1]
         self.model_ = mogp_emulator.GaussianProcess(X, y, nugget=self.nugget)
-        self.model_ = mogp_emulator.fit_GP_MAP(self.model_, n_tries=2)
+        self.model_ = mogp_emulator.fit_GP_MAP(self.model_, n_tries=15)
         self.is_fitted_ = True
         return self
 
diff --git a/autoemulate/emulators/neural_network.py b/autoemulate/emulators/neural_net_sk.py
similarity index 98%
rename from autoemulate/emulators/neural_network.py
rename to autoemulate/emulators/neural_net_sk.py
index 79dfc563..ee4217db 100644
--- a/autoemulate/emulators/neural_network.py
+++ b/autoemulate/emulators/neural_net_sk.py
@@ -7,7 +7,7 @@
 from sklearn.utils.validation import check_X_y, check_array, check_is_fitted
 
 
-class NeuralNetwork(BaseEstimator, RegressorMixin):
+class NeuralNetSk(BaseEstimator, RegressorMixin):
     """Multi-layer perceptron Emulator.
 
     Implements MLPRegressor from scikit-learn.
diff --git a/autoemulate/emulators/neural_net_pt.py b/autoemulate/emulators/neural_net_torch.py
similarity index 91%
rename from autoemulate/emulators/neural_net_pt.py
rename to autoemulate/emulators/neural_net_torch.py
index 8395a760..17e69f7d 100644
--- a/autoemulate/emulators/neural_net_pt.py
+++ b/autoemulate/emulators/neural_net_torch.py
@@ -1,10 +1,12 @@
+# experimental version of a PyTorch neural network emulator wrapped in Skorch
+# to make it compatible with scikit-learn. Works with cross_validate and GridSearchCV,
+# but doesn't pass tests, because we're subclassing
+
 import torch
 import numpy as np
 import skorch
 from torch import nn
 from skorch import NeuralNetRegressor
-from sklearn.base import BaseEstimator, RegressorMixin
-from sklearn.utils.validation import check_X_y, check_array, check_is_fitted
 
 
 class InputShapeSetter(skorch.callbacks.Callback):
@@ -40,7 +42,7 @@ def forward(self, X):
 
 
 # Step 2: Create the Skorch wrapper for the NeuralNetRegressor
-class SkorchMLPRegressor(NeuralNetRegressor):
+class NeuralNetTorch(NeuralNetRegressor):
     def __init__(
         self,
         module=MLPModule,
diff --git a/tests/test_estimators.py b/tests/test_estimators.py
index 22a5326e..431379b2 100644
--- a/tests/test_estimators.py
+++ b/tests/test_estimators.py
@@ -5,6 +5,7 @@
     NeuralNetwork,
     GaussianProcess,
     RadialBasis,
+    SkorchMLPRegressor,
 )
 from functools import partial
 
@@ -12,9 +13,10 @@
 @parametrize_with_checks(
     [  # GaussianProcess(),
         RandomForest(random_state=42),
-        GaussianProcessSk(random_state=1337),
-        NeuralNetwork(random_state=13),
+        # GaussianProcessSk(random_state=1337),
+        # NeuralNetwork(random_state=13),
         # RadialBasis(),
+        SkorchMLPRegressor(random_state=42),
     ]
 )
 def test_check_estimator(estimator, check):