diff --git a/linear_operator/settings.py b/linear_operator/settings.py
index 14c9107b..65ccf628 100644
--- a/linear_operator/settings.py
+++ b/linear_operator/settings.py
@@ -331,6 +331,7 @@ class fast_computations:
     .. _GPyTorch Blackbox Matrix-Matrix Gaussian Process Inference with GPU Acceleration:
         https://arxiv.org/pdf/1809.11165.pdf
     """
+
     covar_root_decomposition = _fast_covar_root_decomposition
     log_prob = _fast_log_prob
     solves = _fast_solves
@@ -394,10 +395,10 @@ class max_cholesky_size(_value_context):
     If the size of of a LinearOperator is less than `max_cholesky_size`,
     then `root_decomposition` and `solve` of LinearOperator will use Cholesky rather than Lanczos/CG.
 
-    (Default: 800)
+    (Default: 50000)
     """
 
-    _global_value = 800
+    _global_value = 50_000
 
 
 class max_lanczos_quadrature_iterations(_value_context):