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.展开更多
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.展开更多
We present a new family of fourth-order splitting methods with positive coefficients especially tailored for the time integration of linear parabolic problems and,in particular,for the time dependent Schrodinger equat...We present a new family of fourth-order splitting methods with positive coefficients especially tailored for the time integration of linear parabolic problems and,in particular,for the time dependent Schrodinger equation,both in real and imaginary time.They are based on the use of a double commutator and a modified processor,and are more efficient than other widely used schemes found in the literature.Moreover,for certain potentials,they achieve order six.Several examples in one,two and three dimensions clearly illustrate the computational advantages of the new schemes.展开更多
The method of recovering a low-rank matrix with an unknown fraction whose entries are arbitrarily corrupted is known as the robust principal component analysis （RPCA）. This RPCA problem, under some conditions, can b...The method of recovering a low-rank matrix with an unknown fraction whose entries are arbitrarily corrupted is known as the robust principal component analysis （RPCA）. This RPCA problem, under some conditions, can be exactly solved via convex optimization by minimizing a combination of the nuclear norm and the 11 norm. In this paper, an algorithm based on the Douglas-Rachford splitting method is proposed for solving the RPCA problem. First, the convex optimization problem is solved by canceling the constraint of the variables, and ~hen the proximity operators of the objective function are computed alternately. The new algorithm can exactly recover the low-rank and sparse components simultaneously, and it is proved to be convergent. Numerical simulations demonstrate the practical utility of the proposed algorithm.展开更多
It has been shown that the alternating direction method of multipliers(ADMM)is not necessarily convergent when it is directly extended to a multiple-block linearly constrained convex minimization model with an objecti...It has been shown that the alternating direction method of multipliers(ADMM)is not necessarily convergent when it is directly extended to a multiple-block linearly constrained convex minimization model with an objective function that is in the sum of more than two functions without coupled variables.Recently,we pro-posed the block-wise ADMM,which was obtained by regrouping the variables and functions of such a model as two blocks and then applying the original ADMM in block-wise.This note is a further study on this topic with the purpose of showing that a well-known relaxation factor proposed by Fortin and Glowinski for iteratively updat-ing the Lagrangian multiplier of the originalADMMcan also be used in the block-wise ADMM.We thus propose the block-wise ADMM with Fortin and Glowinski’s relax-ation factor for the multiple-block convex minimization model.Like the block-wise ADMM,we also suggest further decomposing the resulting subproblems and regular-izing them by proximal terms to ensure the convergence.For the block-wise ADMM with Fortin and Glowinski's relaxation factor,its convergence and worst-case conver-gence rate measured by the iteration complexity in the ergodic sense are derived.展开更多
In this paper the numerical solution of the two-dimensional sine-Gordon equation is studied.Split local artificial boundary conditions are obtained by the operator splitting method.Then the original problem is reduced...In this paper the numerical solution of the two-dimensional sine-Gordon equation is studied.Split local artificial boundary conditions are obtained by the operator splitting method.Then the original problem is reduced to an initial boundary value problem on a bounded computational domain,which can be solved by the finite differencemethod.Several numerical examples are provided to demonstrate the effectiveness and accuracy of the proposed method,and some interesting propagation and collision behaviors of the solitary wave solutions are observed.展开更多
In this paper the slot method is used in the computation of 2-dimensional flow with transient boundary.The slot located inside the computation zone is placed in each space grid, parallel to both X and Y axes.Combined ...In this paper the slot method is used in the computation of 2-dimensional flow with transient boundary.The slot located inside the computation zone is placed in each space grid, parallel to both X and Y axes.Combined with the operator splitting method,the numerical simu- lation of the flow patterns of Yuqiao Reservior is made.The calculated results are in good agree- ment with observed data.展开更多
In this paper,a modified boundary fitted coordinate method is presented.The transformation equations are solved by finite analytic method.The control equations after trans- formation are solved by splitting operator m...In this paper,a modified boundary fitted coordinate method is presented.The transformation equations are solved by finite analytic method.The control equations after trans- formation are solved by splitting operator method.The results of calculation consist well with experiments.展开更多
文摘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.
文摘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 Ministerio de Ciencia e Innovacion(Spain)through projects PID2019-104927GB-C21 and PID2019-104927GB-C22,MCIN/AEI/10.13039/501100011033,ERDF(“A way of making Europe”)the support of the Conselleria d’Innovacio,Universitats,Ciencia i Societat Digital from the Generalitat Valenciana(Spain)through project CIAICO/2021/180.
文摘We present a new family of fourth-order splitting methods with positive coefficients especially tailored for the time integration of linear parabolic problems and,in particular,for the time dependent Schrodinger equation,both in real and imaginary time.They are based on the use of a double commutator and a modified processor,and are more efficient than other widely used schemes found in the literature.Moreover,for certain potentials,they achieve order six.Several examples in one,two and three dimensions clearly illustrate the computational advantages of the new schemes.
基金supported by the National Natural Science Foundation of China(No.61271014)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20124301110003)the Graduated Students Innovation Fund of Hunan Province(No.CX2012B238)
文摘The method of recovering a low-rank matrix with an unknown fraction whose entries are arbitrarily corrupted is known as the robust principal component analysis （RPCA）. This RPCA problem, under some conditions, can be exactly solved via convex optimization by minimizing a combination of the nuclear norm and the 11 norm. In this paper, an algorithm based on the Douglas-Rachford splitting method is proposed for solving the RPCA problem. First, the convex optimization problem is solved by canceling the constraint of the variables, and ~hen the proximity operators of the objective function are computed alternately. The new algorithm can exactly recover the low-rank and sparse components simultaneously, and it is proved to be convergent. Numerical simulations demonstrate the practical utility of the proposed algorithm.
基金Bing-Sheng He and Ming-Hua Xu were supported by the National Natural Science Foundation of China(No.11471156)Xiao-Ming Yuan was supported by the General Research Fund from Hong Kong Research Grants Council(No.HKBU 12313516).
文摘It has been shown that the alternating direction method of multipliers(ADMM)is not necessarily convergent when it is directly extended to a multiple-block linearly constrained convex minimization model with an objective function that is in the sum of more than two functions without coupled variables.Recently,we pro-posed the block-wise ADMM,which was obtained by regrouping the variables and functions of such a model as two blocks and then applying the original ADMM in block-wise.This note is a further study on this topic with the purpose of showing that a well-known relaxation factor proposed by Fortin and Glowinski for iteratively updat-ing the Lagrangian multiplier of the originalADMMcan also be used in the block-wise ADMM.We thus propose the block-wise ADMM with Fortin and Glowinski’s relax-ation factor for the multiple-block convex minimization model.Like the block-wise ADMM,we also suggest further decomposing the resulting subproblems and regular-izing them by proximal terms to ensure the convergence.For the block-wise ADMM with Fortin and Glowinski's relaxation factor,its convergence and worst-case conver-gence rate measured by the iteration complexity in the ergodic sense are derived.
基金support form the National Natural Science Foundation of China(Grant No.10971116).
文摘In this paper the numerical solution of the two-dimensional sine-Gordon equation is studied.Split local artificial boundary conditions are obtained by the operator splitting method.Then the original problem is reduced to an initial boundary value problem on a bounded computational domain,which can be solved by the finite differencemethod.Several numerical examples are provided to demonstrate the effectiveness and accuracy of the proposed method,and some interesting propagation and collision behaviors of the solitary wave solutions are observed.
文摘In this paper the slot method is used in the computation of 2-dimensional flow with transient boundary.The slot located inside the computation zone is placed in each space grid, parallel to both X and Y axes.Combined with the operator splitting method,the numerical simu- lation of the flow patterns of Yuqiao Reservior is made.The calculated results are in good agree- ment with observed data.
文摘In this paper,a modified boundary fitted coordinate method is presented.The transformation equations are solved by finite analytic method.The control equations after trans- formation are solved by splitting operator method.The results of calculation consist well with experiments.
