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

Stefan Takacs

Stefan Takacs 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

14works
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

14 published item(s)

preprint2022arXiv

A IETI-DP method for discontinuous Galerkin discretizations in Isogeometric Analysis with inexact local solvers

We construct solvers for an isogeometric multi-patch discretization, where the patches are coupled via a discontinuous Galerkin approach, which allows the consideration of discretizations that do not match on the interfaces. We solve the resulting linear system using a Dual-Primal IsogEometric Tearing and Interconnecting (IETI-DP) method. We are interested in solving the arising patch-local problems using iterative solvers since this allows the reduction of the memory footprint. We solve the patch-local problems approximately using the Fast Diagonalization method, which is known to be robust in the grid size and the spline degree. To obtain the tensor structure needed for the application of the Fast Diagonalization method, we introduce an orthogonal splitting of the local function spaces. We present a convergence theory that confirms that the condition number of the preconditioned system only grows poly-logarithmically with the grid size. The numerical experiments confirm this finding. Moreover, they show that the convergence of the overall solver only mildly depends on the spline degree. We observe a mild reduction of the computational times and a significant reduction of the memory requirements in comparison to standard IETI-DP solvers using sparse direct solvers for the local subproblems. Furthermore, the experiments indicate good scaling behavior on distributed memory machines.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

IETI-DP for conforming multi-patch Isogeometric Analysis in three dimensions

We consider dual-primal isogeometric tearing and interconnection (IETI-DP) solvers for multi-patch geometries in Isogeometric Analysis. Recently, the authors have published a convergence analysis for those solvers that is explicit in both the grid size and the spline degree for conforming discretizations of two dimensional computational domains. In the present paper, we shortly revisit these results and provide numerical experiments that indicate that similar results may hold for three dimensional domains.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

IETI-DP methods for discontinuous Galerkin multi-patch Isogeometric Analysis with T-junctions

We study Dual-Primal Isogeometric Tearing and Interconnecting (IETI-DP) solvers for non-conforming multi-patch discretizations of a generalized Poisson problem. We realize the coupling between the patches using a symmetric interior penalty discontinuous Galerkin (SIPG) approach. Previously, we have assumed that the interfaces between patches always consist of whole edges. In this paper, we drop this requirement and allow T-junctions. This extension is vital for the consideration of sliding interfaces, for example between the rotor and the stator of an electrical motor. One critical part for the handling of T-junctions in IETI-DP solvers is the choice of the primal degrees of freedom. We propose to add all basis functions that are non-zero at any of the vertices to the primal space. Since there are several such basis functions at any T-junction, we call this concept ''fat vertices''. For this choice, we show a condition number bound that coincides with the bound for the conforming cas

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Towards a IETI-DP solver on non-matching multi-patch domains

Recently, the authors have proposed and analyzed isogeometric tearing and interconnecting (IETI-DP) solvers for multi-patch discretizations in Isogeometric Analysis. Conforming and discontinuous Galerkin settings have been considered. In both cases, we have assumed that the interfaces between the patches consist of whole edges. In this paper, we present a generalization that allows us to drop this requirement. This means that the patches can meet in T-junctions, which increases the flexibility of the geometric model significantly. We use vertex-based primal degrees of freedom. For the T-junctions, we propose to follow the idea of "fat vertices".

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence theory for IETI-DP solvers for discontinuous Galerkin Isogeometric Analysis that is explicit in h and p

In this paper, we develop a convergence theory for Dual-Primal Isogeometric Tearing and Interconnecting (IETI-DP) solvers for isogeometric multi-patch discretizations of the Poisson problem, where the patches are coupled using discontinuous Galerkin. The presented theory provides condition number bounds that are explicit in the grid sizes h and in the spline degrees p. We give an analysis that holds for various choices for the primal degrees of freedom: vertex values, edge averages, and a combination of both. If only the vertex values or both vertex values and edge averages are taken as primal degrees of freedom, the condition number bound is the same as for the conforming case. If only the edge averages are taken, both the convergence theory and the experiments show that the condition number of the preconditioned system grows with the ratio of the grid sizes on neighboring patches.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Condition number bounds for IETI-DP methods that are explicit in h and p

We study the convergence behavior of Dual-Primal Isogeometric Tearing and Interconnecting (IETI-DP) methods for solving large-scale algebraic systems arising from multi-patch Isogeometric Analysis. We focus on the Poisson problem on two dimensional computational domains. We provide a convergence analysis that covers several choices of the primal degrees of freedom: the vertex values, the edge averages, and the combination of both. We derive condition number bounds that show the expected behavior in the grid size h and that are quasi-linear in the spline degree p.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fast multigrid solvers for conforming and non-conforming multi-patch Isogeometric Analysis

Isogeometric Analysis is a high-order discretization method for boundary value problems that uses a number of degrees of freedom which is as small as for a low-order method. Standard isogeometric discretizations require a global parameterization of the computational domain. In non-trivial cases, the domain is decomposed into patches having separate parameterizations and separate discretization spaces. If the discretization spaces agree on the interfaces between the patches, the coupling can be done in a conforming way. Otherwise, non-conforming discretizations (utilizing discontinuous Galerkin approaches) are required. The author and his coworkers have previously introduced multigrid solvers for Isogeometric Analysis for the conforming case. In the present paper, these results are extended to the non-conforming case. Moreover, it is shown that the multigrid solves get even more powerful if the proposed smoother is combined with a (standard) Gauss-Seidel smoother.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A parallel multigrid solver for multi-patch Isogeometric Analysis

