For solving nonlinear parabolic equation on massive parallel computers, the construction of parallel difference schemes with simple design, high parallelism and unconditional stability and second order global accuracy...For solving nonlinear parabolic equation on massive parallel computers, the construction of parallel difference schemes with simple design, high parallelism and unconditional stability and second order global accuracy in space, has long been desired. In the present work, a new kind of general parallel difference schemes for the nonlinear parabolic system is proposed. The general parallel difference schemes include, among others, two new parallel schemes. In one of them, to obtain the interface values on the interface of sub-domains an explicit scheme of Jacobian type is employed, and then the fully implicit scheme is used in the sub-domains. Here, in the explicit scheme of Jacobian type, the values at the points being adjacent to the interface points are taken as the linear combination of values of previous two time layers at the adjoining points of the inner interface. For the construction of another new parallel difference scheme, the main procedure is as follows. Firstly the linear combination of values of previous two time layers at the interface points among the sub-domains is used as the (Dirichlet) boundary condition for solving the sub-domain problems. Then the values in the sub-domains are calculated by the fully implicit scheme. Finally the interface values are computed by the fully implicit scheme, and in fact these calculations of the last step are explicit since the values adjacent to the interface points have been obtained in the previous step. The existence, uniqueness, unconditional stability and the second order accuracy of the discrete vector solutions for the parallel difference schemes are proved. Numerical results are presented to examine the stability, accuracy and parallelism of the parallel schemes.展开更多
Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite ele...Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method (FEM) is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.展开更多
The sequences defined in Example 3 and Example 4 do not serve our purpose for any λ = (λn). Because this sequences are just the sequences x = (xk) = (k) and x = (xk) = (1) respectively and any term of thes...The sequences defined in Example 3 and Example 4 do not serve our purpose for any λ = (λn). Because this sequences are just the sequences x = (xk) = (k) and x = (xk) = (1) respectively and any term of these sequences can not be 0. In this short not we give Example 3* and Example 4* to show that the inclusions given in Theorem 2.4 and Theorem 2.9 are strict for some λ = (λn) , α and β such that 0 α β ≤ 1.展开更多
In this paper, a new class of triangular summation operators based on the equidistant nodes was constructed. It is proved that this class of operators converges uniformly to arbitrary continuous fimctions with the per...In this paper, a new class of triangular summation operators based on the equidistant nodes was constructed. It is proved that this class of operators converges uniformly to arbitrary continuous fimctions with the period 2π on the whole axis, Fttrthermore, the best approximation order and the highest convergence order are obtained. In contrast to certain operators constructed by Bernstein and Kis in the previous works, the convergence properties of the new operator constructed in this paper are superior.展开更多
An equivalent description of u-uniform convergence is presented first. Then the relations among the order convergence, u-uniform convergence and norm convergence of sequences are discussed in Riesz spaces. An equivale...An equivalent description of u-uniform convergence is presented first. Then the relations among the order convergence, u-uniform convergence and norm convergence of sequences are discussed in Riesz spaces. An equivalence of the three convergences is brought forward; namely, {fn} is a u-uniform Cauchy sequence. Finally the relations among the three convergences of sequences are also extended to the relations among the convergences of nets in Riesz spaces.展开更多
In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order...In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order five. Numerical examples show that the new methods are comparable with the well known existing methods and give better results in many aspects.展开更多
There are several ways that can be used to classify or compare iterative methods for nonlinear equations,for instance;order of convergence,informational efficiency,and efficiency index.In this work,we use another way,...There are several ways that can be used to classify or compare iterative methods for nonlinear equations,for instance;order of convergence,informational efficiency,and efficiency index.In this work,we use another way,namely the basins of attraction of the method.The purpose of this study is to compare several iterative schemes for nonlinear equations.All the selected schemes are of the third-order of convergence and most of them have the same efficiency index.The comparison depends on the basins of attraction of the iterative techniques when applied on several polynomials of different degrees.As a comparison,we determine the CPU time(in seconds)needed by each scheme to obtain the basins of attraction,besides,we illustrate the area of convergence of these schemes by finding the number of convergent and divergent points in a selected range for all methods.Comparisons confirm the fact that basins of attraction differ for iterative methods of different orders,furthermore,they vary for iterative methods of the same order even if they have the same efficiency index.Consequently,this leads to the need for a new index that reflects the real efficiency of the iterative scheme instead of the commonly used efficiency index.展开更多
In this paper,we propose a fifth-order scheme for solving systems of nonlinear equations.The convergence analysis of the proposed technique is discussed.The proposed method is generalized and extended to be of any odd...In this paper,we propose a fifth-order scheme for solving systems of nonlinear equations.The convergence analysis of the proposed technique is discussed.The proposed method is generalized and extended to be of any odd order of the form 2n1.The scheme is composed of three steps,of which the first two steps are based on the two-step Homeier’s method with cubic convergence,and the last is a Newton step with an appropriate approximation for the derivative.Every iteration of the presented method requires the evaluation of two functions,two Fréchet derivatives,and three matrix inversions.A comparison between the efficiency index and the computational efficiency index of the presented scheme with existing methods is performed.The basins of attraction of the proposed scheme illustrated and compared to other schemes of the same order.Different test problems including large systems of equations are considered to compare the performance of the proposed method according to other methods of the same order.As an application,we apply the new scheme to some real-life problems,including the mixed Hammerstein integral equation and Burgers’equation.Comparisons and examples show that the presented method is efficient and comparable to the existing techniques of the same order.展开更多
In this paper, by using the Nevanlinna Theory on angular domain, we establish a theorem which concerns the growth of entire function and his zero. As an application, we survey the location of zero of higher order diff...In this paper, by using the Nevanlinna Theory on angular domain, we establish a theorem which concerns the growth of entire function and his zero. As an application, we survey the location of zero of higher order differential equation, which can be regarded as an alternating but precise version of Wu and Yi.展开更多
In this paper, based on the idea of El-Mistikawy and Werle[1] we construct a difference scheme for a singularly perturbed self-adjoint ordinary differential equation in conservation form. We prove that it is a uniform...In this paper, based on the idea of El-Mistikawy and Werle[1] we construct a difference scheme for a singularly perturbed self-adjoint ordinary differential equation in conservation form. We prove that it is a uniformly convergent second order scheme.展开更多
A fractional step scheme with modified characteristic finite differences run- ning in a parallel arithmetic is presented to simulate a nonlinear percolation system of multilayer dynamics of fluids in a porous medium w...A fractional step scheme with modified characteristic finite differences run- ning in a parallel arithmetic is presented to simulate a nonlinear percolation system of multilayer dynamics of fluids in a porous medium with moving boundary values. With the help of theoretical techniques including the change of regions, piecewise threefold quadratic interpolation, calculus of variations, multiplicative commutation rule of differ- ence operators, multiplicative commutation rule of difference operators, decomposition of high order difference operators, induction hypothesis, and prior estimates, an optimal order in 12 norm is displayed to complete the convergence analysis of the numerical algo- rithm. Some numerical results arising in the actual simulation of migration-accumulation of oil resources by this method are listed in the last section.展开更多
The main purpose of this paper is to introduce some approximation properties of a Kantorovich kind q-Bernstein operators related to B′ezier basis functions with shape parameterλ∈[−1,1].Firstly,we compute some basic...The main purpose of this paper is to introduce some approximation properties of a Kantorovich kind q-Bernstein operators related to B′ezier basis functions with shape parameterλ∈[−1,1].Firstly,we compute some basic results such as moments and central moments,and derive the Korovkin type approximation theorem for these operators.Next,we estimate the order of convergence in terms of the usual modulus of continuity,for the functions belong to Lipschitz-type class and Peetre’s K-functional,respectively.Lastly,with the aid of Maple software,we present the comparison of the convergence of these newly defined operators to the certain function with some graphical illustrations and error estimation table.展开更多
This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound cor...This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound corrector formulas for rectangular rule. Examples of numerical calculation have validated theoretical analysis.展开更多
The truncation error of improved Cotes formula is presented in this paper. It also displays an analysis on convergence order of improved Cotes formula. Examples of numerical calculation is given in the end.
In this paper, an interpolation polynomial operator F n(f; l,x) is constructed based on the zeros of a kind of Jacobi polynomials as the interpolation nodes. For any continuous function f(x)∈C b [-1,1] ...In this paper, an interpolation polynomial operator F n(f; l,x) is constructed based on the zeros of a kind of Jacobi polynomials as the interpolation nodes. For any continuous function f(x)∈C b [-1,1] (0≤b≤l) F n(f; l,x) converges to f(x) uniformly, where l is an odd number.展开更多
In this paper we discuss two-stage Miistein methods for solving Ito stochastic differential equations (SDEs). Six fully explicit methods (TSM 1 -- TSM 6) are given in this paper. Their order of strong convergence ...In this paper we discuss two-stage Miistein methods for solving Ito stochastic differential equations (SDEs). Six fully explicit methods (TSM 1 -- TSM 6) are given in this paper. Their order of strong convergence is proved. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of Ito SDEs.展开更多
In this paper, the regularization with closed linear operators is used to solve an operator equation of the first kind. When all initial data are known approximately, we choose the regular parameter by using the gener...In this paper, the regularization with closed linear operators is used to solve an operator equation of the first kind. When all initial data are known approximately, we choose the regular parameter by using the general Arcangeli's criterion to give the convergence and the asymptotic orders of convergence of the regular solution.展开更多
This paper presents a Modified Formula for Cotes rule with fifths derivatives of endpoint and its truncation error. It also displays an analysis on convergence order of compound formula. Though compound modified formu...This paper presents a Modified Formula for Cotes rule with fifths derivatives of endpoint and its truncation error. It also displays an analysis on convergence order of compound formula. Though compound modified formula for Cotes rule with endpoint derivatives just calculates a newly-added fifths derivative of the two endpoints for each time compared with compound Cotes formula calculation, there are 2 more ranks of the convergence order in this modified formula. Examples of numerical calculation have validated theoretical analysis.展开更多
A fifth-order family of an iterative method for solving systems of nonlinear equations and highly nonlinear boundary value problems has been developed in this paper.Convergence analysis demonstrates that the local ord...A fifth-order family of an iterative method for solving systems of nonlinear equations and highly nonlinear boundary value problems has been developed in this paper.Convergence analysis demonstrates that the local order of convergence of the numerical method is five.The computer algebra system CAS-Maple,Mathematica,or MATLAB was the primary tool for dealing with difficult problems since it allows for the handling and manipulation of complex mathematical equations and other mathematical objects.Several numerical examples are provided to demonstrate the properties of the proposed rapidly convergent algorithms.A dynamic evaluation of the presented methods is also presented utilizing basins of attraction to analyze their convergence behavior.Aside from visualizing iterative processes,this methodology provides useful information on iterations,such as the number of diverging-converging points and the average number of iterations as a function of initial points.Solving numerous highly nonlinear boundary value problems and large nonlinear systems of equations of higher dimensions demonstrate the performance,efficiency,precision,and applicability of a newly presented technique.展开更多
基金The project is supported by the Special Funds for Major State Basic Research Projects 2005CB321703, the National Nature Science Foundation of China (No. 10476002, 60533020).
文摘For solving nonlinear parabolic equation on massive parallel computers, the construction of parallel difference schemes with simple design, high parallelism and unconditional stability and second order global accuracy in space, has long been desired. In the present work, a new kind of general parallel difference schemes for the nonlinear parabolic system is proposed. The general parallel difference schemes include, among others, two new parallel schemes. In one of them, to obtain the interface values on the interface of sub-domains an explicit scheme of Jacobian type is employed, and then the fully implicit scheme is used in the sub-domains. Here, in the explicit scheme of Jacobian type, the values at the points being adjacent to the interface points are taken as the linear combination of values of previous two time layers at the adjoining points of the inner interface. For the construction of another new parallel difference scheme, the main procedure is as follows. Firstly the linear combination of values of previous two time layers at the interface points among the sub-domains is used as the (Dirichlet) boundary condition for solving the sub-domain problems. Then the values in the sub-domains are calculated by the fully implicit scheme. Finally the interface values are computed by the fully implicit scheme, and in fact these calculations of the last step are explicit since the values adjacent to the interface points have been obtained in the previous step. The existence, uniqueness, unconditional stability and the second order accuracy of the discrete vector solutions for the parallel difference schemes are proved. Numerical results are presented to examine the stability, accuracy and parallelism of the parallel schemes.
基金the National Natural Science Foundation of China(No.50678093)Program for Changjiang Scholars and Innovative Research Team in University(No.IRT00736)
文摘Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method (FEM) is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.
文摘The sequences defined in Example 3 and Example 4 do not serve our purpose for any λ = (λn). Because this sequences are just the sequences x = (xk) = (k) and x = (xk) = (1) respectively and any term of these sequences can not be 0. In this short not we give Example 3* and Example 4* to show that the inclusions given in Theorem 2.4 and Theorem 2.9 are strict for some λ = (λn) , α and β such that 0 α β ≤ 1.
文摘In this paper, a new class of triangular summation operators based on the equidistant nodes was constructed. It is proved that this class of operators converges uniformly to arbitrary continuous fimctions with the period 2π on the whole axis, Fttrthermore, the best approximation order and the highest convergence order are obtained. In contrast to certain operators constructed by Bernstein and Kis in the previous works, the convergence properties of the new operator constructed in this paper are superior.
文摘An equivalent description of u-uniform convergence is presented first. Then the relations among the order convergence, u-uniform convergence and norm convergence of sequences are discussed in Riesz spaces. An equivalence of the three convergences is brought forward; namely, {fn} is a u-uniform Cauchy sequence. Finally the relations among the three convergences of sequences are also extended to the relations among the convergences of nets in Riesz spaces.
文摘In this paper, we present and analyze a family of fifth-order iterative methods free from second derivative for solving nonlinear equations. It is established that the family of iterative methods has convergence order five. Numerical examples show that the new methods are comparable with the well known existing methods and give better results in many aspects.
基金We are grateful for the financial support from UKM’s research Grant GUP-2019-033。
文摘There are several ways that can be used to classify or compare iterative methods for nonlinear equations,for instance;order of convergence,informational efficiency,and efficiency index.In this work,we use another way,namely the basins of attraction of the method.The purpose of this study is to compare several iterative schemes for nonlinear equations.All the selected schemes are of the third-order of convergence and most of them have the same efficiency index.The comparison depends on the basins of attraction of the iterative techniques when applied on several polynomials of different degrees.As a comparison,we determine the CPU time(in seconds)needed by each scheme to obtain the basins of attraction,besides,we illustrate the area of convergence of these schemes by finding the number of convergent and divergent points in a selected range for all methods.Comparisons confirm the fact that basins of attraction differ for iterative methods of different orders,furthermore,they vary for iterative methods of the same order even if they have the same efficiency index.Consequently,this leads to the need for a new index that reflects the real efficiency of the iterative scheme instead of the commonly used efficiency index.
基金We are grateful for the financial support from UKM’s research Grant GUP-2019-033.
文摘In this paper,we propose a fifth-order scheme for solving systems of nonlinear equations.The convergence analysis of the proposed technique is discussed.The proposed method is generalized and extended to be of any odd order of the form 2n1.The scheme is composed of three steps,of which the first two steps are based on the two-step Homeier’s method with cubic convergence,and the last is a Newton step with an appropriate approximation for the derivative.Every iteration of the presented method requires the evaluation of two functions,two Fréchet derivatives,and three matrix inversions.A comparison between the efficiency index and the computational efficiency index of the presented scheme with existing methods is performed.The basins of attraction of the proposed scheme illustrated and compared to other schemes of the same order.Different test problems including large systems of equations are considered to compare the performance of the proposed method according to other methods of the same order.As an application,we apply the new scheme to some real-life problems,including the mixed Hammerstein integral equation and Burgers’equation.Comparisons and examples show that the presented method is efficient and comparable to the existing techniques of the same order.
文摘In this paper, by using the Nevanlinna Theory on angular domain, we establish a theorem which concerns the growth of entire function and his zero. As an application, we survey the location of zero of higher order differential equation, which can be regarded as an alternating but precise version of Wu and Yi.
文摘In this paper, based on the idea of El-Mistikawy and Werle[1] we construct a difference scheme for a singularly perturbed self-adjoint ordinary differential equation in conservation form. We prove that it is a uniformly convergent second order scheme.
基金Project supported by the Major State Basic Research Program of China (No. 19990328)the National Tackling Key Problems Program (No. 20050200069)+4 种基金the National Natural Science Foundation of China (Nos. 10771124, 10372052, 11101244, and 11271231)the Doctorate Foundation of the Ministry of Education of China (No. 20030422047)the Shandong Province Natural Science Foundation (No. ZR2009AQ012)the Independent Innovation Foundation of Shandong University(No. 2010TS031)the Scientific Research Award Fund for Excellent Middle-Aged and Young Scientists of Shandong Province (No. BS2009NJ003)
文摘A fractional step scheme with modified characteristic finite differences run- ning in a parallel arithmetic is presented to simulate a nonlinear percolation system of multilayer dynamics of fluids in a porous medium with moving boundary values. With the help of theoretical techniques including the change of regions, piecewise threefold quadratic interpolation, calculus of variations, multiplicative commutation rule of differ- ence operators, multiplicative commutation rule of difference operators, decomposition of high order difference operators, induction hypothesis, and prior estimates, an optimal order in 12 norm is displayed to complete the convergence analysis of the numerical algo- rithm. Some numerical results arising in the actual simulation of migration-accumulation of oil resources by this method are listed in the last section.
基金This work is supported by the Natural Science Foundation of Fujian Province of China(Grant No.2020J01783)the Project for High-Level Talent Innovation and Entrepreneurship of Quanzhou(Grant No.2018C087R)the Program for New Century Excellent Talents in Fujian Province University.
文摘The main purpose of this paper is to introduce some approximation properties of a Kantorovich kind q-Bernstein operators related to B′ezier basis functions with shape parameterλ∈[−1,1].Firstly,we compute some basic results such as moments and central moments,and derive the Korovkin type approximation theorem for these operators.Next,we estimate the order of convergence in terms of the usual modulus of continuity,for the functions belong to Lipschitz-type class and Peetre’s K-functional,respectively.Lastly,with the aid of Maple software,we present the comparison of the convergence of these newly defined operators to the certain function with some graphical illustrations and error estimation table.
文摘This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound corrector formulas for rectangular rule. Examples of numerical calculation have validated theoretical analysis.
文摘The truncation error of improved Cotes formula is presented in this paper. It also displays an analysis on convergence order of improved Cotes formula. Examples of numerical calculation is given in the end.
文摘In this paper, an interpolation polynomial operator F n(f; l,x) is constructed based on the zeros of a kind of Jacobi polynomials as the interpolation nodes. For any continuous function f(x)∈C b [-1,1] (0≤b≤l) F n(f; l,x) converges to f(x) uniformly, where l is an odd number.
文摘In this paper we discuss two-stage Miistein methods for solving Ito stochastic differential equations (SDEs). Six fully explicit methods (TSM 1 -- TSM 6) are given in this paper. Their order of strong convergence is proved. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of Ito SDEs.
文摘In this paper, the regularization with closed linear operators is used to solve an operator equation of the first kind. When all initial data are known approximately, we choose the regular parameter by using the general Arcangeli's criterion to give the convergence and the asymptotic orders of convergence of the regular solution.
文摘This paper presents a Modified Formula for Cotes rule with fifths derivatives of endpoint and its truncation error. It also displays an analysis on convergence order of compound formula. Though compound modified formula for Cotes rule with endpoint derivatives just calculates a newly-added fifths derivative of the two endpoints for each time compared with compound Cotes formula calculation, there are 2 more ranks of the convergence order in this modified formula. Examples of numerical calculation have validated theoretical analysis.
基金The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through the Large Groups Project under grant number RGP.2/235/43.
文摘A fifth-order family of an iterative method for solving systems of nonlinear equations and highly nonlinear boundary value problems has been developed in this paper.Convergence analysis demonstrates that the local order of convergence of the numerical method is five.The computer algebra system CAS-Maple,Mathematica,or MATLAB was the primary tool for dealing with difficult problems since it allows for the handling and manipulation of complex mathematical equations and other mathematical objects.Several numerical examples are provided to demonstrate the properties of the proposed rapidly convergent algorithms.A dynamic evaluation of the presented methods is also presented utilizing basins of attraction to analyze their convergence behavior.Aside from visualizing iterative processes,this methodology provides useful information on iterations,such as the number of diverging-converging points and the average number of iterations as a function of initial points.Solving numerous highly nonlinear boundary value problems and large nonlinear systems of equations of higher dimensions demonstrate the performance,efficiency,precision,and applicability of a newly presented technique.