665 results
Search Results
2. High-order linearly implicit exponential integrators conserving quadratic invariants with application to scalar auxiliary variable approach.
- Author
-
Sato, Shun
- Subjects
MATHEMATICAL analysis ,MATRIX multiplications ,ORDINARY differential equations ,QUADRATIC forms ,MATHEMATICS ,NUMERICAL integration - Abstract
This paper proposes a framework for constructing high-order linearly implicit exponential integrators that conserve a quadratic invariant. This is then applied to the scalar auxiliary variable (SAV) approach. Quadratic invariants are significant objects that are present in various physical equations and also in computationally efficient conservative schemes for general invariants. For instance, the SAV approach converts the invariant into a quadratic form by introducing scalar auxiliary variables, which have been intensively studied in recent years. In this vein, Sato et al. (Appl. Numer. Math. 187, 71-88 2023) proposed high-order linearly implicit schemes that conserve a quadratic invariant. In this study, it is shown that their method can be effectively merged with the Lawson transformation, a technique commonly utilized in the construction of exponential integrators. It is also demonstrated that combining the constructed exponential integrators and the SAV approach yields schemes that are computationally less expensive. Specifically, the main part of the computational cost is the product of several matrix exponentials and vectors, which are parallelizable. Moreover, we conduct some mathematical analyses on the proposed schemes. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
3. Convergence analysis of sample average approximation for a class of stochastic nonlinear complementarity problems: from two-stage to multistage
- Author
-
Hailin Sun, Bin Zhou, and Jie Jiang
- Subjects
Original Paper ,Class (set theory) ,Applied Mathematics ,Numerical analysis ,Two-stage ,010103 numerical & computational mathematics ,90C33 ,Multistage ,Lipschitz continuity ,Stochastic complementarity problems ,01 natural sciences ,90C15 ,010101 applied mathematics ,Convergence analysis ,Sample average approximation ,Theory of computation ,Convergence (routing) ,Nonlinear complementarity ,Applied mathematics ,Stage (hydrology) ,0101 mathematics ,Mathematics - Abstract
In this paper, we consider the sample average approximation (SAA) approach for a class of stochastic nonlinear complementarity problems (SNCPs) and study the corresponding convergence properties. We first investigate the convergence of the SAA counterparts of two-stage SNCPs when the first-stage problem is continuously differentiable and the second-stage problem is locally Lipschitz continuous. After that, we extend the convergence results to a class of multistage SNCPs whose decision variable of each stage is influenced only by the decision variables of adjacent stages. Finally, some preliminary numerical tests are presented to illustrate the convergence results.
- Published
- 2020
4. Numerical methods for static shallow shells lying over an obstacle
- Author
-
Paolo Piersanti, Xiaoqin Shen, City University of Hong Kong (CityU), and Xi'an Jiaotong University (Xjtu)
- Subjects
Original Paper ,Applied Mathematics ,Numerical analysis ,Mathematical analysis ,010103 numerical & computational mathematics ,Bilinear form ,Half-space ,Obstacle problems · Elliptic variational inequalities ,01 natural sciences ,Finite element method ,010101 applied mathematics ,Elliptic variational inequalities ,Non-conforming finite element method ,Enriching operator ,Obstacle problems ,Shallow shell ,Obstacle ,Theory of computation ,Convergence (routing) ,Obstacle problem ,0101 mathematics ,Nonconforming finite element method ,[MATH.MATH-NA]Mathematics [math]/Numerical Analysis [math.NA] ,Mathematics - Abstract
In this paper a finite element analysis to approximate the solution of an obstacle problem for a static shallow shell confined in a half space is presented. First, we rigorously prove some estimates for a suitable enriching operator connecting Morley's triangle to Hsieh-Clough-Tocher triangle. Secondly, we establish an estimate for the approximate bilinear form associated with the problem under consideration. Finally, we conduct an error analysis and we prove the convergence of the proposed numerical scheme.
- Published
- 2020
- Full Text
- View/download PDF
5. Waveform relaxation of partial differential equations
- Author
-
Zhen Miao and Yao-Lin Jiang
- Subjects
Partial differential equation ,Applied Mathematics ,Numerical analysis ,Short paper ,Relaxation (iterative method) ,010103 numerical & computational mathematics ,01 natural sciences ,010101 applied mathematics ,Theory of computation ,Convergence (routing) ,Applied mathematics ,Waveform ,0101 mathematics ,Energy (signal processing) ,Mathematics - Abstract
This short paper concludes a general waveform relaxation (WR) method at the PDE level for semi-linear reaction-diffusion equations. For the case of multiple coupled PDE(s), new Jacobi WR and Gauss-Seidel WR are provided to accelerate the convergence result of classical WR. The convergence conditions are proved based on energy estimate. Numerical experiments are demonstrated with several WR methods in parallel to verify the effectiveness of the general WR method.
- Published
- 2018
6. Nitsche's method for elliptic Dirichlet boundary control problems on curved domains.
- Author
-
Zhang, Qian and Zhang, Zhiyue
- Subjects
ENERGY policy ,A priori ,MATHEMATICS - Abstract
We consider Nitsche's method for solving elliptic Dirichlet boundary control problems on curved domains with control constraints. By using Nitsche's method for the treatment of inhomogeneous Dirichlet boundary conditions, the L
2 boundary control enters in the variational formulation in a natural sense. The idea was first used in Chang, et al. (Math. Anal. Appl.453, 529–557 2017) where the curved boundary was approximated by a broken line and a locally defined mapping was needed to obtain the numerical control on the curved boundary. In this paper, we develop a method defined on curved domains directly. We derive a priori estimates of quasi-optimal order for the control in the L2 norm, and quasi-optimal order for the state and adjoint state in energy norms. Numerical examples are provided to show the performance of the proposed method. [ABSTRACT FROM AUTHOR]- Published
- 2023
- Full Text
- View/download PDF
7. An improved tri-coloured rooted-tree theory and order conditions for ERKN methods for general multi-frequency oscillatory systems.
- Author
-
Zeng, Xianyang, Yang, Hongli, and Wu, Xinyuan
- Subjects
SPECTRUM analysis ,MATHEMATICS ,ALGEBRA ,GEOMETRY ,SYMMETRIES (Quantum mechanics) - Abstract
This paper develops an improved tri-coloured rooted-tree theory for the order conditions for ERKN methods solving general multi-frequency and multidimensional second-order oscillatory systems. The bottleneck of the original tricoloured rooted-tree theory is the existence of numerous redundant trees. In light of the fact that the sum of the products of the symmetries and the elementary differentials is meaningful, this paper naturally introduces the so-called extended elementary differential mappings. Then, the new improved tri-coloured rooted tree theory is established based on a subset of the original tri-coloured rooted-tree set. This new theory makes all redundant trees disappear, and thus, the order conditions of ERKN methods for general multi-frequency and multidimensional second-order oscillatory systems are reduced greatly. Furthermore, with this new theory, we present some new ERKN methods of order up to four. Numerical experiments are implemented and the results show that ERKN methods can be competitive with other existing methods in the scientific literature, especially when comparatively large stepsizes are used. [ABSTRACT FROM AUTHOR]
- Published
- 2017
- Full Text
- View/download PDF
8. Spurious eigenvalue-free algorithms of the method of fundamental solutions for solving the Helmholtz equation in bounded multiply connected domains.
- Author
-
Zhang, Li-Ping, Li, Zi-Cai, Wei, Yimin, and Huang, Hung-Tsai
- Subjects
HANKEL functions ,MATHEMATICS ,HELMHOLTZ equation - Abstract
For the Helmholtz equation Δu + k
2 u = 0 in 2D domain S, there exists a unique solution if k2 is not exactly equal to an eigenvalue λ of the Laplace eigenvalue problem Δu + λu = 0 in S. One important criterion for numerical methods is that there must exist no spurious (i.e., superfluous) eigenvalues so that the above unique solution can be obtained correctly. For exterior problems, the method of fundamental solutions (MFS) using Hankel functions suffers from spurious eigenvalues. New modified Hankel functions have been proposed in Zhang et al. (Appl. Numer. Math. 145, 236–260, 2019) to eliminate all spurious eigenvalues. In this paper, we study bounded multiply connected domains. The MFS algorithms without spurious eigenvalues and their strict analysis are our goals. First, we study bounded simply connected domains by the MFS. The algorithms using Hankel functions are free from spurious eigenvalues. A brief error analysis is provided. Next, we focus on bounded multiply connected domains and choose an annular domain for analysis. The Hankel functions and the modified Hankel functions in [47] are chosen as the exterior fundamental solutions (FS) and the interior FS, respectively. Such combined FS eliminate all spurious eigenvalues, and new error and stability analyses are explored. Bounds of errors involve degeneracy, and those of the condition number involve a gap to eigenvalues. Numerical experiments are carried out to support the analysis made, and better pseudo-boundaries of source nodes are also investigated numerically. The new error and stability analyses in this paper are new and essential, thus providing a solid theoretical basis of the MFS for the Helmholtz equation in 2D bounded simply connected and multiply connected domains. [ABSTRACT FROM AUTHOR]- Published
- 2022
- Full Text
- View/download PDF
9. The recursive quasi-orthogonal polynomial algorithm.
- Author
-
Sidki, S. and Sadaka, R.
- Subjects
POLYNOMIALS ,ALGORITHMS ,ORTHOGONAL polynomials ,SCHUR complement ,MATHEMATICS - Abstract
Theory of quasi-orthogonal polynomials is significantly related to constructing vector Padé approximations. The present paper introduces an efficient procedure to compute adjacent families of quasi-orthogonal polynomials as defined in Sadaka (Appl. Numer. Math. 21:57–70, 1996) and Sadaka (Appl. Numer. Math. 24:483–499, 1997). The derived algorithm uses short recursive relations whose coefficients are written in terms of lower triangular block determinants. The strategy of computing such determinants is given and analyzed. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
10. On equivalence of three-parameter iterative methods for singular symmetric saddle-point problem.
- Author
-
Hadjidimos, Apostolos and Tzoumas, Michael
- Subjects
MATHEMATICS - Abstract
There have been a couple of papers for the solution of the nonsingular symmetric saddle-point problem using three-parameter iterative methods. In most of them, regions of convergence for the parameters are found, while in three of them, optimal parameters are determined, and in one of the latter, many more cases, than in all the others, are distinguished, analyzed, and studied. It turns out that two of the optimal parameters coincide making the optimal three-parameter methods be equivalent to the optimal two-parameter known ones. Our aim in this work is manifold: (i) to show that the iterative methods we present are equivalent, (ii) to slightly change some statements in one of the main papers, (iii) to complete the analysis in another one, (iv) to explain how the transition from any of the methods to the others is made, (v) to extend the iterative method to cover the singular symmetric case, and (vi) to present a number of numerical examples in support of our theory. It would be an omission not to mention that the main material which all researchers in the area have inspired from and used is based on the one of the most cited papers by Bai et al. (Numer. Math. 102:1–38, 2005). [ABSTRACT FROM AUTHOR]
- Published
- 2021
- Full Text
- View/download PDF
11. Correction to: convergence rates for Kaczmarz-type algorithms.
- Author
-
Popa, Constantin
- Subjects
ALGORITHMS ,RATES ,MATHEMATICS ,EVIDENCE ,MATHEMATICAL equivalence - Abstract
We will refer to the paper "Convergence rates for Kaczmarz-type algorithms", published in Numerical Algorithms, 79(1)(2018), 1-17. It was observed by Dr. Mokhtar Abbasi (Department of Mathematics, University of Qom, Iran) that in the proof of Theorem 7, related to the linear convergence rate of the MREK algorithm it exists an inequality which is not always true. [ABSTRACT FROM AUTHOR]
- Published
- 2019
- Full Text
- View/download PDF
12. The use of the sinc-Gaussian sampling formula for approximating the derivatives of analytic functions.
- Author
-
Asharabi, Rashad M.
- Subjects
ANALYTIC functions ,ERROR analysis in mathematics ,SAMPLING errors ,INTEGRAL functions ,EXPONENTIAL functions ,MATHEMATICS - Abstract
The sinc-Gaussian sampling formula is used to approximate an analytic function, which satisfies a growth condition, using only finite samples of the function. The error of the sinc-Gaussian sampling formula decreases exponentially with respect to N, i.e., N
− 1/2 e−αN , where α is a positive number. In this paper, we extend this formula to allow the approximation of derivatives of any order of a function from two classes of analytic functions using only finite samples of the function itself. The theoretical error analysis is established based on a complex analytic approach; the convergence rate is also of exponential type. The estimate of Tanaka et al. (Jpan J. Ind. Appl. Math. 25, 209–231 2008) can be derived from ours as an immediate corollary. Various illustrative examples are presented, which show a good agreement with our theoretical analysis. [ABSTRACT FROM AUTHOR]- Published
- 2019
- Full Text
- View/download PDF
13. A new subtraction-free formula for lower bounds of the minimal singular value of an upper bidiagonal matrix.
- Author
-
Yamashita, Takumi, Kimura, Kinji, and Yamamoto, Yusaku
- Subjects
MATHEMATICAL bounds ,MATRICES (Mathematics) ,LINEAR algebra ,ALGEBRA ,MATHEMATICS - Abstract
Traces of inverse powers of a positive definite symmetric tridiagonal matrix give lower bounds of the minimal singular value of an upper bidiagonal matrix. In a preceding work, a formula for the traces which gives the diagonal entries of the inverse powers is presented. In this paper, we present another formula which gives the traces based on a quite different idea from the one in the preceding work. An efficient implementation of the formula for practice is also presented. [ABSTRACT FROM AUTHOR]
- Published
- 2015
- Full Text
- View/download PDF
14. Proximal point algorithms based on S-iterative technique for nearly asymptotically quasi-nonexpansive mappings and applications.
- Author
-
Sahu, D. R., Kumar, Ajeet, and Kang, Shin Min
- Subjects
NONEXPANSIVE mappings ,CONVEX sets ,ALGORITHMS ,POINT set theory ,MATHEMATICS - Abstract
In this paper, we combine the S-iteration process introduced by Agarwal et al. (J. Nonlinear Convex Anal., 8(1), 61–79 2007) with the proximal point algorithm introduced by Rockafellar (SIAM J. Control Optim., 14, 877–898 1976) to propose a new modified proximal point algorithm based on the S-type iteration process for approximating a common element of the set of solutions of convex minimization problems and the set of fixed points of nearly asymptotically quasi-nonexpansive mappings in the framework of CAT(0) spaces and prove the △-convergence of the proposed algorithm for solving common minimization problem and common fixed point problem. Our result generalizes, extends and unifies the corresponding results of Dhompongsa and Panyanak (Comput. Math. Appl., 56, 2572–2579 2008), Khan and Abbas (Comput. Math. Appl., 61, 109–116 2011), Abbas et al. (Math. Comput. Modelling, 55, 1418–1427 2012) and many more. [ABSTRACT FROM AUTHOR]
- Published
- 2021
- Full Text
- View/download PDF
15. A unified semilocal convergence analysis of a family of iterative algorithms for computing all zeros of a polynomial simultaneously.
- Author
-
Ivanov, Stoil
- Subjects
ALGORITHMIC randomness ,FOUNDATIONS of arithmetic ,ALGEBRA ,MATHEMATICS ,POLYNOMIALS - Abstract
In this paper, we first present a family of iterative algorithms for simultaneous determination of all zeros of a polynomial. This family contains two well-known algorithms: Dochev-Byrnev's method and Ehrlich's method. Second, using Proinov's approach to studying convergence of iterative methods for polynomial zeros, we provide a semilocal convergence theorem that unifies the results of Proinov (Appl. Math. Comput. 284: 102-114, 2016) for Dochev-Byrnev's and Ehrlich's methods. [ABSTRACT FROM AUTHOR]
- Published
- 2017
- Full Text
- View/download PDF
16. Optimal parameter of the SOR-like iteration method for solving absolute value equations.
- Author
-
Chen, Cairong, Huang, Bo, Yu, Dongmei, and Han, Deren
- Subjects
ABSOLUTE value ,EQUATIONS ,MATHEMATICS - Abstract
The SOR-like iteration method for solving the system of absolute value equations of finding a vector x such that A x - | x | - b = 0 with ν = ‖ A - 1 ‖ 2 < 1 is investigated. The convergence conditions of the SOR-like iteration method proposed by Ke and Ma (Appl. Math. Comput., 311:195–202, 2017) are revisited and a new proof is given, which exhibits some insights in determining the convergent region and the optimal iteration parameter. Along this line, the optimal parameter which minimizes ‖ T ν (ω) ‖ 2 with T ν (ω) = | 1 - ω | ω 2 ν | 1 - ω | | 1 - ω | + ω 2 ν and the approximate optimal parameter which minimizes an upper bound of ‖ T ν (ω) ‖ 2 are explored. The optimal and approximate optimal parameters are iteration-independent, and the bigger value of ν is, the smaller convergent region of the iteration parameter ω is. Numerical results are presented to demonstrate that the SOR-like iteration method with the optimal parameter is superior to that with the approximate optimal parameter proposed by Guo et al. (Appl. Math. Lett., 97:107–113, 2019). [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
17. Two non-parameter iterative algorithms for identifying strong ℋ-tensors.
- Author
-
Xu, Yangyang, Zhao, Ruijuan, and Zheng, Bing
- Subjects
HOMOGENEOUS polynomials ,ALGORITHMS ,MATHEMATICS - Abstract
The strong ℋ -tensors have important applications in many areas of science and engineering, e.g., the determination of positive definiteness for an even-order homogeneous polynomial form in the real field. In this paper, we propose two iterative algorithms with non-parameter for identifying strong ℋ -tensors, which overcome the drawback of choosing the best value of parameter 𝜖 in some existing algorithms given by Li et al. and Liu et al. (J. Comput. Appl. Math., 255, 1–14, 2014 and Comput. Appl. Math. 36, 1623–1635, 2017). Some numerical experiments are performed to illustrate the feasibility and effectiveness of our algorithms. [ABSTRACT FROM AUTHOR]
- Published
- 2019
- Full Text
- View/download PDF
18. A new integrable convergence acceleration algorithm for computing Brezinski-Durbin-Redivo-Zaglia’s sequence transformation via pfaffians.
- Author
-
Chang, Xiang-Ke, He, Yi, Hu, Xing-Biao, and Li, Shi-Hao
- Subjects
ALGORITHMS ,ALGEBRA ,INTEGERS ,MATHEMATICS ,APPROXIMATION theory - Abstract
In the literature, most known sequence transformations can be written as a ratio of two determinants. But, it is not always this case. One exception is that the sequence transformation proposed by Brezinski, Durbin, and Redivo-Zaglia cannot be expressed as a ratio of two determinants. Motivated by this, we will introduce a new algebraic tool—pfaffians, instead of determinants in the paper. It turns out that Brezinski-Durbin-Redivo-Zaglia’s transformation can be expressed as a ratio of two pfaffians. To the best of our knowledge, this is the first time to introduce pfaffians in the expressions of sequence transformations. Furthermore, an extended transformation of high order is presented in terms of pfaffians and a new convergence acceleration algorithm for implementing the transformation is constructed. Then, the Lax pair of the recursive algorithm is obtained which implies that the algorithm is integrable. Numerical examples with applications of the algorithm are also presented. [ABSTRACT FROM AUTHOR]
- Published
- 2018
- Full Text
- View/download PDF
19. A generalized variant of the deteriorated PSS preconditioner for nonsymmetric saddle point problems.
- Author
-
Huang, Zheng-Ge, Wang, Li-Gong, Xu, Zhong, and Cui, Jing-Jing
- Subjects
MATHEMATICS ,EIGENVECTORS ,MATRICES (Mathematics) ,VECTOR spaces ,MACHINE theory - Abstract
Based on the variant of the deteriorated positive-definite and skew-Hermitian splitting (VDPSS) preconditioner developed by Zhang and Gu (BIT Numer. Math. 56:587-604, 2016), a generalized VDPSS (GVDPSS) preconditioner is established in this paper by replacing the parameter α in (2,2)-block of the VDPSS preconditioner by another parameter β. This preconditioner can also be viewed as a generalized form of the VDPSS preconditioner and the new relaxed HSS (NRHSS) preconditioner which has been exhibited by Salkuyeh and Masoudi (Numer. Algorithms, 2016). The convergence properties of the GVDPSS iteration method are derived. Meanwhile, the distribution of eigenvalues and the forms of the eigenvectors of the preconditioned matrix are analyzed in detail. We also study the upper bounds on the degree of the minimum polynomial of the preconditioned matrix. Numerical experiments are implemented to illustrate the effectiveness of the GVDPSS preconditioner and verify that the GVDPSS preconditioned generalized minimal residual method is superior to the DPSS, relaxed DPSS, SIMPLE-like, NRHSS, and VDPSS preconditioned ones for solving saddle point problems in terms of the iterations and computational times. [ABSTRACT FROM AUTHOR]
- Published
- 2017
- Full Text
- View/download PDF
20. Two formulae with nodes related to zeros of Bessel functions for semi-infinite integrals: extending Gauss–Jacobi-type rules.
- Author
-
Sugiura, Hiroshi and Hasegawa, Takemitsu
- Subjects
INTEGRAL functions ,INTEGRALS ,MATHEMATICS - Abstract
For approximating integrals ∫ 0 ∞ x α f (x) d x ( α > - 1 ) over a semi-infinite interval [ 0 , ∞) with a given function f (x) , two formulae, one of them new and another associated with an existing formula, are presented. They are constructed in a limiting process to a semi-infinite interval [ 0 , ∞ ] with a linear transformation for a well-known approximation method, the Gauss–Jacobi (GJ) rule and its family rules: the Gauss–Jacobi–Radau (GR) and Gauss–Jacobi–Lobatto (GL) rules on a finite interval. This procedure was used in constructing our previous limit Clenshaw–Curtis-type formulae. The limit GJ formula (LGJ) constructed in this way uses as nodes zeros of the Bessel function J α (x) squared after multiplied by a positive constant a and the limit GR (LGR) (and limit GL (LGL)) formula those with zeros of J α + 1 (x) . The LGJ formula is also shown to be derived from the formula developed by Frappier and Olivier (Math. Comp. 60:303–316, 1993) for an integral on [ 0 , ∞) . The LGR and LGL formulae give the same and new formula. We show that for a function f(z) analytic on a domain in the complex plane z and satisfying some appropriate conditions, there exists a constant d > 0 such that the errors of both formulae are O (e - 2 d / a) as a → + 0 . The average of the LGJ and LGR formulae gives smaller quadrature errors than each formula. Numerical examples confirm these behaviors and show that the LGJ and LGR formulae give asymptotically the same quadrature errors of opposite sign. Consequently, the LGR formula behaves like an anti-LGJ formula in the same way as the Lobatto rule for integrals on [ - 1 , 1 ] behaving like the anti-Gauss rule. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
21. Riemannian linearized proximal algorithms for nonnegative inverse eigenvalue problem.
- Author
-
Kum, Sangho, Li, Chong, Wang, Jinhua, Yao, Jen-Chih, and Zhu, Linglingzhi
- Subjects
INVERSE problems ,RIEMANNIAN manifolds ,SPARSE matrices ,ALGORITHMS ,MATHEMATICS - Abstract
We study the issue of numerically solving the nonnegative inverse eigenvalue problem (NIEP). At first, we reformulate the NIEP as a convex composite optimization problem on Riemannian manifolds. Then we develop a scheme of the Riemannian linearized proximal algorithm (R-LPA) to solve the NIEP. Under some mild conditions, the local and global convergence results of the R-LPA for the NIEP are established, respectively. Moreover, numerical experiments are presented. Compared with the Riemannian Newton-CG method in Z. Zhao et al. (Numer. Math. 140:827–855, 2018), this R-LPA owns better numerical performances for large scale problems and sparse matrix cases, which is due to the smaller dimension of the Riemannian manifold derived from the problem formulation of the NIEP as a convex composite optimization problem. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
22. Some results on certain generalized circulant matrices.
- Author
-
Lu, Chengbo
- Subjects
TOEPLITZ matrices ,CIRCULANT matrices ,EIGENVALUES ,INTEGERS ,MATHEMATICS ,POLYNOMIAL time algorithms - Abstract
In this paper a particular partition on blocks of generalized ( h, r)-circulant matrices is determined. We obtain a characterization of generalized ( h, r)-circulant matrices and get some results on the values of the permanent and also on the determination of the eigenvalues of r-circulant matrices. At last, a lower bound for the permanent of these matrices is achieved. [ABSTRACT FROM AUTHOR]
- Published
- 2015
- Full Text
- View/download PDF
23. Convergence for a class of multi-point modified Chebyshev-Halley methods under the relaxed conditions.
- Author
-
Wang, Xiuhua and Kou, Jisheng
- Subjects
STOCHASTIC convergence ,BANACH spaces ,VECTOR spaces ,INTEGRAL equations ,MATHEMATICS ,MATHEMATICS theorems - Abstract
In this paper, the semilocal convergence for a class of multi-point modified Chebyshev-Halley methods in Banach spaces is studied. Different from the results in reference [11], these methods are more general and the convergence conditions are also relaxed. We derive a system of recurrence relations for these methods and based on this, we prove a convergence theorem to show the existence-uniqueness of the solution. A priori error bounds is also given. The R-order of these methods is proved to be 5+ q with ω−conditioned third-order Fréchet derivative, where ω( μ) is a non-decreasing continuous real function for μ > 0 and satisfies ω(0) ≥ 0, ω( tμ) ≤ t ω( μ) for μ > 0, t ∈ [0,1] and q ∈ [0,1]. Finally, we give some numerical results to show our approach. [ABSTRACT FROM AUTHOR]
- Published
- 2015
- Full Text
- View/download PDF
24. Optimal convergence of three iterative methods based on nonconforming finite element discretization for 2D/3D MHD equations
- Author
-
Haiyan Su, Zhilin Li, and Jiali Xu
- Subjects
Piecewise linear function ,Nonlinear system ,Discretization ,Iterative method ,Applied Mathematics ,Numerical analysis ,Convergence (routing) ,Applied mathematics ,Stability (probability) ,Finite element method ,Mathematics::Numerical Analysis ,Mathematics - Abstract
The main purpose of this paper is to analyze nonconforming iterative finite element methods for 2D/3D stationary incompressible magneto-hydrodynamics equations. First, the Crouzeix-Raviart–type finite element is used to approximate the velocity and the conforming piecewise linear element P1 is used for the pressure. Since the finite element method for the velocity field and the pressure is unstable, a simple locally stabilization term is added to satisfy the weak inf-sup condition. Then, the well-posedness and the optimal error estimates of the continuous and discrete problems are analyzed with the nonlinear terms being iteratively updated. Three effective iterative methods are proposed and their stability and convergence analyses are carried out. Finally, the theoretical analysis presented in this paper is verified by numerical experiments.
- Published
- 2021
25. On the best achievable quality of limit points of augmented Lagrangian schemes
- Author
-
Gabriel Haeser, Roberto Andreani, Leonardo D. Secchin, Alberto Ramos, and Leonardo M. Mito
- Subjects
Constraint (information theory) ,Mathematical optimization ,Augmented Lagrangian method ,Applied Mathematics ,Numerical analysis ,media_common.quotation_subject ,Theory of computation ,Convergence (routing) ,Limit point ,Quality (business) ,Algebra over a field ,Mathematics ,media_common - Abstract
The optimization literature is vast in papers dealing with improvements on the global convergence of augmented Lagrangian schemes. Usually, the results are based on weak constraint qualifications, or, more recently, on sequential optimality conditions obtained via penalization techniques. In this paper, we propose a somewhat different approach, in the sense that the algorithm itself is used in order to formulate a new optimality condition satisfied by its feasible limit points. With this tool at hand, we present several new properties and insights on limit points of augmented Lagrangian schemes, in particular, characterizing the strongest possible global convergence result for the safeguarded augmented Lagrangian method.
- Published
- 2021
26. A Tseng extragradient method for solving variational inequality problems in Banach spaces
- Author
-
H. A. Abass, A. A. Mebawondu, Olawale Kazeem Oyewole, and K. O. Aremu
- Subjects
Sequence ,Applied Mathematics ,Numerical analysis ,Variational inequality ,Convergence (routing) ,Banach space ,Applied mathematics ,Lipschitz continuity ,Constant (mathematics) ,Projection (linear algebra) ,Mathematics - Abstract
This paper presents an inertial Tseng extragradient method for approximating a solution of the variational inequality problem. The proposed method uses a single projection onto a half space which can be easily evaluated. The method considered in this paper does not require the knowledge of the Lipschitz constant as it uses variable stepsizes from step to step which are updated over each iteration by a simple calculation. We prove a strong convergence theorem of the sequence generated by this method to a solution of the variational inequality problem in the framework of a 2-uniformly convex Banach space which is also uniformly smooth. Furthermore, we report some numerical experiments to illustrate the performance of this method. Our result extends and unifies corresponding results in this direction in the literature.
- Published
- 2021
27. A matrix-less method to approximate the spectrum and the spectral function of Toeplitz matrices with real eigenvalues
- Author
-
Sven-Erik Ekström and P. Vassalos
- Subjects
Beräkningsmatematik ,Applied Mathematics ,010102 general mathematics ,Generating function ,Order (ring theory) ,Asymptotic expansion ,Spectral analysis ,Numerical Analysis (math.NA) ,010103 numerical & computational mathematics ,Function (mathematics) ,Type (model theory) ,01 natural sciences ,Toeplitz matrix ,Combinatorics ,Computational Mathematics ,Matrix (mathematics) ,Toeplitz matrices ,FOS: Mathematics ,Mathematics - Numerical Analysis ,0101 mathematics ,Structured matrices ,Eigenvalues and eigenvectors ,Mathematics - Abstract
It is known that the generating function f of a sequence of Toeplitz matrices {Tn(f)}n may not describe the asymptotic distribution of the eigenvalues of Tn(f) if f is not real. In this paper, we assume as a working hypothesis that, if the eigenvalues of Tn(f) are real for all n, then they admit an asymptotic expansion of the same type as considered in previous works, where the first function, called the eigenvalue symbol $\mathfrak {f}$ f , appearing in this expansion is real and describes the asymptotic distribution of the eigenvalues of Tn(f). This eigenvalue symbol $\mathfrak {f}$ f is in general not known in closed form. After validating this working hypothesis through a number of numerical experiments, we propose a matrix-less algorithm in order to approximate the eigenvalue distribution function $\mathfrak {f}$ f . The proposed algorithm, which opposed to previous versions, does not need any information about neither f nor $\mathfrak {f}$ f is tested on a wide range of numerical examples; in some cases, we are even able to find the analytical expression of $\mathfrak {f}$ f . Future research directions are outlined at the end of the paper.
- Published
- 2021
28. Stability analysis of the method of fundamental solutions with smooth closed pseudo-boundaries for Laplace’s equation: better pseudo-boundaries
- Author
-
Li-Ping Zhang, Zi-Cai Li, Ming-Gong Lee, and Hung-Tsai Huang
- Subjects
Laplace's equation ,Polynomial ,Laplace transform ,Applied Mathematics ,Numerical analysis ,010103 numerical & computational mathematics ,01 natural sciences ,010101 applied mathematics ,Bounded function ,Applied mathematics ,Method of fundamental solutions ,0101 mathematics ,Condition number ,Circulant matrix ,Mathematics - Abstract
Consider Laplace’s equation in a bounded simply-connected domain S, and use the method of fundamental solutions (MFS). The error and stability analysis is made for circular/elliptic pseudo-boundaries in Dou et al. (J. Comp. Appl. Math. 377:112861, 2020), and polynomial convergence rates and exponential growth rates of the condition number (Cond) are obtained. General pseudo-boundaries are suggested for more complicated solution domains in Dou et al. (J. Comp. Appl. Math. 377:112861, 2020, Section 5). Since the ill-conditioning is severe, the success in computation by the MFS mainly depends on stability. This paper is devoted to stability analysis for smooth closed pseudo-boundaries of source nodes. Bounds of the Cond are derived, and exponential growth rates are also obtained. This paper is the first time to explore stability analysis of the MFS for non-circular/non-elliptic pseudo-boundaries. Circulant matrices are often employed for stability analysis of the MFS; but the stability analysis in this paper is explored based on new techniques without using circulant matrices as in Dou et al. (J. Comp. Appl. Math. 377:112861, 2020). To pursue better pseudo-boundaries, the sensitivity index is proposed from growth/convergence rates of stability via accuracy. Better pseudo-boundaries in the MFS can be found by trial computations, to develop the study in Dou et al. (J. Comp. Appl. Math. 377:112861, 2020) for the selection of pseudo-boundaries. For highly smooth and singular solutions, better pseudo-boundaries are different; an analysis of the sensitivity index is explored. Circular/elliptic pseudo-boundaries are optimal for highly smooth solutions, but not for singular solutions. In this paper, amoeba-like domains are chosen in computation. Several useful types of pseudo-boundaries are developed and their algorithms are simple without using nonlinear solutions. For singular solutions, numerical comparisons are made for different pseudo-boundaries via the sensitivity index.
- Published
- 2021
29. Numerical study on Moore-Penrose inverse of tensors via Einstein product
- Author
-
Baohua Huang
- Subjects
Applied Mathematics ,Numerical analysis ,Inverse ,Space (mathematics) ,General Relativity and Quantum Cosmology ,symbols.namesake ,Product (mathematics) ,Conjugate gradient method ,symbols ,Tensor ,Einstein ,Moore–Penrose pseudoinverse ,Mathematical physics ,Mathematics - Abstract
The notation of Moore-Penrose inverse of matrices has been extended from matrix space to even-order tensor space with Einstein product. In this paper, we give the numerical study on the Moore-Penrose inverse of tensors via the Einstein product. More precisely, we transform the calculation of Moore-Penrose inverse of tensors via the Einstein product into solving a class of tensor equations via the Einstein product. Then, by means of the conjugate gradient method, we obtain the approximate Moore-Penrose inverse of tensors via the Einstein product. Finally, we report some numerical examples to show the efficiency of the proposed methods and testify the conclusion suggested in this paper.
- Published
- 2021
30. Weak Galerkin finite element methods with or without stabilizers
- Author
-
Xiaoshen Wang, Xiu Ye, and Shangyou Zhang
- Subjects
010101 applied mathematics ,Applied Mathematics ,Numerical analysis ,Convergence (routing) ,Theory of computation ,Applied mathematics ,010103 numerical & computational mathematics ,0101 mathematics ,Algebra over a field ,Galerkin method ,01 natural sciences ,Finite element method ,Mathematics - Abstract
The purpose of this paper is to investigate the connections between the weak Galerkin (WG) methods with and without stabilizers. The choices of stabilizers directly affect the convergence rates of the corresponding WG methods in general. However, we observed that the convergence rates are independent of the choices of stabilizers for these WG elements with stabilizers being optional. In this paper, we will verify such phenomena theoretically as well as numerically.
- Published
- 2021
31. Highly efficient schemes for time-fractional Allen-Cahn equation using extended SAV approach
- Author
-
Chuanju Xu, Hongyi Zhu, Dianming Hou, Institut de Mécanique et d'Ingénierie (I2M), Université de Bordeaux (UB)-Institut Polytechnique de Bordeaux-Centre National de la Recherche Scientifique (CNRS)-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE)-Arts et Métiers Sciences et Technologies, and HESAM Université (HESAM)-HESAM Université (HESAM)
- Subjects
Discretization ,Applied Mathematics ,Numerical analysis ,Scalar (physics) ,Numerical Analysis (math.NA) ,010103 numerical & computational mathematics ,01 natural sciences ,Stability (probability) ,[SPI.MAT]Engineering Sciences [physics]/Materials ,010101 applied mathematics ,Nonlinear system ,Theory of computation ,FOS: Mathematics ,Order (group theory) ,Applied mathematics ,Mathematics - Numerical Analysis ,0101 mathematics ,Allen–Cahn equation ,Mathematics - Abstract
In this paper, we propose and analyze high-order efficient schemes for the time-fractional Allen-Cahn equation. The proposed schemes are based on the L1 discretization for the time-fractional derivative and the extended scalar auxiliary variable (SAV) approach developed very recently to deal with the nonlinear terms in the equation. The main contributions of the paper consist of (1) constructing first- and higher order unconditionally stable schemes for different mesh types, and proving the unconditional stability of the constructed schemes for the uniform mesh; (2) carrying out numerical experiments to verify the efficiency of the schemes and to investigate the coarsening dynamics governed by the time-fractional Allen-Cahn equation. In particular, the influence of the fractional order on the coarsening behavior is carefully examined. Our numerical evidence shows that the proposed schemes are more robust than the existing methods, and their efficiency is less restricted to particular forms of the nonlinear potentials.
- Published
- 2021
32. Centrality measures for node-weighted networks via line graphs and the matrix exponential
- Author
-
Mona Matar, Lothar Reichel, and Omar De la Cruz Cabrera
- Subjects
Discrete mathematics ,Applied Mathematics ,Node (networking) ,Directed graph ,law.invention ,law ,Matrix function ,Theory of computation ,Line graph ,Adjacency matrix ,Matrix exponential ,Centrality ,MathematicsofComputing_DISCRETEMATHEMATICS ,Mathematics - Abstract
This paper is concerned with the identification of important nodes in node-weighted graphs by applying matrix functions, in particular the matrix exponential. Many tools that use an adjacency matrix for a graph have been developed to study the importance of the nodes in unweighted or edge-weighted networks. However, adjacency matrices for node-weighted graphs have not received much attention. The present paper proposes using a line graph associated with a node-weighted graph to construct an edge-weighted graph that can be analyzed with available methods. Both undirected and directed graphs with positive node weights are considered. We show that when the weight of a node increases, the importance of this node in the graph increases as well, provided that the adjacency matrix is suitably scaled. Applications to real-life problems are presented.
- Published
- 2021
33. Tensor extrapolation methods with applications
- Author
-
Rachid Sadaka, Khalide Jbilou, Fatemeh Panjeh Ali Beik, and A. El Ichi
- Subjects
Sequence ,Applied Mathematics ,Numerical analysis ,Theory of computation ,Singular value decomposition ,Extrapolation ,Applied mathematics ,Tensor ,Algebra over a field ,Mathematics ,Matrix polynomial - Abstract
In this paper, we mainly develop the well-known vector and matrix polynomial extrapolation methods in tensor framework. To this end, some new products between tensors are defined and the concept of positive definitiveness is extended for tensors corresponding to T-product. Furthermore, we discuss on the solution of least-squares problem associated with a tensor equation using Tensor Singular Value Decomposition (TSVD). Motivated by the effectiveness of some proposed vector extrapolation methods in earlier papers, we describe how an extrapolation technique can be also implemented on the sequence of tensors produced by truncated TSVD (TTSVD) for solving possibly ill-posed tensor equations.
- Published
- 2020
34. A nonnegativity preserving algorithm for multilinear systems with nonsingular ${\mathcal M}$-tensors
- Author
-
Hongjin He, Guanglu Zhou, Chen Ling, and Xueli Bai
- Subjects
Sequence ,Multilinear map ,Applied Mathematics ,Numerical analysis ,MathematicsofComputing_NUMERICALANALYSIS ,010103 numerical & computational mathematics ,System of linear equations ,01 natural sciences ,law.invention ,010101 applied mathematics ,Invertible matrix ,law ,TheoryofComputation_ANALYSISOFALGORITHMSANDPROBLEMCOMPLEXITY ,ComputingMethodologies_SYMBOLICANDALGEBRAICMANIPULATION ,Theory of computation ,Tensor ,0101 mathematics ,Algorithm ,Mathematics ,Numerical partial differential equations - Abstract
This paper addresses multilinear systems of equations which arise in various applications such as data mining and numerical partial differential equations. When the multilinear system under consideration involves a nonsingular ${\mathscr{M}}$ -tensor and a nonnegative right-hand side vector, it may have multiple nonnegative solutions. In this paper, we propose an algorithm which can always preserve the nonnegativity of solutions. Theoretically, we show that the sequence generated by the proposed algorithm is a nonnegative componentwise nonincreasing sequence and converges to a nonnegative solution of the system. Numerical results further support the novelty of the proposed method.
- Published
- 2020
35. On the split common fixed point problem for strict quasi-ϕ-pseudocontractive mappings in Banach spaces
- Author
-
Xindong Liu, Zili Chen, and Jinxing Liu
- Subjects
Sequence ,Pure mathematics ,Iterative method ,Applied Mathematics ,Numerical analysis ,Theory of computation ,Banach space ,Common fixed point ,Null point ,Algebra over a field ,Mathematics - Abstract
The purpose of this paper is to propose an algorithm for solving the split common fixed point problem for strict quasi-ϕ-pseudocontractive mappings in Banach spaces. It is proved that the sequence generated by the proposed iterative algorithm converges strongly to a solution of the split common fixed point problem. Then, the main result is used to study the split common null point problem and the split quasi-inclusion problem. Finally, a numerical example is provided to illustrate our main result. The results presented in this paper extend and improve some recent corresponding results.
- Published
- 2020
36. An implicit-explicit preconditioned direct method for pricing options under regime-switching tempered fractional partial differential models
- Author
-
Xu Chen, Wenfei Wang, Siu-Long Lei, and Deng Ding
- Subjects
Partial differential equation ,Direct sum ,Applied Mathematics ,Direct method ,Numerical analysis ,010103 numerical & computational mathematics ,01 natural sciences ,Toeplitz matrix ,010101 applied mathematics ,Valuation of options ,Applied mathematics ,Partial derivative ,0101 mathematics ,Coefficient matrix ,Mathematics - Abstract
Recently, fractional partial differential equations have been widely applied in option pricing problems, which better explains many important empirical facts of financial markets, but rare paper considers the multi-state options pricing problem based on fractional diffusion models. Thus, multi-state European option pricing problem under regime-switching tempered fractional partial differential equation is considered in this paper. Due to the expensive computational cost caused by the implicit finite difference scheme, a novel implicit-explicit finite difference scheme has been developed with consistency, stability, and convergence guarantee. Since the resulting coefficient matrix equals to the direct sum of several Toeplitz matrices, a preconditioned direct method has been proposed with ${\mathcal O}(\bar {S}N\log N+\bar {S}^{2} N)$ operation cost on each time level with adaptability analysis, where $\bar {S}$ is the number of states and N is the number of grid points. Related numerical experiments including an empirical example have been presented to demonstrate the effectiveness and accuracy of the proposed numerical method.
- Published
- 2020
37. Analysis of optimal superconvergence of the local discontinuous Galerkin method for nonlinear fourth-order boundary value problems
- Author
-
Mahboub Baccouch
- Subjects
Degree (graph theory) ,Discontinuous Galerkin method ,Applied Mathematics ,Order (ring theory) ,Applied mathematics ,Boundary value problem ,Superconvergence ,Type (model theory) ,Finite element method ,Prime (order theory) ,Mathematics::Numerical Analysis ,Mathematics - Abstract
This paper is concerned with the convergence and superconvergence of the local discontinuous Galerkin (LDG) finite element method for nonlinear fourth-order boundary value problems of the type $u^{(4)}=f(x,u,u^{\prime },u^{\prime \prime },u^{\prime \prime \prime })$ , x ∈ [a,b] with classical boundary conditions at the endpoints. Convergence properties for the solution and for all three auxiliary variables are established. More specifically, we use the duality argument to prove that the errors between the LDG solutions and the exact solutions in the L2 norm achieve optimal (p + 1)th-order convergence, when polynomials of degree p are used. We also prove that the derivatives of the errors between the LDG solutions and Gauss-Radau projections of the exact solutions in the L2 norm are superconvergent with order p + 1. Furthermore, a (2p + 1)th-order superconvergent for the errors of the numerical fluxes at mesh nodes as well as for the cell averages is also obtained under quasi-uniform meshes. Finally, we prove that the LDG solutions are superconvergent with an order of p + 2 toward particular projections of the exact solutions. The error analysis presented in this paper is valid for p ≥ 1. Numerical experiments indicate that our theoretical findings are optimal.
- Published
- 2020
38. Continuous and discrete zeroing dynamics models using JMP function array and design formula for solving time-varying Sylvester-transpose matrix inequality
- Author
-
Huanchang Huang, Xiao Liu, Min Yang, Yunong Zhang, and Yihong Ling
- Subjects
Discretization ,Truncation error (numerical integration) ,Applied Mathematics ,Numerical analysis ,010103 numerical & computational mathematics ,Function (mathematics) ,01 natural sciences ,010101 applied mathematics ,Matrix (mathematics) ,Transpose ,Convergence (routing) ,Applied mathematics ,0101 mathematics ,Mathematics ,Variable (mathematics) - Abstract
Zeroing dynamics (ZD) has shown great performance to solve various time-varying problems. In this paper, the problem of time-varying Sylvester-transpose matrix inequality is first investigated. Since it is difficult to solve a matrix inequality with a matrix variable and its transpose by traditional methods, this paper proposes a continuous ZD (CZD) model by employing ZD design formula and JMP function array to solve this challenging problem. Furthermore, for better implementation on digital computers, three discrete ZD (DZD) models are proposed by using three different discretization formulas with different precision, i.e., the Euler-forward formula, the 6-instant Zhang et al discretization (ZeaD) formula and the 7-instant ZeaD formula. What is more, theoretical truncation error analyses and numerical experiments substantiate the convergence, efficacy and superiority of the DZD models for solving time-varying Sylvester-transpose matrix inequality.
- Published
- 2020
39. Newton’s method with fractional derivatives and various iteration processes via visual analysis
- Author
-
Krzysztof Gdawiec, Agnieszka Lisowska, and Wiesław Kotarski
- Subjects
Polynomial ,Applied Mathematics ,Numerical analysis ,Stability (learning theory) ,Fractional derivative ,01 natural sciences ,Fractional calculus ,010101 applied mathematics ,symbols.namesake ,Newton method ,Fixed-point iteration ,ComputingMethodologies_SYMBOLICANDALGEBRAICMANIPULATION ,0103 physical sciences ,Convergence (routing) ,symbols ,Polynomiography ,Applied mathematics ,Iterations ,0101 mathematics ,010301 acoustics ,Complex plane ,Newton's method ,Mathematics - Abstract
The aim of this paper is to visually investigate the dynamics and stability of the process in which the classic derivative is replaced by the fractional Riemann–Liouville or Caputo derivatives in the standard Newton root-finding method. Additionally, instead of the standard Picard iteration, the Mann, Khan, Ishikawa and S iterations are used. This process when applied to polynomials on complex plane produces images showing basins of attractions for polynomial zeros or images representing the number of iterations required to achieve any polynomial root. The images are called polynomiographs. In this paper, we use the colouring according to the number of iterations which reveals the speed of convergence and dynamic properties of processes visualised by polynomiographs. Moreover, to investigate the stability of the methods, we use basins of attraction. To compare numerically the modified root-finding methods among them, we demonstrate their action for polynomialz3− 1 on a complex plane.
- Published
- 2020
40. On equivalence of three-parameter iterative methods for singular symmetric saddle-point problem
- Author
-
M. Tzoumas and Apostolos Hadjidimos
- Subjects
Iterative method ,Applied Mathematics ,Numerical analysis ,010103 numerical & computational mathematics ,Symmetric case ,01 natural sciences ,law.invention ,010101 applied mathematics ,Invertible matrix ,law ,Saddle point ,Theory of computation ,Applied mathematics ,0101 mathematics ,Equivalence (measure theory) ,Saddle ,Mathematics - Abstract
There have been a couple of papers for the solution of the nonsingular symmetric saddle-point problem using three-parameter iterative methods. In most of them, regions of convergence for the parameters are found, while in three of them, optimal parameters are determined, and in one of the latter, many more cases, than in all the others, are distinguished, analyzed, and studied. It turns out that two of the optimal parameters coincide making the optimal three-parameter methods be equivalent to the optimal two-parameter known ones. Our aim in this work is manifold: (i) to show that the iterative methods we present are equivalent, (ii) to slightly change some statements in one of the main papers, (iii) to complete the analysis in another one, (iv) to explain how the transition from any of the methods to the others is made, (v) to extend the iterative method to cover the singular symmetric case, and (vi) to present a number of numerical examples in support of our theory. It would be an omission not to mention that the main material which all researchers in the area have inspired from and used is based on the one of the most cited papers by Bai et al. (Numer. Math. 102:1–38, 2005).
- Published
- 2020
41. Halpern-type iterative process for solving split common fixed point and monotone variational inclusion problem between Banach spaces
- Author
-
Oluwatosin Temitope Mewomo, Timilehin Opeyemi Alakoya, and Adeolu Taiwo
- Subjects
Iterative and incremental development ,Applied Mathematics ,Regular polygon ,Banach space ,010103 numerical & computational mathematics ,Type (model theory) ,01 natural sciences ,010101 applied mathematics ,Monotone polygon ,Corollary ,Theory of computation ,Convergence (routing) ,Applied mathematics ,0101 mathematics ,Mathematics - Abstract
In this paper, we study the split common fixed point and monotone variational inclusion problem in uniformly convex and 2-uniformly smooth Banach spaces. We propose a Halpern-type algorithm with two self-adaptive stepsizes for obtaining solution of the problem and prove strong convergence theorem for the algorithm. Many existing results in literature are derived as corollary to our main result. In addition, we apply our main result to split common minimization problem and fixed point problem and illustrate the efficiency and performance of our algorithm with a numerical example. The main result in this paper extends and generalizes many recent related results in the literature in this direction.
- Published
- 2020
42. Superconvergence in H1-norm of a difference finite element method for the heat equation in a 3D spatial domain with almost-uniform mesh
- Author
-
Ruijian He, Zhangxin Chen, and Xinlong Feng
- Subjects
Backward differentiation formula ,Computational complexity theory ,Applied Mathematics ,Numerical analysis ,Mathematical analysis ,010103 numerical & computational mathematics ,Superconvergence ,01 natural sciences ,Finite element method ,010101 applied mathematics ,Norm (mathematics) ,Bounded function ,Heat equation ,0101 mathematics ,Mathematics - Abstract
In this paper, we propose a novel difference finite element (DFE) method based on the P1-element for the 3D heat equation on a 3D bounded domain. One of the novel ideas of this paper is to use the second-order backward difference formula (BDF) combining DFE method to overcome the computational complexity of conventional finite element (FE) method for the high-dimensional parabolic problem. First, we design a fully discrete difference FE solution ${u^{n}_{h}}$ by the second-order backward difference formula in the temporal t-direction, the center difference scheme in the spatial z-direction, and the P1-element on a almost-uniform mesh Jh in the spatial (x, y)-direction. Next, the H1-stability of ${u_{h}^{n}}$ and the second-order H1-convergence of the interpolation post-processing function on ${u_{h}^{n}}$ with respect to u(tn) are provided. Finally, numerical tests are presented to show the second-order H1-convergence results of the proposed DFE method for the heat equation in a 3D spatial domain.
- Published
- 2020
43. Analysis and application of the interpolating element-free Galerkin method for extended Fisher–Kolmogorov equation which arises in brain tumor dynamics modeling
- Author
-
Mohammad Ilati
- Subjects
Partial differential equation ,Applied Mathematics ,Numerical analysis ,Finite difference ,010103 numerical & computational mathematics ,01 natural sciences ,010101 applied mathematics ,Nonlinear system ,Convergence (routing) ,Biharmonic equation ,Fisher–Kolmogorov equation ,Applied mathematics ,0101 mathematics ,Galerkin method ,Mathematics - Abstract
In this paper, the interpolating element-free Galerkin method is applied for solving the nonlinear biharmonic extended Fisher–Kolmogorov equation which arises in brain tumor dynamics modeling. At first, a finite difference formula is utilized for obtaining a time-discrete scheme. The unconditional stability and convergence of the time-discrete method are proved by the energy method. Then, we use the interpolating element-free Galerkin method to approximate the spatial derivatives. An error analysis of the interpolating element-free Galerkin method is proposed for this nonlinear equation. Moreover, this method is compared with some other meshless local weak-form techniques. The main aim of this paper is to show that the interpolating element-free Galerkin is a suitable technique for solving the nonlinear fourth-order partial differential equations especially extended Fisher–Kolmogorov equation. The numerical experiments confirm the analytical results and show the good efficiency of the interpolating element-free Galerkin method for solving this nonlinear biharmonic equation.
- Published
- 2019
44. Multiscale radial kernels with high-order generalized Strang-Fix conditions
- Author
-
Wenwu Gao and Xuan Zhou
- Subjects
Polynomial ,Applied Mathematics ,Numerical analysis ,ComputingMethodologies_IMAGEPROCESSINGANDCOMPUTERVISION ,010103 numerical & computational mathematics ,Function (mathematics) ,01 natural sciences ,Convolution ,010101 applied mathematics ,Moment (mathematics) ,Linear form ,Theory of computation ,Applied mathematics ,0101 mathematics ,ComputingMethodologies_COMPUTERGRAPHICS ,Variable (mathematics) ,Mathematics - Abstract
The paper provides a general and simple approach for explicitly constructing multiscale radial kernels with high-order generalized Strang-Fix conditions from a given univariate generator. The resulting kernels are constructed by taking a linear functional to the scaled f -form of the generator with respect to the scale variable. Equivalent divided difference forms of the kernels are also derived; based on which, a pyramid-like algorithm for fast and stable computation of multiscale radial kernels is proposed. In addition, characterizations of the kernels in both the spatial and frequency domains are given, which show that the generalized Strang-Fix condition, the moment condition, and the condition of polynomial reproduction in the convolution sense are equivalent to each other. Hence, as a byproduct, the paper provides a unified view of these three classical concepts. These kernels can be used to construct quasi-interpolation with high approximation accuracy and construct convolution operators with high approximation orders, to name a few. As an example, we construct a quasi-interpolation scheme for irregularly spaced data and derived its error estimates and choices of scale parameters of multiscale radial kernels. Numerical results of approximating a bivariate Franke function using our quasi-interpolation are presented at the end of the paper. Both theoretical and numerical results show that quasi-interpolation with multiscale radial kernels satisfying high-order generalized Strang-Fix conditions usually provides high approximation orders.
- Published
- 2019
45. An adaptive local discontinuous Galerkin method for nonlinear two-point boundary-value problems
- Author
-
Mahboub Baccouch
- Subjects
Discretization ,Adaptive mesh refinement ,Applied Mathematics ,Estimator ,010103 numerical & computational mathematics ,Superconvergence ,01 natural sciences ,Prime (order theory) ,010101 applied mathematics ,Rate of convergence ,Discontinuous Galerkin method ,Applied mathematics ,Boundary value problem ,0101 mathematics ,Mathematics - Abstract
In this paper, we propose an adaptive mesh refinement (AMR) strategy based on a posteriori error estimates for the local discontinuous Galerkin (LDG) method for nonlinear two-point boundary-value problems (BVPs) of the form $u^{\prime \prime }=f(x,u),\ x\in [a,b]$ subject to some suitable boundary conditions at the endpoint of the interval [a, b]. We first use the superconvergence results proved in the first part of this paper as reported by Baccouch (Numer. Algorithm. 79(3), 697–718 2018) to show that the significant parts of the local discretization errors are proportional to (p + 1)-degree Radau polynomials, when polynomials of total degree not exceeding p are used. These new results allow us to construct a residual-based a posteriori error estimators which are obtained by solving a local residual problem with no boundary conditions on each element. The proposed error estimates are efficient, reliable, and asymptotically exact. We prove that, for smooth solutions, the proposed a posteriori error estimates converge to the exact errors in the L2-norm with order of convergence p + 3/2. Finally, we present a local AMR procedure that makes use of our local and global a posteriori error estimates. Our proofs are valid for arbitrary regular meshes and for Pp polynomials with p ≥ 1. Several numerical results are presented to validate the theoretical results and to show the efficiency of the grid refinement strategy.
- Published
- 2019
46. Family weak conjugate gradient algorithms and their convergence analysis for nonconvex functions
- Author
-
Zhou Sheng, Gonglin Yuan, and Xiaoliang Wang
- Subjects
Trust region ,Line search ,Applied Mathematics ,Numerical analysis ,Structure (category theory) ,010103 numerical & computational mathematics ,01 natural sciences ,010101 applied mathematics ,Simple (abstract algebra) ,Conjugate gradient method ,Convergence (routing) ,Theory of computation ,0101 mathematics ,Algorithm ,Mathematics - Abstract
It is well-known that conjugate gradient algorithms are widely applied in many practical fields, for instance, engineering problems and finance models, as they are straightforward and characterized by a simple structure and low storage. However, challenging problems remain, such as the convergence of the PRP algorithms for nonconvexity under an inexact line search, obtaining a sufficient descent for all conjugate gradient methods, and other theory properties regarding global convergence and the trust region feature for nonconvex functions. This paper studies family conjugate gradient formulas based on the six classic formulas, PRP, HS, CD, FR, LS, and DY, where the family conjugate gradient algorithms have better theory properties than those of the formulas by themselves. Furthermore, this technique of the presented conjugate gradient formulas can be extended to any two-term conjugate gradient formula. This paper designs family conjugate gradient algorithms for nonconvex functions, which have the following features without other conditions: (i) the sufficient descent property holds, (ii) the trust region feature is true, and (iii) the global convergence holds under normal assumptions. Numerical results show that the given conjugate gradient algorithms are competitive with those of normal methods.
- Published
- 2019
47. Reconstruction algorithms of an inverse source problem for the Helmholtz equation
- Author
-
Ji-Chuan Liu and Xiao-Chen Li
- Subjects
Nonlinear system ,Helmholtz equation ,Salient ,Applied Mathematics ,Numerical analysis ,Theory of computation ,ComputingMethodologies_IMAGEPROCESSINGANDCOMPUTERVISION ,Boundary (topology) ,Cauchy distribution ,Algorithm ,Regularization (mathematics) ,Mathematics - Abstract
In this paper, we study an inverse source problem for the Helmholtz equation from measurements. The purpose of this paper is to reconstruct the salient features of the hidden sources within a body. We propose three stable reconstruction algorithms to detect the number, the location, the size, and the shape of the hidden sources along with compact support from a single measurement of near-field Cauchy data on the external boundary. This problem is nonlinear and ill-posed; thus, we should consider regularization techniques in reconstruction algorithms. We give several numerical experiments to demonstrate the viability of our proposed reconstruction algorithms.
- Published
- 2019
48. A formal construction of a divergence-free basis in the nonconforming virtual element method for the Stokes problem.
- Author
-
Kwak, Do Y. and Park, Hyeokjoo
- Subjects
POSITIVE systems ,DIVERGENCE theorem ,MATHEMATICS ,GENERALIZATION - Abstract
We develop a formal construction of a pointwise divergence-free basis in the nonconforming virtual element method of arbitrary order for the Stokes problem introduced in Zhao et al. (SIAM J. Numer. Anal. 57(6):2730–2759, 2019). The proposed construction can be seen as a generalization of the divergence-free basis in Crouzeix-Raviart finite element space (Brenner, Math. Comp. 55(192):411–437, 1990; Thomasset, 1981) to the virtual element space. Using the divergence-free basis obtained from our construction, we can eliminate the pressure variable from the mixed system and obtain a symmetric positive definite system. Several numerical tests are presented to confirm the efficiency and the accuracy of our construction. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
49. Interval methods of Adams-Bashforth type with variable step sizes
- Author
-
Andrzej Marciniak and Malgorzata A. Jankowska
- Subjects
Applied Mathematics ,Numerical analysis ,010103 numerical & computational mathematics ,01 natural sciences ,010101 applied mathematics ,Exact solutions in general relativity ,Theory of computation ,Initial value problem ,Interval (graph theory) ,Applied mathematics ,0101 mathematics ,Constant (mathematics) ,Mathematics ,Variable (mathematics) ,Linear multistep method - Abstract
In a number of our previous papers, we have proposed interval versions of multistep methods (explicit and implicit), including interval predictor-corrector methods, in which the step size was constant. In this paper, we present interval versions of Adams-Bashforth methods with a possibility to change step sizes. This possibility can be used to obtain interval enclosures of the exact solution with a width given beforehand.
- Published
- 2019
50. Split-step cubic B-spline collocation methods for nonlinear Schrödinger equations in one, two, and three dimensions with Neumann boundary conditions
- Author
-
Luming Zhang and Shanshan Wang
- Subjects
Nonlinear system ,symbols.namesake ,Collocation ,Applied Mathematics ,B-spline ,Numerical analysis ,Neumann boundary condition ,symbols ,Finite difference ,Applied mathematics ,Nonlinear Schrödinger equation ,Schrödinger equation ,Mathematics - Abstract
In this paper, split-step cubic B-spline collocation (SS3BC) schemes are constructed by combining the split-step approach with the cubic B-spline collocation (3BC) method for the nonlinear Schrodinger (NLS) equation in one, two, and three dimensions with Neumann boundary conditions. Unfortunately, neither of the advantages of the two methods can be maintained for the multi-dimensional problems, if one combines them in the usual manner. For overcoming the difficulty, new medium quantities are introduced in this paper to successfully reduce the multi-dimensional problems into one-dimensional ones, which are essential for the SS3BC methods. Numerical tests are carried out, and the schemes are verified to be convergent with second-order both in time and space. The proposed method is also compared with the split-step finite difference (SSFD) scheme. Finally, the present method is applied to two problems of the Bose-Einstein condensate. The proposed SS3BC method is numerically verified to be effective and feasible.
- Published
- 2019
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.