]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add a testcase.
authorWolfgang Bangerth <bangerth@colostate.edu>
Tue, 6 May 2025 20:38:52 +0000 (14:38 -0600)
committerWolfgang Bangerth <bangerth@colostate.edu>
Wed, 7 May 2025 16:33:20 +0000 (10:33 -0600)
tests/sundials/kinsol_06_v3.cc [new file with mode: 0644]
tests/sundials/kinsol_06_v3.output [new file with mode: 0644]
tests/sundials/kinsol_06_v3_input.prm [new file with mode: 0644]

diff --git a/tests/sundials/kinsol_06_v3.cc b/tests/sundials/kinsol_06_v3.cc
new file mode 100644 (file)
index 0000000..97e0afd
--- /dev/null
@@ -0,0 +1,102 @@
+// ------------------------------------------------------------------------
+//
+// SPDX-License-Identifier: LGPL-2.1-or-later
+// Copyright (C) 2017 - 2023 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// Part of the source code is dual licensed under Apache-2.0 WITH
+// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
+// governing the source code and code contributions can be found in
+// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
+//
+// ------------------------------------------------------------------------
+
+
+#include <deal.II/base/parameter_handler.h>
+
+#include <deal.II/lac/vector.h>
+
+#include <deal.II/sundials/kinsol.h>
+
+#include "../tests.h"
+
+
+// Like the kinsol_06 test, but with a case where KINSOL, after a few
+// tries with the given residual/Jacobian, gives up and terminates
+// because it can't find a solution. We need to make sure that we
+// throw a catchable exception.
+//
+// This testcase is a variation of the kinsol_06 test, modified by
+// Simon Wiesheier, and posted on the mailing list. Then further
+// adapted by Wolfgang Bangerth.
+
+
+int
+main()
+{
+  initlog();
+
+  using VectorType = Vector<double>;
+
+  // Size of the problem
+  const unsigned int N = 1;
+
+  SUNDIALS::KINSOL<VectorType>::AdditionalData data;
+  ParameterHandler                             prm;
+  data.add_parameters(prm);
+
+  std::ifstream ifile(SOURCE_DIR "/kinsol_06_v3_input.prm");
+  prm.parse_input(ifile);
+
+  SUNDIALS::KINSOL<VectorType> kinsol(data);
+
+  kinsol.reinit_vector = [N](VectorType &v) { v.reinit(N); };
+
+  kinsol.residual = [&](const VectorType &u, VectorType &F) {
+    // Count how many times this function has been called. Let it call
+    // *every* time after the first attempt:
+    static int count = 0;
+    deallog << "Computing residual for the " << count + 1
+            << "th time, at u=" << u[0] << std::endl;
+    if ((u[0] < -10) || (u[0] > 20))
+      {
+        deallog << "Reporting recoverable failure." << std::endl;
+        throw RecoverableUserCallbackError();
+      }
+    count++;
+
+
+    F.reinit(u);
+    F[0] = std::atan(u[0]) - 0.5;
+  };
+
+  double J_inverse;
+
+  kinsol.setup_jacobian = [&J_inverse](const VectorType &u,
+                                       const VectorType &F) {
+    deallog << "Setting up Jacobian system at u=" << u[0] << std::endl;
+
+    const double J = 1. / (1 + u[0] * u[0]);
+    J_inverse      = 1. / J;
+  };
+
+
+  kinsol.solve_with_jacobian = [&](const VectorType &rhs,
+                                   VectorType       &dst,
+                                   double) { dst[0] = J_inverse * rhs[0]; };
+
+  VectorType v(N);
+  v[0] = 10;
+
+  try
+    {
+      auto niter = kinsol.solve(v);
+    }
+  catch (const ExceptionBase &e)
+    {
+      deallog << "KINSOL threw an exception with the following message:"
+              << std::endl;
+      e.print_info(deallog.get_file_stream());
+    }
+}
diff --git a/tests/sundials/kinsol_06_v3.output b/tests/sundials/kinsol_06_v3.output
new file mode 100644 (file)
index 0000000..dda6d57
--- /dev/null
@@ -0,0 +1,23 @@
+
+DEAL::Computing residual for the 1th time, at u=10.0000
+DEAL::Setting up Jacobian system at u=10.0000
+DEAL::Computing residual for the 2th time, at u=10.0000
+DEAL::Computing residual for the 3th time, at u=-88.0839
+DEAL::Reporting recoverable failure.
+DEAL::Computing residual for the 3th time, at u=-39.0419
+DEAL::Reporting recoverable failure.
+DEAL::Computing residual for the 3th time, at u=-14.5210
+DEAL::Reporting recoverable failure.
+DEAL::Computing residual for the 3th time, at u=-2.26049
+DEAL::Setting up Jacobian system at u=-2.26049
+DEAL::Computing residual for the 4th time, at u=-2.26049
+DEAL::Computing residual for the 5th time, at u=7.84693
+DEAL::Computing residual for the 6th time, at u=7.84693
+DEAL::Computing residual for the 7th time, at u=2.07902
+DEAL::Setting up Jacobian system at u=2.07902
+DEAL::Computing residual for the 8th time, at u=2.07902
+DEAL::Computing residual for the 9th time, at u=-1.23396
+DEAL::Computing residual for the 10th time, at u=-1.23396
+DEAL::Computing residual for the 11th time, at u=6.16274
+DEAL::KINSOL threw an exception with the following message:
+    One of the SUNDIALS KINSOL internal functions returned a negative error code: -6. Please consult SUNDIALS manual.
diff --git a/tests/sundials/kinsol_06_v3_input.prm b/tests/sundials/kinsol_06_v3_input.prm
new file mode 100644 (file)
index 0000000..fd9f762
--- /dev/null
@@ -0,0 +1,19 @@
+set Function norm stopping tolerance       = 0.000000
+set Maximum number of nonlinear iterations = 5
+set Scaled step stopping tolerance         = 0.000000
+set Solution strategy                      = newton
+
+subsection Fixed point and Picard parameters
+  set Anderson acceleration subspace size = 5
+end
+
+subsection Linesearch parameters
+  set Maximum number of beta-condition failures = 0
+end
+
+subsection Newton parameters
+  set Maximum allowable scaled length of the Newton step = 0.000000
+  set Maximum iterations without matrix setup            = 0
+  set No initial matrix setup                            = false
+  set Relative error for different quotient computation  = 0.000000
+end

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.