-
The Beijing-Arizona Sky Survey (BASS) is a new wide-field legacy imaging
survey in the northern Galactic cap using the 2.3m Bok telescope. The survey
will cover about 5400 deg$^2$ in the $g$ and $r$ bands, and the expected
5$\sigma$ depths (corrected for the Galactic extinction) in the two bands are
24.0 and 23.4 mag, respectively. BASS started observations in January 2015, and
has completed about 41% of the whole area as of July 2016. The first data
release contains both calibrated images and photometric catalogs obtained in
2015 and 2016. The depths of single-epoch images in the two bands are 23.4 and
22.9 mag, and the full depths of three epochs are about 24.1 and 23.5 mag,
respectively.
-
The characteristics of IR properties and MIR variability of RSGs in the LMC
are analyzed based on 12 bands of NIR to MIR co-added data from 2MASS, Spitzer
and WISE, and $\sim$6.6 years of MIR time-series data collected by the ALLWISE
and NEOWISE-R projects. 773 RSGs candidates were compiled from the literature
and verified by using the CMD, SED and MIR variability. About 15\% of valid
targets in the $IRAC1-IRAC2$/$IRAC2-IRAC3$ diagram may show PAH emission. We
show that arbitrary dereddening Q parameters related to the IRAC4, S9W, WISE3,
WISE4, and MIPS24 bands could be constructed based on a precise measurement of
MIR interstellar extinction law. Several peculiar outliers in our sample are
discussed, in which one outlier might be a RSG right before the explosion or an
x-AGB star in the very late evolutionary stage based on the MIR spectrum and
photometry. There are 744 identified RSGs in the final sample having both the
WISE1- and WISE2-band time-series data. The results show that the MIR
variability is increasing along with the increasing of brightness. There is a
relatively tight correlation between the MIR variability, MLR, and the warm
dust or continuum, where the MIR variability is evident for the targets with
$K_S-WISE3>1.0~mag$ and $WISE4<6.5~mag$, while the rest of the targets show
much smaller MIR variability. The MIR variability is also correlated with the
MLR for which targets with larger variability also show larger MLR with an
approximate upper limit of $-6.1~M_\odot/yr^{-1}$. Both the variability and the
luminosity may be important for the MLR since the WISE4-band flux is increasing
exponentially along with the degeneracy of luminosity and variability. The
identified RSG sample has been compared with the theoretical evolutionary
models and shown that the discrepancy between observation and evolutionary
models can be mitigated by considering both variability and extinction.
-
Generative Adversarial Network (GAN) and its variants have recently attracted
intensive research interests due to their elegant theoretical foundation and
excellent empirical performance as generative models. These tools provide a
promising direction in the studies where data availability is limited. One
common issue in GANs is that the density of the learned generative distribution
could concentrate on the training data points, meaning that they can easily
remember training samples due to the high model complexity of deep networks.
This becomes a major concern when GANs are applied to private or sensitive data
such as patient medical records, and the concentration of distribution may
divulge critical patient information. To address this issue, in this paper we
propose a differentially private GAN (DPGAN) model, in which we achieve
differential privacy in GANs by adding carefully designed noise to gradients
during the learning procedure. We provide rigorous proof for the privacy
guarantee, as well as comprehensive empirical evidence to support our analysis,
where we demonstrate that our method can generate high quality data points at a
reasonable privacy level.
-
Classical Cepheids are well-known and widely used distance indicators. As
distance and extinction are usually degenerate, it is important to develop
suitable methods to robustly anchor the distance scale. Here, we introduce a
near-infrared (near-IR) optimal distance method to determine both the
extinction values of and distances to a large sample of 288 Galactic classical
Cepheids. The overall uncertainty in the derived distances is less than 4.9%.
We compare our newly determined distances to the Cepheids in our sample with
previously published distances to the same Cepheids with Hubble Space Telescope
parallax measurements and distances based on the IR surface brightness method,
Wesenheit functions, and the main-sequence fitting method. The systematic
deviations in the distances determined here with respect to those of previous
publications is less than 1-2%. Hence, we constructed Galactic mid-IR
period-luminosity (PL) relations for classical Cepheids in the four Wide-Field
Infrared Survey Explorer (WISE) bands (W1, W2, W3, and W4) and the four Spitzer
Space Telescope bands ([3.6], [4.5], [5.8] and [8.0]). Based on our sample of
hundreds of Cepheids, the WISE PL relations have been determined for the first
time; their dispersion is approximately 0.10 mag. Using the currently most
complete sample, our Spitzer PL relations represent a significant improvement
in accuracy, especially in the [3.6] band which has the smallest dispersion
(0.066 mag). In addition, the average mid-IR extinction curve for Cepheids has
been obtained: A_W1/A_Ks=0.560, A_W2/A_Ks=0.479, A_W3/A_Ks=0.507,
A_W4/A_Ks=0.406, A_[3.6]/A_Ks=0.481, A_[4.5]/A_Ks=0.469, A_[5.8]/A_Ks=0.427,
and A_[8.0]/A_Ks=0.427 mag.
-
Modern software systems are often equipped with hundreds to thousands of
configuration options, many of which greatly affect performance. Unfortunately,
properly setting these configurations is challenging for developers due to the
complex and dynamic nature of system workload and environment. In this paper,
we first conduct an empirical study to understand performance-related
configurations and the challenges of setting them in the real-world. Guided by
our study, we design a systematic and general control-theoretic framework,
SmartConf, to automatically set and dynamically adjust performance-related
configurations to meet required operating constraints while optimizing other
performance metrics. Evaluation shows that SmartConf is effective in solving
real-world configuration problems, often providing better performance than even
the best static configuration developers can choose under existing
configuration systems.
-
Understanding the effects of dust extinction is important to properly
interpret observations. The optical total-to-selective extinction ratio, Rv =
Av/E(B-V), is widely used to describe extinction variations in ultraviolet and
optical bands. Since the Rv=3.1 extinction curve adequately represents the
average extinction law of diffuse regions in the Milky Way, it is commonly used
to correct observational measurements along sightlines toward diffuse regions
in the interstellar medium. However, the Rv value may vary even along different
diffuse interstellar medium sightlines. In this paper, we investigate the
optical--mid-infrared (mid-IR) extinction law toward a very diffuse region at l
= 165 in the Galactic plane, which was selected based on a CO emission map.
Adopting red clump stars as extinction tracers, we determine the
optical-to-mid-IR extinction law for our diffuse region in the two APASS bands
(B, V), the three XSTPS-GAC bands (g, r, i), the three 2MASS bands (J, H, Ks,
and the two WISE bands (W1, W2). Specifically, 18 red clump stars were selected
from the APOGEE--RC catalog based on spectroscopic data in order to explore the
diversity of the extinction law. We find that the optical extinction curves
exhibit appreciable diversity. The corresponding Rv ranges from 1.7 to 3.8,
while the mean Rv value of 2.8 is consistent with the widely adopted average
value of 3.1 for Galactic diffuse clouds. There is no apparent correlation
between Rv value and color excess E(B-V) in the range of interest, from 0.2 to
0.6 mag, or with specific visual extinction per kiloparsec, AV/d.
-
The existing molecular relaxation models based on both parallel relaxation
theory and series relaxation theory cannot extract the contributions of gas
compositions to acoustic relaxation absorption in mixtures. In this paper, we
propose an analytical model to predict acoustic relaxation absorption and
clarify composition relaxation contributions based on the rate-determining
energy transfer processes in molecular relaxation in excitable gases. By
combining parallel and series relaxation theory, the proposed model suggests
that the vibration-translation process of the lowest vibrational mode in each
composition provides the primary deexcitation path to the relaxation energy,
and the rate-determining vibration-vibration processes between the lowest modes
and others dominate the coupling energy transfer between different modes. Thus,
each gas composition contributes directly one single relaxation process to the
molecular relaxation in mixture, which can be illustrated by the decomposed
acoustic relaxation absorption spectrum of the single relaxation process. The
proposed model is validated by simulation results in good agreement with
experimental data such as $\mathrm{N_2}$, $\mathrm{O_2}$, $\mathrm{CO_2}$,
$\mathrm{CH_4}$ and their mixtures.
-
We consider the periodic problem for two-fluid non-isentropic Euler-Maxwell
systems in plasmas. By means of suitable choices of symmetrizers and an
induction argument on the order of the time-space derivatives of solutions in
energy estimates, the global smooth solution with small amplitude is
established near a non-constant equilibrium solution with asymptotic stability
properties. This improves the results obtained in \cite{LWF16a} for models with
temperature diffusion terms by using the pressure functions $p^\nu$ in place of
the unknown variables densities $n^\nu$.
-
This paper concerns the initial-boundary value problem to 2D micropolar
equations without angular viscosity in a smooth bounded domain. It is shown
that such a system admits a unique and global weak solution. The main idea of
this paper is to fully exploit the structure of this system and establish high
order estimates via introducing an auxiliary field which is at the energy level
of one order lower than micro-rotation.
-
We present initial results from the first systematic survey of luminous
$z\sim 5.5$ quasars. Quasars at $z \sim$ 5.5, the post-reionization epoch, are
crucial tools to explore the evolution of intergalactic medium, quasar
evolution and the early super-massive black hole growth. However, it has been
very challenging to select quasars at redshifts 5.3 $\le z \le$ 5.7 using
conventional color selections, due to their similar optical colors to late-type
stars, especially M dwarfs, resulting in a glaring redshift gap in quasar
redshift distributions. We develop a new selection technique for $z \sim$ 5.5
quasars based on optical, near-IR and mid-IR photometric data from Sloan
Digital Sky Survey (SDSS), UKIRT InfraRed Deep Sky Surveys - Large Area Survey
(ULAS), VISTA Hemisphere Survey (VHS) and Wide field Infrared Survey Explorer
(WISE). From our pilot observations in SDSS-ULAS/VHS area, we have discovered
15 new quasars at 5.3 $\le z \le$ 5.7 and 6 new lower redshift quasars, with
SDSS z band magnitude brighter than 20.5. Including other two $z \sim$ 5.5
quasars already published in our previous work, we now construct an uniform
quasar sample at 5.3 $\le z \le$ 5.7 with 17 quasars in a $\sim$ 4800 square
degree survey area. For further application in a larger survey area, we apply
our selection pipeline to do a test selection by using the new wide field J
band photometric data from a preliminary version of the UKIRT Hemisphere Survey
(UHS). We successfully discover the first UHS selected $z \sim$ 5.5 quasar.
-
Besides independent learning, human learning process is highly improved by
summarizing what has been learned, communicating it with peers, and
subsequently fusing knowledge from different sources to assist the current
learning goal. This collaborative learning procedure ensures that the knowledge
is shared, continuously refined, and concluded from different perspectives to
construct a more profound understanding. The idea of knowledge transfer has led
to many advances in machine learning and data mining, but significant
challenges remain, especially when it comes to reinforcement learning,
heterogeneous model structures, and different learning tasks. Motivated by
human collaborative learning, in this paper we propose a collaborative deep
reinforcement learning (CDRL) framework that performs adaptive knowledge
transfer among heterogeneous learning agents. Specifically, the proposed CDRL
conducts a novel deep knowledge distillation method to address the
heterogeneity among different learning tasks with a deep alignment network.
Furthermore, we present an efficient collaborative Asynchronous Advantage
Actor-Critic (cA3C) algorithm to incorporate deep knowledge distillation into
the online training of agents, and demonstrate the effectiveness of the CDRL
framework using extensive empirical evaluation on OpenAI gym.
-
Complex activity recognition is challenging due to the inherent uncertainty
and diversity of performing a complex activity. Normally, each instance of a
complex activity has its own configuration of atomic actions and their temporal
dependencies. We propose in this paper an atomic action-based Bayesian model
that constructs Allen's interval relation networks to characterize complex
activities with structural varieties in a probabilistic generative way: By
introducing latent variables from the Chinese restaurant process, our approach
is able to capture all possible styles of a particular complex activity as a
unique set of distributions over atomic actions and relations. We also show
that local temporal dependencies can be retained and are globally consistent in
the resulting interval network. Moreover, network structure can be learned from
empirical data. A new dataset of complex hand activities has been constructed
and made publicly available, which is much larger in size than any existing
datasets. Empirical evaluations on benchmark datasets as well as our in-house
dataset demonstrate the competitiveness of our approach.
-
Multispectral pedestrian detection is essential for around-the-clock
applications, e.g., surveillance and autonomous driving. We deeply analyze
Faster R-CNN for multispectral pedestrian detection task and then model it into
a convolutional network (ConvNet) fusion problem. Further, we discover that
ConvNet-based pedestrian detectors trained by color or thermal images
separately provide complementary information in discriminating human instances.
Thus there is a large potential to improve pedestrian detection by using color
and thermal images in DNNs simultaneously. We carefully design four ConvNet
fusion architectures that integrate two-branch ConvNets on different DNNs
stages, all of which yield better performance compared with the baseline
detector. Our experimental results on KAIST pedestrian benchmark show that the
Halfway Fusion model that performs fusion on the middle-level convolutional
features outperforms the baseline method by 11% and yields a missing rate 3.5%
lower than the other proposed architectures.
-
We present the discovery of nine quasars at $z\sim6$ identified in the Sloan
Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$
quasars in the SDSS footprint. Our final sample consists of 52 quasars at
$5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from
11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10
quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS
overlap regions (regions with two or more imaging scans), and 13 quasars down
to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a
wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample
is used to derive the quasar luminosity function (QLF) at $z\sim6$. After
combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from
the literature, we obtain the parameters for a double power-law fit to the QLF.
The bright-end slope $\beta$ of the QLF is well constrained to be
$\beta=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the
faint-end slope $\alpha$ and the characteristic magnitude $M_{1450}^{\ast}$ are
less well constrained, with $\alpha=-1.90_{-0.44}^{+0.58}$ and
$M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars,
parametrized as $\rho(M_{1450}<-26,z)=\rho(z=6)\,10^{k(z-6)}$, drops rapidly
from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming
an IGM clumping factor of $C=3$, we find that the observed quasar population
cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\%
confidence. Quasars may still provide a significant fraction of the required
photons, although much larger samples of faint quasars are needed for more
stringent constraints on the quasar contribution to reionization.
-
In this paper, we study the quasineutral limit and asymptotic behaviors for
the quantum Navier-Stokes-Possion equation. We apply a formal expansion
according to Debye length and derive the neutral incompressible Navier-Stokes
equation. To establish this limit mathematically rigorously, we derive uniform
(in Debye length) estimates for the remainders, for well-prepared initial data.
It is demonstrated that the quantum effect do play important roles in the
estimates and the norm introduced depends on the Planck constant $\hbar>0$.
-
This is the second paper in a series on a new luminous z ~ 5 quasar survey
using optical and near-infrared colors. Here we present a new determination of
the bright end of the quasar luminosity function (QLF) at z ~ 5. Combined our
45 new quasars with previously known quasars that satisfy our selections, we
construct the largest uniform luminous z ~ 5 quasar sample to date, with 99
quasars in the range 4.7 <= z < 5.4 and -29 < M1450 <= -26.8, within the Sloan
Digital Sky Survey (SDSS) footprint. We use a modified 1/Va method including
flux limit correction to derive a binned QLF, and we model the parametric QLF
using maximum likelihood estimation. With the faint-end slope of the QLF fixed
as alpha = -2.03 from previous deeper samples, the best fit of our QLF gives a
flatter bright end slope beta = -3.58+/-0.24 and a fainter break magnitude
M*1450 = -26.98+/-0.23 than previous studies at similar redshift. Combined with
previous work at lower and higher redshifts, our result is consistent with a
luminosity evolution and density evolution (LEDE) model. Using the best fit
QLF, the contribution of quasars to the ionizing background at z ~ 5 is found
to be 18% - 45% with a clumping factor C of 2 - 5. Our sample suggests an
evolution of radio loud fraction with optical luminosity but no obvious
evolution with redshift.
-
Consider the problem of finding a point in a unit $n$-dimensional
$\ell_p$-ball ($p\ge 2$) such that the minimum of the weighted Euclidean
distance from given $m$ points is maximized. We show in this paper that the
recent SDP-relaxation-based approximation algorithm [SIAM J. Optim. 23(4),
2264-2294, 2013] will not only provide the first theoretical approximation
bound of $\frac{1-O\left(\sqrt{ \ln(m)/n}\right)}{2}$, but also perform much
better in practice, if the SDP relaxation is removed and the optimal solution
of the SDP relaxation is replaced by a simple scalar matrix.
-
The ball-constrained weighted maximin dispersion problem $(\rm P_{ball})$ is
to find a point in an $n$-dimensional Euclidean ball such that the minimum of
the weighted Euclidean distance from given $m$ points is maximized. We propose
a new second-order cone programming relaxation for $(\rm P_{ball})$. Under the
condition $m\le n$, $(\rm P_{ball})$ is polynomial-time solvable since the new
relaxation is shown to be tight. In general, we prove that $({\rm P_{ball}})$
is NP-hard. Then, we propose a new randomized approximation algorithm for
solving $({\rm P_{ball}})$, which provides a new approximation bound of
$\frac{1-O(\sqrt{\ln(m)/n})}{2}$.
-
High-redshift quasars are important tracers of structure and evolution in the
early universe. However, they are very rare and difficult to find when using
color selection because of contamination from late-type dwarfs. High-redshift
quasar surveys based on only optical colors suffer from incompleteness and low
identification efficiency, especially at $z\gtrsim4.5$. We have developed a new
method to select $4.7\lesssim z \lesssim 5.4$ quasars with both high efficiency
and completeness by combining optical and mid-IR Wide-field Infrared Survey
Explorer (WISE) photometric data, and are conducting a luminous $z\sim5$ quasar
survey in the whole Sloan Digital Sky Survey (SDSS) footprint. We have
spectroscopically observed 99 out of 110 candidates with $z$-band magnitudes
brighter than 19.5 and 64 (64.6\%) of them are quasars with redshifts of
$4.4\lesssim z \lesssim 5.5$ and absolute magnitudes of $-29\lesssim M_{1450}
\lesssim -26.4$. In addition, we also observed 14 fainter candidates selected
with the same criteria and identified 8 (57.1\%) of them as quasars with
$4.7<z<5.4$ . Among 72 newly identified quasars, 12 of them are at $5.2 < z <
5.7$, which leads to an increase of $\sim$36\% of the number of known quasars
at this redshift range. More importantly, our identifications doubled the
number of quasars with $M_{1450}<-27.5$ at $z>4.5$, which will set strong
constraints on the bright end of the quasar luminosity function. We also expand
our method to select quasars at $z\gtrsim5.7$. In this paper we report the
discovery of four new luminous $z\gtrsim5.7$ quasars based on SDSS-WISE
selection.
-
In this paper, we propose a novel visual tracking framework that
intelligently discovers reliable patterns from a wide range of video to resist
drift error for long-term tracking tasks. First, we design a Discrete Fourier
Transform (DFT) based tracker which is able to exploit a large number of
tracked samples while still ensures real-time performance. Second, we propose a
clustering method with temporal constraints to explore and memorize consistent
patterns from previous frames, named as reliable memories. By virtue of this
method, our tracker can utilize uncontaminated information to alleviate
drifting issues. Experimental results show that our tracker performs favorably
against other state of-the-art methods on benchmark datasets. Furthermore, it
is significantly competent in handling drifts and able to robustly track
challenging long videos over 4000 frames, while most of others lose track at
early frames.
-
The sizes of interstellar grains are widely distributed, ranging from a few
angstroms to a few micrometers. The ultraviolet (UV) and optical extinction
constrains the dust in the size range of a couple hundredth micrometers to
several submicrometers. The near and mid infrared (IR) emission constrains the
nanometer-sized grains and angstrom-sized very large molecules. However, the
quantity and size distribution of micrometer-sized grains remain unknown as
they are gray in the UV/optical extinction and they are too cold and emit too
little in the IR to be detected by IRAS, Spitzer, or Herschel. In this work, we
employ the ~3-8 micron mid-IR extinction which is flat in both diffuse and
dense regions to constrain the quantity, size, and composition of the
micron-sized grain component. We find that, together with nano- and
submicron-sized silicate and graphite (as well as PAHs), micron-sized graphite
grains with C/H=137 ppm and a mean size of ~1.2 micron closely fit the observed
interstellar extinction of the Galactic diffuse interstellar medium from the
far-UV to the mid-IR as well as the near-IR to millimeter thermal emission
obtained by COBE/DIRBE, COBE/FIRAS, and Planck up to lambda < 1000 micron. The
micron-sized graphite component accounts for ~14.6% of the total dust mass and
~2.5% of the total IR emission.
-
The interstellar medium (ISM) seems to have a significant surplus of oxygen
which was dubbed as the "O crisis": independent of the adopted interstellar
reference abundance, the total number of O atoms depleted from the gas phase
far exceeds that tied up in solids by as much as ~160ppm of O/H. Recently, it
has been hypothesized that the missing O could be hidden in micrometer-sized
H2O ice grains. We examine this hypothesis by comparing the infrared (IR)
extinction and far-IR emission arising from these grains with that observed in
the Galactic diffuse ISM. We find that it is possible for the diffuse ISM to
accommodate ~160ppm of O/H in micron-sized H2O ice grains without violating the
observational constraints including the absence of the 3.1micron O-H absorption
feature. More specifically, H2O ice grains of radii ~4micron and O/H = 160 ppm
are capable of accounting for the observed flat extinction at ~ 3-8 micron and
produce no excessive emission in the far-IR. These grains could be present in
the diffuse ISM through rapid exchange of material between dense molecular
clouds where they form and diffuse clouds where they are destroyed by
photosputtering.
-
The uniform quadratic optimizatin problem (UQ) is a nonconvex quadratic
constrained quadratic programming (QCQP) sharing the same Hessian matrix. Based
on the second-order cone programming (SOCP) relaxation, we establish a new
sufficient condition to guarantee strong duality for (UQ) and then extend it to
(QCQP), which not only covers several well-known results in literature but also
partially gives answers to a few open questions. For convex constrained
nonconvex (UQ), we propose an improved approximation algorithm based on (SOCP).
Our approximation bound is dimensional independent. As an application, we
establish the first approximation bound for the problem of finding the
Chebyshev center of the intersection of several balls.
-
How dust absorbs and scatters starlight as a function of wavelength (known as
the interstellar extinction curve) is crucial for correcting for the effects of
dust extinction in inferring the true luminosity and colors of reddened
astrophysical objects. Together with the extinction spectral features, the
extinction curve contains important information about the dust size
distribution and composition. This review summarizes our current knowledge of
the dust extinction of the Milky Way, three Local Group galaxies (i.e., the
Small and Large Magellanic Clouds, and M31), and galaxies beyond the Local
Group.
-
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 $\mu$ such that $f(x) + \mu 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) + \mu
h(x)\ge0, ~\forall x\in \mathbb{R}^n$ for some $\mu\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.