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

Manuel Radons

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

ResearcherUnclaimed source record
math.NAmath.COmath.OCNumerical Analysis

Catalog footprint

What is connected

5works
4topics
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

5 published item(s)

preprint2021arXiv

Sign controlled solvers for the absolute value equation with an application to support vector machines

Let $A$ be a real $n\times n$ matrix and $z,b\in \mathbb R^n$. The piecewise linear equation system $z-A\vert z\vert = b$ is called an absolute value equation. It is equivalent to the general linear complementarity problem, and thus NP hard in general. Concerning the latter problem, three solvers are presented: One direct, one semi-iterative and one discrete variant of damped Newton. Their previously proved ranges of correctness and convergence, respectively, are extended. Their performance is compared on instances of the XOR separation problem for support vector machines which can be reformulated as an absolute value equation.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Integrating Lipschitzian Dynamical Systems using Piecewise Algorithmic Differentiation

In this article we analyze a generalized trapezoidal rule for initial value problems with piecewise smooth right hand side \(F:\R^n\to\R^n\). When applied to such a problem the classical trapezoidal rule suffers from a loss of accuracy if the solution trajectory intersects a nondifferentiability of \(F\). The advantage of the proposed generalized trapezoidal rule is threefold: Firstly we can achieve a higher convergence order than with the classical method. Moreover, the method is energy preserving for piecewise linear Hamiltonian systems. Finally, in analogy to the classical case we derive a third order interpolation polynomial for the numerical trajectory. In the smooth case the generalized rule reduces to the classical one. Hence, it is a proper extension of the classical theory. An error estimator is given and numerical results are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Solving piecewise linear equations in abs-normal form

With the ultimate goal of iteratively solving piecewise smooth (PS) systems, we consider the solution of piecewise linear (PL) equations. PL models can be derived in the fashion of automatic or algorithmic differentiation as local approximations of PS functions with a second order error in the distance to a given reference point. The resulting PL functions are obtained quite naturally in what we call the abs-normal form, a variant of the state representation proposed by Bokhoven in his dissertation. Apart from the tradition of PL modeling by electrical engineers, which dates back to the Master thesis of Thomas Stern in 1956, we take into account more recent results on linear complementarity problems and semi-smooth equations originating in the optimization community. We analyze simultaneously the original PL problem (OPL) in abs-normal form and a corresponding complementary system (CPL), which is closely related to the absolute value equation (AVE) studied by Mangasarian et al and a corresponding linear complementarity problem (LCP). We show that the CPL, like KKT conditions and other simply switched systems, cannot be open without being injective. Hence some of the intriguing PL structure described by Scholtes is lost in the transformation from OPL to CPL. To both problems one may apply Newton variants with appropriate generalized Jacobians directly computable from the abs-normal representation. Alternatively, the CPL can be solved by Bokhoven's modulus method and related fixed point iterations. We compile the properties of the various schemes and highlight the connection to the properties of the Schur complement matrix, in particular its signed real spectral radius as analyzed by Rump. Numerical experiments and suitable combinations of the fixed point solvers and stabilized generalized Newton variants remain to be realized.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$\mathcal O(n)$ working precision inverses for symmetric tridiagonal Toeplitz matrices with $\mathcal O(1)$ floating point calculations

A well known numerical task is the inversion of large symmetric tridiagonal Toeplitz matrices, i.e., matrices whose entries equal $a$ on the diagonal and $b$ on the extra diagonals ($a, b\in \mathbb R$). The inverses of such matrices are dense and there exist well known explicit formulas by which they can be calculated in $\mathcal O(n^2)$. In this note we present a simplification of the problem that has proven to be rather useful in everyday practice: If $\vert a\vert > 2\vert b\vert$, that is, if the matrix is strictly diagonally dominant, its inverse is a band matrix to working precision and the bandwidth is independent of $n$ for sufficiently large $n$. Employing this observation, we construct a linear time algorithm for an explicit tridiagonal inversion that only uses $\mathcal O(1)$ floating point operations. On the basis of this simplified inversion algorithm we outline the cornerstones for an efficient parallelizable approximative equation solver.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Direct solution of piecewise linear systems

Let $S$ be a real $n\times n$ matrix, $z,\hat c\in \mathbb R^n$, and $| z|$ the componentwise modulus of $z$. Then the piecewise linear equation system $$z-S| z| = \hat c$$ is called an \textit{absolute value equation} (AVE). It has been proven to be equivalent to the general \textit{linear complementarity problem}, which means that it is NP hard in general. We will show that for several system classes the AVE essentially retains the good natured solvability properties of regular linear systems. I.e., it can be solved directly by a slightly modified Gaussian elimination that we call the signed Gaussian elimination. For dense matrices $S$ this algorithm has the same operations count as the classical Gaussian elimination with symmetric pivoting. For tridiagonal systems in $n$ variables its computational cost is roughly that of sorting $n$ floating point numbers. The sharpness of the proposed restrictions on $S$ will be established.

0 citations0 reviewsNo code linkedOpen access