اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدSolving nonlinear equations by a derivative-free form of the Kings family with memory
36
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
In this paper, we present an iterative three-point method with memory based on the family of Kings methods to solve nonlinear equations. This proposed method has eighth order convergence and costs only four function evaluations per iteration which supports the Kung-Traub conjecture on the optimal order of convergence. An acceleration of the convergence speed is achieved by an appropriate variation of a free parameter in each step. This self accelerator parameter is estimated using Newtons interpolation polynomial of fourth degree. The order of convergence is increased from 8 to 12 without any extra function evaluation. Consequently, this method, possesses a high computational efficiency. Finally, a numerical comparison of the proposed method with related methods shows its effectiveness and performance in high precision computations.
قيم البحث
اقرأ أيضاً
Spectral residual methods are powerful tools for solving nonlinear systems of equations without derivatives. In a recent paper, it was shown that an acceleration technique based on the Sequential Secant Method can greatly improve its efficiency and r
obustness. In the present work, an R implementation of the method is presented. Numerical experiments with a widely used test bed compares the presented approach with its plain (i.e. non-accelerated) version that makes part of the R package BB. Additional numerical experiments compare the proposed method with NITSOL, a state-of-the-art solver for nonlinear systems. The comparison shows that the acceleration process greatly improves the robustness of its counterpart included in the existent R package. As a by-product, an interface is provided between R and the consolidated CUTEst collection, which contains over a thousand nonlinear programming problems of all types and represents a standard for evaluating the performance of optimization methods.
In this paper, we propose a $mu$-mode integrator for computing the solution of stiff evolution equations. The integrator is based on a d-dimensional splitting approach and uses exact (usually precomputed) one-dimensional matrix exponentials. We show
that the action of the exponentials, i.e. the corresponding batched matrix-vector products, can be implemented efficiently on modern computer systems. We further explain how $mu$-mode products can be used to compute spectral transformations efficiently even if no fast transform is available. We illustrate the performance of the new integrator by solving three-dimensional linear and nonlinear Schrodinger equations, and we show that the $mu$-mode integrator can significantly outperform numerical methods well established in the field. We also discuss how to efficiently implement this integrator on both multi-core CPUs and GPUs. Finally, the numerical experiments show that using GPUs results in performance improvements between a factor of 10 and 20, depending on the problem.
Differential algebraic Riccati equations are at the heart of many applications in control theory. They are time-depent, matrix-valued, and in particular nonlinear equations that require special methods for their solution. Low-rank methods have been u
sed heavily computing a low-rank solution at every step of a time-discretization. We propose the use of an all-at-once space-time solution leading to a large nonlinear space-time problem for which we propose the use of a Newton-Kleinman iteration. Approximating the space-time problem in low-rank form requires fewer applications of the discretized differential operator and gives a low-rank approximation to the overall solution.
In this paper, several two-grid finite element algorithms for solving parabolic integro-differential equations (PIDEs) with nonlinear memory are presented. Analysis of these algorithms is given assuming a fully implicit time discretization. It is sho
wn that these algorithms are as stable as the standard fully discrete finite element algorithm, and can achieve the same accuracy as the standard algorithm if the coarse grid size $H$ and the fine grid size $h$ satisfy $H=O(h^{frac{r-1}{r}})$. Especially for PIDEs with nonlinear memory defined by a lower order nonlinear operator, our two-grid algorithm can save significant storage and computing time. Numerical experiments are given to confirm the theoretical results.
We introduce a novel spatio-temporal discretization for nonlinear Fokker-Planck equations on the multi-dimensional unit cube. This discretization is based on two structural properties of these equations: the first is the representation as a gradient
flow of an entropy functional in the $L^2$-Wasserstein metric, the second is the Lagrangian nature, meaning that solutions can be written as the push forward transformation of the initial density under suitable flow maps. The resulting numerical scheme is entropy diminishing and mass conserving. Further, the scheme is weakly stable, which allows us to prove convergence under certain regularity assumptions. Finally, we present results from numerical experiments in space dimension $d=2$.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