Isogeometric Analysis (IgA) is a framework for setting up spline-based discretizations of partial differential equations, which has been introduced around a decade ago and has gained much attention since then. If large spline degrees are considered, one obtains the approximation power of a high-order method, but the number of degrees of freedom behaves like for a low-order method. One important ingredient to use a discretization with large spline degree, is a robust and preferably parallelizable solver. While numerical evidence shows that multigrid solvers with standard smoothers (like Gauss Seidel) does not perform well if the spline degree is increased, the multigrid solvers proposed by the authors and their co-workers proved to behave optimal both in the grid size and the spline degree. In the present paper, the authors want to show that those solvers are parallelizable and that they scale well in a parallel environment.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Robust approximation error estimates and multigrid solvers for isogeometric multi-patch discretizations

In recent publications, the author and his coworkers have shown robust approximation error estimates for B-splines of maximum smoothness and have proposed multigrid methods based on them. These methods allow to solve the linear system arizing from the discretization of a partial differential equation in Isogeometric Analysis in a single-patch setting with convergence rates that are provably robust both in the grid size and the spline degree. In real-world problems, the computational domain cannot be nicely represented by just one patch. In computer aided design, such domains are typically represented as a union of multiple patches. In the present paper, we extend the approximation error estimates and the multigrid solver to this multi-patch case.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Robust multigrid methods for isogeometric discretizations of the Stokes equations

In recent publications, the author and his coworkers have proposed a multigrid method for solving linear systems arizing from the discretization of partial differential equations in isogeometric analysis and have proven that the convergence rates are robust in both the grid size and the polynomial degree. So, far the method has only been discussed for the Poisson problem. In the present paper, we want to face the question if it is possible to extend the method to the Stokes equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A robust multigrid method for the time-dependent Stokes problem

In the present paper we propose a coupled multigrid method for generalized Stokes flow problems. Such problems occur as subproblems in implicit time-stepping approaches for time-dependent Stokes problems. The discretized Stokes system is a large-scale linear system whose condition number depends on the grid size of the spatial discretization and of the length of the time step. Recently, for this problem a coupled multigrid method has been proposed, where in each smoothing step a Poisson problem has to be solved (approximately) for the pressure field. In the present paper, we propose a coupled multigrid method where the solution of such sub-problems is not needed. We prove that the proposed method shows robust convergence behavior in the grid size of the spatial discretization and of the length of the time-step.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems

In the present paper we concentrate on an important issue in constructing a good multigrid solver: the choice of an efficient smoother. We will introduce all-at-once multigrid solvers for optimal control problems which show robust convergence in the grid size and in the regularization parameter. We will refer to recent publications that guarantee such a convergence behavior. These publications do not pay much attention to the construction of the smoother and suggest to use a normal equation smoother. We will see that using a Gauss Seidel like variant of this smoother, the overall multigrid solver is speeded up by a factor of about two with no additional work. The author will give a proof which indicates that also the Gauss Seidel like variant of the smoother is covered by the convergence theory. Numerical experiments suggest that the proposed method are competitive with Vanka type methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A robust all-at-once multigrid method for the Stokes control problem

In this paper we present an all-at-once multigrid method for a distributed Stokes control problem (velocity tracking problem). For solving such a problem, we use the fact that the solution is characterized by the optimality system (Karush-Kuhn-Tucker-system). The discretized optimality system is a large-scale linear system whose condition number depends on the grid size and on the choice of the regularization parameter forming a part of the problem. Recently, block-diagonal preconditioners have been proposed, which allow to solve the problem using a Krylov space method with convergence rates that are robust in both, the grid size and the regularization parameter or cost parameter. In the present paper, we develop an all-at-once multigrid method for a Stokes control problem and show robust convergence, more precisely, we show that the method converges with rates which are bounded away from one by a constant which is independent of the grid size and the choice of the regularization or cost parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Using cylindrical algebraic decomposition and local Fourier analysis to study numerical methods: two examples

Local Fourier analysis is a strong and well-established tool for analyzing the convergence of numerical methods for partial differential equations. The key idea of local Fourier analysis is to represent the occurring functions in terms of a Fourier series and to use this representation to study certain properties of the particular numerical method, like the convergence rate or an error estimate. In the process of applying a local Fourier analysis, it is typically necessary to determine the supremum of a more or less complicated term with respect to all frequencies and, potentially, other variables. The problem of computing such a supremum can be rewritten as a quantifier elimination problem, which can be solved with cylindrical algebraic decomposition, a well-known tool from symbolic computation. The combination of local Fourier analysis and cylindrical algebraic decomposition is a machinery that can be applied to a wide class of problems. In the present paper, we will discuss two examples. The first example is to compute the convergence rate of a multigrid method. As second example we will see that the machinery can also be used to do something rather different: We will compare approximation error estimates for different kinds of discretizations.

0 citations0 reviewsNo code linkedOpen access