Researcher profile

Ruey-Lin Sheu

Ruey-Lin Sheu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

8 published item(s)

preprint2023arXiv

Simultaneous diagonalization via congruence of $m$ real symmetric matrices and its implications in optimization

Let $\{C_1, C_2, \ldots, C_m\},~m\ge2$ be a collection of $n\times n$ real symmetric matrices. The objective of the paper is to offer an algorithm that finds a common congruence matrix $R$ such that $R^TC_iR$ is real diagonal for every $C_i;$ or reports none of such kind. The problem, referred to as the simultaneously diagonalization via congruence (SDC in short), seems to be of pure linear algebra at first glance. However, for quadratically constrained quadratic programming (QCQP), if the quadratic forms are SDC, their joint range set is a closed convex polyhedral cone, which opens the possibility to extend the classical $\mathcal{S}$-lemma for more than two symmetric matrices. In addition, under the SDC assumption of quadratic forms, QCQP can be recast in separable forms which is usually easier to tackle. It is thus important to have a standard procedure for determining whether or not the SDC property holds for the underlined quadratic optimization problem. Our result solves a long standing problem posed by Hiriart-Urruty in 2007.

preprint2021arXiv

On Separation of level sets for a pair of quadratic functions

Given a quadratic function $f(x)=x^TAx+2a^Tx+a_0,$ it is possible that its level set $\{x\in\mathbb{R}^n: f(x)=0\}$ has two connected components and thus can be separated by the level set $\{x\in\mathbb{R}^n: g(x)=0\}$ of another quadratic function $g(x)=x^TBx+2b^Tx+b_0.$ It turns out that the separation property of such kind has great implication in quadratic optimization problems and thus deserves careful studies. In this paper, we characterize the separation property analytically by necessary and sufficient conditions as a new tool to solving optimization problems.

preprint2015arXiv

S-Lemma with Equality and Its Applications

Let $f(x)=x^TAx+2a^Tx+c$ and $h(x)=x^TBx+2b^Tx+d$ be two quadratic functions having symmetric matrices $A$ and $B$. The S-lemma with equality asks when the unsolvability of the system $f(x)<0, h(x)=0$ implies the existence of a real number $μ$ such that $f(x) + μh(x)\ge0, ~\forall x\in \mathbb{R}^n$. The problem is much harder than the inequality version which asserts that, under Slater condition, $f(x)<0, h(x)\le0$ is unsolvable if and only if $f(x) + μh(x)\ge0, ~\forall x\in \mathbb{R}^n$ for some $μ\ge0$. In this paper, we show that the S-lemma with equality does not hold only when the matrix $A$ has exactly one negative eigenvalue and $h(x)$ is a non-constant linear function ($B=0, b\not=0$). As an application, we can globally solve $\inf\{f(x)\vert h(x)=0\}$ as well as the two-sided generalized trust region subproblem $\inf\{f(x)\vert l\le h(x)\le u\}$ without any condition. Moreover, the convexity of the joint numerical range $\{(f(x), h_1(x),\ldots, h_p(x)):~x\in\Bbb R^n\}$ where $f$ is a (possibly non-convex) quadratic function and $h_1(x),\ldots,h_p(x)$ are affine functions can be characterized using the newly developed S-lemma with equality.

preprint2014arXiv

An SDP Approach For Solving Quadratic Fractional Programming Problems

This paper considers a fractional programming problem (P) which minimizes a ratio of quadratic functions subject to a two-sided quadratic constraint. As is well-known, the fractional objective function can be replaced by a parametric family of quadratic functions, which makes (P) highly related to, but more difficult than a single quadratic programming problem subject to a similar constraint set. The task is to find the optimal parameter $λ^*$ and then look for the optimal solution if $λ^*$ is attained. Contrasted with the classical Dinkelbach method that iterates over the parameter, we propose a suitable constraint qualification under which a new version of the S-lemma with an equality can be proved so as to compute $λ^*$ directly via an exact SDP relaxation. When the constraint set of (P) is degenerated to become an one-sided inequality, the same SDP approach can be applied to solve (P) {\it without any condition}. We observe that the difference between a two-sided problem and an one-sided problem lies in the fact that the S-lemma with an equality does not have a natural Slater point to hold, which makes the former essentially more difficult than the latter. This work does not, either, assume the existence of a positive-definite linear combination of the quadratic terms (also known as the dual Slater condition, or a positive-definite matrix pencil), our result thus provides a novel extension to the so-called &#34;hard case&#34; of the generalized trust region subproblem subject to the upper and the lower level set of a quadratic function.

