Related papers: New rigorous perturbation bounds for the LU and QR…
For linear recurrence systems, the problem of finding rational solutions is reduced to the problem of computing polynomial solutions by computing a content bound or a denominator bound. There are several bounds in the literature. The…
We take advantage of the fact that in lambda phi ^4 problems a large field cutoff phi_max makes perturbative series converge toward values exponentially close to the exact values, to make optimal choices of phi_max. For perturbative series…
We propose a new risk-constrained reformulation of the standard Linear Quadratic Regulator (LQR) problem. Our framework is motivated by the fact that the classical (risk-neutral) LQR controller, although optimal in expectation, might be…
Factor-revealing linear programs (LPs) and policy-revealing LPs arise in various contexts of algorithm design and analysis. They are commonly used techniques for analyzing the performance of approximation and online algorithms, especially…
The Cholesky QR algorithm is an efficient communication-minimizing algorithm for computing the QR factorization of a tall-skinny matrix. Unfortunately it has the inherent numerical instability and breakdown when the matrix is…
We present a novel method to compute componentwise transient bounds, ultimate bounds, and invariant regions for a class of switching continuous-time linear systems with perturbation bounds that may depend nonlinearly on a delayed state. The…
This paper provides a comprehensive estimation framework via nuclear norm plus $l_1$ norm penalization for high-dimensional approximate factor models with a sparse residual covariance. The underlying assumptions allow for non-pervasive…
Floating-point round-off errors are ubiquitous in numerically intensive programs arising in fields such as scientific computing and optimization. As floating-point errors potentially lead to unexpected and catastrophic program failures, one…
Perturbative QCD results in the MSbar scheme can be dramatically improved by switching to a scheme that accounts for the dominant power law dependence on the factorization scale in the operator product expansion. We introduce the ``MSR…
We propose a new stopping criterion for Krylov subspace iterative regularization of large-scale ill-posed inverse problems. Our stopping criterion accurately filters the data using a generalization of the Picard parameter that was…
Two new sum rules were recently discovered by Le Yaouanc et al. by applying the operator product expansion to the nonforward matrix element of a time-ordered product of $b \to c$ currents in the heavy-quark limit of QCD. They lead to the…
Quantum Annealing (QA) is a computational framework where a quantum system's continuous evolution is used to find the global minimum of an objective function over an unstructured search space. It can be seen as a general metaheuristic for…
We present the analytic expressions for the color-planar contributions to the heavy-light form factors at three loops in perturbative QCD. These form factors play an important role in the precision predictions of various observables in top…
A finite element methodology for large classes of variational boundary value problems is defined which involves discretizing two linear operators: (1) the differential operator defining the spatial boundary value problem; and (2) a Riesz…
We revisit the problem of computing submatrices of the Cram\'er-Rao bound (CRB), which lower bounds the variance of any unbiased estimator of a vector parameter $\vth$. We explore iterative methods that avoid direct inversion of the Fisher…
We present a high order perturbative computation of the renormalization constants Z_V, Z_A and of the ratio Z_P/Z_S for Wilson fermions. The computational setup is the one provided by the RI'-MOM scheme. Three- and four-loop expansions are…
We calculate the branching ratio and $CP$ violations for $B\to K\pi$ decays in a modified perturbative QCD approach based on $k_{T}$ factorization. The resummation effect of the transverse momentum regulates the endpoint singularity. Using…
We consider the solution of large linear systems of equations that arise when two-dimensional singularly perturbed reaction-diffusion equations are discretized. Standard methods for these problems, such as central finite differences, lead…
We establish new results of first-order necessary conditions of optimality for finite-dimensional problems with inequality constraints and for problems with equality and inequality constraints, in the form of John's theorem and in the form…
Quadratic programmingis a class of constrained optimization problem with quadratic objective functions and linear constraints. It has applications in many areas and is also used to solve nonlinear optimization problems. This article focuses…