diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion.prj b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion.prj
index 47ce3430c0dbfc37d99c947c311934c834c0ea9f..a8810ce04a7ea06f8c181a4f3da6cf2133ef11aa 100644
--- a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion.prj
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion.prj
@@ -37,10 +37,10 @@
                             <type>Constant</type>
                             <value>1e-3</value>
                         </property>
-                         <property>
-                             <name>density</name>
-                             <type>Constant</type>
-                             <value>1000</value>
+                        <property>
+                            <name>density</name>
+                            <type>Constant</type>
+                            <value>1000</value>
                         </property>
                         <property>
                             <name>specific_heat_capacity</name>
@@ -335,6 +335,15 @@
             <eigen>
                 <solver_type>SparseLU</solver_type>
             </eigen>
+            <petsc>
+                <parameters>
+                    -ksp_type bcgs
+                    -pc_type lu
+                    -ksp_rtol 1e-8
+                    -ksp_atol 1e-12
+                    -ksp_max_it 4000
+                </parameters>
+            </petsc>
         </linear_solver>
     </linear_solvers>
 </OpenGeoSysProject>
diff --git a/Tests/Python/test_matrix_debug_output.py b/Tests/Python/test_matrix_debug_output.py
index dad4970179de3d0a19ba49c917e80b819e63d01d..9a218ecb21da6eddddf9a62f5bc3872e685fccb7 100644
--- a/Tests/Python/test_matrix_debug_output.py
+++ b/Tests/Python/test_matrix_debug_output.py
@@ -3,6 +3,7 @@ import platform
 import tempfile
 from pathlib import Path
 
+import ogs
 import pytest
 from ogs import cli
 
@@ -154,6 +155,10 @@ def test_local_matrix_debug_output(monkeypatch, prefix_parameter, elements_param
         ("", True, True),  # set env var to an empty string
     ],
 )
+@pytest.mark.skipif(
+    ogs.OGS_USE_PETSC == "ON",
+    reason="global matrix output is not implemented for PETSc matrices",
+)
 def test_global_matrix_debug_output(monkeypatch, prefix_parameter):
     srcdir = Path(__file__).parent.parent.parent
     prjpath = srcdir / "Tests/Data/Mechanics/Linear/square_1e2.prj"