preprint2014arXiv

Double Well Potential Function and Its Optimization in the n-dimensional Real Space - Part II

In contrast to taking the dual approach for finding a global minimum solution of a double well potential function, in Part II of the paper, we characterize a local minimizer, local maximizer, and global minimizer directly from the primal side. It is proven that, for a ``nonsingular&#34; double well function, there exists at most one local, but non-global, minimizer and at most one local maximizer. Moreover, when it exists, the local maximizer is ``surrounded&#34; by local minimizers in the sense that the norm of the local maximizer is strictly less than that of any local minimizer. We also establish some necessary and sufficient optimality conditions for the global minimizer, local non-global minimizer and local maximizer by studying a convex secular function over specific intervals. These conditions lead to three algorithms for identifying different types of critical points of a given double well function.

preprint2014arXiv

Double Well Potential Function and Its Optimization in The n-dimensional Real Space -- Part I

A special type of multi-variate polynomial of degree 4, called the double well potential function, is studied. When the function is bounded from below, it has a very unique property that two or more local minimum solutions are separated by one local maximum solution, or one saddle point. Our intension in this paper is to categorize all possible configurations of the double well potential functions mathematically. In part I, we begin the study with deriving the double well potential function from a numerical estimation of the generalized Ginzburg-Landau functional. Then, we solve the global minimum solution from the dual side by introducing a geometrically nonlinear measure which is a type of Cauchy-Green strain. We show that the dual of the dual problem is a linearly constrained convex minimization problem, which is mapped equivalently to a portion of the original double well problem subject to additional linear constraints. Numerical examples are provided to illustrate the important features of the problem and the mapping in between.

preprint2013arXiv

On RIC bounds of Compressed Sensing Matrices for Approximating Sparse Solutions Using $\ell_q$ Quasi Norms

This paper follows the recent discussion on the sparse solution recovery with quasi-norms $\ell_q,~q\in(0,1)$ when the sensing matrix possesses a Restricted Isometry Constant $δ_{2k}$ (RIC). Our key tool is an improvement on a version of &#34;the converse of a generalized Cauchy-Schwarz inequality&#34; extended to the setting of quasi-norm. We show that, if $δ_{2k}\le 1/2$, any minimizer of the $l_q$ minimization, at least for those $q\in(0,0.9181]$, is the sparse solution of the corresponding underdetermined linear system. Moreover, if $δ_{2k}\le0.4931$, the sparse solution can be recovered by any $l_q, q\in(0,1)$ minimization. The values $0.9181$ and $0.4931$ improves those reported previously in the literature.

preprint2013arXiv

Trust Region Subproblem with a Fixed Number of Additional Linear Inequality Constraints has Polynomial Complexity

The trust region subproblem with a fixed number m additional linear inequality constraints, denoted by (Tm), have drawn much attention recently. The question as to whether Problem (Tm) is in Class P or Class NP remains open. So far, the only affirmative general result is that (T1) has an exact SOCP/SDP reformulation and thus is polynomially solvable. By adopting an early result of Martinez on local non-global minimum of the trust region subproblem, we can inductively reduce any instance in (Tm) to a sequence of trust region subproblems (T0). Although the total number of (T0) to be solved takes an exponential order of m, the reduction scheme still provides an argument that the class (Tm) has polynomial complexity for each fixed m. In contrast, we show by a simple example that, solving the class of extended trust region subproblems which contains more linear inequality constraints than the problem dimension; or the class of instances consisting of an arbitrarily number of linear constraints is NP-hard. When m is small such as m = 1,2, our inductive algorithm should be more efficient than the SOCP/SDP reformulation since at most 2 or 5 subproblems of (T0), respectively, are to be handled. In the end of the paper, we improve a very recent dimension condition by Jeyakumar and Li under which (Tm) admits an exact SDP relaxation. Examples show that such an improvement can be strict indeed.