For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of calculus of variations, multiplicative commutation r...For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L 2 norm are derived to determine the error, in the approximate solution.展开更多
After discretization by the finite volume method,the numerical solution of fractional diffusion equations leads to a linear system with the Toeplitz-like structure.The theoretical analysis gives sufficient conditions ...After discretization by the finite volume method,the numerical solution of fractional diffusion equations leads to a linear system with the Toeplitz-like structure.The theoretical analysis gives sufficient conditions to guarantee the positive-definite property of the discretized matrix.Moreover,we develop a class of positive-definite operator splitting iteration methods for the numerical solution of fractional diffusion equations,which is unconditionally convergent for any positive constant.Meanwhile,the iteration methods introduce a new preconditioner for Krylov subspace methods.Numerical experiments verify the convergence of the positive-definite operator splitting iteration methods and show the efficiency of the proposed preconditioner,compared with the existing approaches.展开更多
In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor depos...In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.展开更多
The operator-splitting methods for the mathematic model of one kind of oin reactions for the problem of groundwater are considered.Optimal error estimates in L 2 and H 1 norm are obtained for the approximation solut...The operator-splitting methods for the mathematic model of one kind of oin reactions for the problem of groundwater are considered.Optimal error estimates in L 2 and H 1 norm are obtained for the approximation solution.展开更多
The Operator Splitting method is applied to differential equations occurring as mathematical models in financial models. This paper provides various operator splitting methods to obtain an effective and accurate solut...The Operator Splitting method is applied to differential equations occurring as mathematical models in financial models. This paper provides various operator splitting methods to obtain an effective and accurate solution to the Black-Scholes equation with appropriate boundary conditions for a European option pricing problem. Finally brief comparisons of option prices are given by different models.展开更多
Let H;, H;, H;be real Hilbert spaces, let A : H;→ H;, B : H;→ H;be two bounded linear operators. The split equality common fixed point problem(SECFP) in the infinite-dimensional Hilbert spaces introduced by Moudaf...Let H;, H;, H;be real Hilbert spaces, let A : H;→ H;, B : H;→ H;be two bounded linear operators. The split equality common fixed point problem(SECFP) in the infinite-dimensional Hilbert spaces introduced by Moudafi(Alternating CQ-algorithm for convex feasibility and split fixed-point problems. Journal of Nonlinear and Convex Analysis)is to find x ∈ F(U), y ∈ F(T) such that Ax = By,(1)where U : H;→ H;and T : H;→ H;are two nonlinear operators with nonempty fixed point sets F(U) = {x ∈ H;: Ux = x} and F(T) = {x ∈ H;: Tx = x}. Note that,by taking B = I and H;= H;in(1), we recover the split fixed point problem originally introduced in Censor and Segal. Recently, Moudafi introduced alternating CQ-algorithms and simultaneous iterative algorithms with weak convergence for the SECFP(1) of firmly quasi-nonexpansive operators. In this paper, we introduce two viscosity iterative algorithms for the SECFP(1) governed by the general class of quasi-nonexpansive operators. We prove the strong convergence of algorithms. Our results improve and extend previously discussed related problems and algorithms.展开更多
Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by t...Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.展开更多
Based on non-maximally entangled four-particle cluster states, we propose a new hierarchical information splitting protocol to probabilistically realize the quantum state sharing of an arbitrary unknown two-qubit stat...Based on non-maximally entangled four-particle cluster states, we propose a new hierarchical information splitting protocol to probabilistically realize the quantum state sharing of an arbitrary unknown two-qubit state. In this scheme, the sender transmits the two-qubit secret state to three agents who are divided into two grades with two Bell-state measurements,and broadcasts the measurement results via a classical channel. One agent is in the upper grade and two agents are in the lower grade. The agent in the upper grade only needs to cooperate with one of the other two agents to recover the secret state but both of the agents in the lower grade need help from all of the agents. Every agent who wants to recover the secret state needs to introduce two ancillary qubits and performs a positive operator-valued measurement(POVM) instead of the usual projective measurement. Moreover, due to the symmetry of the cluster state, we extend this protocol to multiparty agents.展开更多
This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we stu...This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.展开更多
Deep neural network is a powerful tool for many tasks.Understanding why it is so successful and providing a mathematical explanation is an important problem and has been one popular research direction in past years.In...Deep neural network is a powerful tool for many tasks.Understanding why it is so successful and providing a mathematical explanation is an important problem and has been one popular research direction in past years.In the literature of mathematical analysis of deep neural networks,a lot of works is dedicated to establishing representation theories.How to make connections between deep neural networks and mathematical algorithms is still under development.In this paper,we give an algorithmic explanation for deep neural networks,especially in their connections with operator splitting.We show that with certain splitting strategies,operator-splitting methods have the same structure as networks.Utilizing this connection and the Potts model for image segmentation,two networks inspired by operator-splitting methods are proposed.The two networks are essentially two operator-splitting algorithms solving the Potts model.Numerical experiments are presented to demonstrate the effectiveness of the proposed networks.展开更多
In this paper, we study the propagation and its failure to propagate (pinning) of a travelling wave in a Nagumo type equation, an equation that describes impulse propagation in nerve axons that also models population ...In this paper, we study the propagation and its failure to propagate (pinning) of a travelling wave in a Nagumo type equation, an equation that describes impulse propagation in nerve axons that also models population growth with Allee effect. An analytical solution is derived for the traveling wave and the work is extended to a discrete formulation with a piecewise linear reaction function. We propose an operator splitting numerical scheme to solve the equation and demonstrate that the wave either propagates or gets pinned based on how the spatial mesh is chosen.展开更多
The Peaceman-Rachford splitting method is efficient for minimizing a convex optimization problem with a separable objective function and linear constraints.However,its convergence was not guaranteed without extra requ...The Peaceman-Rachford splitting method is efficient for minimizing a convex optimization problem with a separable objective function and linear constraints.However,its convergence was not guaranteed without extra requirements.He et al.(SIAM J.Optim.24:1011-1040,2014)proved the convergence of a strictly contractive Peaceman-Rachford splitting method by employing a suitable underdetermined relaxation factor.In this paper,we further extend the so-called strictly contractive Peaceman-Rachford splitting method by using two different relaxation factors.Besides,motivated by the recent advances on the ADMM type method with indefinite proximal terms,we employ the indefinite proximal term in the strictly contractive Peaceman-Rachford splitting method.We show that the proposed indefinite-proximal strictly contractive Peaceman-Rachford splitting method is convergent and also prove the o(1/t)convergence rate in the nonergodic sense.The numerical tests on the l 1 regularized least square problem demonstrate the efficiency of the proposed method.展开更多
基金the Major State Basic Research Program of China(19990328)NNSF of China(19871051,19972039) the Doctorate Foundation of the State Education Commission
文摘For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L 2 norm are derived to determine the error, in the approximate solution.
基金This work was supported by the National Natural Science Foundation of China(No.11971354)The author Yi-Shu Du acknowledges the financial support from the China Scholarship Council(File No.201906260146).
文摘After discretization by the finite volume method,the numerical solution of fractional diffusion equations leads to a linear system with the Toeplitz-like structure.The theoretical analysis gives sufficient conditions to guarantee the positive-definite property of the discretized matrix.Moreover,we develop a class of positive-definite operator splitting iteration methods for the numerical solution of fractional diffusion equations,which is unconditionally convergent for any positive constant.Meanwhile,the iteration methods introduce a new preconditioner for Krylov subspace methods.Numerical experiments verify the convergence of the positive-definite operator splitting iteration methods and show the efficiency of the proposed preconditioner,compared with the existing approaches.
文摘In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.
基金the National Sciences Foundation of China and the Doctorial Program of Higher Edua-tion
文摘The operator-splitting methods for the mathematic model of one kind of oin reactions for the problem of groundwater are considered.Optimal error estimates in L 2 and H 1 norm are obtained for the approximation solution.
文摘The Operator Splitting method is applied to differential equations occurring as mathematical models in financial models. This paper provides various operator splitting methods to obtain an effective and accurate solution to the Black-Scholes equation with appropriate boundary conditions for a European option pricing problem. Finally brief comparisons of option prices are given by different models.
基金supported by National Natural Science Foundation of China(61503385)Fundamental Research Funds for the Central Universities of China(3122016L002)
文摘Let H;, H;, H;be real Hilbert spaces, let A : H;→ H;, B : H;→ H;be two bounded linear operators. The split equality common fixed point problem(SECFP) in the infinite-dimensional Hilbert spaces introduced by Moudafi(Alternating CQ-algorithm for convex feasibility and split fixed-point problems. Journal of Nonlinear and Convex Analysis)is to find x ∈ F(U), y ∈ F(T) such that Ax = By,(1)where U : H;→ H;and T : H;→ H;are two nonlinear operators with nonempty fixed point sets F(U) = {x ∈ H;: Ux = x} and F(T) = {x ∈ H;: Tx = x}. Note that,by taking B = I and H;= H;in(1), we recover the split fixed point problem originally introduced in Censor and Segal. Recently, Moudafi introduced alternating CQ-algorithms and simultaneous iterative algorithms with weak convergence for the SECFP(1) of firmly quasi-nonexpansive operators. In this paper, we introduce two viscosity iterative algorithms for the SECFP(1) governed by the general class of quasi-nonexpansive operators. We prove the strong convergence of algorithms. Our results improve and extend previously discussed related problems and algorithms.
文摘Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.
基金Supported by the Major Project of National Natural Science Foundation of China (No.20409205) and National High Technology Research and Development Program of China (No.G20070040).
基金Project supported by the National Natural Science Foundation of China(Grant No.61671087)
文摘Based on non-maximally entangled four-particle cluster states, we propose a new hierarchical information splitting protocol to probabilistically realize the quantum state sharing of an arbitrary unknown two-qubit state. In this scheme, the sender transmits the two-qubit secret state to three agents who are divided into two grades with two Bell-state measurements,and broadcasts the measurement results via a classical channel. One agent is in the upper grade and two agents are in the lower grade. The agent in the upper grade only needs to cooperate with one of the other two agents to recover the secret state but both of the agents in the lower grade need help from all of the agents. Every agent who wants to recover the secret state needs to introduce two ancillary qubits and performs a positive operator-valued measurement(POVM) instead of the usual projective measurement. Moreover, due to the symmetry of the cluster state, we extend this protocol to multiparty agents.
文摘This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.
基金Supported by the Program for New Century Excellent Talents at the University of China under Grant No.NCET-06-0554the National Natural Science Foundation of China under Grant Nos.10975001,60677001,10747146,and 10874122+3 种基金the Science-technology Fund of Anhui Province for Outstanding Youth under Grant No.06042087the Key Fund of the Ministry of Education of China under Grant No.206063 the General Fund of the Educational Committee of Anhui Province under Grant No.2006KJ260Bthe Natural Science Foundation of Guangdong Province under Grant Nos.06300345 and 7007806
基金supported by HKBU 179356,NSFC 12201530 and HKRGC ECS 22302123supported by NSFC/RGC grant N-HKBU214-19 and NORCE Kompetanseoppbygging programsupported by HKRGC GRF grants CityU1101120,CityU11309922,CRF grant C1013-21GF,and HKRGC-NSFC Grant N CityU214/19.
文摘Deep neural network is a powerful tool for many tasks.Understanding why it is so successful and providing a mathematical explanation is an important problem and has been one popular research direction in past years.In the literature of mathematical analysis of deep neural networks,a lot of works is dedicated to establishing representation theories.How to make connections between deep neural networks and mathematical algorithms is still under development.In this paper,we give an algorithmic explanation for deep neural networks,especially in their connections with operator splitting.We show that with certain splitting strategies,operator-splitting methods have the same structure as networks.Utilizing this connection and the Potts model for image segmentation,two networks inspired by operator-splitting methods are proposed.The two networks are essentially two operator-splitting algorithms solving the Potts model.Numerical experiments are presented to demonstrate the effectiveness of the proposed networks.
文摘In this paper, we study the propagation and its failure to propagate (pinning) of a travelling wave in a Nagumo type equation, an equation that describes impulse propagation in nerve axons that also models population growth with Allee effect. An analytical solution is derived for the traveling wave and the work is extended to a discrete formulation with a piecewise linear reaction function. We propose an operator splitting numerical scheme to solve the equation and demonstrate that the wave either propagates or gets pinned based on how the spatial mesh is chosen.
基金supported by the Natural Science Foundation of Jiangsu Province (Grant No.BK20210267)supported by the National Natural Science Foundation of China (Grant No.11971239)+1 种基金the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province (Grant No.21KJA110002)supported by the National Natural Science Foundation of China (Grant Nos.12131004,11625105).
文摘The Peaceman-Rachford splitting method is efficient for minimizing a convex optimization problem with a separable objective function and linear constraints.However,its convergence was not guaranteed without extra requirements.He et al.(SIAM J.Optim.24:1011-1040,2014)proved the convergence of a strictly contractive Peaceman-Rachford splitting method by employing a suitable underdetermined relaxation factor.In this paper,we further extend the so-called strictly contractive Peaceman-Rachford splitting method by using two different relaxation factors.Besides,motivated by the recent advances on the ADMM type method with indefinite proximal terms,we employ the indefinite proximal term in the strictly contractive Peaceman-Rachford splitting method.We show that the proposed indefinite-proximal strictly contractive Peaceman-Rachford splitting method is convergent and also prove the o(1/t)convergence rate in the nonergodic sense.The numerical tests on the l 1 regularized least square problem demonstrate the efficiency of the proposed method.