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

Matthias Bolten

Matthias Bolten contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical AnalysisComputational Engineering, Finance, and Sciencemath.OCMathematical SoftwareNeurons and Cognition

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
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

9 published item(s)

preprint2022arXiv

Symbol based convergence analysis in multigrid methods for saddle point problems

Saddle point problems arise in a variety of applications, e.g., when solving the Stokes equations. They can be formulated such that the system matrix is symmetric, but indefinite, so the variational convergence theory that is usually used to prove multigrid convergence cannot be applied. In a 2016 paper in Numerische Mathematik Notay has presented a different algebraic approach that analyzes properly preconditioned saddle point problems, proving convergence of the Two-Grid method. In the present paper we analyze saddle point problems where the blocks are circulant within this framework. We are able to derive sufficient conditions for convergence and provide optimal parameters for the preconditioning of the saddle point problem and for the point smoother that is used. The analysis is based on the generating symbols of the circulant blocks. Further, we show that the structure can be kept on the coarse level, allowing for a recursive application of the approach in a W- or V-cycle and proving the "level independency" property. Numerical results demonstrate the efficiency of the proposed method in the circulant and the Toeplitz case.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Toeplitz Momentary Symbols: definition, results, and limitations in the spectral analysis of Structured Matrices

A powerful tool for analyzing and approximating the singular values and eigenvalues of structured matrices is the theory of GLT sequences. By the GLT theory one can derive a function, which describes the singular value or the eigenvalue distribution of the sequence, the latter under precise assumptions. However, for small values of the matrix size of the considered sequence, the approximations may not be as good as it is desirable, since in the construction of the GLT symbol one disregards small norm and low-rank perturbations. On the other hand, LFA can be used to construct polynomial symbols in a similar manner for discretizations, where the geometric information is present, but the small norm perturbations are retained. The main focus of this paper is the introduction of the concept of sequence of "Toeplitz momentary symbols", associated with a given sequence of truncated Toeplitz-like matrices. We construct the symbol in the same way as in the GLT theory, but we keep the information of the small norm contributions. The low-rank contributions are still disregarded, and we give an idea on the reason why this is negligible in certain cases and why it is not in other cases, being aware that in presence of high nonnormality the same low-rank perturbation can produce a dramatic change in the eigenvalue distribution. Moreover, a difference with respect to the LFA symbols is that GLT symbols and Toeplitz momentary symbols are more general and are applicable to a larger class of matrices. We show the applicability of the approach which leads to higher accuracy in some cases when compared with the GLT symbol. Finally, since for many applications and their analysis it is often necessary to consider non-square Toeplitz matrices, we formalize and provide some useful definitions, applicable for non-square Toeplitz momentary symbols.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Event-based update of synapses in voltage-based learning rules

Due to the point-like nature of neuronal spiking, efficient neural network simulators often employ event-based simulation schemes for synapses. Yet many types of synaptic plasticity rely on the membrane potential of the postsynaptic cell as a third factor in addition to pre- and postsynaptic spike times. Synapses therefore require continuous information to update their strength which a priori necessitates a continuous update in a time-driven manner. The latter hinders scaling of simulations to realistic cortical network sizes and relevant time scales for learning. Here, we derive two efficient algorithms for archiving postsynaptic membrane potentials, both compatible with modern simulation engines based on event-based synapse updates. We theoretically contrast the two algorithms with a time-driven synapse update scheme to analyze advantages in terms of memory and computations. We further present a reference implementation in the spiking neural network simulator NEST for two prototypical voltage-based plasticity rules: the Clopath rule and the Urbanczik-Senn rule. For both rules, the two event-based algorithms significantly outperform the time-driven scheme. Depending on the amount of data to be stored for plasticity, which heavily differs between the rules, a strong performance increase can be achieved by compressing or sampling of information on membrane potentials. Our results on computational efficiency related to archiving of information provide guidelines for the design of learning rules in order to make them practically usable in large-scale networks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GivEn -- Shape Optimization for Gas Turbines in Volatile Energy Networks

This paper describes the project GivEn that develops a novel multicriteria optimization process for gas turbine blades and vanes using modern "adjoint" shape optimization algorithms. Given the many start and shut-down processes of gas power plants in volatile energy grids, besides optimizing gas turbine geometries for efficiency, the durability understood as minimization of the probability of failure is a design objective of increasing importance. We also describe the underlying coupling structure of the multiphysical simulations and use modern, gradient based multicriteria optimization procedures to enhance the exploration of Pareto-optimal solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-overlapping block smoothers for the Stokes equations

