• 대한수학회논문집
• /
• 제34권2호
• /
• pp.657-670
• /
• 2019

In this paper, we consider the preconditioned iterative methods for solving linear complementarity problem associated with an M-matrix. Based on the generalized Gunawardena's preconditioner, two preconditioned SSOR methods for solving the linear complementarity problem are proposed. The convergence of the proposed methods are analyzed, and the comparison results are derived. The comparison results showed that preconditioned SSOR methods accelerate the convergent rate of the original SSOR method. Numerical examples are used to illustrate the theoretical results.

• Journal of applied mathematics & informatics
• /
• 제31권5_6호
• /
• pp.769-784
• /
• 2013

We consider the modulus-based successive overrelaxation method for the linear complementarity problems from the discretization of Black-Scholes American options model. The $H_+$-matrix property of the system matrix discretized from American option pricing which guarantees the convergence of the proposed method for the linear complementarity problem is analyzed. Numerical experiments confirm the theoretical analysis, and further show that the modulus-based successive overrelaxation method is superior to the classical projected successive overrelaxation method with optimal parameter.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 2009년도 춘계 학술대회 제21권1호
• /
• pp.579-580
• /
• 2009

이 논문은 초고강도 I 형보의 quasi-brittle 파괴역학적 수치해석을 수행하였다. Non-holonomic rate 형태로 변형경로에 의존하는 경화-연화거동 관계 방정식을 구성하기 위해 linear Complementarity 방법을 사용하였으며, PATH solver 를 사용하여 LCP 방법의 해를 구하였다.

• 호남수학학술지
• /
• 제19권1호
• /
• pp.157-164
• /
• 1997

We are concerned with three classes of matrices that are relevant to the linear complementary problem. We prove that within the class of $P_0$-matrices, the Q-matrices are precisely the regular matrices and we show that the same characterizations hold for an L-matrix as well, and that the symmetric copositive-plus Q-matrices are precisely those which are strictly copositive.

• 대한수학회논문집
• /
• 제25권4호
• /
• pp.655-669
• /
• 2010

In this paper we propose new primal-dual interior point methods (IPMs) for $P_*(\kappa)$ linear complementarity problems (LCPs) and analyze the iteration complexity of the algorithm. New search directions and proximity measures are defined based on a class of kernel functions, $\psi(t)=\frac{t^2-1}{2}-{\int}^t_1e{^{q(\frac{1}{\xi}-1)}d{\xi}$, $q\;{\geq}\;1$. If a strictly feasible starting point is available and the parameter $q\;=\;\log\;\(1+a{\sqrt{\frac{2{\tau}+2{\sqrt{2n{\tau}}+{\theta}n}}{1-{\theta}}\)$, where $a\;=\;1\;+\;\frac{1}{\sqrt{1+2{\kappa}}}$, then new large-update primal-dual interior point algorithms have $O((1\;+\;2{\kappa})\sqrt{n}log\;n\;log\;{\frac{n}{\varepsilon}})$ iteration complexity which is the best known result for this method. For small-update methods, we have $O((1\;+\;2{\kappa})q{\sqrt{qn}}log\;{\frac{n}{\varepsilon}})$ iteration complexity.

• 대한기계학회논문집
• /
• 제18권8호
• /
• pp.2026-2038
• /
• 1994

A numerical scheme is developed for the analysis of incipient sliding contact with orthotropic friction condition subjected to tangential load and twisting moment. The inherent nonlinearity in the orthotropic friction law has been treated by a polyhedral friction law. Then, a three-dimensional linear complementarity problem(LCP) formulation in an incremental form is obtained, and the existence of a solution is investigated. A Lemke's complementary pivoting algorithm is used for solving the LCP. The scheme is illustrated by spherical contact problems, and the effects of eccentricity of elliptical friction domain on the traction and stick region are discussed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제13권3호
• /
• pp.189-202
• /
• 2009

In this paper we propose new primal-dual interior point algorithms for $P_*(\kappa)$ linear complementarity problems based on a new class of kernel functions which contains the kernel function in [8] as a special case. We show that the iteration bounds are $O((1+2\kappa)n^{\frac{9}{14}}\;log\;\frac{n{\mu}^0}{\epsilon}$) for large-update and $O((1+2\kappa)\sqrt{n}log\frac{n{\mu}^0}{\epsilon}$) for small-update methods, respectively. This iteration complexity for large-update methods improves the iteration complexity with a factor $n^{\frac{5}{14}}$ when compared with the method based on the classical logarithmic kernel function. For small-update, the iteration complexity is the best known bound for such methods.

• Journal of applied mathematics & informatics
• /
• 제29권5_6호
• /
• pp.1511-1523
• /
• 2011

In this paper, we consider the smoothing Newton method for the nonlinear complementarity problems with $P_0$-function. The proposed algorithm is based on a new smoothing function and it needs only to solve one linear system of equations and perform one line search per iteration. Under the condition that the solution set is nonempty and bounded, the proposed algorithm is proved to be convergent globally. Furthermore, the local superlinearly(quadratic) convergence is established under suitable conditions. Preliminary numerical results show that the proposed algorithm is very promising.

• Journal of applied mathematics & informatics
• /
• 제32권1_2호
• /
• pp.211-225
• /
• 2014

A new class of smoothing functions is introduced in this paper, which includes some important smoothing complementarity functions as its special cases. Based on this new smoothing function, we proposed a smoothing Newton method. Our algorithm needs only to solve one linear system of equations. Without requiring the nonemptyness and boundedness of the solution set, the proposed algorithm is proved to be globally convergent. Numerical results indicate that the smoothing Newton method based on the new proposed class of smoothing functions with ${\theta}{\in}(0,1)$ seems to have better numerical performance than those based on some other important smoothing functions, which also demonstrate that our algorithm is promising.

• Journal of applied mathematics & informatics
• /
• 제29권3_4호
• /
• pp.653-668
• /
• 2011

This paper presents a new hybrid method for solving nonlinear complementarity problems with $P_0$-functions. It can be regarded as a combination of smoothing trust region method with ODE-based method and line search technique. A feature of the proposed method is that at each iteration, a linear system is only solved once to obtain a trial step, thus avoiding solving a trust region subproblem. Another is that when a trial step is not accepted, the method does not resolve the linear system but generates an iterative point whose step-length is defined by a line search. Under some conditions, the method is proven to be globally and superlinearly convergent. Preliminary numerical results indicate that the proposed method is promising.