-
Logarithmic regret in the ergodic Avellaneda-Stoikov market making model
Authors:
Jialun Cao,
David Šiška,
Lukasz Szpruch,
Tanut Treetanthiploet
Abstract:
We analyse the regret arising from learning the price sensitivity parameter $κ$ of liquidity takers in the ergodic version of the Avellaneda-Stoikov market making model. We show that a learning algorithm based on a regularised maximum-likelihood estimator for the parameter achieves the regret upper bound of order $\ln^2 T$ in expectation. To obtain the result we need two key ingredients. The first…
▽ More
We analyse the regret arising from learning the price sensitivity parameter $κ$ of liquidity takers in the ergodic version of the Avellaneda-Stoikov market making model. We show that a learning algorithm based on a regularised maximum-likelihood estimator for the parameter achieves the regret upper bound of order $\ln^2 T$ in expectation. To obtain the result we need two key ingredients. The first are tight upper bounds on the derivative of the ergodic constant in the Hamilton-Jacobi-Bellman (HJB) equation with respect to $κ$. The second is the learning rate of the maximum-likelihood estimator which is obtained from concentration inequalities for Bernoulli signals. Numerical experiment confirms the convergence and the robustness of the proposed algorithm.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
Towards Automatic Linearization via SMT Solving
Authors:
Jian Cao,
Liyong Lin,
Lele Li
Abstract:
Mathematical optimization is ubiquitous in modern applications. However, in practice, we often need to use nonlinear optimization models, for which the existing optimization tools such as Cplex or Gurobi may not be directly applicable and an (error-prone) manual transformation often has to be done. Thus, to address this issue, in this paper we investigate the problem of automatically verifying and…
▽ More
Mathematical optimization is ubiquitous in modern applications. However, in practice, we often need to use nonlinear optimization models, for which the existing optimization tools such as Cplex or Gurobi may not be directly applicable and an (error-prone) manual transformation often has to be done. Thus, to address this issue, in this paper we investigate the problem of automatically verifying and synthesizing reductions, the solution of which may allow an automatic linearization of nonlinear models. We show that the synthesis of reductions can be formulated as an $\exists^* \forall^*$ synthesis problem, which can be solved by an SMT solver via the counter-example guided inductive synthesis approach (CEGIS).
△ Less
Submitted 24 August, 2024;
originally announced August 2024.
-
The motivic fundamental group of a punctured elliptic curve and algebraic cycles
Authors:
Jin Cao,
Tomohide Terasoma
Abstract:
In this paper, we consider the motivic fundamental group of the punctured elliptic curves as a DG complex in the DG category of elliptic motives and describe its resolution via Schur complexes. During this process, we find the algebraic cycles analogous to the Bloch-Totaro cycles.
In this paper, we consider the motivic fundamental group of the punctured elliptic curves as a DG complex in the DG category of elliptic motives and describe its resolution via Schur complexes. During this process, we find the algebraic cycles analogous to the Bloch-Totaro cycles.
△ Less
Submitted 30 June, 2024;
originally announced July 2024.
-
A posteriori error estimates for the exponential midpoint method for linear and semilinear parabolic equations
Authors:
Xianfa Hu,
Wansheng Wang,
Mengli Mao,
Jiliang Cao
Abstract:
In this paper, the a posteriori error estimates of the exponential midpoint method for time discretization are studied for linear and semilinear parabolic equations. Using the exponential midpoint approximation defined by a continuous and piecewise linear interpolation of nodal values yields the suboptimal order estimates. Based on the property of the entire function, we introduce a continuous and…
▽ More
In this paper, the a posteriori error estimates of the exponential midpoint method for time discretization are studied for linear and semilinear parabolic equations. Using the exponential midpoint approximation defined by a continuous and piecewise linear interpolation of nodal values yields the suboptimal order estimates. Based on the property of the entire function, we introduce a continuous and piecewise quadratic time reconstruction of the exponential midpoint method to derive the optimal order estimates, and the error bounds are solely dependent on the discretization parameters, the data of the problem and the approximation of the entire function. Several numerical examples are implemented to illustrate the theoretical results.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Mixed Finite Element Method for Multi-layer Elastic Contact Systems
Authors:
Zhizhuo Zhang,
Mikaël Barboteu,
Xiaobing Nie,
Serge Dumont,
Mahmoud Abdel-Aty,
Jinde Cao
Abstract:
With the development of multi-layer elastic systems in the field of engineering mechanics, the corresponding variational inequality theory and algorithm design have received more attention and research. In this study, a class of equivalent saddle point problems with interlayer Tresca friction conditions and the mixed finite element method are proposed and analyzed. Then, the convergence of the num…
▽ More
With the development of multi-layer elastic systems in the field of engineering mechanics, the corresponding variational inequality theory and algorithm design have received more attention and research. In this study, a class of equivalent saddle point problems with interlayer Tresca friction conditions and the mixed finite element method are proposed and analyzed. Then, the convergence of the numerical solution of the mixed finite element method is theoretically proven, and the corresponding algebraic dual algorithm is given. Finally, through numerical experiments, the mixed finite element method is not only compared with the layer decomposition method, but also its convergence relationship with respect to the spatial discretization parameter $H$ is verified.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
A layer decomposition method for multi-layer elastic contact systems with interlayer Tresca friction
Authors:
Zhizhuo Zhang,
Xiaobing Nie,
Mikaël Barboteu,
Jinde Cao
Abstract:
With the increasing demand for the accuracy of numerical simulation of pavement mechanics, the variational inequality model and its induced finite element method which can simulate the interlayer contact state becomes a potential solution. In this paper, a layer decomposition algorithm for solving variational inequality models of multi-layer elastic contact systems with interlayer Tresca friction…
▽ More
With the increasing demand for the accuracy of numerical simulation of pavement mechanics, the variational inequality model and its induced finite element method which can simulate the interlayer contact state becomes a potential solution. In this paper, a layer decomposition algorithm for solving variational inequality models of multi-layer elastic contact systems with interlayer Tresca friction conditions is studied. Continuous and discrete versions of the algorithm and their convergence theorems have been proposed and proved successively. Then, the algebraic form of the executable optimization algorithm and the numerical experimental results verify the practicability of the variational inequality model and its algorithm in the pavement mechanics modeling.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
An Accelerated Gradient Method for Convex Smooth Simple Bilevel Optimization
Authors:
Jincheng Cao,
Ruichen Jiang,
Erfan Yazdandoost Hamedani,
Aryan Mokhtari
Abstract:
In this paper, we focus on simple bilevel optimization problems, where we minimize a convex smooth objective function over the optimal solution set of another convex smooth constrained optimization problem. We present a novel bilevel optimization method that locally approximates the solution set of the lower-level problem using a cutting plane approach and employs an accelerated gradient-based upd…
▽ More
In this paper, we focus on simple bilevel optimization problems, where we minimize a convex smooth objective function over the optimal solution set of another convex smooth constrained optimization problem. We present a novel bilevel optimization method that locally approximates the solution set of the lower-level problem using a cutting plane approach and employs an accelerated gradient-based update to reduce the upper-level objective function over the approximated solution set. We measure the performance of our method in terms of suboptimality and infeasibility errors and provide non-asymptotic convergence guarantees for both error criteria. Specifically, when the feasible set is compact, we show that our method requires at most $\mathcal{O}(\max\{1/\sqrt{ε_{f}}, 1/ε_g\})$ iterations to find a solution that is $ε_f$-suboptimal and $ε_g$-infeasible. Moreover, under the additional assumption that the lower-level objective satisfies the $r$-th Hölderian error bound, we show that our method achieves an iteration complexity of $\mathcal{O}(\max\{ε_{f}^{-\frac{2r-1}{2r}},ε_{g}^{-\frac{2r-1}{2r}}\})$, which matches the optimal complexity of single-level convex constrained optimization when $r=1$.
△ Less
Submitted 31 May, 2024; v1 submitted 12 February, 2024;
originally announced February 2024.
-
Optimal Insurance to Maximize Exponential Utility when Premium is Computed by a Convex Functional
Authors:
Jingyi Cao,
Dongchen Li,
Virginia R. Young,
Bin Zou
Abstract:
We find the optimal indemnity to maximize the expected utility of terminal wealth of a buyer of insurance whose preferences are modeled by an exponential utility. The insurance premium is computed by a convex functional. We obtain a necessary condition for the optimal indemnity; then, because the candidate optimal indemnity is given implicitly, we use that necessary condition to develop a numerica…
▽ More
We find the optimal indemnity to maximize the expected utility of terminal wealth of a buyer of insurance whose preferences are modeled by an exponential utility. The insurance premium is computed by a convex functional. We obtain a necessary condition for the optimal indemnity; then, because the candidate optimal indemnity is given implicitly, we use that necessary condition to develop a numerical algorithm to compute it. We prove that the numerical algorithm converges to a unique indemnity that, indeed, equals the optimal policy. We also illustrate our results with numerical examples.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
Scattering from time-modulated subwavelength resonators
Authors:
Habib Ammari,
Jinghao Cao,
Erik Orvehed Hiltunen,
Liora Rueff
Abstract:
We consider wave scattering from a system of highly contrasting resonators with time-modulated material parameters. In this setting, the wave equation reduces to a system of coupled Helmholtz equations that models the scattering problem. We consider the one-dimensional setting. In order to understand the energy of the system, we prove a novel higher-order discrete, capacitance matrix approximation…
▽ More
We consider wave scattering from a system of highly contrasting resonators with time-modulated material parameters. In this setting, the wave equation reduces to a system of coupled Helmholtz equations that models the scattering problem. We consider the one-dimensional setting. In order to understand the energy of the system, we prove a novel higher-order discrete, capacitance matrix approximation of the subwavelength resonant quasifrequencies. Further, we perform numerical experiments to support and illustrate our analytical results and show how periodically time-dependent material parameters affect the scattered wave field.
△ Less
Submitted 2 January, 2024;
originally announced January 2024.
-
Universal bifurcation scenarios in delay-differential equations with one delay
Authors:
Yu Wang,
Jinde Cao,
Jürgen Kurths,
Serhiy Yanchuk
Abstract:
We show that delay-differential equations (DDE) exhibit universal bifurcation scenarios, which are observed in large classes of DDEs with a single delay. Each such universality class has the same sequence of stabilizing or destabilizing Hopf bifurcations. These bifurcation sequences and universality classes can be explicitly described by using the asymptotic continuous spectrum for DDEs with large…
▽ More
We show that delay-differential equations (DDE) exhibit universal bifurcation scenarios, which are observed in large classes of DDEs with a single delay. Each such universality class has the same sequence of stabilizing or destabilizing Hopf bifurcations. These bifurcation sequences and universality classes can be explicitly described by using the asymptotic continuous spectrum for DDEs with large delays. Here, we mainly study linear DDEs, provide a general transversality result for the delay-induced bifurcations, and consider three most common universality classes. For each of them, we explicitly describe the sequence of stabilizing and destabilizing bifurcations. We also illustrate the implications for a nonlinear Stuart-Landau oscillator with time-delayed feedback.
△ Less
Submitted 29 December, 2023;
originally announced December 2023.
-
MG-Skip: Random Multi-Gossip Skipping Method for Nonsmooth Distributed Optimization
Authors:
Luyao Guo,
Luqing Wang,
Xinli Shi,
Jinde Cao
Abstract:
Distributed optimization methods with probabilistic local updates have recently gained attention for their provable ability to communication acceleration. Nevertheless, this capability is effective only when the loss function is smooth and the network is sufficiently well-connected. In this paper, we propose the first linear convergent method MG-Skip with probabilistic local updates for nonsmooth…
▽ More
Distributed optimization methods with probabilistic local updates have recently gained attention for their provable ability to communication acceleration. Nevertheless, this capability is effective only when the loss function is smooth and the network is sufficiently well-connected. In this paper, we propose the first linear convergent method MG-Skip with probabilistic local updates for nonsmooth distributed optimization. Without any extra condition for the network connectivity, MG-Skip allows for the multiple-round gossip communication to be skipped in most iterations, while its iteration complexity is $\mathcal{O}\left(κ\log \frac{1}ε\right)$ and communication complexity is only $\mathcal{O}\left(\sqrt{\fracκ{(1-ρ)}} \log \frac{1}ε\right)$, where $κ$ is the condition number of the loss function and $ρ$ reflects the connectivity of the network topology. To the best of our knowledge, MG-Skip achieves the best communication complexity when the loss function has the smooth (strongly convex)+nonsmooth (convex) composite form.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
Numerical approximation of discontinuous solutions of the semilinear wave equation
Authors:
Jiachuan Cao,
Buyang Li,
Yanping Lin,
Fangyan Yao
Abstract:
A fully discrete low-regularity integrator with high-frequency recovery techniques is constructed to approximate rough and possibly discontinuous solutions of the semilinear wave equation. The proposed method can capture the discontinuities of the solutions correctly without spurious oscillations and can approximate rough and discontinuous solutions with a higher convergence rate than pre-existing…
▽ More
A fully discrete low-regularity integrator with high-frequency recovery techniques is constructed to approximate rough and possibly discontinuous solutions of the semilinear wave equation. The proposed method can capture the discontinuities of the solutions correctly without spurious oscillations and can approximate rough and discontinuous solutions with a higher convergence rate than pre-existing methods. Rigorous analysis is presented for the convergence rates of the proposed method in approximating solutions such that $(u,\partial_{t}u)\in C([0,T];H^γ\times H^{γ-1})$ for $γ\in(0,1]$. For discontinuous solutions of bounded variation in one dimension (which allow jump discontinuities), the proposed method is proved to have almost first-order convergence under the step size condition $τ\sim N^{-1}$, where $τ$ and $N$ denote the time step size and the number of Fourier terms in the space discretization, respectively. Extensive numerical examples are presented in both one and two dimensions to illustrate the advantages of the proposed method in improving the accuracy in approximating rough and discontinuous solutions of the semilinear wave equation. The numerical results are consistent with the theoretical results and show the efficiency of the proposed method.
△ Less
Submitted 16 December, 2023;
originally announced December 2023.
-
The Non-Hermitian Skin Effect With Three-Dimensional Long-Range Coupling
Authors:
Habib Ammari,
Silvio Barandun,
Jinghao Cao,
Bryn Davies,
Erik Orvehed Hiltunen,
Ping Liu
Abstract:
We study the non-Hermitian skin effect in a three-dimensional system of finitely many subwavelength resonators with an imaginary gauge potential. We introduce a discrete approximation of the eigenmodes and eigenfrequencies of the system in terms of the eigenvectors and eigenvalues of the so-called gauge capacitance matrix $\mathcal{C}_N^γ$, which is a dense matrix due to long-range interactions in…
▽ More
We study the non-Hermitian skin effect in a three-dimensional system of finitely many subwavelength resonators with an imaginary gauge potential. We introduce a discrete approximation of the eigenmodes and eigenfrequencies of the system in terms of the eigenvectors and eigenvalues of the so-called gauge capacitance matrix $\mathcal{C}_N^γ$, which is a dense matrix due to long-range interactions in the system. Based on translational invariance of this matrix and the decay of its off-diagonal entries, we prove the condensation of the eigenmodes at one edge of the structure by showing the exponential decay of its pseudo-eigenvectors. In particular, we consider a range-k approximation to keep the long-range interaction to a certain extent, thus obtaining a k-banded gauge capacitance matrix $\mathcal{C}_{N,k}^γ$ . Using techniques for Toeplitz matrices and operators, we establish the exponential decay of the pseudo-eigenvectors of $\mathcal{C}_{N,k}^γ$ and demonstrate that they approximate those of the gauge capacitance matrix $\mathcal{C}_N^γ$ well. Our results are numerically verified. In particular, we show that long-range interactions affect only the first eigenmodes in the system. As a result, a tridiagonal approximation of the gauge capacitance matrix, similar to the nearest-neighbour approximation in quantum mechanics, provides a good approximation for the higher modes. Moreover, we also illustrate numerically the behaviour of the eigenmodes and the stability of the non-Hermitian skin effect with respect to disorder in a variety of three-dimensional structures.
△ Less
Submitted 17 November, 2023;
originally announced November 2023.
-
Statistical Parameterized Physics-Based Machine Learning Digital Twin Models for Laser Powder Bed Fusion Process
Authors:
Yangfan Li,
Satyajit Mojumder,
Ye Lu,
Abdullah Al Amin,
Jiachen Guo,
Xiaoyu Xie,
Wei Chen,
Gregory J. Wagner,
Jian Cao,
Wing Kam Liu
Abstract:
A digital twin (DT) is a virtual representation of physical process, products and/or systems that requires a high-fidelity computational model for continuous update through the integration of sensor data and user input. In the context of laser powder bed fusion (LPBF) additive manufacturing, a digital twin of the manufacturing process can offer predictions for the produced parts, diagnostics for m…
▽ More
A digital twin (DT) is a virtual representation of physical process, products and/or systems that requires a high-fidelity computational model for continuous update through the integration of sensor data and user input. In the context of laser powder bed fusion (LPBF) additive manufacturing, a digital twin of the manufacturing process can offer predictions for the produced parts, diagnostics for manufacturing defects, as well as control capabilities. This paper introduces a parameterized physics-based digital twin (PPB-DT) for the statistical predictions of LPBF metal additive manufacturing process. We accomplish this by creating a high-fidelity computational model that accurately represents the melt pool phenomena and subsequently calibrating and validating it through controlled experiments. In PPB-DT, a mechanistic reduced-order method-driven stochastic calibration process is introduced, which enables the statistical predictions of the melt pool geometries and the identification of defects such as lack-of-fusion porosity and surface roughness, specifically for diagnostic applications. Leveraging data derived from this physics-based model and experiments, we have trained a machine learning-based digital twin (PPB-ML-DT) model for predicting, monitoring, and controlling melt pool geometries. These proposed digital twin models can be employed for predictions, control, optimization, and quality assurance within the LPBF process, ultimately expediting product development and certification in LPBF-based metal additive manufacturing.
△ Less
Submitted 13 November, 2023;
originally announced November 2023.
-
Fixed-Time Gradient Flows for Solving Constrained Optimization: A Unified Approach
Authors:
Xinli Shi,
Xiangping Xu,
Guanghui Wen,
Jinde Cao
Abstract:
The accelerated method in solving optimization problems has always been an absorbing topic. Based on the fixed-time (FxT) stability of nonlinear dynamical systems, we provide a unified approach for designing FxT gradient flows (FxTGFs). First, a general class of nonlinear functions in designing FxTGFs is provided. A unified method for designing first-order FxTGFs is shown under PolyakL jasiewicz i…
▽ More
The accelerated method in solving optimization problems has always been an absorbing topic. Based on the fixed-time (FxT) stability of nonlinear dynamical systems, we provide a unified approach for designing FxT gradient flows (FxTGFs). First, a general class of nonlinear functions in designing FxTGFs is provided. A unified method for designing first-order FxTGFs is shown under PolyakL jasiewicz inequality assumption, a weaker condition than strong convexity. When there exist both bounded and vanishing disturbances in the gradient flow, a specific class of nonsmooth robust FxTGFs with disturbance rejection is presented. Under the strict convexity assumption, Newton-based FxTGFs is given and further extended to solve time-varying optimization. Besides, the proposed FxTGFs are further used for solving equation-constrained optimization. Moreover, an FxT proximal gradient flow with a wide range of parameters is provided for solving nonsmooth composite optimization. To show the effectiveness of various FxTGFs, the static regret analysis for several typical FxTGFs are also provided in detail. Finally, the proposed FxTGFs are applied to solve two network problems, i.e., the network consensus problem and solving a system linear equations, respectively, from the respective of optimization. Particularly, by choosing component-wisely sign-preserving functions, these problems can be solved in a distributed way, which extends the existing results. The accelerated convergence and robustness of the proposed FxTGFs are validated in several numerical examples stemming from practical applications.
△ Less
Submitted 12 November, 2023;
originally announced November 2023.
-
Optimal $L^\infty(L^2)$ and $L^1(L^2)$ a posteriori error estimates for the fully discrete approximations of time fractional parabolic differential equations
Authors:
Jiliang Cao,
Wansheng Wang,
Aiguo Xiao
Abstract:
We derive optimal order a posteriori error estimates in the $L^\infty(L^2)$ and $L^1(L^2)$-norms for the fully discrete approximations of time fractional parabolic differential equations. For the discretization in time, we use the $L1$ methods, while for the spatial discretization, we use standard conforming finite element methods. The linear and quadratic space-time reconstructions are introduced…
▽ More
We derive optimal order a posteriori error estimates in the $L^\infty(L^2)$ and $L^1(L^2)$-norms for the fully discrete approximations of time fractional parabolic differential equations. For the discretization in time, we use the $L1$ methods, while for the spatial discretization, we use standard conforming finite element methods. The linear and quadratic space-time reconstructions are introduced, which are generalizations of the elliptic space reconstruction. Then the related a posteriori error estimates for the linear and quadratic space-time reconstructions play key roles in deriving global and pointwise final error estimates. Numerical experiments verify and complement our theoretical results.
△ Less
Submitted 11 November, 2023;
originally announced November 2023.
-
PID-inspired Continuous-time Distributed Optimization
Authors:
Meng Tao,
Dongdong Yue,
Jinde Cao
Abstract:
This paper proposes two novel distributed continuous-time algorithms inspired by PID control to solve distributed optimization problems. The algorithms are referred to as first-order and second-order, respectively, depend on the intrinsic dynamics of the agents in the network. Sufficient conditions are derived so that both algorithms converge exponentially over undirected connected graphs. Finally…
▽ More
This paper proposes two novel distributed continuous-time algorithms inspired by PID control to solve distributed optimization problems. The algorithms are referred to as first-order and second-order, respectively, depend on the intrinsic dynamics of the agents in the network. Sufficient conditions are derived so that both algorithms converge exponentially over undirected connected graphs. Finally, numerical simulations illustrate the effectiveness and efficiency of the proposed algorithms.
△ Less
Submitted 8 November, 2023;
originally announced November 2023.
-
Multi/Single-stage structured zero-gradient-sum approach for prescribed-time optimization
Authors:
Shuaiyu Zhou,
Yiheng Wei,
Jinde Cao,
Yang Liu
Abstract:
Prescribed-time convergence mechanism has become a prominent research focus in the current field of optimization and control due to its ability to precisely control the target completion time. The recently arisen prescribed-time algorithms for distributed optimization, currently necessitate multi-stage structures to achieve global convergence. This paper introduces two modified zero-gradient-sum a…
▽ More
Prescribed-time convergence mechanism has become a prominent research focus in the current field of optimization and control due to its ability to precisely control the target completion time. The recently arisen prescribed-time algorithms for distributed optimization, currently necessitate multi-stage structures to achieve global convergence. This paper introduces two modified zero-gradient-sum algorithms, each based on a multi-stage and a single-stage structural frameworks established in this work. These algorithms are designed to achieve prescribed-time convergence and relax two common yet stringent conditions. This work also bridges the gap in current research on single-stage structured PTDO algorithm. The excellent convergence performance of the proposed algorithms is validated through a case study.
△ Less
Submitted 29 October, 2023;
originally announced October 2023.
-
Revisiting Decentralized ProxSkip: Achieving Linear Speedup
Authors:
Luyao Guo,
Sulaiman A. Alghunaim,
Kun Yuan,
Laurent Condat,
Jinde Cao
Abstract:
The ProxSkip algorithm for decentralized and federated learning is gaining increasing attention due to its proven benefits in accelerating communication complexity while maintaining robustness against data heterogeneity. However, existing analyses of ProxSkip are limited to the strongly convex setting and do not achieve linear speedup, where convergence performance increases linearly with respect…
▽ More
The ProxSkip algorithm for decentralized and federated learning is gaining increasing attention due to its proven benefits in accelerating communication complexity while maintaining robustness against data heterogeneity. However, existing analyses of ProxSkip are limited to the strongly convex setting and do not achieve linear speedup, where convergence performance increases linearly with respect to the number of nodes. So far, questions remain open about how ProxSkip behaves in the non-convex setting and whether linear speedup is achievable.
In this paper, we revisit decentralized ProxSkip and address both questions. We demonstrate that the leading communication complexity of ProxSkip is $\mathcal{O}\left(\frac{pσ^2}{nε^2}\right)$ for non-convex and convex settings, and $\mathcal{O}\left(\frac{pσ^2}{nε}\right)$ for the strongly convex setting, where $n$ represents the number of nodes, $p$ denotes the probability of communication, $σ^2$ signifies the level of stochastic noise, and $ε$ denotes the desired accuracy level. This result illustrates that ProxSkip achieves linear speedup and can asymptotically reduce communication overhead proportional to the probability of communication. Additionally, for the strongly convex setting, we further prove that ProxSkip can achieve linear speedup with network-independent stepsizes.
△ Less
Submitted 19 April, 2024; v1 submitted 11 October, 2023;
originally announced October 2023.
-
Almost sure dimensional properties for the spectrum and the density of states of Sturmian Hamiltonians
Authors:
Jie Cao,
Yanhui Qu
Abstract:
In this paper, we find a full Lebesgue measure set of frequencies $\check \II\subset [0,1]\setminus \Q$ such that for any $(α,λ)\in \check \II\times [24,\infty)$, the Hausdorff and box dimensions of the spectrum of the Sturmian Hamiltonian $H_{α,λ,θ}$ coincide and are independent of $α$. Denote the common value by $D(λ)$, we show that $D(λ)$ satisfies a Bowen type formula, and is locally Lipschitz…
▽ More
In this paper, we find a full Lebesgue measure set of frequencies $\check \II\subset [0,1]\setminus \Q$ such that for any $(α,λ)\in \check \II\times [24,\infty)$, the Hausdorff and box dimensions of the spectrum of the Sturmian Hamiltonian $H_{α,λ,θ}$ coincide and are independent of $α$. Denote the common value by $D(λ)$, we show that $D(λ)$ satisfies a Bowen type formula, and is locally Lipschitz. We obtain the exact asymptotic behavior of $D(λ)$ as $λ$ tends to $ \infty.$ This considerably improves the result of Damanik and Gorodetski (Comm. Math. Phys. 337, 2015). We also show that for any $(α,λ)\in \check \II\times [24,\infty)$, the density of states measure of $H_{α,λ,θ}$ is exact-dimensional; its Hausdorff and packing dimensions coincide and are independent of $α$. Denote the common value by $d(λ)$, we show that $d(λ)$ satisfies a Young type formula, and is Lipschitz. We obtain the exact asymptotic behavior of $d(λ)$ as $λ$ tends to $ \infty.$ During the course of study, we also answer several questions in the same paper of Damanik and Gorodetski.
△ Less
Submitted 11 October, 2023;
originally announced October 2023.
-
$pκ$-Curves: Interpolatory curves with curvature approximating a parabola
Authors:
Zhihao Wang,
Juan Cao,
Tuan Guan,
Zhonggui Chen,
Yongjie Jessica Zhang
Abstract:
This paper introduces a novel class of fair and interpolatory curves called $pκ$-curves. These curves are comprised of smoothly stitched Bézier curve segments, where the curvature distribution of each segment is made to closely resemble a parabola, resulting in an aesthetically pleasing shape. Moreover, each segment passes through an interpolated point at a parameter where the parabola has an extr…
▽ More
This paper introduces a novel class of fair and interpolatory curves called $pκ$-curves. These curves are comprised of smoothly stitched Bézier curve segments, where the curvature distribution of each segment is made to closely resemble a parabola, resulting in an aesthetically pleasing shape. Moreover, each segment passes through an interpolated point at a parameter where the parabola has an extremum, encouraging the alignment of interpolated points with curvature extrema. To achieve these properties, we tailor an energy function that guides the optimization process to obtain the desired curve characteristics. Additionally, we develop an efficient algorithm and an initialization method, enabling interactive modeling of the $pκ$-curves without the need for global optimization. We provide various examples and comparisons with existing state-of-the-art methods to demonstrate the curve modeling capabilities and visually pleasing appearance of $pκ$-curves.
△ Less
Submitted 11 October, 2023;
originally announced October 2023.
-
High-order BDF convolution quadrature for fractional evolution equations with hyper-singular source term
Authors:
Jiankang Shi,
Minghua Chen,
Jianxiong Cao
Abstract:
Anomalous diffusion in the presence or absence of an external force field is often modelled in terms of the fractional evolution equations, which can involve the hyper-singular source term.
For this case, conventional time stepping methods may exhibit a severe order reduction. Although a second-order numerical algorithm is provided for the subdiffusion model with a simple hyper-singular source t…
▽ More
Anomalous diffusion in the presence or absence of an external force field is often modelled in terms of the fractional evolution equations, which can involve the hyper-singular source term.
For this case, conventional time stepping methods may exhibit a severe order reduction. Although a second-order numerical algorithm is provided for the subdiffusion model with a simple hyper-singular source term $t^μ$, $-2<μ<-1$ in [arXiv:2207.08447], the convergence analysis remain to be proved.
To fill in these gaps, we present a simple and robust smoothing method for the hyper-singular source term, where the Hadamard finite-part integral is introduced. This method is based on the smoothing/ID$m$-BDF$k$ method proposed by the authors [Shi and Chen, SIAM J. Numer. Anal., to appear] for subdiffusion equation with a weakly singular source term.
We prove that the $k$th-order convergence rate can be restored for the diffusion-wave case $γ\in (1,2)$ and sketch the proof for the subdiffusion case $γ\in (0,1)$, even if the source term is hyper-singular and the initial data is not compatible.
Numerical experiments are provided to confirm the theoretical results.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
Projection-Free Methods for Stochastic Simple Bilevel Optimization with Convex Lower-level Problem
Authors:
Jincheng Cao,
Ruichen Jiang,
Nazanin Abolfazli,
Erfan Yazdandoost Hamedani,
Aryan Mokhtari
Abstract:
In this paper, we study a class of stochastic bilevel optimization problems, also known as stochastic simple bilevel optimization, where we minimize a smooth stochastic objective function over the optimal solution set of another stochastic convex optimization problem. We introduce novel stochastic bilevel optimization methods that locally approximate the solution set of the lower-level problem via…
▽ More
In this paper, we study a class of stochastic bilevel optimization problems, also known as stochastic simple bilevel optimization, where we minimize a smooth stochastic objective function over the optimal solution set of another stochastic convex optimization problem. We introduce novel stochastic bilevel optimization methods that locally approximate the solution set of the lower-level problem via a stochastic cutting plane, and then run a conditional gradient update with variance reduction techniques to control the error induced by using stochastic gradients. For the case that the upper-level function is convex, our method requires $\tilde{\mathcal{O}}(\max\{1/ε_f^{2},1/ε_g^{2}\}) $ stochastic oracle queries to obtain a solution that is $ε_f$-optimal for the upper-level and $ε_g$-optimal for the lower-level. This guarantee improves the previous best-known complexity of $\mathcal{O}(\max\{1/ε_f^{4},1/ε_g^{4}\})$. Moreover, for the case that the upper-level function is non-convex, our method requires at most $\tilde{\mathcal{O}}(\max\{1/ε_f^{3},1/ε_g^{3}\}) $ stochastic oracle queries to find an $(ε_f, ε_g)$-stationary point. In the finite-sum setting, we show that the number of stochastic oracle calls required by our method are $\tilde{\mathcal{O}}(\sqrt{n}/ε)$ and $\tilde{\mathcal{O}}(\sqrt{n}/ε^{2})$ for the convex and non-convex settings, respectively, where $ε=\min \{ε_f,ε_g\}$.
△ Less
Submitted 14 August, 2023;
originally announced August 2023.
-
Resilient Controller Synthesis Against DoS Attacks for Vehicular Platooning in Spatial Domain
Authors:
Jian Gong,
Carlos Murguia,
Anggera Bayuwindra,
Jinde Cao
Abstract:
This paper proposes a vehicular platoon control approach under Denial-of-Service (DoS) attacks and external disturbances. DoS attacks increase the service time on the communication network and cause additional transmission delays, which consequently increase the risk of rear-end collisions of vehicles in the platoon. To counter DoS attacks, we propose a resilient control scheme that exploits polyt…
▽ More
This paper proposes a vehicular platoon control approach under Denial-of-Service (DoS) attacks and external disturbances. DoS attacks increase the service time on the communication network and cause additional transmission delays, which consequently increase the risk of rear-end collisions of vehicles in the platoon. To counter DoS attacks, we propose a resilient control scheme that exploits polytopic overapproximations of the closed-loop dynamics under DoS attacks. This scheme allows synthesizing robust controllers that guarantee tracking of both the desired spacing policy and spatially varying reference velocity for all space-varying DoS attacks satisfying a hard upper bound on the attack duration. In addition, L2 string stability conditions are derived to ensure that external perturbations do not grow as they propagate through the platoon, thus ensuring the string stability. Numerical simulations illustrate the effectiveness of the proposed control method.
△ Less
Submitted 28 July, 2023;
originally announced July 2023.
-
Geometry of $K$-trivial Moishezon manifolds : decomposition theorem and holomorphic geometric structures
Authors:
Indranil Biswas,
Junyan Cao,
Sorin Dumitrescu,
Henri Guenancia
Abstract:
Let $X$ be a compact complex manifold such that its canonical bundle $K_X$ is numerically trivial. Assume additionally that $X$ is Moishezon or $X$ is Fujiki with dimension at most four. Using the MMP and classical results in foliation theory, we prove a Beauville-Bogomolov type decomposition theorem for $X$. We deduce that holomorphic geometric structures of affine type on $X$ are in fact locally…
▽ More
Let $X$ be a compact complex manifold such that its canonical bundle $K_X$ is numerically trivial. Assume additionally that $X$ is Moishezon or $X$ is Fujiki with dimension at most four. Using the MMP and classical results in foliation theory, we prove a Beauville-Bogomolov type decomposition theorem for $X$. We deduce that holomorphic geometric structures of affine type on $X$ are in fact locally homogeneous away from an analytic subset of complex codimension at least two, and that they cannot be rigid unless $X$ is an étale quotient of a compact complex torus. Moreover, we establish a characterization of torus quotients using the vanishing of the first two Chern classes which is valid for any compact complex $n$-folds of algebraic dimension at least $n-1$. Finally, we show that a compact complex manifold with trivial canonical bundle bearing a rigid geometric structure must have infinite fundamental group if either $X$ is Fujiki, $X$ is a threefold, or $X$ is of algebraic dimension at most one.
△ Less
Submitted 10 July, 2023; v1 submitted 29 June, 2023;
originally announced June 2023.
-
Mathematical foundations of the non-Hermitian skin effect
Authors:
Habib Ammari,
Silvio Barandun,
Jinghao Cao,
Bryn Davies,
Erik Orvehed Hiltunen
Abstract:
We study the skin effect in a one-dimensional system of finitely many subwavelength resonators with a non-Hermitian imaginary gauge potential. Using Toeplitz matrix theory, we prove the condensation of bulk eigenmodes at one of the edges of the system. By introducing a generalised (complex) Brillouin zone, we can compute spectral bands of the associated infinitely periodic structure and prove that…
▽ More
We study the skin effect in a one-dimensional system of finitely many subwavelength resonators with a non-Hermitian imaginary gauge potential. Using Toeplitz matrix theory, we prove the condensation of bulk eigenmodes at one of the edges of the system. By introducing a generalised (complex) Brillouin zone, we can compute spectral bands of the associated infinitely periodic structure and prove that this is the limit of the spectra of the finite structures with arbitrarily large size. Finally, we contrast the non-Hermitian systems with imaginary gauge potentials considered here with systems where the non-Hermiticity arises due to complex material parameters, showing that the two systems are fundamentally distinct.
△ Less
Submitted 27 June, 2023;
originally announced June 2023.
-
Functional Adaptive Double-Sparsity Estimator for Functional Linear Regression Model with Multiple Functional Covariates
Authors:
Cheng Cao,
Jiguo Cao,
Hailiang Wang,
Kwok-Leung Tsui,
Xinyue Li
Abstract:
Sensor devices have been increasingly used in engineering and health studies recently, and the captured multi-dimensional activity and vital sign signals can be studied in association with health outcomes to inform public health. The common approach is the scalar-on-function regression model, in which health outcomes are the scalar responses while high-dimensional sensor signals are the functional…
▽ More
Sensor devices have been increasingly used in engineering and health studies recently, and the captured multi-dimensional activity and vital sign signals can be studied in association with health outcomes to inform public health. The common approach is the scalar-on-function regression model, in which health outcomes are the scalar responses while high-dimensional sensor signals are the functional covariates, but how to effectively interpret results becomes difficult. In this study, we propose a new Functional Adaptive Double-Sparsity (FadDoS) estimator based on functional regularization of sparse group lasso with multiple functional predictors, which can achieve global sparsity via functional variable selection and local sparsity via zero-subinterval identification within coefficient functions. We prove that the FadDoS estimator converges at a bounded rate and satisfies the oracle property under mild conditions. Extensive simulation studies confirm the theoretical properties and exhibit excellent performances compared to existing approaches. Application to a Kinect sensor study that utilized an advanced motion sensing device tracking human multiple joint movements and conducted among community-dwelling elderly demonstrates how the FadDoS estimator can effectively characterize the detailed association between joint movements and physical health assessments. The proposed method is not only effective in Kinect sensor analysis but also applicable to broader fields, where multi-dimensional sensor signals are collected simultaneously, to expand the use of sensor devices in health studies and facilitate sensor data analysis.
△ Less
Submitted 27 June, 2024; v1 submitted 17 May, 2023;
originally announced May 2023.
-
Kirchhoff-Love shell representation and analysis using triangle configuration B-splines
Authors:
Zhihao Wang,
Juan Cao,
Xiaodong Wei,
Zhonggui Chen,
Hugo Casquero,
Yongjie Jessica Zhang
Abstract:
This paper presents the application of triangle configuration B-splines (TCB-splines) for representing and analyzing the Kirchhoff-Love shell in the context of isogeometric analysis (IGA). The Kirchhoff-Love shell formulation requires global $C^1$-continuous basis functions. The nonuniform rational B-spline (NURBS)-based IGA has been extensively used for developing Kirchhoff-Love shell elements. H…
▽ More
This paper presents the application of triangle configuration B-splines (TCB-splines) for representing and analyzing the Kirchhoff-Love shell in the context of isogeometric analysis (IGA). The Kirchhoff-Love shell formulation requires global $C^1$-continuous basis functions. The nonuniform rational B-spline (NURBS)-based IGA has been extensively used for developing Kirchhoff-Love shell elements. However, shells with complex geometries inevitably need multiple patches and trimming techniques, where stitching patches with high continuity is a challenge. On the other hand, due to their unstructured nature, TCB-splines can accommodate general polygonal domains, have local refinement, and are flexible to model complex geometries with $C^1$ continuity, which naturally fit into the Kirchhoff-Love shell formulation with complex geometries. Therefore, we propose to use TCB-splines as basis functions for geometric representation and solution approximation. We apply our method to both linear and nonlinear benchmark shell problems, where the accuracy and robustness are validated. The applicability of the proposed approach to shell analysis is further exemplified by performing geometrically nonlinear Kirchhoff-Love shell simulations of a pipe junction and a front bumper represented by a single patch of TCB-splines.
△ Less
Submitted 20 April, 2023;
originally announced April 2023.
-
Sensor Fault Detection and Isolation in Autonomous Nonlinear Systems Using Neural Network-Based Observers
Authors:
John Cao,
Muhammad Umar B. Niazi,
Matthieu Barreau,
Karl Henrik Johansson
Abstract:
This paper presents a novel observer-based approach to detect and isolate faulty sensors in nonlinear systems. The proposed sensor fault detection and isolation (s-FDI) method applies to a general class of nonlinear systems. Our focus is on s-FDI for two types of faults: complete failure and sensor degradation. The key aspect of this approach lies in the utilization of a neural network-based Kazan…
▽ More
This paper presents a novel observer-based approach to detect and isolate faulty sensors in nonlinear systems. The proposed sensor fault detection and isolation (s-FDI) method applies to a general class of nonlinear systems. Our focus is on s-FDI for two types of faults: complete failure and sensor degradation. The key aspect of this approach lies in the utilization of a neural network-based Kazantzis-Kravaris/Luenberger (KKL) observer. The neural network is trained to learn the dynamics of the observer, enabling accurate output predictions of the system. Sensor faults are detected by comparing the actual output measurements with the predicted values. If the difference surpasses a theoretical threshold, a sensor fault is detected. To identify and isolate which sensor is faulty, we compare the numerical difference of each sensor meassurement with an empirically derived threshold. We derive both theoretical and empirical thresholds for detection and isolation, respectively. Notably, the proposed approach is robust to measurement noise and system uncertainties. Its effectiveness is demonstrated through numerical simulations of sensor faults in a network of Kuramoto oscillators.
△ Less
Submitted 22 November, 2023; v1 submitted 18 April, 2023;
originally announced April 2023.
-
$\partial\bar\partial$-lemmas and a conjecture of O. Fujino
Authors:
Junyan Cao,
Mihai Păun
Abstract:
It contains the proof of a very general $\partial\bar\partial$-lemma, together with a decomposition theorem for currents with values in a (singular) Hermitian line bundle. As a corollary, we establish the Kähler version on an injectivity theorem due to O. Fujino in the projective case.
It contains the proof of a very general $\partial\bar\partial$-lemma, together with a decomposition theorem for currents with values in a (singular) Hermitian line bundle. As a corollary, we establish the Kähler version on an injectivity theorem due to O. Fujino in the projective case.
△ Less
Submitted 28 March, 2023;
originally announced March 2023.
-
Hermite--Einstein metrics in singular settings
Authors:
Junyan Cao,
Patrick Graf,
Philipp Naumann,
Mihai Paun,
Thomas Peternell,
Xiaojun Wu
Abstract:
In this article we pursue the following main goals. In the first place, we establish the existence of "estimable" Hermite--Einstein metrics for stable reflexive coherent sheaves on compact normal Kähler spaces. If moreover the background variety has klt singularities, we obtain a much more precise result.
In this article we pursue the following main goals. In the first place, we establish the existence of "estimable" Hermite--Einstein metrics for stable reflexive coherent sheaves on compact normal Kähler spaces. If moreover the background variety has klt singularities, we obtain a much more precise result.
△ Less
Submitted 21 January, 2024; v1 submitted 15 March, 2023;
originally announced March 2023.
-
Direct images of pseudoeffective cotangent bundles
Authors:
Junyan Cao,
Andreas Höring
Abstract:
Let A be an elliptic curve, and let $V_A$ be the Serre vector bundle on A. A famous example of Demailly-Peternell-Schneider shows that the tautological class of $V_A$ contains a unique closed positive current. In this survey we start by generalising this statement to arbitrary compact Kähler manifolds. We then give an application to abelian fibrations $X \rightarrow Y$ where the total space X has…
▽ More
Let A be an elliptic curve, and let $V_A$ be the Serre vector bundle on A. A famous example of Demailly-Peternell-Schneider shows that the tautological class of $V_A$ contains a unique closed positive current. In this survey we start by generalising this statement to arbitrary compact Kähler manifolds. We then give an application to abelian fibrations $X \rightarrow Y$ where the total space X has pseudoeffective cotangent bundle and raise some questions about nonvanishing properties of these bundles.
△ Less
Submitted 24 February, 2023;
originally announced February 2023.
-
A novel class of fractional adams method for solving uncertain fractional differential equation
Authors:
Chenlei Tian,
Jing Cao,
Yifu Song,
Ting Jin
Abstract:
Uncertain fractional differential equation (UFDE) is a kind of differential equation about uncertain process. As an significant mathematical tool to describe the evolution process of dynamic system, UFDE is better than the ordinary differential equation with integer derivatives because of its hereditability and memorability characteristics. However, in most instances, the precise analytical soluti…
▽ More
Uncertain fractional differential equation (UFDE) is a kind of differential equation about uncertain process. As an significant mathematical tool to describe the evolution process of dynamic system, UFDE is better than the ordinary differential equation with integer derivatives because of its hereditability and memorability characteristics. However, in most instances, the precise analytical solutions of UFDE is difficult to obtain due to the complex form of the UFDE itself. Up to now, there is not plenty of researches about the numerical method of UFDE, as for the existing numerical algorithms, their accuracy is also not high. In this research, derive from the interval weighting method, a class of fractional adams method is innovatively proposed to solve UFDE. Meanwhile, such fractional adams method extends the traditional predictor-corrector method to higher order cases. The stability and truncation error limit of the improved algorithm are analyzed and deduced. As the application, several numerical simulations (including $α$-path, extreme value and the first hitting time of the UFDE) are provided to manifest the higher accuracy and efficiency of the proposed numerical method.
△ Less
Submitted 24 February, 2023;
originally announced February 2023.
-
On some analytic properties of nabla tempered fractional calculus
Authors:
Yiheng Wei,
Linlin Zhao,
Xuan Zhao,
Jinde Cao
Abstract:
Despite many applications regarding fractional calculus have been reported in literature, it is still unknown how to model some practical process. One major challenge in solving such a problem is that, the nonlocal property is needed while the infinite memory is undesired. Under this context, a new kind nabla fractional calculus accompanied by a tempered function is formulated. However, many prope…
▽ More
Despite many applications regarding fractional calculus have been reported in literature, it is still unknown how to model some practical process. One major challenge in solving such a problem is that, the nonlocal property is needed while the infinite memory is undesired. Under this context, a new kind nabla fractional calculus accompanied by a tempered function is formulated. However, many properties of such fractional calculus needed to be discovered. From this, this paper gives particular emphasis to the topic. Some remarkable properties like the equivalence relation, the nabla Taylor formula, and the nabla Laplace transform for such nabla fractional calculus are developed and analyzed. It is believed that this work greatly enriches the mathematical theory of nabla tempered fractional calculus and provides high value and huge potential for further applications.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
Decentralized Inexact Proximal Gradient Method With Network-Independent Stepsizes for Convex Composite Optimization
Authors:
Luyao Guo,
Xinli Shi,
Jinde Cao,
Zihao Wang
Abstract:
This paper proposes a novel CTA (Combine-Then-Adapt)-based decentralized algorithm for solving convex composite optimization problems over undirected and connected networks. The local loss function in these problems contains both smooth and nonsmooth terms. The proposed algorithm uses uncoordinated network-independent constant stepsizes and only needs to approximately solve a sequence of proximal…
▽ More
This paper proposes a novel CTA (Combine-Then-Adapt)-based decentralized algorithm for solving convex composite optimization problems over undirected and connected networks. The local loss function in these problems contains both smooth and nonsmooth terms. The proposed algorithm uses uncoordinated network-independent constant stepsizes and only needs to approximately solve a sequence of proximal mappings, which is advantageous for solving decentralized composite optimization problems where the proximal mappings of the nonsmooth loss functions may not have analytical solutions. For the general convex case, we prove an O(1/k) convergence rate of the proposed algorithm, which can be improved to o(1/k) if the proximal mappings are solved exactly. Furthermore, with metric subregularity, we establish a linear convergence rate for the proposed algorithm. Numerical experiments demonstrate the efficiency of the algorithm.
△ Less
Submitted 4 March, 2023; v1 submitted 6 February, 2023;
originally announced February 2023.
-
Transmission properties of time-dependent one-dimensional metamaterials
Authors:
Habib Ammari,
Jinghao Cao,
Erik Orvehed Hiltunen,
Liora Rueff
Abstract:
We solve the wave equation with periodically time-modulated material parameters in a one-dimensional high-contrast resonator structure in the subwavelength regime exactly, for which we compute the subwavelength quasifrequencies numerically using Muller's method. We prove a formula in the form of an ODE using a capacitance matrix approximation. Comparison of the exact results with the approximation…
▽ More
We solve the wave equation with periodically time-modulated material parameters in a one-dimensional high-contrast resonator structure in the subwavelength regime exactly, for which we compute the subwavelength quasifrequencies numerically using Muller's method. We prove a formula in the form of an ODE using a capacitance matrix approximation. Comparison of the exact results with the approximations reveals that the method of capacitance matrix approximation is accurate and significantly more efficient. We prove various transmission properties in the aforementioned structure and illustrate them with numerical simulations. In particular, we investigate the effect of time-modulated material parameters on the formation of degenerate points, band gaps and k-gaps.
△ Less
Submitted 31 May, 2023; v1 submitted 24 January, 2023;
originally announced January 2023.
-
Edge modes in subwavelength resonators in one dimension
Authors:
Habib Ammari,
Silvio Barandun,
Jinghao Cao,
Florian Feppon
Abstract:
We present the mathematical theory of one-dimensional infinitely periodic chains of subwavelength resonators. We analyse both Hermitian and non-Hermitian systems. Subwavelength resonances and associated modes can be accurately predicted by a finite dimensional eigenvalue problem involving a capacitance matrix. We are able to compute the Hermitian and non-Hermitian Zak phases, showing that the form…
▽ More
We present the mathematical theory of one-dimensional infinitely periodic chains of subwavelength resonators. We analyse both Hermitian and non-Hermitian systems. Subwavelength resonances and associated modes can be accurately predicted by a finite dimensional eigenvalue problem involving a capacitance matrix. We are able to compute the Hermitian and non-Hermitian Zak phases, showing that the former is quantised and the latter is not. Furthermore, we show the existence of localised edge modes arising from defects in the periodicity in both the Hermitian and non-Hermitian cases. In the non-Hermitian case, we provide a complete characterisation of the edge modes.
△ Less
Submitted 20 April, 2023; v1 submitted 17 January, 2023;
originally announced January 2023.
-
Some useful inequalities for nabla tempered fractional calculus
Authors:
Yiheng Wei,
Linlin Zhao,
Kai Cao,
Jinde Cao
Abstract:
This paper gives particular emphasis to the nabla tempered fractional calculus, which involves the multiplication of the rising function kernel with tempered functions, and provides a more flexible alternative with considerable promise for practical applications. Some remarkable inequalities for such nabla fractional calculus are developed and analyzed, which greatly enrich the mathematical theory…
▽ More
This paper gives particular emphasis to the nabla tempered fractional calculus, which involves the multiplication of the rising function kernel with tempered functions, and provides a more flexible alternative with considerable promise for practical applications. Some remarkable inequalities for such nabla fractional calculus are developed and analyzed, which greatly enrich the mathematical theory of nabla tempered fractional calculus. Numerical results confirm the validity of the developed properties once again, which also reveals that the introduction of tempered function provides high value and huge potential.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
BALPA: A Balanced Primal-Dual Algorithm for Nonsmooth Optimization with Application to Distributed Optimization
Authors:
Luyao Guo,
Jinde Cao,
Xinli Shi,
Shaofu Yang
Abstract:
In this paper, we propose a novel primal-dual proximal splitting algorithm (PD-PSA), named BALPA, for the composite optimization problem with equality constraints, where the loss function consists of a smooth term and a nonsmooth term composed with a linear mapping. In BALPA, the dual update is designed as a proximal point for a time-varying quadratic function, which balances the implementation of…
▽ More
In this paper, we propose a novel primal-dual proximal splitting algorithm (PD-PSA), named BALPA, for the composite optimization problem with equality constraints, where the loss function consists of a smooth term and a nonsmooth term composed with a linear mapping. In BALPA, the dual update is designed as a proximal point for a time-varying quadratic function, which balances the implementation of primal and dual update and retains the proximity-induced feature of classic PD-PSAs. In addition, by this balance, BALPA eliminates the inefficiency of classic PD-PSAs for composite optimization problems in which the Euclidean norm of the linear mapping or the equality constraint mapping is large. Therefore, BALPA not only inherits the advantages of simple structure and easy implementation of classic PD-PSAs but also ensures a fast convergence when these norms are large. Moreover, we propose a stochastic version of BALPA (S-BALPA) and apply the developed BALPA to distributed optimization to devise a new distributed optimization algorithm. Furthermore, a comprehensive convergence analysis for BALPA and S-BALPA is conducted, respectively. Finally, numerical experiments demonstrate the efficiency of the proposed algorithms.
△ Less
Submitted 6 December, 2022;
originally announced December 2022.
-
Frobenius integrability of certain $p$-forms on singular spaces
Authors:
Junyan Cao,
Andreas Höring
Abstract:
Demailly proved that on a smooth compact Kähler manifold the distribution defined by a holomorphic $p$-form with values in an anti-pseudoeffective line bundle is always integrable. We generalise his result to compact Kähler spaces with klt singularities.
Demailly proved that on a smooth compact Kähler manifold the distribution defined by a holomorphic $p$-form with values in an anti-pseudoeffective line bundle is always integrable. We generalise his result to compact Kähler spaces with klt singularities.
△ Less
Submitted 6 October, 2022;
originally announced October 2022.
-
Learning-based Design of Luenberger Observers for Autonomous Nonlinear Systems
Authors:
Muhammad Umar B. Niazi,
John Cao,
Xudong Sun,
Amritam Das,
Karl Henrik Johansson
Abstract:
Designing Luenberger observers for nonlinear systems involves the challenging task of transforming the state to an alternate coordinate system, possibly of higher dimensions, where the system is asymptotically stable and linear up to output injection. The observer then estimates the system's state in the original coordinates by inverting the transformation map. However, finding a suitable injectiv…
▽ More
Designing Luenberger observers for nonlinear systems involves the challenging task of transforming the state to an alternate coordinate system, possibly of higher dimensions, where the system is asymptotically stable and linear up to output injection. The observer then estimates the system's state in the original coordinates by inverting the transformation map. However, finding a suitable injective transformation whose inverse can be derived remains a primary challenge for general nonlinear systems. We propose a novel approach that uses supervised physics-informed neural networks to approximate both the transformation and its inverse. Our method exhibits superior generalization capabilities to contemporary methods and demonstrates robustness to both neural network's approximation errors and system uncertainties.
△ Less
Submitted 5 April, 2023; v1 submitted 4 October, 2022;
originally announced October 2022.
-
DISA: A Dual Inexact Splitting Algorithm for Distributed Convex Composite Optimization
Authors:
Luyao Guo,
Xinli Shi,
Shaofu Yang,
Jinde Cao
Abstract:
In this paper, we propose a novel Dual Inexact Splitting Algorithm (DISA) for distributed convex composite optimization problems, where the local loss function consists of a smooth term and a possibly nonsmooth term composed with a linear mapping. DISA, for the first time, eliminates the dependence of the convergent step-size range on the Euclidean norm of the linear mapping, while inheriting the…
▽ More
In this paper, we propose a novel Dual Inexact Splitting Algorithm (DISA) for distributed convex composite optimization problems, where the local loss function consists of a smooth term and a possibly nonsmooth term composed with a linear mapping. DISA, for the first time, eliminates the dependence of the convergent step-size range on the Euclidean norm of the linear mapping, while inheriting the advantages of the classic Primal-Dual Proximal Splitting Algorithm (PD-PSA): simple structure and easy implementation. This indicates that DISA can be executed without prior knowledge of the norm, and tiny step-sizes can be avoided when the norm is large. Additionally, we prove sublinear and linear convergence rates of DISA under general convexity and metric subregularity, respectively. Moreover, we provide a variant of DISA with approximate proximal mapping and prove its global convergence and sublinear convergence rate. Numerical experiments corroborate our theoretical analyses and demonstrate a significant acceleration of DISA compared to existing PD-PSAs.
△ Less
Submitted 22 April, 2023; v1 submitted 5 September, 2022;
originally announced September 2022.
-
Variational inequalities of multilayer elastic systems with interlayer friction: existence and uniqueness of solution and convergence of numerical solution
Authors:
Zhizhuo Zhang,
Xiaobing Nie,
Jinde Cao
Abstract:
Based on the mathematical-physical model of pavement mechanics, a multilayer elastic system with interlayer friction conditions is constructed. Given the complex boundary conditions, the corresponding variational inequalities of the partial differential equations are derived, so that the problem can be analyzed under the variational framework. First, the existence and uniqueness of the solution of…
▽ More
Based on the mathematical-physical model of pavement mechanics, a multilayer elastic system with interlayer friction conditions is constructed. Given the complex boundary conditions, the corresponding variational inequalities of the partial differential equations are derived, so that the problem can be analyzed under the variational framework. First, the existence and uniqueness of the solution of the variational inequality is proved; then the approximation error of the numerical solution based on the finite element method is analyzed, and when the finite element space satisfies certain approximation conditions, the convergence of the numerical solution is proved; finally, in the trivial finite element space, the convergence order of the numerical solution is derived. The above conclusions provide basic theoretical support for solving the displacement-strain problem of multilayer elastic systems under the framework of variational inequalities.
△ Less
Submitted 3 June, 2024; v1 submitted 26 July, 2022;
originally announced July 2022.
-
Transmission properties of space-time modulated metamaterials
Authors:
Habib Ammari,
Jinghao Cao,
Xinmeng Zeng
Abstract:
We prove the possibility of achieving exponentially growing wave propagation in space-time modulated media and give an asymptotic analysis of the quasifrequencies in terms of the amplitude of the time modulation at the degenerate points of the folded band structure. Our analysis provides the first proof of existence of k-gaps in the band structures of space-time modulated systems of subwavelength…
▽ More
We prove the possibility of achieving exponentially growing wave propagation in space-time modulated media and give an asymptotic analysis of the quasifrequencies in terms of the amplitude of the time modulation at the degenerate points of the folded band structure. Our analysis provides the first proof of existence of k-gaps in the band structures of space-time modulated systems of subwavelength resonators.
△ Less
Submitted 29 June, 2022;
originally announced June 2022.
-
Boundaries of capture hyperbolic components
Authors:
Jie Cao,
Xiaoguang Wang,
Yongcheng Yin
Abstract:
In complex dynamics, the boundaries of higher dimensional hyperbolic components in holomorphic families of polynomials or rational maps are mysterious objects, whose topological and analytic properties are fundamental problems.
In this paper, we show that in some typical families of polynomials (i.e. algebraic varieties defined by periodic critical relations), the boundary of a capture hyperboli…
▽ More
In complex dynamics, the boundaries of higher dimensional hyperbolic components in holomorphic families of polynomials or rational maps are mysterious objects, whose topological and analytic properties are fundamental problems.
In this paper, we show that in some typical families of polynomials (i.e. algebraic varieties defined by periodic critical relations), the boundary of a capture hyperbolic component $\mathcal H$ is homeomorphic to the sphere $S^{2\dim_\mathbb{C}(\mathcal{H})-1}$. Furthermore, we establish an unexpected identity for the Hausdorff dimension of $\partial \mathcal H$: $$\operatorname{H{.}dim}(\partial\mathcal{H}) = 2 \dim_\mathbb{C}(\mathcal{H})-2+\max_{f\in\partial\mathcal{H}} \operatorname{H{.}dim}(\partial A^J(f)),$$ where $A^J(f)$ is the union of the bounded attracting Fatou components of $f$ associated with the free critical points in the Julia set $J(f)$.
In the proof, some new results with independent interests are discovered.
△ Less
Submitted 15 June, 2022;
originally announced June 2022.
-
Unidirectional edge modes in time-modulated metamaterials
Authors:
Habib Ammari,
Jinghao Cao
Abstract:
We prove the possibility of achieving unidirectional edge modes in time-modulated supercell structures. Such finite structures consist of two trimers repeated periodically. Because of their symmetry, they admit degenerate edge eigenspaces. When the trimers are time-modulated with two opposite orientations, the degenerate eigenspace splits into two one-dimensional eigenspaces described by an analyt…
▽ More
We prove the possibility of achieving unidirectional edge modes in time-modulated supercell structures. Such finite structures consist of two trimers repeated periodically. Because of their symmetry, they admit degenerate edge eigenspaces. When the trimers are time-modulated with two opposite orientations, the degenerate eigenspace splits into two one-dimensional eigenspaces described by an analytical formula, each corresponds to a mode which is localized at one edge of the structure. Our results on the localization and stability of these edge modes with respect to fluctuations in the time-modulation amplitude are illustrated by several numerical simulations.
△ Less
Submitted 18 May, 2022;
originally announced May 2022.
-
Spectral Analysis and its applications for a class of scale-free network based on the weighted m-clique annex operation
Authors:
Zhizhuo Zhang,
Jinde Cao,
Bo Wu
Abstract:
The spectrum of network is an important tool to study the function and dynamic properties of network, and graph operation and product is an effective mechanism to construct a specific local and global topological structure. In this study, a class of weighted $m-$clique annex operation $τ_m^r(\cdot)$ controlled by scale factor $m$ and weight factor $r$ is defined, through which an iterative weighte…
▽ More
The spectrum of network is an important tool to study the function and dynamic properties of network, and graph operation and product is an effective mechanism to construct a specific local and global topological structure. In this study, a class of weighted $m-$clique annex operation $τ_m^r(\cdot)$ controlled by scale factor $m$ and weight factor $r$ is defined, through which an iterative weighted network model $G_t$ with small-world and scale-free properties is constructed. In particular, when the number of iterations $t$ tends to infinity, the network has transfinite fractal property. Then, through the iterative features of the network structure, the iterative relationship of the eigenvalues of the normalized Laplacian matrix corresponding to the network is studied. Accordingly, some applications of the spectrum of the network, including the Kenemy constant, Multiplicative Degree-Kirchhoff index and the number of weighted spanning trees, are further given. In addition, we also study the effect of the two factors controlling network operation on the structure and function of the iterative weighted network $G_t$, so that the network operation can better simulate the real network and have more application potential in the field of artificial network.
△ Less
Submitted 10 April, 2022;
originally announced April 2022.
-
Correction of high-order $L_k$ approximation for subdiffusion
Authors:
Jiankang Shi,
Minghua Chen,
Yubin Yan,
Jianxiong Cao
Abstract:
The subdiffusion equations with a Caputo fractional derivative of order $α\in (0,1)$ arise in a wide variety of practical problems, which is describing the transport processes, in the force-free limit, slower than Brownian diffusion. In this work, we derive the correction schemes of the Lagrange interpolation with degree $k$ ($k\leq 6$) convolution quadrature, called $L_k$ approximation, for the s…
▽ More
The subdiffusion equations with a Caputo fractional derivative of order $α\in (0,1)$ arise in a wide variety of practical problems, which is describing the transport processes, in the force-free limit, slower than Brownian diffusion. In this work, we derive the correction schemes of the Lagrange interpolation with degree $k$ ($k\leq 6$) convolution quadrature, called $L_k$ approximation, for the subdiffusion, which are easy to implement on variable grids. The key step of designing correction algorithm is to calculate the explicit form of the coefficients of $L_k$ approximation by the polylogarithm function or Bose-Einstein integral. To construct a $τ_8$ approximation of Bose-Einstein integral, the desired $(k+1-α)$th-order convergence rate can be proved for the correction $L_k$ scheme with nonsmooth data, which is higher than $k$th-order BDF$k$ method in [Jin, Li, and Zhou, SIAM J. Sci. Comput., 39 (2017), A3129--A3152; Shi and Chen, J. Sci. Comput., (2020) 85:28]. The numerical experiments with spectral method are given to illustrate theoretical results.
△ Less
Submitted 27 December, 2021;
originally announced December 2021.
-
Multiphysics mixed finite element method with Nitsche's technique for Stokes poroelasticity problem
Authors:
Zhihao Ge,
Jin'ge Pang,
Jiwei Cao
Abstract:
In this paper, we propose a multiphysics mixed finite element method with Nitsche's technique for Stokes-poroelasticity problem. Firstly, we present a multiphysics reformulation of poroelasticity part of the original problem by introducing two pseudo-pressures to reveal the underlying deformation and diffusion multi physical processes in the Stokes-poroelasticity problem. Then, we prove the existe…
▽ More
In this paper, we propose a multiphysics mixed finite element method with Nitsche's technique for Stokes-poroelasticity problem. Firstly, we present a multiphysics reformulation of poroelasticity part of the original problem by introducing two pseudo-pressures to reveal the underlying deformation and diffusion multi physical processes in the Stokes-poroelasticity problem. Then, we prove the existence and uniqueness of weak solution of the reformulated and original problem. And we use Nitsche's technique to approximate the coupling condition at the interface to propose a loosely-coupled time-stepping method -- multiphysics mixed finite element method for space variables, and we decouple the reformulated problem into three sub-problems at each time step -- a Stokes problem, a generalized Stokes problem and a mixed diffusion problem. Also, we give the stability analysis and error estimates of the loosely-coupled time-stepping method.
△ Less
Submitted 23 December, 2021;
originally announced December 2021.
-
Generalized q-difference equations for general q-polynomials with double q-binomial coefficients
Authors:
Jian Cao,
Sama Arjika,
Mahouton Norbert Hounkonnou
Abstract:
In this paper, we use the generalized q-polynomials with double q-binomial coefficients and homogeneous q-operators [J. Difference Equ. Appl. 20 (2014), 837--851.] to construct q-difference equations with seven variables, which generalize recent works of Jia et al [Symmetry 2021, 13, 1222.]. In addition, we derive Rogers formulas, extended Rogers formulas and Srivastava--Agarwal type bilinear gene…
▽ More
In this paper, we use the generalized q-polynomials with double q-binomial coefficients and homogeneous q-operators [J. Difference Equ. Appl. 20 (2014), 837--851.] to construct q-difference equations with seven variables, which generalize recent works of Jia et al [Symmetry 2021, 13, 1222.]. In addition, we derive Rogers formulas, extended Rogers formulas and Srivastava--Agarwal type bilinear generating functions for generalized q-polynomials, which generalize generating functions for Cigler's polynomials [J. Difference Equ. Appl. 24 (2018), 479--502.]. Finally, we also derive mixed generating functions using q-difference equations.
△ Less
Submitted 17 December, 2021;
originally announced December 2021.