The influence of isolated largest eigenvalues on the numerical convergence of the CG method

2000 ◽  
Vol 101 (4) ◽  
pp. 3231-3236
Author(s):  
A. Y. Yeremin ◽  
I. E. Kaporin
Mathematics ◽  
2021 ◽  
Vol 9 (2) ◽  
pp. 160
Author(s):  
Rafael Company ◽  
Vera N. Egorova ◽  
Lucas Jódar

In this paper, we consider random hyperbolic partial differential equation (PDE) problems following the mean square approach and Laplace transform technique. Randomness requires not only the computation of the approximating stochastic processes, but also its statistical moments. Hence, appropriate numerical methods should allow for the efficient computation of the expectation and variance. Here, we analyse different numerical methods around the inverse Laplace transform and its evaluation by using several integration techniques, including midpoint quadrature rule, Gauss–Laguerre quadrature and its extensions, and the Talbot algorithm. Simulations, numerical convergence, and computational process time with experiments are shown.


2019 ◽  
Vol 2019 ◽  
pp. 1-10
Author(s):  
Darae Jeong ◽  
Yibao Li ◽  
Chaeyoung Lee ◽  
Junxiang Yang ◽  
Yongho Choi ◽  
...  

In this paper, we propose a verification method for the convergence rates of the numerical solutions for parabolic equations. Specifically, we consider the numerical convergence rates of the heat equation, the Allen–Cahn equation, and the Cahn–Hilliard equation. Convergence test results show that if we refine the spatial and temporal steps at the same time, then we have the second-order convergence rate for the second-order scheme. However, in the case of the first-order in time and the second-order in space scheme, we may have the first-order or the second-order convergence rates depending on starting spatial and temporal step sizes. Therefore, for a rigorous numerical convergence test, we need to perform the spatial and the temporal convergence tests separately.


2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
S. M. Nuugulu ◽  
F. Gideon ◽  
K. C. Patidar

AbstractDividend paying European stock options are modeled using a time-fractional Black–Scholes (tfBS) partial differential equation (PDE). The underlying fractional stochastic dynamics explored in this work are appropriate for capturing market fluctuations in which random fractional white noise has the potential to accurately estimate European put option premiums while providing a good numerical convergence. The aim of this paper is two fold: firstly, to construct a time-fractional (tfBS) PDE for pricing European options on continuous dividend paying stocks, and, secondly, to propose an implicit finite difference method for solving the constructed tfBS PDE. Through rigorous mathematical analysis it is established that the implicit finite difference scheme is unconditionally stable. To support these theoretical observations, two numerical examples are presented under the proposed fractional framework. Results indicate that the tfBS and its proposed numerical method are very effective mathematical tools for pricing European options.


Author(s):  
Никита Сергеевич Олейник ◽  
Владислав Юрьевич Щеколдин

Рассмотрена задача выявления аномальных наблюдений в данных больших размерностей на основе метода многомерного шкалирования с учетом возможности построения качественной визуализации данных. Предложен алгоритм модифицированного метода главных проекций Торгерсона, основанный на построении подпространства проектирования исходных данных путем изменения способа факторизации матрицы скалярных произведений при помощи метода анализа кумулятивных кривых. Построено и проанализировано эмпирическое распределение F -меры для разных вариантов проектирования исходных данных Purpose. Purpose of the article. The paper aims at the development of methods for multidimensional data presentation for solving classification problems based on the cumulative curves analysis. The paper considers the outlier detection problem for high-dimensional data based on the multidimensional scaling, in order to construct high-quality data visualization. An abnormal observation (or outlier), according to D. Hawkins, is an observation that is so different from others that it may be assumed as appeared in the sample in a fundamentally different way. Methods. One of the conceptual approaches that allow providing the classification of sample observations is multidimensional scaling, representing by the classical Orlochi method, the Torgerson main projections and others. The Torgerson method assumes that when converting data to construct the most convenient classification, the origin must be placed at the gravity center of the analyzed data, after which the matrix of scalar products of vectors with the origin at the gravity center is calculated, the two largest eigenvalues and corresponding eigenvectors are chosen and projection matrix is evaluated. Moreover, the method assumes the linear partitioning of regular and anomalous observations, which arises rarely. Therefore, it is logical to choose among the possible axes for designing those that allow obtaining more effective results for solving the problem of detecting outlier observations. A procedure of modified CC-ABOD (Cumulative Curves for Angle Based Outlier Detection) to estimate the visualization quality has been applied. It is based on the estimation of the variances of angles assumed by particular observation and remaining observations in multidimensional space. Further the cumulative curves analysis is implemented, which allows partitioning out groups of closely localized observations (in accordance with the chosen metric) and form classes of regular, intermediate, and anomalous observations. Results. A proposed modification of the Torgerson method is developed. The F1-measure distribution is constructed and analyzed for different design options in the source data. An analysis of the empirical distribution showed that in a number of cases the best axes are corresponding to the second, third, or even fourth largest eigenvalues. Findings. The multidimensional scaling methods for constructing visualizations of multi-dimensional data and solving problems of outlier detection have been considered. It was found out that the determination of design is an ambiguous problem.


