Related papers: An iterative method to compute the overlap Dirac o…
One of the most computationally expensive steps of the low-rank ADI method for large-scale Lyapunov equations is the solution of a shifted linear system at each iteration. We propose the use of the extended Krylov subspace method for this…
This paper presents a new algorithm KIOPS for computing linear combinations of $\varphi$-functions that appear in exponential integrators. This algorithm is suitable for large-scale problems in computational physics where little or no…
We consider the effect of discretization errors on the microscopic spectrum of the Wilson Dirac operator using both chiral Perturbation Theory and chiral Random Matrix Theory. A graded chiral Lagrangian is used to evaluate the microscopic…
This paper establishes the iteration-complexity of an inner accelerated inexact proximal augmented Lagrangian (IAIPAL) method for solving linearly-constrained smooth nonconvex composite optimization problems that is based on the classical…
We present non-perturbative renormalization constants of fermionic bilinears on the lattice in the quenched approximation at beta=6.1 using an overlap fermion action with hypercubic(HYP)-blocked links. We consider the effects of the exact…
We determine the free field hypercubic Dirac operator which is optimally close to satisfying the Ginsparg-Wilson relation. Inserting this operator into the overlap formula, we show that the analytic locality bound on the resulting overlap…
Computing quark propagators with overlap fermions requires the solution of a shifted unitary linear system. Jagels and Reichel have shown that for such systems it is possible to construct a minimal residual algorithm by short recurrences.…
In this paper, we develop algorithms for computing the recurrence coefficients corresponding to multiple orthogonal polynomials on the step-line. We reformulate the problem as an inverse eigenvalue problem, which can be solved using…
Since being analyzed by Rokhlin, Szlam, and Tygert and popularized by Halko, Martinsson, and Tropp, randomized Simultaneous Power Iteration has become the method of choice for approximate singular value decomposition. It is more accurate…
We propose a nonlocal extension of the generalized Dirac oscillator (GDO) in $(1+1)$ dimensions by replacing the multiplicative interaction $f(x)$ with an integral operator $\hat F$ with kernel $f(x,x')$. The resulting Dirac equation…
In quantum chromodynamics (QCD) at nonzero chemical potential, the eigenvalues of the Dirac operator are scattered in the complex plane. Can the fluctuation properties of the Dirac spectrum be described by universal predictions of…
Understanding dynamics in complex systems is challenging because there are many degrees of freedom, and those that are most important for describing events of interest are often not obvious. The leading eigenfunctions of the transition…
We consider perturbations of Dirac type operators on complete, connected metric spaces equipped with a doubling measure. Under a suitable set of assumptions, we prove quadratic estimates for such operators and hence deduce that these…
We present algorithmic improvements to the overlap Hybrid Monte Carlo algorithm, including preconditioning techniques and improvements to the correction step, used when one of the eigenvalues of the Kernel operator changes sign, which is…
An algorithm is proposed, analyzed, and tested experimentally for solving stochastic optimization problems in which the decision variables are constrained to satisfy equations defined by deterministic, smooth, and nonlinear functions. It is…
We summarize our recent investigations of lattice QCD with dynamical overlap fermions. We sketch algorithmic issues and our approach to solving them. We show our measurement of the topological susceptibility. We describe a computation of…
In many statistical learning problems, it is desired that the optimal solution conforms to an a priori known sparsity structure represented by a directed acyclic graph. Inducing such structures by means of convex regularizers requires…
In this paper, two new efficient algorithms for calculating the sign function of the large-scale sparse matrix are proposed by combining filtering algorithm with Newton method and Newton Schultz method respectively. Through the theoretical…
We propose an efficient numerical method for a non-selfadjoint Steklov eigenvalue problem. The Lagrange finite element is used for discretization. The convergence is proved using the spectral perturbation theory for compact operators. The…
We consider the two-dimensional Dirac operator with Lorentz-scalar $\delta$-shell interactions on each edge of a star-graph. An orthogonal decomposition is performed which shows such an operator is unitarily equivalent to an orthogonal sum…