From 0edfbafa8286a60a3ad9efad2a7fbc5f33ed59f2 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Tue, 9 May 2023 10:11:03 -0600 Subject: [PATCH] Add a test with an irrecoverable failure of KINSOL. --- tests/sundials/kinsol_07.cc | 99 +++++++++++++++++++++++++++++++++ tests/sundials/kinsol_07.output | 17 ++++++ 2 files changed, 116 insertions(+) create mode 100644 tests/sundials/kinsol_07.cc create mode 100644 tests/sundials/kinsol_07.output diff --git a/tests/sundials/kinsol_07.cc b/tests/sundials/kinsol_07.cc new file mode 100644 index 0000000000..b45c69da18 --- /dev/null +++ b/tests/sundials/kinsol_07.cc @@ -0,0 +1,99 @@ +//----------------------------------------------------------- +// +// Copyright (C) 2017 - 2021 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE.md at +// the top level directory of deal.II. +// +//----------------------------------------------------------- + +#include + +#include +#include + +#include + +#include "../tests.h" + +// This is a variation of the kinsol_06 test, which throws +// "recoverable" errors and forces KINSOL to backtrack. This test +// instead throws an irrecoverable error and forces KINSOL to +// abort. The test checks that the thrown exception is re-thrown after +// KINSOL returns. + + +int +main() +{ + initlog(); + + using VectorType = Vector; + + // Size of the problem + const unsigned int N = 1; + + SUNDIALS::KINSOL::AdditionalData data; + ParameterHandler prm; + data.add_parameters(prm); + + std::ifstream ifile(SOURCE_DIR "/kinsol_newton.prm"); + prm.parse_input(ifile); + + SUNDIALS::KINSOL kinsol(data); + + kinsol.reinit_vector = [N](VectorType &v) { v.reinit(N); }; + + kinsol.residual = [&](const VectorType &u, VectorType &F) { + deallog << "Computing residual at " << u[0] << std::endl; + + if ((u[0] < -10) || (u[0] > 20)) + { + deallog << "Reporting irrecoverable failure." << std::endl; + throw ExcMessage("Irrecoverable failure."); + } + + + 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); + + deallog << "Found solution at " << std::flush; + v.print(deallog.get_file_stream()); + deallog << "Converged in " << niter << " iterations." << std::endl; + } + + catch (const std::exception &exc) + { + deallog << "Caught an irrecoverable exception in a callback:" << std::endl + << exc.what() << std::endl; + } +} diff --git a/tests/sundials/kinsol_07.output b/tests/sundials/kinsol_07.output new file mode 100644 index 0000000000..ba030d4cdf --- /dev/null +++ b/tests/sundials/kinsol_07.output @@ -0,0 +1,17 @@ + +DEAL::Computing residual at 10.0000 +DEAL::Setting up Jacobian system at u=10.0000 +DEAL::Computing residual at 10.0000 +DEAL::Computing residual at -88.0839 +DEAL::Reporting irrecoverable failure. +DEAL::Caught an irrecoverable exception in a callback: +DEAL:: +-------------------------------------------------------- +An error occurred in line <0> of file <> in function + +The violated condition was: + +Additional information: + Irrecoverable failure. +-------------------------------------------------------- + -- 2.39.5