Researcher profile

Stefan Takacs

Stefan Takacs contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
2topics
4close 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

10 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.

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.

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

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".

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.

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.

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.

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.

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.

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.