BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Rasmus Tamstorf

Rasmus Tamstorf contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysis

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
2topics
3close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

3 published item(s)

preprint2022arXiv

Constructing Nitsche's method for variational problems

Nitsche's method is a well-established approach for weak enforcement of boundary conditions for partial differential equations (PDEs). It has many desirable properties, including the preservation of variational consistency and the fact that it yields symmetric, positive-definite discrete linear systems that are not overly ill-conditioned. In recent years, the method has gained in popularity in a number of areas, including isogeometric analysis, immersed methods, and contact mechanics. However, arriving at a formulation based on Nitsche's method can be a mathematically arduous process, especially for high-order PDEs. Fortunately, the derivation is conceptually straightforward in the context of variational problems. To facilitate the process, we devised an abstract framework for constructing Nitsche's method for these types of problems in [J. Benzaken, J. A. Evans, S. McCormick, and R. Tamstorf, Nitsche's method for linear Kirchhoff-Love shells: Formulation, error analysis, and verification, Comput. Methods Appl. Mech. Eng., 374 (2021), p. 113544]. The goal of this paper is to elucidate the process through a sequence of didactic examples. First, we show the derivation of Nitsche's method for Poisson's equation to gain an intuition for the various steps. Next, we present the abstract framework and then revisit the derivation for Poisson's equation to use the framework and add mathematical rigor. In the process, we extend our derivation to cover the vector-valued setting. Armed with a basic recipe, we then show how to handle a higher-order problem by considering the vector-valued biharmonic equation and the linearized Kirchhoff-Love plate. In the end, the hope is that the reader will be able to apply Nitsche's method to any problem that arises from variational principles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Algebraic error analysis for mixed-precision multigrid solvers

This paper establishes the first theoretical framework for analyzing the rounding-error effects on multigrid methods using mixed-precision iterative-refinement solvers. While motivated by the sparse symmetric positive definite (SPD) matrix equations that arise from discretizing linear elliptic PDEs, the framework is purely algebraic such that it applies to matrices that do not necessarily come from the continuum. Based on the so-called energy or $A$ norm, which is the natural norm for many problems involving SPD matrices, we provide a normwise forward error analysis, and introduce the notion of progressive precision for multigrid solvers. Each level of the multigrid hierarchy uses three different precisions that each increase with the fineness of the level, but at different rates, thereby ensuring that the bulk of the computation uses the lowest possible precision. The theoretical results developed here in the energy norm differ notably from previous theory based on the Euclidean norm in important ways. In particular, we show that simply rounding an exact result to finite precision causes an error in the energy norm that is proportional to the square root of $κ$, the associated matrix condition number. (By contrast, this error is of order $1$ when measured in the Euclidean norm.) Given this observation, we show that the limiting accuracy for both V-cycles and full multigrid is optimal in the sense that it is also proportional to $κ^{1/2}$ in energy. Additionally, we show that the loss of convergence rate due to rounding grows in proportion to $κ^{1/2}$, but argue that this loss is insignificant in practice. The theory presented here is the first forward error analysis in the energy norm of iterative refinement and the first rounding error analysis of multigrid in general.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discretization-error-accurate mixed-precision multigrid solvers

This paper builds on the algebraic theory in the companion paper [Algebraic Error Analysis for Mixed-Precision Multigrid Solvers] to obtain discretization-error-accurate solutions for linear elliptic partial differential equations (PDEs) by mixed-precision multigrid solvers. It is often assumed that the achievable accuracy is limited by discretization or algebraic errors. On the contrary, we show that the quantization error incurred by simply storing the matrix in any fixed precision quickly begins to dominate the total error as the discretization is refined. We extend the existing theory to account for these quantization errors and use the resulting bounds to guide the choice of four different precision levels in order to balance quantization, algebraic, and discretization errors in the progressive-precision scheme proposed in the companion paper. A remarkable result is that while iterative refinement is susceptible to quantization errors during the residual and update computation, the V-cycle used to compute the correction in each iteration is much more resilient, and continues to work if the system matrices in the hierarchy become indefinite due to quantization. As a result, the V-cycle only requires relatively few bits of precision per level. Based on our findings, we outline a simple way to implement a progressive precision FMG solver with minimal overhead, and demonstrate as an example that the one dimensional biharmonic equation can be solved reliably to any desired accuracy using just a few V-cycles when the underlying smoother works well. In the process we also confirm many of the theoretical results numerically.

0 citations0 reviewsNo code linkedOpen access