-
Energetic Spectral-Element Time Marching Methods for Phase-Field Nonlinear Gradient Systems
Authors:
Shiqin Liu,
Haijun Yu
Abstract:
We propose two efficient energetic spectral-element methods in time for marching nonlinear gradient systems with the phase-field Allen--Cahn equation as an example: one fully implicit nonlinear method and one semi-implicit linear method. Different from other spectral methods in time using spectral Petrov-Galerkin or weighted Galerkin approximations, the presented implicit method employs an energet…
▽ More
We propose two efficient energetic spectral-element methods in time for marching nonlinear gradient systems with the phase-field Allen--Cahn equation as an example: one fully implicit nonlinear method and one semi-implicit linear method. Different from other spectral methods in time using spectral Petrov-Galerkin or weighted Galerkin approximations, the presented implicit method employs an energetic variational Galerkin form that can maintain the mass conservation and energy dissipation property of the continuous dynamical system. Another advantage of this method is its superconvergence. A high-order extrapolation is adopted for the nonlinear term to get the semi-implicit method. The semi-implicit method does not have superconvergence, but can be improved by a few Picard-like iterations to recover the superconvergence of the implicit method. Numerical experiments verify that the method using Legendre elements of degree three outperforms the 4th-order implicit-explicit backward differentiation formula and the 4th-order exponential time difference Runge-Kutta method, which were known to have best performances in solving phase-field equations. In addition to the standard Allen--Cahn equation, we also apply the method to a conservative Allen--Cahn equation, in which the conservation of discrete total mass is verified. The applications of the proposed methods are not limited to phase-field Allen--Cahn equations. They are suitable for solving general, large-scale nonlinear dynamical systems.
△ Less
Submitted 23 June, 2024;
originally announced June 2024.
-
Foliation of area minimizing hypersurfaces in asymptotically flat manifolds and Schoen's conjecture
Authors:
Shihang He,
Yuguang Shi,
Haobin Yu
Abstract:
In this paper, we demonstrate that any asymptotically flat manifold $(M^n, g)$ with $4\leq n\leq 7$ can be foliated by a family of area-minimizing hypersurfaces, each of which is asymptotic to Cartesian coordinate hyperplanes defined at an end of $(M^n, g)$. As an application of this foliation, we show that for any asymptotically flat manifold $(M^n, g)$ with $4\leq n\leq 7$, nonnegative scalar cu…
▽ More
In this paper, we demonstrate that any asymptotically flat manifold $(M^n, g)$ with $4\leq n\leq 7$ can be foliated by a family of area-minimizing hypersurfaces, each of which is asymptotic to Cartesian coordinate hyperplanes defined at an end of $(M^n, g)$. As an application of this foliation, we show that for any asymptotically flat manifold $(M^n, g)$ with $4\leq n\leq 7$, nonnegative scalar curvature and positive mass, the solution of free boundary problem for area-minimizing hypersurface in coordinate cylinder $C_{R_i}$ in $(M^n, g)$ either does not exist or drifts to infinity of $(M^n, g)$ as $R_i$ tends to infinity. Additionally, we introduce a concept of globally minimizing hypersurface in $(M^n, g)$, and verify a version of the Schoen Conjecture.
△ Less
Submitted 23 June, 2024;
originally announced June 2024.
-
Non-geodesically-convex optimization in the Wasserstein space
Authors:
Hoang Phuc Hau Luu,
Hanlin Yu,
Bernardo Williams,
Petrus Mikkola,
Marcelo Hartmann,
Kai Puolamäki,
Arto Klami
Abstract:
We study a class of optimization problems in the Wasserstein space (the space of probability measures) where the objective function is \emph{nonconvex} along generalized geodesics. When the regularization term is the negative entropy, the optimization problem becomes a sampling problem where it minimizes the Kullback-Leibler divergence between a probability measure (optimization variable) and a ta…
▽ More
We study a class of optimization problems in the Wasserstein space (the space of probability measures) where the objective function is \emph{nonconvex} along generalized geodesics. When the regularization term is the negative entropy, the optimization problem becomes a sampling problem where it minimizes the Kullback-Leibler divergence between a probability measure (optimization variable) and a target probability measure whose logarithmic probability density is a nonconvex function. We derive multiple convergence insights for a novel {\em semi Forward-Backward Euler scheme} under several nonconvex (and possibly nonsmooth) regimes. Notably, the semi Forward-Backward Euler is just a slight modification of the Forward-Backward Euler whose convergence is -- to our knowledge -- still unknown in our very general non-geodesically-convex setting.
△ Less
Submitted 1 June, 2024;
originally announced June 2024.
-
Dihedral-product groups of nonabelian simple groups
Authors:
Hao Yu
Abstract:
Given a finite group $G$, a group $X$ is called a dihedral-product group of $G$ if $X=GD$ where $D$ is a finite dihedral group. In particular, it is called a dihedral-skew product group of $G$ if $G\cap D=1$ and $D$ has a trivial core in $X$. In this paper, dihedral-skew product groups of all nonabelian simple groups are classified and dihedral-product groups of all nonabelian simple groups are ch…
▽ More
Given a finite group $G$, a group $X$ is called a dihedral-product group of $G$ if $X=GD$ where $D$ is a finite dihedral group. In particular, it is called a dihedral-skew product group of $G$ if $G\cap D=1$ and $D$ has a trivial core in $X$. In this paper, dihedral-skew product groups of all nonabelian simple groups are classified and dihedral-product groups of all nonabelian simple groups are characterized.
△ Less
Submitted 26 April, 2024; v1 submitted 23 April, 2024;
originally announced April 2024.
-
Fast OMP for Exact Recovery and Sparse Approximation
Authors:
Huiyuan Yu,
Jia He,
Maggie Cheng
Abstract:
Orthogonal Matching Pursuit (OMP) has been a powerful method in sparse signal recovery and approximation. However OMP suffers computational issue when the signal has large number of non-zeros. This paper advances OMP in two fronts: it offers a fast algorithm for the orthogonal projection of the input signal at each iteration, and a new selection criterion for making the greedy choice, which reduce…
▽ More
Orthogonal Matching Pursuit (OMP) has been a powerful method in sparse signal recovery and approximation. However OMP suffers computational issue when the signal has large number of non-zeros. This paper advances OMP in two fronts: it offers a fast algorithm for the orthogonal projection of the input signal at each iteration, and a new selection criterion for making the greedy choice, which reduces the number of iterations it takes to recover the signal. The proposed modifications to OMP directly reduce the computational complexity. Experiment results show significant improvement over the classical OMP in computation time. The paper also provided a sufficient condition for exact recovery under the new greedy choice criterion. For general signals that may not have sparse representations, the paper provides a bound for the approximation error. The approximation error is at the same order as OMP but is obtained within fewer iterations and less time.
△ Less
Submitted 29 March, 2024;
originally announced April 2024.
-
Optimizing Vaccine Site Locations While Considering Travel Inconvenience and Public Health Outcomes
Authors:
Suyanpeng Zhang,
Sze-chuan Suen,
Han Yu,
Maged Dessouky,
Fernando Ordonez
Abstract:
During the COVID-19 pandemic, there were over three million infections in Los Angeles County (LAC). To facilitate distribution when vaccines first became available, LAC set up six mega-sites for dispensing a large number of vaccines to the public. To understand if another choice of mega-site location would have improved accessibility and health outcomes, and to provide insight into future vaccine…
▽ More
During the COVID-19 pandemic, there were over three million infections in Los Angeles County (LAC). To facilitate distribution when vaccines first became available, LAC set up six mega-sites for dispensing a large number of vaccines to the public. To understand if another choice of mega-site location would have improved accessibility and health outcomes, and to provide insight into future vaccine allocation problems, we propose a multi-objective mixed integer linear programming model that balances travel convenience, infection reduction, and equitable distribution. We provide a tractable objective formulation that effectively proxies real-world public health goals of reducing infections while considering travel inconvenience and equitable distribution of resources. Compared with the solution empirically used in LAC in 2020, we recommend more dispersed mega-site locations that result in a 28% reduction in travel inconvenience and avert an additional 1,000 infections.
△ Less
Submitted 26 March, 2024;
originally announced March 2024.
-
Positivity-preserving and energy-dissipating discontinuous Galerkin methods for nonlinear nonlocal Fokker-Planck equations
Authors:
José A. Carrillo,
Hailiang Liu,
Hui Yu
Abstract:
This paper is concerned with structure-preserving numerical approximations for a class of nonlinear nonlocal Fokker-Planck equations, which admit a gradient flow structure and find application in diverse contexts. The solutions, representing density distributions, must be non-negative and satisfy a specific energy dissipation law. We design an arbitrary high-order discontinuous Galerkin (DG) metho…
▽ More
This paper is concerned with structure-preserving numerical approximations for a class of nonlinear nonlocal Fokker-Planck equations, which admit a gradient flow structure and find application in diverse contexts. The solutions, representing density distributions, must be non-negative and satisfy a specific energy dissipation law. We design an arbitrary high-order discontinuous Galerkin (DG) method tailored for these model problems. Both semi-discrete and fully discrete schemes are shown to admit the energy dissipation law for non-negative numerical solutions. To ensure the preservation of positivity in cell averages at all time steps, we introduce a local flux correction applied to the DDG diffusive flux. Subsequently, a hybrid algorithm is presented, utilizing a positivity-preserving limiter, to generate positive and energy-dissipating solutions. Numerical examples are provided to showcase the high resolution of the numerical solutions and the verified properties of the DG schemes.
△ Less
Submitted 22 March, 2024;
originally announced March 2024.
-
Event-triggered Boundary Control of Mixed-autonomy Traffic
Authors:
Yihuai Zhang,
Huan Yu
Abstract:
Control problems of mixed-autonomy traffic system consisting of both Human-driven Vehicles (HV) and Autonomous Vehicles (AV) have gained increasing attention. This paper is focused on suppressing traffic oscillations of the mixed-autonomy traffic system using boundary control design. The mixed traffic dynamics are described by a 4 x 4 hyperbolic partial differential equations (PDE) which governs p…
▽ More
Control problems of mixed-autonomy traffic system consisting of both Human-driven Vehicles (HV) and Autonomous Vehicles (AV) have gained increasing attention. This paper is focused on suppressing traffic oscillations of the mixed-autonomy traffic system using boundary control design. The mixed traffic dynamics are described by a 4 x 4 hyperbolic partial differential equations (PDE) which governs propagation of four properties in traffic including density of HV, density of AV, friction between two classes of vehicles from driving interactions, and averaged velocity. We propose event-triggered boundary control design since control signal of traffic light on ramp or varying speed limit cannot be updated in a continuous time fashion. We apply event-triggered mechanism for a PDE backstepping controller and obtain dynamic triggering condition. Lyapunov analysis is conducted to prove the exponential stability of the closed loop system with the event-triggered controller. Numerical simulation demonstrates how car-following spacing of AV affects event-triggering mechanism of control input in mixed-autonomy traffic.
△ Less
Submitted 21 March, 2024;
originally announced March 2024.
-
Frequency-Reactive Power Optimization Strategy of Grid-forming Offshore Wind Farm Using DRU-HVDC Transmission
Authors:
Zhekai Li,
Kun Han,
Xu Cai,
Renxin Yang,
Haotian Yu,
Kepeng Xia,
Lulu Liu
Abstract:
The diode rectifier unit-based high voltage direct current (DRU-HVDC) transmission with grid-forming (GFM) wind turbine is becoming a promising scheme for offshore wind farm(OWF) integration due to its high reliability and low cost. In this scheme, the AC network of the OWF and the DRU has completely different synchronization mechanisms and power flow characteristics from the traditional power sys…
▽ More
The diode rectifier unit-based high voltage direct current (DRU-HVDC) transmission with grid-forming (GFM) wind turbine is becoming a promising scheme for offshore wind farm(OWF) integration due to its high reliability and low cost. In this scheme, the AC network of the OWF and the DRU has completely different synchronization mechanisms and power flow characteristics from the traditional power system. To optimize the power flow and reduce the net loss, this paper carries out the power flow modeling and optimization analysis for the DRU-HVDC transmission system with grid-forming OWFs. The influence of the DRU and the GFM wind turbines on the power flow of the system is analyzed. On this basis, improved constraint conditions are proposed and an optimal power flow (OPF) method is established. This method can minimize the power loss by adjusting the reactive power output of each wind turbine and internal network frequency. Finally, based on MATLAB, this paper uses YALMIP toolkit and CPLEX mathematical solver to realize the programming solution of the OPF model proposed in this paper. The results show that the proposed optimization strategy can effectively reduce the power loss of the entire OWF and the transmission system with an optimization ratio of network losses exceeding 25.3%.
△ Less
Submitted 16 March, 2024;
originally announced March 2024.
-
Counting rationals and diophantine approximation in missing-digit Cantor sets
Authors:
Sam Chow,
Peter Varju,
Han Yu
Abstract:
We establish a new upper bound for the number of rationals up to a given height in a missing-digit set, making progress towards a conjecture of Broderick, Fishman, and Reich. This enables us to make novel progress towards another conjecture of those authors about the corresponding intrinsic diophantine approximation problem. Moreover, we make further progress towards conjectures of Bugeaud--Durand…
▽ More
We establish a new upper bound for the number of rationals up to a given height in a missing-digit set, making progress towards a conjecture of Broderick, Fishman, and Reich. This enables us to make novel progress towards another conjecture of those authors about the corresponding intrinsic diophantine approximation problem. Moreover, we make further progress towards conjectures of Bugeaud--Durand and Levesley--Salp--Velani on the distribution of diophantine exponents in missing-digit sets.
A key tool in our study is Fourier $\ell^1$ dimension introduced by the last named author in [H. Yu, Rational points near self-similar sets, arXiv:2101.05910]. An important technical contribution of the paper is a method to compute this quantity.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
Sharp variance estimator and causal bootstrap in stratified randomized experiments
Authors:
Haoyang Yu,
Ke Zhu,
Hanzhong Liu
Abstract:
The design-based finite-population asymptotic theory provides a normal approximation for the sampling distribution of the average treatment effect estimator in stratified randomized experiments. The asymptotic variance could be estimated by a Neyman-type conservative variance estimator. However, the variance estimator can be overly conservative, and the asymptotic theory may fail in small samples.…
▽ More
The design-based finite-population asymptotic theory provides a normal approximation for the sampling distribution of the average treatment effect estimator in stratified randomized experiments. The asymptotic variance could be estimated by a Neyman-type conservative variance estimator. However, the variance estimator can be overly conservative, and the asymptotic theory may fail in small samples. To solve these issues, we propose a sharp variance estimator for the weighted difference-in-means in stratified randomized experiments. Furthermore, we propose two causal bootstrap procedures to more accurately approximate the sampling distribution of the weighted difference-in-means estimator. The first causal bootstrap procedure is based on rank-preserving imputation and we prove its second-order refinement over normal approximation. The second causal bootstrap procedure is based on constant-treatment-effect imputation and is applicable in paired experiments. We prove its validity even when the assumption of constant treatment effect is violated for the true potential outcomes. Our analysis is randomization-based or design-based by conditioning on the potential outcomes, with treatment assignment being the sole source of randomness. Numerical studies and two real data applications demonstrate advantages of our proposed methods in finite samples.
△ Less
Submitted 26 June, 2024; v1 submitted 29 January, 2024;
originally announced January 2024.
-
$L_x^p\rightarrow L^q_{x,u}$ estimates for dilated averages over planar curves
Authors:
Junfeng Li,
Naijia Liu,
Zengjian Lou,
Haixia Yu
Abstract:
In this paper, we consider the $L_x^p(\mathbb{R}^2)\rightarrow L_{x,u}^q(\mathbb{R}^2\times [1,2])$ estimate for the operator $T$ along a dilated plane curve $(ut,uγ(t))$, where
$$Tf(x,u):=\int_{0}^{1}f(x_1-ut,x_2-u γ(t))\,\textrm{d}t,$$
$x:=(x_1,x_2)$ and $γ$ is a general plane curve satisfying some suitable smoothness and curvature conditions. We show that $T$ is $L_x^p(\mathbb{R}^2)$ to…
▽ More
In this paper, we consider the $L_x^p(\mathbb{R}^2)\rightarrow L_{x,u}^q(\mathbb{R}^2\times [1,2])$ estimate for the operator $T$ along a dilated plane curve $(ut,uγ(t))$, where
$$Tf(x,u):=\int_{0}^{1}f(x_1-ut,x_2-u γ(t))\,\textrm{d}t,$$
$x:=(x_1,x_2)$ and $γ$ is a general plane curve satisfying some suitable smoothness and curvature conditions. We show that $T$ is $L_x^p(\mathbb{R}^2)$ to $L_{x,u}^q(\mathbb{R}^2\times [1,2])$ bounded whenever $(\frac{1}{p},\frac{1}{q})\in \square \cup \{(0,0)\}\cup \{(\frac{2}{3},\frac{1}{3})\}$ and $1+(1 +ω)(\frac{1}{q}-\frac{1}{p})>0$, where the trapezium $\square:=\{(\frac{1}{p},\frac{1}{q}):\ \frac{2}{p}-1\leq\frac{1}{q}\leq \frac{1}{p}, \frac{1}{q}>\frac{1}{3p}, \frac{1}{q}>\frac{1}{p}-\frac{1}{3}\}$ and $ω:=\limsup_{t\rightarrow 0^{+}}\frac{\ln|γ(t)|}{\ln t}$. This result is sharp except for some borderline cases. On the other hand, in a smaller $(\frac{1}{p},\frac{1}{q})$ region, we also obtain the almost sharp estimate $T : L_x^p(\mathbb{R}^2)\rightarrow L_{x}^q(\mathbb{R}^2)$ uniformly for $u\in [1,2]$. These results imply that the operator $T$ has the so called local smoothing phenomenon, i.e., the $L^q$ integral about $u$ on $[1,2]$ extends the region of $(\frac{1}{p},\frac{1}{q})$ in uniform estimate $T : L_x^p(\mathbb{R}^2)\rightarrow L_{x}^q(\mathbb{R}^2)$.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
Extending Dynamic Origin-Destination Estimation to Understand Traffic Patterns During COVID-19
Authors:
Han Yu,
Suyanpeng Zhang,
Sze-chuan Suen,
Maged Dessouky,
Fernando Ordonez
Abstract:
Estimating dynamic Origin-Destination (OD) traffic flow is crucial for understanding traffic patterns and the traffic network. While dynamic origin-destination estimation (DODE) has been studied for decades as a useful tool for estimating traffic flow, few existing models have considered its potential in evaluating the influence of policy on travel activity. This paper proposes a data-driven appro…
▽ More
Estimating dynamic Origin-Destination (OD) traffic flow is crucial for understanding traffic patterns and the traffic network. While dynamic origin-destination estimation (DODE) has been studied for decades as a useful tool for estimating traffic flow, few existing models have considered its potential in evaluating the influence of policy on travel activity. This paper proposes a data-driven approach to estimate OD traffic flow using sensor data on highways and local roads. We extend prior DODE models to improve accuracy and realism in order to estimate how policies affect OD traffic flow in large urban networks. We applied our approach to a case study in Los Angeles County, where we developed a traffic network, estimated OD traffic flow between health districts during COVID-19, and analyzed the relationship between OD traffic flow and demographic characteristics such as income. Our findings demonstrate that the proposed approach provides valuable insights into traffic flow patterns and their underlying demographic factors for a large-scale traffic network. Specifically, our approach allows for evaluating the impact of policy changes on travel activity. The approach has practical applications for transportation planning and traffic management, enabling a better understanding of traffic flow patterns and the impact of policy changes on travel activity.
△ Less
Submitted 18 January, 2024;
originally announced January 2024.
-
Robust Boundary Stabilization of Stochastic Hyperbolic PDEs
Authors:
Yihuai Zhang,
Jean Auriol,
Huan Yu
Abstract:
This paper proposes a backstepping boundary control design for robust stabilization of linear first-order coupled hyperbolic partial differential equations (PDEs) with Markov-jumping parameters. The PDE system consists of 4 X 4 coupled hyperbolic PDEs whose first three characteristic speeds are positive and the last one is negative. We first design a full-state feedback boundary control law for a…
▽ More
This paper proposes a backstepping boundary control design for robust stabilization of linear first-order coupled hyperbolic partial differential equations (PDEs) with Markov-jumping parameters. The PDE system consists of 4 X 4 coupled hyperbolic PDEs whose first three characteristic speeds are positive and the last one is negative. We first design a full-state feedback boundary control law for a nominal, deterministic system using the backstepping method. Then, by applying Lyapunov analysis methods, we prove that the nominal backstepping control law can stabilize the PDE system with Markov jumping parameters if the nominal parameters are sufficiently close to the stochastic ones on average. The mean-square exponential stability conditions are theoretically derived and then validated via numerical simulations.
△ Less
Submitted 27 December, 2023;
originally announced December 2023.
-
A Note on Stability in Asynchronous Stochastic Approximation without Communication Delays
Authors:
Huizhen Yu,
Yi Wan,
Richard S. Sutton
Abstract:
In this paper, we study asynchronous stochastic approximation algorithms without communication delays. Our main contribution is a stability proof for these algorithms that extends a method of Borkar and Meyn by accommodating more general noise conditions. We also derive convergence results from this stability result and discuss their application in important average-reward reinforcement learning p…
▽ More
In this paper, we study asynchronous stochastic approximation algorithms without communication delays. Our main contribution is a stability proof for these algorithms that extends a method of Borkar and Meyn by accommodating more general noise conditions. We also derive convergence results from this stability result and discuss their application in important average-reward reinforcement learning problems.
△ Less
Submitted 22 December, 2023;
originally announced December 2023.
-
Deciding Foot-sortability and Minimal 2-bounded Non-foot-sortable Sock Orderings
Authors:
Hung-Hsun Hans Yu
Abstract:
A sock ordering is a sequence of socks with different colors. A sock ordering is foot-sortable if the sequence of socks can be sorted by a stack so that socks with the same color form a contiguous block. The problem of deciding whether a given sock ordering is foot-sortable was first considered by Defant and Kravitz, who resolved the case for alignment-free 2-uniform sock orderings. In this paper,…
▽ More
A sock ordering is a sequence of socks with different colors. A sock ordering is foot-sortable if the sequence of socks can be sorted by a stack so that socks with the same color form a contiguous block. The problem of deciding whether a given sock ordering is foot-sortable was first considered by Defant and Kravitz, who resolved the case for alignment-free 2-uniform sock orderings. In this paper, we resolve the problem in a more general setting, where each color appears in the sock ordering at most twice. A key component of the argument is a fast algorithm that determines the foot-sortability of a sock ordering of length $N$ in time $O(N\log N)$, which is also an interesting result on its own.
△ Less
Submitted 21 December, 2023;
originally announced December 2023.
-
Subsonic time-periodic solution to damped compressible Euler equations with large entropy
Authors:
Peng Qu,
Huimin Yu,
Xiaomin Zhang
Abstract:
In this paper, one-dimensional nonisentropic compressible Euler equations with linear damping $α(x)ρu$ are analyzed.~We want to explore the conditions under which a subsonic temporal periodic boundary can trigger a time-periodic $C^{1}$ solution. To achieve this aim, we use a technically constructed iteration scheme and give the sufficient conditions to guarantee the existence, uniqueness and stab…
▽ More
In this paper, one-dimensional nonisentropic compressible Euler equations with linear damping $α(x)ρu$ are analyzed.~We want to explore the conditions under which a subsonic temporal periodic boundary can trigger a time-periodic $C^{1}$ solution. To achieve this aim, we use a technically constructed iteration scheme and give the sufficient conditions to guarantee the existence, uniqueness and stability of the $C^{1}$ time-periodic solutions on the perturbation of a subsonic Fanno flow.~It is worthy to be pointed out that the entropy exhibits large amplitude under the assumption that the inflow sound speed is small.~However, it is crucial to assume that the boundary conditions possess a kind of dissipative structure at least on one side, which is used to cancel the nonlinear accelerating effect in the system.~The results indicate that the time-periodic feedback boundary control with dissipation can stabilize the nonisentropic compressible Euler equations around the Fanno flows.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
High Probability Guarantees for Random Reshuffling
Authors:
Hengxu Yu,
Xiao Li
Abstract:
We consider the stochastic gradient method with random reshuffling ($\mathsf{RR}$) for tackling smooth nonconvex optimization problems. $\mathsf{RR}$ finds broad applications in practice, notably in training neural networks. In this work, we first investigate the concentration property of $\mathsf{RR}$'s sampling procedure and establish a new high probability sample complexity guarantee for drivin…
▽ More
We consider the stochastic gradient method with random reshuffling ($\mathsf{RR}$) for tackling smooth nonconvex optimization problems. $\mathsf{RR}$ finds broad applications in practice, notably in training neural networks. In this work, we first investigate the concentration property of $\mathsf{RR}$'s sampling procedure and establish a new high probability sample complexity guarantee for driving the gradient (without expectation) below $\varepsilon$, which effectively characterizes the efficiency of a single $\mathsf{RR}$ execution. Our derived complexity matches the best existing in-expectation one up to a logarithmic term while imposing no additional assumptions nor changing $\mathsf{RR}$'s updating rule. Furthermore, by leveraging our derived high probability descent property and bound on the stochastic error, we propose a simple and computable stopping criterion for $\mathsf{RR}$ (denoted as $\mathsf{RR}$-$\mathsf{sc}$). This criterion is guaranteed to be triggered after a finite number of iterations, and then $\mathsf{RR}$-$\mathsf{sc}$ returns an iterate with its gradient below $\varepsilon$ with high probability. Moreover, building on the proposed stopping criterion, we design a perturbed random reshuffling method ($\mathsf{p}$-$\mathsf{RR}$) that involves an additional randomized perturbation procedure near stationary points. We derive that $\mathsf{p}$-$\mathsf{RR}$ provably escapes strict saddle points and efficiently returns a second-order stationary point with high probability, without making any sub-Gaussian tail-type assumptions on the stochastic gradient errors. Finally, we conduct numerical experiments on neural network training to support our theoretical findings.
△ Less
Submitted 7 December, 2023; v1 submitted 20 November, 2023;
originally announced November 2023.
-
Standard solutions of complex differential equations
Authors:
Janne Heittokangas,
Samu Pulkkinen,
Hui Yu,
Amine Zemirni
Abstract:
A meromorphic solution of a complex linear differential equation (with meromorphic coefficients) for which the value zero is the only possible finite deficient/deviated value is called a standard solution. Conditions for the existence and the number of standard solutions are discussed for various types of deficient and deviated values.
A meromorphic solution of a complex linear differential equation (with meromorphic coefficients) for which the value zero is the only possible finite deficient/deviated value is called a standard solution. Conditions for the existence and the number of standard solutions are discussed for various types of deficient and deviated values.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
Mean-Square Exponential Stabilization of Mixed-Autonomy Traffic PDE System
Authors:
Yihuai Zhang,
Huan Yu,
Jean Auriol,
Mike Pereira
Abstract:
Control of mixed-autonomy traffic where Human-driven Vehicles (HVs) and Autonomous Vehicles (AVs) coexist on the road has gained increasing attention over the recent decades. This paper addresses the boundary stabilization problem for mixed traffic on freeways. The traffic dynamics are described by uncertain coupled hyperbolic partial differential equations (PDEs) with Markov jumping parameters, w…
▽ More
Control of mixed-autonomy traffic where Human-driven Vehicles (HVs) and Autonomous Vehicles (AVs) coexist on the road has gained increasing attention over the recent decades. This paper addresses the boundary stabilization problem for mixed traffic on freeways. The traffic dynamics are described by uncertain coupled hyperbolic partial differential equations (PDEs) with Markov jumping parameters, which aim to address the distinctive driving strategies between AVs and HVs. Considering that the spacing policies of AVs vary in mixed traffic, the stochastic impact area of AVs is governed by a continuous Markov chain. The interactions between HVs and AVs such as overtaking or lane changing are mainly induced by impact areas. Using backstepping design, we develop a full-state feedback boundary control law to stabilize the deterministic system (nominal system). Applying Lyapunov analysis, we demonstrate that the nominal backstepping control law is able to stabilize the traffic system with Markov jumping parameters, provided the nominal parameters are sufficiently close to the stochastic ones on average. The mean-square exponential stability conditions are derived, and the results are validated by numerical simulations.
△ Less
Submitted 25 June, 2024; v1 submitted 24 October, 2023;
originally announced October 2023.
-
The Product of a Semi Dihedral Group And a Cyclic Group
Authors:
Hao Yu
Abstract:
Let $X(G)=GC$ be a group, where $G$ is a semi dihedral group and $C$ is a cyclic group such that $G\cap C=1$. In this paper, $X(G)$ will be characterized.
Let $X(G)=GC$ be a group, where $G$ is a semi dihedral group and $C$ is a cyclic group such that $G\cap C=1$. In this paper, $X(G)$ will be characterized.
△ Less
Submitted 30 September, 2023;
originally announced October 2023.
-
An Isomorphism Theorem for Arithmetic Complexes
Authors:
Luca Fiorindo,
Ethan Reed,
Shahriyar Roshan-Zamir,
Hongmiao Yu
Abstract:
We consider generalizations of certain arithmetic complexes appearing in work of Raicu and VandeBogert in connection with the study of stable sheaf cohomology on flag varieties. Defined over the ring of integer valued polynomials, we prove an isomorphism of these complexes as conjectured by Gao, Raicu, and VandeBogert. In particular, this gives a more conceptual proof of an identification between…
▽ More
We consider generalizations of certain arithmetic complexes appearing in work of Raicu and VandeBogert in connection with the study of stable sheaf cohomology on flag varieties. Defined over the ring of integer valued polynomials, we prove an isomorphism of these complexes as conjectured by Gao, Raicu, and VandeBogert. In particular, this gives a more conceptual proof of an identification between the stable sheaf cohomology of hook and two column partition Schur functors applied to the cotangent sheaf of projective space.
△ Less
Submitted 29 September, 2023;
originally announced September 2023.
-
On Schur Rings Over a Semi Dihedral Group
Authors:
Hao Yu
Abstract:
In this paper, we shall show that the only primitive Schur ring over a semi dihedral group is the trivial one and every semi dihedral subgroup is Burnside group, that is a primitive group containing a regular subgroup isomorphic to the semi dihedral subgroup is necessarily 2-transitive.
In this paper, we shall show that the only primitive Schur ring over a semi dihedral group is the trivial one and every semi dihedral subgroup is Burnside group, that is a primitive group containing a regular subgroup isomorphic to the semi dihedral subgroup is necessarily 2-transitive.
△ Less
Submitted 4 May, 2024; v1 submitted 13 September, 2023;
originally announced September 2023.
-
Skew Product Groups for 2-Groups of Maximal Class
Authors:
Wenjuan Luo,
Hao Yu
Abstract:
Skew morphisms, which generalise automorphisms for groups, provide a fundamental tool for the study of regular Cayley maps and, more generally, for finite groups with a complementary factorisation $X = GY$, where $Y$ is cyclic and core-free in $X$. $X$ is called the skew product group associated with $G$ and $Y$.
In this paper, we classify skew product groups for the maximal class 2-groups.
Skew morphisms, which generalise automorphisms for groups, provide a fundamental tool for the study of regular Cayley maps and, more generally, for finite groups with a complementary factorisation $X = GY$, where $Y$ is cyclic and core-free in $X$. $X$ is called the skew product group associated with $G$ and $Y$.
In this paper, we classify skew product groups for the maximal class 2-groups.
△ Less
Submitted 6 September, 2023;
originally announced September 2023.
-
$L^p$-improving bounds of maximal functions along planar curves
Authors:
Naijia Liu,
Haixia Yu
Abstract:
In this paper, we study the $L^p(\mathbb{R}^2)$-improving bounds, i.e., $L^p(\mathbb{R}^2)\rightarrow L^q(\mathbb{R}^2)$ estimates, of the maximal function $M_γ$ along a plane curve $(t,γ(t))$, where
$$M_γf(x_1,x_2):=\sup_{u\in [1,2]}\left|\int_{0}^{1}f(x_1-ut,x_2-u γ(t))\,\textrm{d}t\right|,$$
and $γ$ is a general plane curve satisfying some suitable smoothness and curvature conditions. We ob…
▽ More
In this paper, we study the $L^p(\mathbb{R}^2)$-improving bounds, i.e., $L^p(\mathbb{R}^2)\rightarrow L^q(\mathbb{R}^2)$ estimates, of the maximal function $M_γ$ along a plane curve $(t,γ(t))$, where
$$M_γf(x_1,x_2):=\sup_{u\in [1,2]}\left|\int_{0}^{1}f(x_1-ut,x_2-u γ(t))\,\textrm{d}t\right|,$$
and $γ$ is a general plane curve satisfying some suitable smoothness and curvature conditions. We obtain $M_γ : L^p(\mathbb{R}^2)\rightarrow L^q(\mathbb{R}^2)$ if $(\frac{1}{p},\frac{1}{q})\in Δ\cup \{(0,0)\}$ and $(\frac{1}{p},\frac{1}{q})$ satisfying $1+(1 +ω)(\frac{1}{q}-\frac{1}{p})>0$, where $Δ:=\{(\frac{1}{p},\frac{1}{q}):\ \frac{1}{2p}<\frac{1}{q}\leq \frac{1}{p}, \frac{1}{q}>\frac{3}{p}-1 \}$ and $ω:=\limsup_{t\rightarrow 0^{+}}\frac{\ln|γ(t)|}{\ln t}$. This result is sharp except for some borderline cases. As Hickman stated in [J. Funct. Anal. 270 (2016), pp. 560--608], this is a very different situation.
△ Less
Submitted 5 September, 2023;
originally announced September 2023.
-
On the absolute continuity of radial and linear projections of missing digits measures
Authors:
Han Yu
Abstract:
In this paper, we study the absolute continuity of radial projections of missing digits measures. We show that for large enough missing digits measures $λ$ on $\mathbb{R}^n,n\geq 2,$ for all $x\in\mathbb{R}^n\setminus \mathrm{supp}(λ),$ $Π_x(λ)$ is absolutely continuous with a density function in $L^2(S^{n-1}).$ Our method applies to linear projections as well. In particular, we show that for $λ$…
▽ More
In this paper, we study the absolute continuity of radial projections of missing digits measures. We show that for large enough missing digits measures $λ$ on $\mathbb{R}^n,n\geq 2,$ for all $x\in\mathbb{R}^n\setminus \mathrm{supp}(λ),$ $Π_x(λ)$ is absolutely continuous with a density function in $L^2(S^{n-1}).$ Our method applies to linear projections as well. In particular, we show that for $λ$ as above, the linearly projected measure $P_θ(λ)$ is absolutely continuous with a continuous density function for almost all directions $θ\in S^{n-1}.$ This implies a version of Palis' conjecture for missing digits sets.
△ Less
Submitted 3 September, 2023;
originally announced September 2023.
-
Missing digits points near manifolds
Authors:
Han Yu
Abstract:
We consider a problem concerning the distribution of points with missing digits coordinates that are close to non-degenerate analytic submanifolds. We show that large enough (to be specified in the paper) sets of points with missing digits coordinates distribute 'equally' around non-degenerate submanifolds. As a consequence, we show that intersecting those missing digits sets with non-degenerate s…
▽ More
We consider a problem concerning the distribution of points with missing digits coordinates that are close to non-degenerate analytic submanifolds. We show that large enough (to be specified in the paper) sets of points with missing digits coordinates distribute 'equally' around non-degenerate submanifolds. As a consequence, we show that intersecting those missing digits sets with non-degenerate submanifolds always achieve the optimal dimension reduction. On the other hand, we also prove that there is no lack of points with missing digits that are contained in non-degenerate submanifolds. Among the other results,
1. we prove that the pinned distance sets of those missing digits sets contain non-trivial intervals regardless of where the pin is.
2. we prove that for each $ε>0,$ for missing digits sets $K$ with large bases, simple digit sets (to be specified in the paper), and $\dim_{H} K>3/4+ε,$ the arithmetic product sets $K\cdot K$ contains non-trivial intervals.
△ Less
Submitted 31 August, 2023;
originally announced September 2023.
-
Frequency-domain criterion on the stabilizability for infinite-dimensional linear control systems
Authors:
Karl Kunisch,
Gengsheng Wang,
Huaiqiang Yu
Abstract:
A quantitative frequency-domain condition related to the exponential stabilizability for infinite-dimensional linear control systems is presented. It is proven that this condition is necessary and sufficient for the stabilizability of special systems, while it is a necessary condition for the stabilizability in general. Applications are provided.
A quantitative frequency-domain condition related to the exponential stabilizability for infinite-dimensional linear control systems is presented. It is proven that this condition is necessary and sufficient for the stabilizability of special systems, while it is a necessary condition for the stabilizability in general. Applications are provided.
△ Less
Submitted 6 March, 2024; v1 submitted 29 August, 2023;
originally announced August 2023.
-
Generic properties in free boundary problems
Authors:
Xavier Fernández-Real,
Hui Yu
Abstract:
In this work, we show the generic uniqueness of minimizers for a large class of energies, including the Alt-Caffarelli and Alt-Phillips functionals.
We then prove the generic regularity of free boundaries for minimizers of the one-phase Alt-Caffarelli and Alt-Phillips functionals, for a monotone family of boundary data $\{\varphi_t\}_{t\in(-1,1)}$. More precisely, we show that for a co-countable…
▽ More
In this work, we show the generic uniqueness of minimizers for a large class of energies, including the Alt-Caffarelli and Alt-Phillips functionals.
We then prove the generic regularity of free boundaries for minimizers of the one-phase Alt-Caffarelli and Alt-Phillips functionals, for a monotone family of boundary data $\{\varphi_t\}_{t\in(-1,1)}$. More precisely, we show that for a co-countable subset of $\{\varphi_t\}_{t\in(-1,1)}$, minimizers have smooth free boundaries in $\mathbb{R}^5$ for the Alt-Caffarelli and in $\mathbb{R}^3$ for the Alt-Phillips functional. In general dimensions, we show that the singular set is one dimension smaller than expected for almost every boundary datum in $\{\varphi_t\}_{t\in(-1,1)}$.
△ Less
Submitted 25 August, 2023;
originally announced August 2023.
-
KdV limit for the Vlasov-Poisson-Landau system
Authors:
Renjun Duan,
Dongcheng Yang,
Hongjun Yu
Abstract:
Two fundamental models in plasma physics are given by the Vlasov-Poisson-Landau system and the compressible Euler-Poisson system which both capture the complex dynamics of plasmas under the self-consistent electric field interactions at the kinetic and fluid levels, respectively. Although there have been extensive studies on the long wave limit of the Euler-Poisson system towards Korteweg-de Vries…
▽ More
Two fundamental models in plasma physics are given by the Vlasov-Poisson-Landau system and the compressible Euler-Poisson system which both capture the complex dynamics of plasmas under the self-consistent electric field interactions at the kinetic and fluid levels, respectively. Although there have been extensive studies on the long wave limit of the Euler-Poisson system towards Korteweg-de Vries equations, few results on this topic are known for the Vlasov-Poisson-Landau system due to the complexity of the system and its underlying multiscale feature. In this article, we derive and justify the Korteweg-de Vries equations from the Vlasov-Poisson-Landau system modelling the motion of ions under the Maxwell-Boltzmann relation. Specifically, under the Gardner-Morikawa transformation $$ (t,x,v)\to (δ^{\frac{3}{2}}t,δ^{\frac{1}{2}}(x-\sqrt{\frac{8}{3}}t),v) $$ with $ \varepsilon^{\frac{2}{3}}\leq δ\leq \varepsilon^{\frac{2}{5}}$ and $\varepsilon>0$ being the Knudsen number, we construct smooth solutions of the rescaled Vlasov-Poisson-Landau system over an arbitrary finite time interval that can converge uniformly to smooth solutions of the Korteweg-de Vries equations as $δ\to 0$. Moreover, the explicit rate of convergence in $δ$ is also obtained. The proof is obtained by an appropriately chosen scaling and the intricate weighted energy method through the micro-macro decomposition around local Maxwellians.
△ Less
Submitted 17 August, 2023;
originally announced August 2023.
-
MSAT: Matrix stability analysis tool for shock-capturing schemes
Authors:
Weijie Ren,
Wenjia Xie,
Ye Zhang,
Hang Yu,
Zhengyu Tian
Abstract:
The simulation of supersonic or hypersonic flows often suffers from numerical shock instabilities if the flow field contains strong shocks, limiting the further application of shock-capturing schemes. In this paper, we develop the unified matrix stability analysis method for schemes with three-point stencils and present MSAT, an open-source tool to quantitatively analyze the shock instability prob…
▽ More
The simulation of supersonic or hypersonic flows often suffers from numerical shock instabilities if the flow field contains strong shocks, limiting the further application of shock-capturing schemes. In this paper, we develop the unified matrix stability analysis method for schemes with three-point stencils and present MSAT, an open-source tool to quantitatively analyze the shock instability problem. Based on the finite-volume approach on the structured grid, MSAT can be employed to investigate the mechanism of the shock instability problem, evaluate the robustness of numerical schemes, and then help to develop robust schemes. Also, MSAT has the ability to analyze the practical simulation of supersonic or hypersonic flows, evaluate whether it will suffer from shock instabilities, and then assist in selecting appropriate numerical schemes accordingly. As a result, MSAT is a helpful tool that can investigate the shock instability problem and help to cure it.
△ Less
Submitted 15 August, 2023;
originally announced August 2023.
-
Traceability of Water Pollution: An Inversion Scheme Via Dynamic Complex Geometrical Optics Solutions
Authors:
Lingyun Qiu,
Zhongjing Wang,
Hui Yu,
Shenwen Yu
Abstract:
We investigate the identification of the time-dependent source term in the diffusion equation using boundary measurements. This facilitates tracing back the origins of environmental pollutants. Employing the concept of dynamic complex geometrical optics (CGO) solutions, a variational formulation of the inverse source problem is analyzed, leading to a proof of uniqueness result. Our proposed two-st…
▽ More
We investigate the identification of the time-dependent source term in the diffusion equation using boundary measurements. This facilitates tracing back the origins of environmental pollutants. Employing the concept of dynamic complex geometrical optics (CGO) solutions, a variational formulation of the inverse source problem is analyzed, leading to a proof of uniqueness result. Our proposed two-step reconstruction algorithm first determines the point source locations and subsequently reconstructs the Fourier components of the emission concentration functions. Numerical experiments on simulated data are conducted. The results demonstrate that the proposed two-step reconstruction algorithm can reliably reconstruct multiple point sources and accurately reconstruct the emission concentration functions. Additionally, by partitioning the algorithm into online and offline computations, and concentrating computational demand offline, real-time pollutant traceability becomes feasible. This method, applicable in various fields - especially those related to water pollution, can identify the source of a contaminant in the environment, thus serving as a valuable tool in environmental protection.
△ Less
Submitted 17 November, 2023; v1 submitted 11 August, 2023;
originally announced August 2023.
-
Numerical stability analysis of shock-capturing methods for strong shocks II: high-order finite-volume schemes
Authors:
Weijie Ren,
Wenjia Xie,
Ye Zhang,
Hang Yu,
Zhengyu Tian
Abstract:
The shock instability problem commonly arises in flow simulations involving strong shocks, particularly when employing high-order schemes, limiting their applications in hypersonic flow simulations. This study focuses on exploring the numerical characteristics and underlying mechanisms of shock instabilities in fifth-order finite-volume WENO schemes. To this end, for the first time, we have establ…
▽ More
The shock instability problem commonly arises in flow simulations involving strong shocks, particularly when employing high-order schemes, limiting their applications in hypersonic flow simulations. This study focuses on exploring the numerical characteristics and underlying mechanisms of shock instabilities in fifth-order finite-volume WENO schemes. To this end, for the first time, we have established the matrix stability analysis method for the fifth-order scheme. By predicting the evolution of perturbation errors in the exponential growth stage, this method provides quantitative insights into the behavior of shock-capturing and helps elucidate the mechanisms that cause shock instabilities. Results reveal that even dissipative solvers also suffer from shock instabilities when the spatial accuracy is increased to fifth-order. Further investigation indicates that this is due to the excessively high spatial accuracy of the WENO scheme near the numerical shock structure. Moreover, the shock instability problem of fifth-order schemes is demonstrated to be a multidimensional coupling problem. To stably capture strong shocks, it is crucial to have sufficient dissipation on transverse faces and ensure at least two points within the numerical shock structure in the direction perpendicular to the shock. The source location of instability is also clarified by the matrix stability analysis method, revealing that the instability arises from the numerical shock structure. Additionally, stability analysis demonstrates that local characteristic decomposition helps mitigate shock instabilities in high-order schemes, although the instability still persists. These conclusions pave the way for a better understanding of the shock instability in fifth-order schemes and provide guidance for the development of more reliable high-order shock-capturing methods for compressible flows with high Mach numbers.
△ Less
Submitted 7 August, 2023;
originally announced August 2023.
-
Tight Bound and Structural Theorem for Joints
Authors:
Ting-Wei Chao,
Hung-Hsun Hans Yu
Abstract:
A joint of a set of lines $\mathcal{L}$ in $\mathbb{F}^d$ is a point that is contained in $d$ lines with linearly independent directions. The joints problem asks for the maximum number of joints that are formed by $L$ lines. Guth and Katz showed that the number of joints is at most $O(L^{3/2})$ in $\mathbb{R}^3$ using polynomial method. This upper bound is met by the construction given by taking t…
▽ More
A joint of a set of lines $\mathcal{L}$ in $\mathbb{F}^d$ is a point that is contained in $d$ lines with linearly independent directions. The joints problem asks for the maximum number of joints that are formed by $L$ lines. Guth and Katz showed that the number of joints is at most $O(L^{3/2})$ in $\mathbb{R}^3$ using polynomial method. This upper bound is met by the construction given by taking the joints and the lines to be all the $d$-wise intersections and all the $(d-1)$-wise intersections of $M$ hyperplanes in general position. Furthermore, this construction is conjectured to be optimal.
In this paper, we verify the conjecture and show that this is the only optimal construction by using a more sophisticated polynomial method argument. This is the first tight bound and structural theorem obtained using this method. We also give a new definition of multiplicity that strengthens the main result of a previous work by Tidor, Zhao and the second author. Lastly, we relate the joints problem to some set-theoretic problems and prove conjectures of Bollobás and Eccles regarding partial shadows.
△ Less
Submitted 21 December, 2023; v1 submitted 28 July, 2023;
originally announced July 2023.
-
Kruskal--Katona-Type Problems via the Entropy Method
Authors:
Ting-Wei Chao,
Hung-Hsun Hans Yu
Abstract:
In this paper, we investigate several extremal combinatorics problems that ask for the maximum number of copies of a fixed subgraph given the number of edges. We call problems of this type Kruskal--Katona-type problems. Most of the problems that will be discussed in this paper are related to the joints problem. There are two main results in this paper. First, we prove that, in a $3$-edge-colored g…
▽ More
In this paper, we investigate several extremal combinatorics problems that ask for the maximum number of copies of a fixed subgraph given the number of edges. We call problems of this type Kruskal--Katona-type problems. Most of the problems that will be discussed in this paper are related to the joints problem. There are two main results in this paper. First, we prove that, in a $3$-edge-colored graph with $R$ red, $G$ green, $B$ blue edges, the number of rainbow triangles is at most $\sqrt{2RGB}$, which is sharp. Second, we give a generalization of the Kruskal--Katona theorem that implies many other previous generalizations. Both arguments use the entropy method, and the main innovation lies in a more clever argument that improves bounds given by Shearer's inequality.
△ Less
Submitted 2 June, 2024; v1 submitted 28 July, 2023;
originally announced July 2023.
-
Temporal periodic solutions to nonhomogeneous quasilinear hyperbolic equations driven by time-periodic boundary conditions
Authors:
Xixi Fang,
Peng Qu,
Huimin Yu
Abstract:
We consider the temporal periodic solutions to general nonhomogeneous quasilinear hyperbolic equations with a kind of weak diagonal dominant structure. Under the temporal periodic boundary conditions, the existence, stability and uniqueness of the time-periodic classical solutions are obtained.Moreover, the W2, regularity and stability around the time-periodic solutions are discussed. Our results…
▽ More
We consider the temporal periodic solutions to general nonhomogeneous quasilinear hyperbolic equations with a kind of weak diagonal dominant structure. Under the temporal periodic boundary conditions, the existence, stability and uniqueness of the time-periodic classical solutions are obtained.Moreover, the W2, regularity and stability around the time-periodic solutions are discussed. Our results reveal that the feedback boundary control with dissipative structure can stabilize the K-weakly diagonally dominant nonhomogeneous quasilinear hyperbolic system around the temporal periodic solution.
△ Less
Submitted 20 June, 2023; v1 submitted 16 June, 2023;
originally announced June 2023.
-
Data-driven optimal control under safety constraints using sparse Koopman approximation
Authors:
Hongzhe Yu,
Joseph Moyalan,
Umesh Vaidya,
Yongxin Chen
Abstract:
In this work we approach the dual optimal reach-safe control problem using sparse approximations of Koopman operator. Matrix approximation of Koopman operator needs to solve a least-squares (LS) problem in the lifted function space, which is computationally intractable for fine discretizations and high dimensions. The state transitional physical meaning of the Koopman operator leads to a sparse LS…
▽ More
In this work we approach the dual optimal reach-safe control problem using sparse approximations of Koopman operator. Matrix approximation of Koopman operator needs to solve a least-squares (LS) problem in the lifted function space, which is computationally intractable for fine discretizations and high dimensions. The state transitional physical meaning of the Koopman operator leads to a sparse LS problem in this space. Leveraging this sparsity, we propose an efficient method to solve the sparse LS problem where we reduce the problem dimension dramatically by formulating the problem using only the non-zero elements in the approximation matrix with known sparsity pattern. The obtained matrix approximation of the operators is then used in a dual optimal reach-safe problem formulation where a linear program with sparse linear constraints naturally appears. We validate our proposed method on various dynamical systems and show that the computation time for operator approximation is greatly reduced with high precision in the solutions.
△ Less
Submitted 1 June, 2023;
originally announced June 2023.
-
Solutions to the nonlinear obstacle problem with compact contact sets
Authors:
Simon Eberle,
Hui Yu
Abstract:
For the obstacle problem with a nonlinear operator, we characterize the space of global solutions with compact contact sets. This is achieved by constructing a bijection onto a class of quadratic polynomials describing the asymptotic behavior of solutions.
For the obstacle problem with a nonlinear operator, we characterize the space of global solutions with compact contact sets. This is achieved by constructing a bijection onto a class of quadratic polynomials describing the asymptotic behavior of solutions.
△ Less
Submitted 31 May, 2023;
originally announced May 2023.
-
Optimal subsampling for large scale Elastic-net regression
Authors:
Hang Yu,
Zhenxing Dou,
Zhiwei Chen,
Xiaomeng Yan
Abstract:
Datasets with sheer volume have been generated from fields including computer vision, medical imageology, and astronomy whose large-scale and high-dimensional properties hamper the implementation of classical statistical models. To tackle the computational challenges, one of the efficient approaches is subsampling which draws subsamples from the original large datasets according to a carefully-des…
▽ More
Datasets with sheer volume have been generated from fields including computer vision, medical imageology, and astronomy whose large-scale and high-dimensional properties hamper the implementation of classical statistical models. To tackle the computational challenges, one of the efficient approaches is subsampling which draws subsamples from the original large datasets according to a carefully-design task-specific probability distribution to form an informative sketch. The computation cost is reduced by applying the original algorithm to the substantially smaller sketch. Previous studies associated with subsampling focused on non-regularized regression from the computational efficiency and theoretical guarantee perspectives, such as ordinary least square regression and logistic regression. In this article, we introduce a randomized algorithm under the subsampling scheme for the Elastic-net regression which gives novel insights into L1-norm regularized regression problem. To effectively conduct consistency analysis, a smooth approximation technique based on alpha absolute function is firstly employed and theoretically verified. The concentration bounds and asymptotic normality for the proposed randomized algorithm are then established under mild conditions. Moreover, an optimal subsampling probability is constructed according to A-optimality. The effectiveness of the proposed algorithm is demonstrated upon synthetic and real data datasets.
△ Less
Submitted 29 May, 2023; v1 submitted 24 May, 2023;
originally announced May 2023.
-
The Product of a Generalized Quaternion Group And a Cyclic Group
Authors:
Shaofei Du,
Hao Yu,
Wenjuan Luo
Abstract:
Let $X(Q)=QC$ be a group, where $Q$ is a generalized quaternion group and $C$ is a cyclic group such that $Q\cap C=1$. In this paper, $X(Q)$ will be characterized and moreover, a complete classification for that will be given, provided $C$ is core-free. For the reason of self-constraint, in this paper a classification of the group $X(D)=DC$ is also given, where $D$ is a dihedral group and $C$ is a…
▽ More
Let $X(Q)=QC$ be a group, where $Q$ is a generalized quaternion group and $C$ is a cyclic group such that $Q\cap C=1$. In this paper, $X(Q)$ will be characterized and moreover, a complete classification for that will be given, provided $C$ is core-free. For the reason of self-constraint, in this paper a classification of the group $X(D)=DC$ is also given, where $D$ is a dihedral group and $C$ is a cyclic group such that $D\cap C=1$ and $C$ is core-free. Remind that the group $X(D)$ was recently classified in [12], based on a number of papers on skew-morphisms of dihedral groups. In this paper, a different approach from that in [12] will be used.
△ Less
Submitted 16 May, 2023; v1 submitted 15 May, 2023;
originally announced May 2023.
-
Numerical stability analysis of shock-capturing methods for strong shocks I: second-order MUSCL schemes
Authors:
Weijie Ren,
Wenjia Xie,
Ye Zhang,
Hang Yu,
Zhengyu Tian
Abstract:
Modern shock-capturing schemes often suffer from numerical shock anomalies if the flow field contains strong shocks, which may limit their further application in hypersonic flow computations. In the current study, we devote our efforts to exploring the primary numerical characteristics and the underlying mechanism of shock instability for second-order finite-volume schemes. To this end, we, for th…
▽ More
Modern shock-capturing schemes often suffer from numerical shock anomalies if the flow field contains strong shocks, which may limit their further application in hypersonic flow computations. In the current study, we devote our efforts to exploring the primary numerical characteristics and the underlying mechanism of shock instability for second-order finite-volume schemes. To this end, we, for the first time, develop the matrix stability analysis method for the finite-volume MUSCL approach. Such a linearized analysis method allows to investigate the shock instability problem of the finite-volume shock-capturing schemes in a quantitative and efficient manner. Results of the stability analysis demonstrate that the shock stability of second-order scheme is strongly related to the Riemann solver, Mach number, limiter function, numerical shock structure, and computational grid. Unique stability characteristics associated with these factors for second-order methods are revealed quantitatively with the established method. Source location of instability is also clarified by the matrix stability analysis method. Results show that the shock instability originates from the numerical shock structure. Such conclusions pave the way to better understand the shock instability problem and may shed new light on developing more reliable shock-capturing methods for compressible flows with high Mach number.
△ Less
Submitted 5 May, 2023;
originally announced May 2023.
-
$\ell$-adic local systems and Higgs bundles: the generic case
Authors:
Hongjie Yu
Abstract:
Let $X$ be a projective smooth geometrically connected curve defined over a finite field $\mathbb{F}_q$ of cardinality $q$. Let $S$ be a finite set of closed points of $X$. Let $\bar{X}$ and $\bar{S}$ be the base change of $X$, $S$ to an algebraic closure. We consider the set of $\ell$-adic ($\ell\nmid q$) local systems of rank $n$ over $\bar{X}-\bar{S}$ with prescribed tame regular semisimple and…
▽ More
Let $X$ be a projective smooth geometrically connected curve defined over a finite field $\mathbb{F}_q$ of cardinality $q$. Let $S$ be a finite set of closed points of $X$. Let $\bar{X}$ and $\bar{S}$ be the base change of $X$, $S$ to an algebraic closure. We consider the set of $\ell$-adic ($\ell\nmid q$) local systems of rank $n$ over $\bar{X}-\bar{S}$ with prescribed tame regular semisimple and generic ramifications in $\bar{S}$. The genericity ensures that such an $\ell$-adic local system is automatically irreducible. We show that the number of these $\ell$-adic local systems fixed by Frobenius endomorphism equals the number of stable logarithmic Higgs bundles of rank $n$ and degree $e$ coprime to $n$, with a fixed residue, up to a power of $q$. In the split case, this number is equal to the number of stable parabolic Higgs bundles (with full flag structures) fixed by $\mathbb{G}_m$-action with generic parabolic weights.
△ Less
Submitted 13 April, 2023;
originally announced April 2023.
-
Compact contact sets of sub-quadratic solutions to the thin obstacle problem
Authors:
Simon Eberle,
Hui Yu
Abstract:
We study global solutions to the thin obstacle problem with at most quadratic growth at infinity. We show that every ellipsoid can be realized as the contact set of such a solution. On the other hand, if such a solution has a compact contact set, we show that it must be an ellipsoid.
We study global solutions to the thin obstacle problem with at most quadratic growth at infinity. We show that every ellipsoid can be realized as the contact set of such a solution. On the other hand, if such a solution has a compact contact set, we show that it must be an ellipsoid.
△ Less
Submitted 31 March, 2024; v1 submitted 8 April, 2023;
originally announced April 2023.
-
Sign-twisted generating functions of the odd length in Weyl groups of type $D$
Authors:
Haihang Gu,
Houyi Yu
Abstract:
The odd length in Weyl groups is a new statistic analogous to the classical Coxeter length, and features combinatorial and parity conditions. We establish an explicit closed product formula for the sign-twisted generating functions of the odd length for parabolic quotients of Weyl groups of type $D$. As a consequence, we verify three conjectures of Brenti and Carnevale on evaluating closed forms f…
▽ More
The odd length in Weyl groups is a new statistic analogous to the classical Coxeter length, and features combinatorial and parity conditions. We establish an explicit closed product formula for the sign-twisted generating functions of the odd length for parabolic quotients of Weyl groups of type $D$. As a consequence, we verify three conjectures of Brenti and Carnevale on evaluating closed forms for these generating functions. We then give an equivalent condition for the sign-twisted generating functions to be expressible as products of cyclotomic polynomials, settling a conjecture of Stembridge.
△ Less
Submitted 20 February, 2024; v1 submitted 14 March, 2023;
originally announced March 2023.
-
Componentwise Linearity Under Square-Free Gröbner Degenerations
Authors:
Hongmiao Yu
Abstract:
Using the recent results on square-free Gröbner degenerations by Conca and Varbaro, we proved that if a homogeneous ideal $I$ of a polynomial ring is such that its initial ideal $\mathrm{in}_<(I)$ is square-free and $β_0(I) = β_0(\mathrm{in}_<(I))$, then $I$ is a componentwise linear ideal if and only if $\mathrm{in}_<(I)$ is a componentwise linear ideal. In particular, if furthermore one of $I$ a…
▽ More
Using the recent results on square-free Gröbner degenerations by Conca and Varbaro, we proved that if a homogeneous ideal $I$ of a polynomial ring is such that its initial ideal $\mathrm{in}_<(I)$ is square-free and $β_0(I) = β_0(\mathrm{in}_<(I))$, then $I$ is a componentwise linear ideal if and only if $\mathrm{in}_<(I)$ is a componentwise linear ideal. In particular, if furthermore one of $I$ and $\mathrm{in}_<(I)$ is componentwise linear, then their graded Betti numbers coincide.
△ Less
Submitted 30 March, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
Rank 2 $\ell$-adic local systems and Higgs bundles over a curve
Authors:
Hongjie Yu
Abstract:
Let $X$ be a smooth, projective, and geometrically connected curve defined over a finite field $\mathbb{F}_q$ and $S\subseteq X$ a subset of closed points. Let $\bar{X}$ and $\bar{S}$ be their base changes to an algebraic closure of $\mathbb{F}_q$. We study the number of $\ell$-adic local systems in rank $2$ over $\bar{X}-\bar{S}$ with prescribed tame local monodromies fixed by $k$-fold iterated a…
▽ More
Let $X$ be a smooth, projective, and geometrically connected curve defined over a finite field $\mathbb{F}_q$ and $S\subseteq X$ a subset of closed points. Let $\bar{X}$ and $\bar{S}$ be their base changes to an algebraic closure of $\mathbb{F}_q$. We study the number of $\ell$-adic local systems in rank $2$ over $\bar{X}-\bar{S}$ with prescribed tame local monodromies fixed by $k$-fold iterated action of Frobenius endomorphism for every $k\geq 1$. We confirm some conjectures of Deligne predicting that these numbers behave as if they were obtained from a Lefschetz fixed point formula. In fact, in all cases, our counting results are expressed in terms of the numbers of some Higgs bundles.
△ Less
Submitted 30 January, 2023;
originally announced January 2023.
-
Gaps between consecutive eigenvalues for compact metric graphs
Authors:
David Borthwick,
Evans M. Harrell II,
Haozhe Yu
Abstract:
On a compact metric graph, we consider the spectrum of the Laplacian defined with a mix of standard and Dirichlet vertex conditions. A Cheeger-type lower bound on the gap $λ_2 - λ_1$ is established, with a constant that depends only on the total length of the graph and minimum edge length. We also prove some improvements of known upper bounds for eigenvalue gaps and ratios for metric trees and ext…
▽ More
On a compact metric graph, we consider the spectrum of the Laplacian defined with a mix of standard and Dirichlet vertex conditions. A Cheeger-type lower bound on the gap $λ_2 - λ_1$ is established, with a constant that depends only on the total length of the graph and minimum edge length. We also prove some improvements of known upper bounds for eigenvalue gaps and ratios for metric trees and extensions to certain other types of graphs.
△ Less
Submitted 17 January, 2023;
originally announced January 2023.
-
Improved Laguerre Spectral Methods with Less Round-off Errors and Better Stability
Authors:
Shenghe Huang,
Haijun Yu
Abstract:
Laguerre polynomials are orthogonal polynomials defined on positive half line with respect to weight $e^{-x}$. They have wide applications in scientific and engineering computations. However, the exponential growth of Laguerre polynomials of high degree makes it hard to apply them to complicated systems that need to use large numbers of Laguerre bases. In this paper, we introduce modified three-te…
▽ More
Laguerre polynomials are orthogonal polynomials defined on positive half line with respect to weight $e^{-x}$. They have wide applications in scientific and engineering computations. However, the exponential growth of Laguerre polynomials of high degree makes it hard to apply them to complicated systems that need to use large numbers of Laguerre bases. In this paper, we introduce modified three-term recurrence formula to reduce the round-off error and to avoid overflow and underflow issues in generating generalized Laguerre polynomials and Laguerre functions. We apply the improved Laguerre methods to solve an elliptic equation defined on the half line. More than one thousand Laguerre bases are used in this application and meanwhile accuracy close to machine precision is achieved. The optimal scaling factor of Laguerre methods are studied and found to be independent of number of quadrature points in two cases that Laguerre methods have better convergence speeds than mapped Jacobi methods.
△ Less
Submitted 26 December, 2022;
originally announced December 2022.
-
Graphical solutions to one-phase free boundary problems
Authors:
Max Engelstein,
Xavier Fernández-Real,
Hui Yu
Abstract:
We study viscosity solutions to the classical one-phase problem and its thin counterpart. In low dimensions, we show that when the free boundary is the graph of a continuous function, the solution is the half-plane solution. This answers, in the salient dimensions, a one-phase free boundary analogue of Bernstein's problem for minimal surfaces.
As an application, we also classify monotone solutio…
▽ More
We study viscosity solutions to the classical one-phase problem and its thin counterpart. In low dimensions, we show that when the free boundary is the graph of a continuous function, the solution is the half-plane solution. This answers, in the salient dimensions, a one-phase free boundary analogue of Bernstein's problem for minimal surfaces.
As an application, we also classify monotone solutions of semilinear equations with a bump-type nonlinearity.
△ Less
Submitted 9 November, 2023; v1 submitted 17 December, 2022;
originally announced December 2022.
-
Global quasineutral Euler limit for the Vlasov-Poisson-Landau system with rarefaction waves
Authors:
Renjun Duan,
Dongcheng Yang,
Hongjun Yu
Abstract:
In the paper, we consider the Cauchy problem on the spatially one-dimensional Vlasov-Poisson-Landau system modelling the motion of ions under a generalized Boltzmann relation. Let the Knudsen number and the Debye length be given as $\varepsilon>0$ and $\varepsilon^{b}$ with $\frac{3}{5}\leq b\leq 1$, respectively. As $\varepsilon\to 0$ the formal Hilbert expansion gives the fluid limit to the quas…
▽ More
In the paper, we consider the Cauchy problem on the spatially one-dimensional Vlasov-Poisson-Landau system modelling the motion of ions under a generalized Boltzmann relation. Let the Knudsen number and the Debye length be given as $\varepsilon>0$ and $\varepsilon^{b}$ with $\frac{3}{5}\leq b\leq 1$, respectively. As $\varepsilon\to 0$ the formal Hilbert expansion gives the fluid limit to the quasineutral compressible Euler system. We start from the small-amplitude rarefaction wave of the Euler system that admits a smooth approximation with a parameter $δ\sim\varepsilon^{\frac{3}{5}-\frac{2}{5}a}$, where the wave strength is independent of $\varepsilon$ and we take $\frac{2}{3}\leq a\leq 1$ if $\frac{2}{3}\leq b\leq 1$ and $4-5b\leq a\leq 1$ if $\frac{3}{5}\leq b< \frac{2}{3}$. Under the scaling $(t,x)\to (\varepsilon^{-a}t,\varepsilon^{-a}x)$, for well-prepared initial data we construct the unique global classical solution to the Vlasov-Poisson-Landau system around the rarefaction wave in the vanishing limit $\varepsilon\to 0$ and also obtain the global-in-time convergence of solutions toward the rarefaction wave with rate $\varepsilon^{\frac{3}{5}-\frac{2}{5}a}|\ln\varepsilon|$ in the $L^{\infty}_xL^2_v$ norm. The best rate is $\varepsilon^{\frac{1}{3}}|\ln\varepsilon|$ with the choice of $a=\frac{2}{3}$ and $\frac{2}{3}\leq b\leq 1$. Note that the nontrivial electric potential in the solution connects two fixed distinct states at the far fields $x=\pm\infty$ for all $t\geq 0$ and tends asymptotically as $\varepsilon\to 0$ toward a profile determined by the macro density function under the quasineutral assumption. Our strategy is based on an intricate weighted energy method capturing the quartic dissipation to give uniform bounds of the nonlinear dynamics around rarefaction waves.
△ Less
Submitted 15 December, 2022;
originally announced December 2022.