From 4565bbf9df7c8b5e5f3499202dfa366a38274f08 Mon Sep 17 00:00:00 2001
From: Quang Hoang <qxhmorrow@gmail.com>
Date: Mon, 14 Oct 2024 21:23:27 -0400
Subject: [PATCH] add burstedde weighted bfbt test

---
 tests/burstedde_bfbt.cc            | 21 ++++++++++++++++++++
 tests/burstedde_bfbt.prm           | 32 ++++++++++++++++++++++++++++++
 tests/burstedde_bfbt/screen-output | 25 +++++++++++++++++++++++
 3 files changed, 78 insertions(+)
 create mode 100644 tests/burstedde_bfbt.cc
 create mode 100644 tests/burstedde_bfbt.prm
 create mode 100644 tests/burstedde_bfbt/screen-output

diff --git a/tests/burstedde_bfbt.cc b/tests/burstedde_bfbt.cc
new file mode 100644
index 00000000000..f43d3026353
--- /dev/null
+++ b/tests/burstedde_bfbt.cc
@@ -0,0 +1,21 @@
+/*
+  Copyright (C) 2022 by the authors of the ASPECT code.
+
+  This file is part of ASPECT.
+
+  ASPECT is free software; you can redistribute it and/or modify
+  it under the terms of the GNU General Public License as published by
+  the Free Software Foundation; either version 2, or (at your option)
+  any later version.
+
+  ASPECT is distributed in the hope that it will be useful,
+  but WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+  GNU General Public License for more details.
+
+  You should have received a copy of the GNU General Public License
+  along with ASPECT; see the file LICENSE.  If not see
+  <http://www.gnu.org/licenses/>.
+*/
+
+#include "../benchmarks/burstedde/burstedde.cc"
diff --git a/tests/burstedde_bfbt.prm b/tests/burstedde_bfbt.prm
new file mode 100644
index 00000000000..6c05cdbef72
--- /dev/null
+++ b/tests/burstedde_bfbt.prm
@@ -0,0 +1,32 @@
+# Test burstedde benchmark using BFBT preconditioner
+
+include $ASPECT_SOURCE_DIR/benchmarks/burstedde/burstedde.prm
+
+set Dimension=3
+
+subsection Mesh refinement
+  set Initial adaptive refinement        = 0
+  set Initial global refinement          = 2
+  set Refinement fraction                = 0.2
+  set Strategy                           = velocity
+end
+
+
+subsection Solver parameters
+subsection Stokes solver parameters
+    set Stokes solver type= block AMG
+    set Use full A block as preconditioner              = true
+    set Number of cheap Stokes solver steps             = 500
+    set Maximum number of expensive Stokes solver steps = 0
+    set Linear solver tolerance                         = 1e-6
+    set GMRES solver restart length                     = 100
+    set Use weighted BFBT for Schur complement          = true
+    set Krylov method for cheap solver steps            = GMRES
+  end
+
+  subsection AMG parameters
+    set AMG aggregation threshold = 0.02
+  end
+
+
+end
diff --git a/tests/burstedde_bfbt/screen-output b/tests/burstedde_bfbt/screen-output
new file mode 100644
index 00000000000..54ab302ef19
--- /dev/null
+++ b/tests/burstedde_bfbt/screen-output
@@ -0,0 +1,25 @@
+
+Loading shared library <./libburstedde_bfbt.debug.so>
+
+Number of active cells: 64 (on 3 levels)
+Number of degrees of freedom: 3,041 (2,187+125+729)
+
+*** Timestep 0:  t=0 seconds, dt=0 seconds
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system (AMG-BFBT)... 159+0 iterations.
+      Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1
+
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system (AMG-BFBT)... 0+0 iterations.
+      Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 8.87564e-07
+
+
+   Postprocessing:
+     Writing graphical output:      output-burstedde_bfbt/solution/solution-00000
+     RMS, max velocity:             4.3 m/s, 13.2 m/s
+     Errors u_L1, p_L1, u_L2, p_L2: 6.030147e-01, 1.447579e-03, 1.059175e+00, 3.017047e-03
+
+Termination requested by criterion: end time
+
+
+