Overlapping block smoothers efficiently damp the error contributions from highly oscillatory components within multigrid methods for the Stokes equations but they are computationally expensive. This paper is concentrated on the development and analysis of new block smoothers for the Stokes equations that are discretized on staggered grids. These smoothers are non-overlapping and therefore desirable due to reduced computational costs. Traditional geometric multigrid methods are based on simple pointwise smoothers. However, the efficiency of multigrid methods for solving more difficult problems such as the Stokes equations lead to computationally more expensive smoothers, e.g., overlapping block smoothers. Non-overlapping smoothers are less expensive, but have been considered less efficient in the literature. In this paper, we develop new non-overlapping smoothers, the so-called triad-wise smoothers, and show their efficiency within multigrid methods to solve the Stokes equations. In addition, we compare overlapping and non-overlapping smoothers by measuring their computational costs and analyzing their behavior by the use of local Fourier analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Parallel-in-Time Simulation of an Electrical Machine using MGRIT

We apply the multigrid-reduction-in-time (MGRIT) algorithm to an eddy current simulation of a two-dimensional induction machine supplied by a pulse-width-modulation signal. To resolve the fast-switching excitations, small time steps are needed, such that parallelization in time becomes highly relevant for reducing the simulation time. The MGRIT algorithm is well suited for introducing time parallelism in the simulation of electrical machines using existing application codes, as MGRIT is a non-intrusive approach that essentially uses the same time integrator as a traditional time-stepping algorithm. We investigate effects of spatial coarsening on MGRIT convergence when applied to two numerical models of an induction machine, one with linear material laws and a full nonlinear model. Parallel results demonstrate significant speedup in the simulation time compared to sequential time stepping, even for moderate numbers of processors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PyMGRIT: A Python Package for the parallel-in-time method MGRIT

In this paper, we introduce the Python framework PyMGRIT, which implements the multigrid-reduction-in-time (MGRIT) algorithm for solving the (non-)linear systems arising from the discretization of time-dependent problems. The MGRIT algorithm is a reduction-based iterative method that allows parallel-in-time simulations, i. e., calculating multiple time steps simultaneously in a simulation, by using a time-grid hierarchy. The PyMGRIT framework features many different variants of the MGRIT algorithm, ranging from different multigrid cycle types and relaxation schemes, as well as various coarsening strategies, including time-only and space-time coarsening, to using different time integrators on different levels in the multigrid hierachy. PyMGRIT allows serial runs for prototyping and testing of new approaches, as well as parallel runs using the Message Passing Interface (MPI). Here, we describe the implementation of the MGRIT algorithm in PyMGRIT and present the usage from both user and developer point of views. Three examples illustrate different aspects of the package, including pure time parallelism as well as space-time parallelism by coupling PyMGRIT with PETSc or Firedrake, which enable spatial parallelism through MPI.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards Real-Time Magnetic Dosimetry Simulations for Inductive Charging Systems

The exposure of a human by magneto-quasistatic fields from wireless charging systems is to be determined from magnetic field measurements in near real-time. This requires a fast linear equations solver for the discrete Poisson system of the Co-Simulation Scalar Potential Finite Difference (Co-Sim. SPFD) scheme. Here, the use of the AmgX library on NVIDIA GPUs is presented for this task. It enables solving the equation system resulting from an ICNIRP recommended human voxel model resolution of 2 mm in less than 0.5 seconds on a single NVIDIA Tesla V100 GPU.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tracing locally Pareto optimal points by numerical integration

We suggest a novel approach for the efficient and reliable approximation of the Pareto front of sufficiently smooth unconstrained bi-criteria optimization problems. Optimality conditions formulated for weighted sum scalarizations of the problem yield a description of (parts of) the Pareto front as a parametric curve, parameterized by the scalarization parameter (i.e., the weight in the weighted sum scalarization). Its sensitivity w.r.t. parameter variations can be described by an ordinary differential equation (ODE). Starting from an arbitrary initial Pareto optimal solution, the Pareto front can then be traced by numerical integration. We provide an error analysis based on Lipschitz properties and suggest an explicit Runge-Kutta method for the numerical solution of the ODE. The method is validated on bi-criteria convex quadratic programming problems for which the exact solution is explicitly known, and numerically tested on complex bi-criteria shape optimization problems involving finite element discretizations of the state equation.

0 citations0 reviewsNo code linkedOpen access