diff --git a/test/multidomain/facet/1p_1p/analytical/main.cc b/test/multidomain/facet/1p_1p/analytical/main.cc
index 403393e9f7126aab78b3ba797ee6ce7db31a1df7..749f84704831ee40db1b64cea6945de57492d33e 100644
--- a/test/multidomain/facet/1p_1p/analytical/main.cc
+++ b/test/multidomain/facet/1p_1p/analytical/main.cc
@@ -334,7 +334,7 @@ int main(int argc, char** argv) try
                                                   couplingManager);
 
     // the linear solver
-    using LinearSolver = UMFPackBackend;
+    using LinearSolver = BlockDiagAMGBiCGSTABSolver;
     auto linearSolver = std::make_shared<LinearSolver>();
 
     // the non-linear solver
diff --git a/test/multidomain/facet/1p_1p/analytical/params.input b/test/multidomain/facet/1p_1p/analytical/params.input
index 2e045e5d4a049292efcdf822a537ed73db2c8189..81c8acbd5bb99a6cd6f7da7c027ffb5109a63020 100644
--- a/test/multidomain/facet/1p_1p/analytical/params.input
+++ b/test/multidomain/facet/1p_1p/analytical/params.input
@@ -19,6 +19,10 @@ Problem.Name = lowdim
 LiquidDensity = 1.0
 LiquidKinematicViscosity = 1.0
 
+[LinearSolver]
+ResidualReduction = 1e-20
+MaxIterations = 2000
+
 [L2Error]
 QuadratureOrder = 1