BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thomas P. Wihler

Thomas P. Wihler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysis

Catalog footprint

What is connected

16works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2023arXiv

A note on the Sassenfeld criterion and its relation to H-matrices

The starting point of this note is a decades-old yet little-noticed sufficient condition, presented by Sassenfeld in 1951, for the convergence of the classical Gauss-Seidel method. The purpose of the present paper is to shed new light on Sassenfeld's criterion and to demonstrate that it is directly related to H-matrices. In particular, our results yield a new characterization of H-matrices. In addition, the convergence of iterative linear solvers that involve H-matrix preconditioners is briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A numerical energy minimisation approach for semilinear diffusion-reaction boundary value problems based on steady state iterations

We present a novel energy-based numerical analysis of semilinear diffusion-reaction boundary value problems. Based on a suitable variational setting, the proposed computational scheme can be seen as an energy minimisation approach. More specifically, this procedure aims to generate a sequence of numerical approximations, which results from the iterative solution of related (stabilised) linearised discrete problems, and tends to a local minimum of the underlying energy functional. Simultaneously, the finite-dimensional approximation spaces are adaptively refined; this is implemented in terms of a new mesh refinement strategy in the context of finite element discretisations, which again relies on the energy structure of the problem under consideration, and does not involve any a posteriori error indicators. In combination, the resulting adaptive algorithm consists of an iterative linearisation procedure on a sequence of hierarchically refined discrete spaces, which we prove to converge towards a solution of the continuous problem in an appropriate sense. Numerical experiments demonstrate the robustness and reliability of our approach for a series of examples.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Energy contraction and optimal convergence of adaptive iterative linearized finite element methods

We revisit a unified methodology for the iterative solution of nonlinear equations in Hilbert spaces. Our key observation is that the general approach from [Heid & Wihler, Math. Comp. 89 (2020), Calcolo 57 (2020)] satisfies an energy contraction property in the context of (abstract) strongly monotone problems. This property, in turn, is the crucial ingredient in the recent convergence analysis in [Gantner et al., arXiv:2003.10785]. In particular, we deduce that adaptive iterative linearized finite element methods (AILFEMs) lead to full linear convergence with optimal algebraic rates with respect to the degrees of freedom as well as the total computational time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive local minimax Galerkin methods for variational problems

In many applications of practical interest, solutions of partial differential equation models arise as critical points of an underlying (energy) functional. If such solutions are saddle points, rather than being maxima or minima, then the theoretical framework is non-standard, and the development of suitable numerical approximation procedures turns out to be highly challenging. In this paper, our aim is to present an iterative discretization methodology for the numerical solution of nonlinear variational problems with multiple (saddle point) solutions. In contrast to traditional numerical approximation schemes, which typically fail in such situations, the key idea of the current work is to employ a simultaneous interplay of a previously developed local minimax approach and adaptive Galerkin discretizations. We thereby derive an adaptive local minimax Galerkin (LMMG) method, which combines the search for saddle point solutions and their approximation in finite-dimensional spaces in a highly effective way. Under certain assumptions, we will prove that the generated sequence of approximate solutions converges to the solution set of the variational problem. This general framework will be applied to the specific context of finite element discretizations of (singularly perturbed) semilinear elliptic boundary value problems, and a series of numerical experiments will be presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gradient Flow Finite Element Discretizations with Energy-Based Adaptivity for the Gross-Pitaevskii Equation

We present an effective adaptive procedure for the numerical approximation of the steady-state Gross-Pitaevskii equation. Our approach is solely based on energy minimization, and consists of a combination of gradient flow iterations and adaptive finite element mesh refinements. Numerical tests show that this strategy is able to provide highly accurate results, with optimal convergence rates with respect to the number of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive Pseudo-Transient-Continuation-Galerkin Methods for Semilinear Elliptic Partial Differential Equations

In this paper we investigate the application of pseudo-transient-continuation (PTC) schemes for the numerical solution of semilinear elliptic partial differential equations, with possible singular perturbations. We will outline a residual reduction analysis within the framework of general Hilbert spaces, and, subsequently, employ the PTC-methodology in the context of finite element discretizations of semilinear boundary value problems. Our approach combines both a prediction-type PTC-method (for infinite dimensional problems) and an adaptive finite element discretization (based on a robust a posteriori residual analysis), thereby leading to a fully adaptive PTC-Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for different examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An $hp$-Adaptive Newton-Discontinuous-Galerkin Finite Element Approach for Semilinear Elliptic Boundary Value Problems