2018 ◽  
Vol 7 (1) ◽  
pp. 1 ◽  
Author(s):  
F. MUHAMMAD ZAIN ◽  
M. GARDA KHADAFI ◽  
P. H. GUNAWAN

The diffusion equation or known as heat equation is a parabolic and linear type of partial differential equation. One of the numerical method to approximate the solution of diffusion equations is Finite Difference Method (FDM). In this study, the analysis of numerical convergence of FDM to the solution of diffusion equation is discussed. The analytical solution of diffusion equation is given by the separation of variables approach. Here, the result show the convergence of rate the numerical method is approximately approach 2. This result is in a good agreement with the spatial error from Taylor expansion of spatial second derivative.


2017 ◽  
Vol 21 (3) ◽  
pp. 835-866 ◽  
Author(s):  
Meng Wu ◽  
Bernard Mourrain ◽  
André Galligo ◽  
Boniface Nkonga

AbstractMotivated by the magneto hydrodynamic (MHD) simulation for Tokamaks with Isogeometric analysis, we present splines defined over a rectangular mesh with a complex topological structure, i.e., with extraordinary vertices. These splines are piecewise polynomial functions of bi-degree (d,d) and parameter continuity. And we compute their dimension and exhibit basis functions called Hermite bases for bicubic spline spaces. We investigate their potential applications for solving partial differential equations (PDEs) over a physical domain in the framework of Isogeometric analysis. For instance, we analyze the property of approximation of these spline spaces for the L2-norm; we show that the optimal approximation order and numerical convergence rates are reached by setting a proper parameterization, although the fact that the basis functions are singular at extraordinary vertices.


2018 ◽  
Author(s):  
Jesse A Sharp ◽  
Alexander P Browning ◽  
Tarunendu Mapder ◽  
Kevin Burrage ◽  
Matthew J Simpson

AbstractAcute myeloid leukaemia (AML) is a blood cancer affecting haematopoietic stem cells. AML is routinely treated with chemotherapy, and so it is of great interest to develop optimal chemotherapy treatment strategies. In this work, we incorporate an immune response into a stem cell model of AML, since we find that previous models lacking an immune response are inappropriate for deriving optimal control strategies. Using optimal control theory, we produce continuous controls and bang-bang controls, corresponding to a range of objectives and parameter choices. Through example calculations, we provide a practical approach to applying optimal control using Pontryagin’s Maximum Principle. In particular, we describe and explore factors that have a profound influence on numerical convergence. We find that the convergence behaviour is sensitive to the method of control updating, the nature of the control, and to the relative weighting of terms in the objective function. All codes we use to implement optimal control are made available.


2012 ◽  
Vol 170-173 ◽  
pp. 1238-1242
Author(s):  
Xue Wei Li ◽  
Xin Yuan ◽  
Xiao Wei Li

Abstract. Combined the strength reduction method with ABAQUS, the development of the slope plastic strain of different reduction coefficient is obtained by constantly adjusting reduction coefficient to change the strength index of the soil. The reduction coefficient is obtained from the criterion of numerical convergence and displacement mutation and plastic zone breakthrough. Through the analysis and comparison with the results, the reduction coefficient by the criterion of displacement mutation is consistent with the result of the criterion of plastic zone breakthrough. The reduction coefficient is the safety coefficient of the slope, and compared and analyzed with the slope factor of limit equilibrium method Bishop. The result shows that the displacement mutation and the plastic zone breakthrough as criterions to judge the slope instability are reasonable.


Sign in / Sign up

Export Citation Format

Share Document