In this paper we develop an $hp$-adaptive procedure for the numerical solution of general second-order semilinear elliptic boundary value problems, with possible singular perturbation. Our approach combines both adaptive Newton schemes and an $hp$-version adaptive discontinuous Galerkin finite element discretisation, which, in turn, is based on a robust $hp$-version a posteriori residual analysis. Numerical experiments underline the robustness and reliability of the proposed approach for various examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Adaptive Newton-Method Based on a Dynamical Systems Approach

The traditional Newton method for solving nonlinear operator equations in Banach spaces is discussed within the context of the continuous Newton method. This setting makes it possible to interpret the Newton method as a discrete dynamical system and thereby to cast it in the framework of an adaptive step size control procedure. In so doing, our goal is to reduce the chaotic behavior of the original method without losing its quadratic convergence property close to the roots. The performance of the modified scheme is illustrated with various examples from algebraic and differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Continuous and Discontinuous Galerkin Time Stepping Methods for Nonlinear Initial Value Problems with Application to Finite Time Blow-Up

We consider continuous and discontinuous Galerkin time stepping methods of arbitrary order as applied to nonlinear initial value problems in real Hilbert spaces. Our only assumption is that the nonlinearities are continuous; in particular, we include the case of unbounded nonlinear operators. Specifically, we develop new techniques to prove general Peano-type existence results for discrete solutions. In particular, our results show that the existence of solutions is independent of the local approximation order, and only requires the local time steps to be sufficiently small (independent of the polynomial degree). The uniqueness of (local) solutions is addressed as well. In addition, our theory is applied to finite time blow-up problems with nonlinearities of algebraic growth. For such problems we develop a time step selection algorithm for the purpose of numerically computing the blow-up time, and provide a convergence result.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Posteriori Error Analysis of $hp$-FEM for singularly perturbed problems

We consider the approximation of singularly perturbed linear second-order boundary value problems by $hp$-finite element methods. In particular, we include the case where the associated differential operator may not be coercive. Within this setting we derive an a posteriori error estimate for a natural residual norm. The error bound is robust with respect to the perturbation parameter and fully explicit with respect to both the local mesh size $h$ and the polynomial degree $p$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Adaptive Energy Minimisation for $hp$-Finite Element Methods

This article is concerned with the numerical solution of convex variational problems. More precisely, we develop an iterative minimisation technique which allows for the successive enrichment of an underlying discrete approximation space in an adaptive manner. Specifically, we outline a new approach in the context of $hp$-adaptive finite element methods employed for the efficient numerical solution of linear and nonlinear second-order boundary value problems. Numerical experiments are presented which highlight the practical performance of this new $hp$-refinement technique for both one- and two-dimensional problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Adaptive Variable Order Quadrature Strategy

In this article we propose a new adaptive numerical quadrature procedure which includes both local subdivision of the integration domain, as well as local variation of the number of quadrature points employed on each subinterval. In this way we aim to account for local smoothness properties of the function to be integrated as effectively as possible, and thereby achieve highly accurate results in a very efficient manner. Indeed, this idea originates from so-called hp-version finite element methods which are known to deliver high-order convergence rates, even for nonsmooth functions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fully Adaptive Newton-Galerkin Time Stepping Methods for Singularly Perturbed Parabolic Evolution Equations

In this paper we develop an adaptive procedure for the numerical solution of semilinear parabolic problems, with possible singular perturbations. Our approach combines a linearization technique using Newton's method with an adaptive discretization-which is based on a spatial finite element method and the backward Euler time stepping scheme-of the resulting sequence of linear problems. Upon deriving a robust a posteriori error analysis, we design a fully adaptive Newton-Galerkin time stepping algorithm. Numerical experiments underline the robustness and reliability of the proposed approach for various examples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Iterative Galerkin Discretizations for Strongly Monotone Problems

In this article we investigate a finite element formulation of strongly monotone quasi-linear elliptic PDEs in the context of fixed-point iterations. As opposed to Newton's method, which requires information from the previous iteration in order to linearise the iteration matrix (and thereby to recompute it) in each step, the alternative method used in this article exploits the monotonicity properties of the problem, and only needs the iteration matrix calculated once for all iterations of the fixed-point method. We outline the a priori and a posteriori error estimates for iteratively obtained solutions, and show both theoretically as well as numerically how the number of iterations of the fixed-point method can be restricted in dependence of the mesh size, or of the polynomial degree, to obtain optimal convergence.

0 citations0 reviewsNo code linkedOpen access